Метод минимизации энергии. Основы классической молекулярной динамики

Содержание

Слайд 2

Минимизация энергии (молекулярная статика) Назначение: определение равновесной структуры систем многих атомов

Минимизация энергии (молекулярная статика)

Назначение: определение равновесной структуры систем многих атомов

Используемые методы:

математические методы минимизации функций многих переменных; метод МД релаксации

Исх. структура Минимизация Равнов. структура

Слайд 3

Трудности минимизации энергии У сложной функции много локальных минимумов, невозможно однозначно найти глобальный минимум

Трудности минимизации энергии

У сложной функции много локальных минимумов,
невозможно однозначно найти глобальный

минимум
Слайд 4

Слайд 5

Суть метода молекулярной динамики МД – наиболее универсальный, мощный метод моделирования

Суть метода молекулярной динамики

МД – наиболее универсальный, мощный метод моделирования атомной

структуры материалов и процессов, происходящих в материалах

Исх. состояние «Путешествие по фазо- Анализ, выводы
ri, vi вому пространству»

Структура,
Т/д свойства
(энергия, энтропия,
теплоемкость,…)
Кинетические свойства
(коэф-т диффузии,
теплопроводность,…)
Мех-змы деформации,
Фазовые переходы,

Слайд 6

Вехи развития МД Alder B.J., Weinwright T.E. 1957. Phase transition for

Вехи развития МД

Alder B.J., Weinwright T.E. 1957. Phase transition for a

hard sphere system. IBM-704 (4 Кфлоп)
Gibson J.B., Goland A.N., Milgram M., Vineyard G.-H. 1960. Dynamics of radiation damage. 500 атомов. IBM-704, 1 мин. на шаг МД
Rahman A. 1964. Correlations in the motion of atoms in liquid argon. 864 атома.
Parinello M., Rahman A. 1981. Polymorphic transitions in single crystals: a new MD method.
Nosé S. 1984. A MD method for simulations in the canonical ensemble.
Roth J., Gähler F., Trebin H.-R. 2000. A molecular dynamics run with 5.180.116.000 particles. 5×109 атомов. (Мощности компьютеров ≈1014-1015 флоп)
Слайд 7

Современные возможности МД T.C. Germann, K. Kadau. Trillion-atom molecular dynamics becomes

Современные возможности МД

T.C. Germann, K. Kadau. Trillion-atom molecular dynamics
becomes a reality.

Int. J. Modern Phys. 2008.
Los Alamos National Laboratory

Суперкомпьютер: Blue Gene/L (212992 процессора IBM 700 МГц) в Lawrence Livermore Nat. Lab.
Общий объем памяти: 72 ТБ

Требуемая память на 1 атом
3 вектора (радиус-вектор, скорость, сила) – 9 чисел по 4 байта
2 целых числа (тип атома и номер атома) – 2 числа по 4 байта
Итого 44 байт
44 ТБ занимают 1012 (1 триллион) атомов

Система занимает куб со стороной 2,5 мкм
Проведено моделирование поведения в течение 10 пс

Слайд 8

Основные задачи, решаемые с помощью МД Жидкости: равновесные, неравновесные, простые, многокомпонентные,

Основные задачи, решаемые с помощью МД

Жидкости: равновесные, неравновесные, простые, многокомпонентные,

вязкость, теплопроводность, кипение,…
Дефекты в кристаллах: атомная структура, энергия, напряжения вакансий, межузельных атомов, дислокаций, дефектов упаковки, границ зерен…
Процессы в твердых телах: пластическая деформация, разрушение, диффузия, трения
Фазовые превращения, в том числе между агрегатными состояниями одного и того же вещества, построение фазовых диаграмм
Процессы нанотехнологии: процессы на поверхности твердых тел (перестройка поверхности, осаждение…), структура и свойства кластеров и наночастиц, больших молекул, в том числе биологических…
Слайд 9

Ограничения классической МД ħ=1.05×10-34 Дж⋅с, b=3 ×10-10 м Ограничения, связанные с возможностями интегрирования уравнений движения:

Ограничения классической МД

ħ=1.05×10-34 Дж⋅с, b=3 ×10-10 м

Ограничения, связанные с возможностями интегрирования

уравнений движения:
Слайд 10

Инициализация систем для моделирования в МД Описание потенциала межатомного взаимодействия; Задание

Инициализация систем для моделирования в МД

Описание потенциала межатомного взаимодействия;
Задание

исходного состояния, то есть координат и скоростей частиц;
Задание граничных условий
Слайд 11

Роль поверхности в свойствах атомных систем R Ns N С уменьшением

Роль поверхности в свойствах атомных систем

R

Ns

N

С уменьшением R влияние поверхностных атомов

возрастает.
Для моделирования поведения макроскопических систем или дефектов в макросистемах необходимо накладывать специальные условия на атомы на границе моделируемой системы, называемые граничными условиями.
Слайд 12

Периодические граничные условия

Периодические граничные условия

Слайд 13

Правило ближайшей частицы H>2rc j j' i rc- радиус обрезания потенциала

Правило ближайшей частицы

H>2rc

j

j'

i

rc- радиус обрезания потенциала

H>2rc

из всех пар, которые составляет

частица i в ячейке и все образы другой частицы j, взаимодействует не более чем одна пара
Действительно, одно из расстояний между i и j или j’ будет больше rc

Правило ближайшей частицы: из всех возможных образов частицы j мы оставляем только ближайший, выбрасывая все остальные. Только ближайшая частица является кандидатом для взаимодействия, все остальные не взаимодействуют. Правило сильно упрощает программу МД и используется повсеместно. Однако размер расчетной ячейки во всех направлениях, в которых наложены ПГУ, должен превышать удвоенный радиус обрезания потенциала, 2rc .

Слайд 14

Методы интегрирования уравнений движения Ошибки при решении уравнений движения 1. Ошибки

Методы интегрирования уравнений движения

Ошибки при решении уравнений движения
1. Ошибки

отбрасывания (усечения), связанные с неточностью метода конечных разностей по сравнению с истинным решением. Методы конечных разностей основаны разложении в ряд Тейлора, усеченный на некотором члене, откуда и происходит название ошибок. Присущи алгоритму решения.
2. Ошибки округления, связанные с реализацией алгоритма. Например, они связаны с конечным числом цифр в представлении числа в компьютере.
Слайд 15

Алгоритм Верле

Алгоритм Верле

Слайд 16

Список соседей При расчете взаимодействий атома i учитываются только атомы, находящиеся

Список соседей

При расчете взаимодействий атома i учитываются только атомы, находящиеся

в сфере радиуса r1, которые вносятся в список соседей этого атома; через определенное число шагов список обновляется.
Слайд 17

Расчет термодинамических величин Средняя потенциальная энергия Средняя кинетическая энергия Полная энергия Температура

Расчет термодинамических величин

Средняя потенциальная энергия
Средняя кинетическая энергия
Полная энергия
Температура

Слайд 18

Калорическая кривая

Калорическая кривая

Слайд 19

Определение температуры плавления твердого тела По определению, температура плавления – это

Определение температуры плавления твердого тела

По определению, температура плавления – это температура, при

которой твердая и жидкая фазы сосуществуют, имея одинаковую свободную энергию.

Точка мех. неуст-сти кр-ла
Tкр≈1.4Tпл

Т

Ж