Quantidade permissível de forças sobre os íons em uma estrutura relaxada

Question

Uma estrutura de equilíbrio estaria em um mínimo local, mas isso não acontece com todos os cálculos de relaxamento. Estruturas relaxadas tendem a ter algumas forças resultantes sobre os íons. Que quantidade dessas forças pode ser aceita? É 0,05 eV /$\unicode{xC5}$ aceitável?

Answer 1

Isso depende completamente do que você está procurando modelar e do tipo de material que está modelando. Uma regra prática padrão geralmente não é superior a 0,05 eV / Å para a força líquida máxima em qualquer átomo, mas mesmo isso pode ser muito grande, particularmente (mas não exclusivamente) no caso de materiais altamente flexíveis ou cristais moleculares . Pessoalmente, gosto de 0,03 eV / Å, mas isso é novamente baseado na experiência de uma classe muito específica de materiais que estudo (MOFs). Para outros materiais, isso pode não ser suficiente. Também dependerá muito da propriedade de interesse, com as energias sendo menos sensíveis do que os modos vibracionais como um exemplo. Infelizmente, não há uma resposta adequada para essa pergunta. Você simplesmente tem que testá-lo e ver.

Answer 2

Como outros mencionaram, não existe uma 'regra prática', mas eu faço cálculos do phonon (PHONOPY e DFPT) e fiz algumas verificações há muito tempo para convergir as energias do phonon. Descobri que ~ 0,001 eV / A é mais do que baixo o suficiente para todos os casos que verifiquei. Observe que, embora o relaxamento para alta precisão seja caro, uma vez que se aproxima do mínimo, as posições geralmente convergem mais rapidamente. Observe que, em alguns casos, o corte de energia (para códigos de onda plana que eu uso) pode ter que ser suficientemente alto para remover o ruído nas forças se você quiser um valor tão baixo.