atomsk --density <fichier> <propriété> <1d|2d|3d> [<x|y|z>] <σ> [options]
Ce mode calcule la densité d'une <propriété> donnée, en employant une grille d'éléments finis.
Les valeurs discrètes (atomiques) de la <propriété> donnée sont remplacées par des fonctions gaussiennes de varianceσ2 qui sont évaluées sur une grille d'éléments finis. La density de la <property> dans un élément de la grille est définie comme la somme des contributions des gaussiennes de tous les atomes.
Les paramètres associés à ce mode sont :
Un profil de densité 1-D remplacera les atomes par des fonctions gaussiennes et calculera leur densité suivant l'<axe> donné. Ceci peut être employé pour déterminer la densité d'un élément ou d'une propriété le long d'une direction.
Une carte de densité 2-D remplacera les atomes par des fonctions gaussiennes en deux dimensions, et calculera leur densité dans le plan normal à l'<axe> donné. Ceci peut être utilisé pour appliquer un flou gaussien aux positions atomiques et produire une image comparable à un cliché de microscopie électronique, ou bien pour déterminer la densité d'une propriété dans un plan.
Une carte de densité 3-D remplacera les atomes par des fonctions gaussiennes 3-D et calculera la densité dans tout le volume. In addition, the positions where the density is significantly higher or lower than the average will be written into file(s) (see below). This can be useful to find the positions of interstitial atoms, or the positions of vacancies.
Le système est toujours considéré comme périodique, et les répliques périodiques des atomes contribueront à la densité.
Les fichiers de sortie sont les suivants :
Dans le cas d'un calcul de densité 2-D, le fichier suivant sera aussi produit :
Enfin, dans le cas d'un calcul de densité 3-D, le fichier suivant sera aussi produit :
Si ce mode est utilisé avec des options alors ces options seront appliquées avant que la densité ne soit calculée.
Si une sélection a été définie (avec l'option -select) alors seuls les atomes sélectionnés contribueront à la densité.
Notez que ce mode suppose que tous les atomes se trouvent dans la boîte. Si ce n'est pas le cas alors le calcul peut donner de mauvais résultats. Les atomes peuvent être replacés dans la boîte grâce à l'option -wrap.
atomsk --density monfichier.cfg Fe 1d X 0.6
Ceci va lire le fichier monfichier.cfg, et calculera la densité linéique d'atomes de fer (Fe) le long de l'axe X. Chaque atome de fer sera remplacé par une fonction gaussienne de variance σ2=0.62. Le résultat final sera écrit dans le fichier monfichier_density.dat.
atomsk --density Al.xsf mass 2d Z 2.0
Ceci va lire le fichier Al.xsf, et calculera la densité de masse dans le plan normal à l'axe Z, en utilisant des fonctions gaussiennes 2-D telles que σ=2. Le résultat sera écrit dans le fichier Al_density.dat. Un script Gnuplot nommé Al_density.gp sera aussi produit pour faciliter la visualisation des données.
atomsk --density MD_5000.lmc mass 3d 2.0
Ceci va lire le fichier MD_5000.lmc, et calculera la densité de masse en 3-D dans tout le système, en utilisant des fonctions gaussiennes 3-D telles que σ=2. Notez qu'aucun <axe> ne doit apparaître dans les paramètres. Les résultats seront écrits dans le fichier MD_5000_density.dat. In addition, if the variations of mass are significant then the positions where the density is high will be written into MD_5000_interstitials.xsf, and the positions where it is low will be written into MD_5000_vacancies.xsf.