|
|
Потенциалы для ковалентных кристалловЭлементы четвертой группы (C, Ge, Si) образуют решетку алмаза с четырьмя ковалентными связями на каждый атом, которые направлены к вершинам тетраэдра и образуют угол 109,5° друг с другом. Главным требованием к потенциалам для этих элементов является описание устойчивости этой решетки. Один из первых потенциалов для ковалентных кристаллов был предложен Стиллингером и Вебером в 1985 г. Направленность связей в этом потенциале вводится с помощью трехчастичного члена в разложении потенциала:
Первый множитель в выражении для парной части потенциала (1.38) имеет форму потенциала ЛД, а второй, экспоненциальный, представляет собой “обрезающий” член, который плавно сводит потенциал к нулю на расстоянии
где
При Потенциал Стиллингера-Вебера завоевал большую популярность при моделировании кремния и был использован во многих исследованиях. Этому способствовали относительная простота потенциала и довольно реалистическое описание кристаллического кремния. Однако этот потенциал имеет несколько серьезных недостатков. Трехчастичный член определяет только одну равновесную конфигурацию, поэтому его трудно распространить на углерод, для которого существует множество равновесных углов: 180°, 120° и 109.47°, благодаря чему углерод имеет множество модификаций, как графит, алмаз, фуллерены, нанотрубки и т.д. Кроме того, слишком жесткая установка тетраэдрического расположения связей приводит к неправильному описанию релаксации на поверхности и около дефектов, где координация нарушена. Поэтому для элементов IV группы были разработаны так называемые потенциалы кратной связи (bond order potentials), которые основаны на учете зависимости прочности связи от локального окружения (потенциалы Терсоффа-Абеля для Si и Ge, потенциал Бреннера для C). Эти потенциалы являются в настоящее время наиболее употребительными при моделировании ковалентных кристаллов.
Метод минимизации энергии. Построенная для моделирования атомная система в исходном состоянии, как правило, имеет конфигурацию, далеко не соответствующую минимуму потенциальной энергии взаимодействия атомов. Если эту систему предоставить самой себе, она будет релаксировать к состоянию с минимумом энергии. Во многих исследованиях, направленных, например, не на выяснение динамического поведения, а на определение структуры дефектов, нахождение такого состояния является основной задачей. Такая задача, например, возникает, при исследовании структуры точечных дефектов или дислокаций или границ зерен. Метод компьютерного моделирования, реализующий нахождение конфигурации атомной системы, соответствующей минимуму потенциальной энергии, называется методом минимизации энергии или, очень часто, применительно к атомному моделированию он называется молекулярной статикой. Физически структура, которая определяется методом минимизации энергии, представляет собой равновесную структуру, которую имеет атомная система при температуре Т = 0 К. Математически задача нахождения равновесной структуры состоит в минимизации полной потенциальной энергии взаимодействия системы Таких алгоритмов имеется много. Практически все они основаны на расчете градиента минимизируемой функции (в данном случае Одним из широко используемых в программах атомного компьютерного моделирования методов является метод сопряженных градиентов.
Рис. 4. Схематический график функции многих переменных с локальными и глобальным минимумами Наиболее трудной проблемой при минимизации энергии является нахождение глобального минимума, то есть состояния с действительно наименьшей энергией, которая возможна среди всех конфигураций. Сложные многомерные функции имеют множество локальных минимумов (см. рис. 4), и тот из них, в котором функция обладает наименьшим значением, называется глобальным. Если исходная конфигурация оказалась вблизи одного из локальных минимумов, не совпадающих с глобальным, то релаксация любым из методов приводит к этому минимуму. Для атомной системы такое состояние является метастабильным, из которого при нулевой температуре она никогда не выйдет. Таким образом, стабильная структура с минимальной энергией оказывается доступной не из любого начального состояния. Найдя то или иное метастабильное состояние, никогда нельзя определить, является ли оно стабильным, если не перебраны все возможные локальные минимумы. Более того, нет способов определения, сколько локальных минимумов существует для данной системы. Для решения этой проблемы не существует общих методов. Частично преодолеть ее помогает использование метода молекулярной динамики (МД). Ниже, при изучении методов моделирования границ зерен, мы ознакомимся с несколькими методами, которые позволяют преодолевать проблему нахождения глобального минимума именно для этих дефектов.
  Что будет с Землей, если ось ее сместится на 6666 км? Что будет с Землей? - задался я вопросом...  Что делает отдел по эксплуатации и сопровождению ИС? Отвечает за сохранность данных (расписания копирования, копирование и пр.)...  ЧТО И КАК ПИСАЛИ О МОДЕ В ЖУРНАЛАХ НАЧАЛА XX ВЕКА Первый номер журнала «Аполлон» за 1909 г. начинался, по сути, с программного заявления редакции журнала...  ЧТО ПРОИСХОДИТ, КОГДА МЫ ССОРИМСЯ Не понимая различий, существующих между мужчинами и женщинами, очень легко довести дело до ссоры... Не нашли то, что искали? Воспользуйтесь поиском гугл на сайте:
|
, (1.37)
. (1.38)
. Трехчастичный член потенциала имеет вид
, (1.39)
‑ угол с вершиной на атоме i и h – следующая функция с тремя параметрами 
:
. (1.40)
зависящий от угла член (1.39) проводит к энергетически выгодной кристаллической структуре алмаза с тетраэдрическими углами для Si.
, которая представляет собой сложную функцию многих переменных. Поэтому для решения этой задачи используются наиболее эффективные методы и алгоритмы, разработанные в математике для минимизации функций.
) на каждом шаге и составлении некоторой итерационной схемы, приводящей к многомерной точке, в которой этот градиент равен нулю (необходимое условие минимума), а вторые производные положительны (достаточное условие минимума). Алгоритмы различаются именно реализацией этой итерационной схемы.