Допустимое количество сил на ионы в релаксированной структуре
Равновесная структура будет находиться в локальном минимуме, но этого не происходит при всех расчетах релаксации. Расслабленные структуры, как правило, имеют некоторые результирующие силы на ионах. Какое количество этих сил можно принять? 0,05 эВ /$\unicode{xC5}$ приемлемо?
Ответы
Это полностью зависит от того, что вы хотите моделировать, и от типа моделируемого материала. Стандартное практическое правило обычно не превышает 0,05 эВ / Å для максимальной чистой силы, действующей на любой атом, но даже это часто может быть слишком большим, особенно (но не исключительно) в случае очень гибких материалов или молекулярных кристаллов. . Лично мне нравится 0,03 эВ / Å, но это опять же основано на опыте работы с очень специфическим классом материалов, которые я изучаю (MOF). Для других материалов этого может быть недостаточно. Это также будет сильно зависеть от интересующего свойства, например, энергии будут менее чувствительны, чем колебательные моды. К сожалению, на этот вопрос нет правильного ответа. Вам просто нужно проверить это и увидеть.
Как уже упоминали другие, не существует «практического правила», но я провожу расчеты фононов (PHONOPY и DFPT) и много лет назад провела некоторые проверки, чтобы сходиться энергии фононов. Я обнаружил, что ~ 0,001 эВ / А более чем достаточно для всех проверенных мной случаев. Обратите внимание, что, хотя снижение точности до высокой является дорогостоящим, как только оно приближается к минимуму, позиции обычно сходятся быстрее. Обратите внимание, что в некоторых случаях ограничение энергии (для кодов плоских волн, которые я использую), возможно, должно быть достаточно высоким, чтобы удалить шум в силах, если вы хотите такое низкое значение.