Nom : fichier Quantum Espresso PWscf
Extension : pw
Spécification : docuentation de Quantum Espresso
Programmes de visualisation : VESTA, XCrySDen
Atomsk peut lire les fichiers au format Quantum Espresso PWscf (*.pw). Si les coordonnées sont réduites (c-à-d si elles sont dans la section ATOMIC_POSITIONS crystal) alors Atomsk les convertira immédiatement en coordonnées cartésiennes. Pour retouver des coordonnées réduites il est possible d'utiliser l'option -frac.
Par défaut dans le format Quantum Espresso PWscf, les quantités dont les dimensions ne sont pas définies explicitement sont en unités atomiques, et donc les vecteurs de boîte et positions atomiques sont en Bohrs. Ainsi, si les vecteurs de boîte et les coordonnées atomiques sont exprimées dans le même système d'unités (par ex. tous en Bohrs, ou tous en Å) alors Atomsk conserve cette unité (autrement dit les unités ne sont pas automatiquement converties en Å ni en une autre unité). Cependant, si aucune unité n'est spécifiée pour les vecteurs de boîte (ce qui signifie qu'ils sont en Bohrs), tandis que les positions atomiques sont exprimées en Å, alors les vecteurs de boîte seront convertis en angströms pour rétablir la cohérence. Pour convertir les unités il est aussi possible d'utiliser l'option -unit.
Lorsqu'il écrit un fichier au format PW, Atomsk utiliser un format similaire au suivant :
&CONTROL
title = '<commentaire>'
pseudo_dir = '/votre/chemin/vers/les/pseudopotentiels'
calculation = 'scf'
/
&SYSTEM
nat= <NP>
ntyp= <Nspecies>
ibrav= 0
ecutwfc= 20.0
/
&ELECTRONS
mixing_beta = 0.7
conv_thr = 1.0d-8
/
&IONS
/
&CELL
/
ATOMIC_SPECIES
<sp1> <mass1> sp1.fixme.upf
...
<spN> <massN> spN.fixme.upf
CELL_PARAMETERS angstrom
<H(1,1)> <H(1,2)> <H(1,3)>
<H(2,1)> <H(2,2)> <H(2,3)>
<H(3,1)> <H(3,2)> <H(3,3)>
ATOMIC_POSITIONS angstrom
<atom1> <x1> <y1> <z1> [<fixx1> <fixy1> <fixz1>]
<atom2> <x2> <y2> <z2> [<fixx2> <fixy2> <fixz2>]
... ... ...
<atomN> <xN> <yN> <zN> [<fixxN> <fixyN> <fixzN>]
K_POINTS automatic
2 2 2 0 0 0
ATOMIC_FORCES
<atom1> <fx1> <fy1> <fz1>
<atom2> <fx2> <fy2> <fz2>
... ... ...
<atomN> <fxN> <fyN> <fzN>
Un tel fichier est juste un brouillon et devra être édité, il n'est PAS suffisant pour une simulation.
Premièrement, le chemin vers les pseudopotentiels (suivant le mot-clé "pseudo_dir") correspond à la variable d'environnement "$ESPRESSO_PSEUDO". Si celle-ci n'est pas définie, alors "pseudo_dir" a la valeur par défaut "$HOME/espresso/pseudo". Il appartient à l'utilisateur de vérifier que ce chemin est correct.
Deuxièmement, pour remplir la section "ATOMIC_SPECIES", Atomsk cherche des fichiers de pseudopotentiels dans le dossier "pseudo_dir" défini auparavant. Pour chaque élément présent dans le système, Atomsk utilise le nom du premier fichier correspondant à cet élément. Cela signifie que, si le chemin "pseudo_dir" est bien défini et qu'il ne contient qu'un seul fichier de pseudopotentiel pour chaque élément (le pseudopotentiel que vous souhaitez utiliser), alors les noms de fichiers de pseudopotentiels seront déjà correctement définis dans le fichier *.pw écrit par Atomsk. Si le dossier contient plusieurs fichiers de pseudopotentiels pour le même élément, alors le nom du premier fichier sera utilisé (ce qui n'est peut-être pas ce que vous souhaitez). Si aucun fichier n'est trouvé dans ce dossier, un nom de fichier fictif sera utilisé, contenant la chaîne de caractère "fixme". Dans tous les cas, il appartient à l'utilisateur de vérifier que les noms de fichiers de pseudopotentiels sont corrects et correspondent à ses besoins.
Troisièmement, les paramètres du calcul sont définis à certaines valeurs par défaut : calculation='scf', ecutwfc=20, grille de points k 2x2x2, et ainsi de suite. Ces valeurs sont arbitraires et ne doivent pas être utilisées telles quelles ! Il appartient à l'utilisateur d'utiliser des valeurs qui correspondent à ses besoins.
Enfin, si les forces sur les atomes sont définies, alors elles seront écrites dans la section "ATOMIC_FORCES". Sinon cette section n'apparaîtra pas.
Par défaut Atomsk utilise ATOMIC_POSITIONS crystal si les positions des atomes sont exprimées en coordonnées réduites (par ex. si l'option -frac a été utilisée), ou ATOMIC_POSITIONS angstrom si elles sont en coordonnées cartésiennes. De fait, si les coordonnées atomiques et les vecteurs de boîte sont en réalité en Bohrs, les lignes ATOMIC_POSITIONS ansgtrom et CELL_PARAMETERS angstrom seront incorrectes et devront être modifiées à la main; le script qepw_bohr.sh fourni avec Atomsk permet de réaliser ce changement.
i des propriétés auxiliaires existent, alors seules les directions dans lesquelles les atomes sont fixés (fixx, fixy, fixz) sont utilisées et écrites pour chaque atome.