Densité d'état du Nickel
La précision de la densité d'état dépend fortement de la grille de points k
utilisée. Elle doit être suffisamment dense pour paver finement la première zone
de Brillouin. Des méthodes spéciales d'intégration de la zone de
Brillouin comme la tetahedron method (ISMEAR = -5) peuvent être utilisée.
L'utilisation d'un grand nombre de points k étant coûteuse en temps de calcul, on peut procéder en plusieurs étapes :
1 - Relaxer la densité électronique sans changer la grille de points k en
choisissant une valeur de MAGMOM.
2 - Recopier dans un nouveau dossier les fichiers INCAR, CONTCAR, POTCAR,
KPOINTS et CHGCAR. Par exemple :
[vasptp@pyrene] > mkdir DOS
[vasptp@pyrene] > cp INCAR POSCAR POTCAR KPOINTS CHGCAR DOS/
[vasptp@pyrene] > cd DOS/
Le fichier CHGCAR contient la densité électronique en chaque point de la maille.
3 - Choisir une grille de points k de 21 x 21 x 21 et relancer le calcul en
ajoutant ICHARG = 11 dans le fichier INCAR.
ICHARG = 11 signifie que le calcul suit un processus SCF contraint qui ne modifiera
pas la densité électronique qui est lue dans le fichier CHGCAR. Pour chaque
point k de la nouvelle grille demandée on calcule donc les orbitales cristallines
et les valeurs propres associées de telle sorte qu'elle reproduisent la densité
électronique fixée au départ. Cette technique permet d'augmenter le nombre
de points k à moindre coût dans le but d'obtenir une DOS plus précise.
Pour visualiser la DOS, utiliser la commande v :
[vasptp@pyrene] v tdos
Observer la dissymétrie entre la DOS des électrons et , caractéristique d'un système ferromagnétique.
4 - Pour avoir plus d'information, on peut visualiser les contributions des différentes OA du nickel à la DOS totale. Il faut alors tracer les densités d'états projetées.
[vasptp@pyrene] v pdos
5 - Refaire le calcul avec une autre valeur de MAGMOM pour voir l'impact sur
la densité d'état obtenues.
Voir les résultats de cette partie.