Quantité admissible de forces sur les ions dans une structure détendue
Une structure d'équilibre serait à un minimum local, mais cela ne se produit pas avec tous les calculs de relaxation. Les structures relâchées ont tendance à avoir des forces résultantes sur les ions. Quelle quantité de ces forces peut être acceptée? Est de 0,05 eV /$\unicode{xC5}$ acceptable?
Réponses
Cela dépend complètement de ce que vous cherchez à modéliser et du type de matériau que vous modélisez. Une règle empirique standard n'est généralement pas supérieure à 0,05 eV / Å pour la force nette maximale sur n'importe quel atome, mais même celle-ci peut souvent être trop grande, en particulier (mais pas exclusivement) dans le cas de matériaux hautement flexibles ou de cristaux moléculaires . Personnellement, j'aime 0,03 eV / Å, mais c'est encore une fois basé sur l'expérience d'une classe très spécifique de matériaux que j'étudie (MOF). Pour d'autres matériaux, cela peut très bien ne pas être suffisant. Cela dépendra aussi fortement de la propriété d'intérêt, les énergies étant moins sensibles que les modes vibrationnels à titre d'exemple. Il n'y a malheureusement pas de réponse appropriée à cette question. Vous devez simplement le tester et voir.
Comme d'autres l'ont mentionné, il n'y a pas de «règle empirique», mais je fais des calculs de phonons (PHONOPIE et DFPT) et j'ai fait quelques vérifications il y a longtemps pour faire converger les énergies des phonons. J'ai trouvé ~ 0,001 eV / A plus que suffisamment bas pour tous les cas que j'ai vérifiés. Notez que, si la relaxation à haute précision est coûteuse, une fois qu'elle se rapproche du minimum, les positions convergent généralement plus rapidement. Notez que dans certains cas, la coupure d'énergie (pour les codes d'ondes planes que j'utilise) peut devoir être suffisamment élevée pour supprimer le bruit dans les forces si vous voulez une valeur aussi basse.