Продолжительность молекулярно-динамического моделирования, необходимая для вычисления коэффициента самодиффузии при миграции различных точечных дефектов в никеле
В работе проведена оценка продолжительности молекулярно-динамического эксперимента, необходимой для расчета коэффициента самодиффузии при миграции различных точечных дефектов: вакансии, бивакансии, собственного межузельного атома, атома водорода. Перечисленные дефекты обладают различной подвижностью, в результате чего смещения атомов, возникающие вследствие миграции дефекта, имеют разную интенсивность. Погрешность определения коэффициента диффузии связана с точностью определения среднеквадратических изменений координат атомов, которая в свою очередь повышается с ростом продолжительности молекулярно-динамического эксперимента t, температуры Т и подвижности дефекта, инициирующего диффузию. Для описания межатомных взаимодействий в работе использовался многочастичный потенциал Клери - Розато. Показано, что для расчета коэффициента диффузии при миграции вакансии, бивакансии и межузель-ного атома при температурах выше 0,6 от температуры плавления достаточно моделирования в течение 100 пс. При расчете коэффициента диффузии примеси в кристалле металла, например водорода, одного примесного атома, как показано в настоящей работе, недостаточно. В данном случае снизить погрешность определения среднеквадратических смещений примесных атомов можно путем введения большого числа атомов примеси.