Потенциал Леннарда-Джонса — это простая модель, которая все еще может описать основные черты взаимодействия между простыми атомами и молекулами: две взаимодействующие частицы отталкиваются друг от друга на очень близком расстоянии, притягиваются друг к другу на умеренном расстоянии и в конечном итоге прекращают взаимодействовать на бесконечном расстоянии, как показано на рисунке. Потенциал Леннарда-Джонса — это парный потенциал, т. е. никакие трех- или многочастичные взаимодействия не покрываются потенциалом. [3] [4]
Общий потенциал Леннарда-Джонса объединяет отталкивающий потенциал, , с притягивающим потенциалом, , используя эмпирически определенные коэффициенты и : [5] [6]
В своем обзоре 1931 года [5] Леннард-Джонс предложил использовать для согласования силу дисперсии Лондона и основанные на согласовании экспериментальные данные. [1] Установка и дает широко используемый потенциал Леннарда-Джонса 12-6: [7]
где r - расстояние между двумя взаимодействующими частицами, ε - глубина потенциальной ямы , а σ - расстояние, на котором потенциальная энергия частица-частица V равна нулю. Потенциал Леннарда-Джонса 12-6 имеет свой минимум на расстоянии, где потенциальная энергия имеет значение
Потенциал Леннарда-Джонса обычно является стандартным выбором для разработки теорий материи (особенно мягкой материи), а также для разработки и тестирования вычислительных методов и алгоритмов.
В прошлом было предложено множество межмолекулярных потенциалов для моделирования простых мягких отталкивающих и притягивающих взаимодействий между сферически симметричными частицами, т. е. общей формы, показанной на рисунке. Примерами других потенциалов являются потенциал Морзе , потенциал Ми , [8] потенциал Букингема и потенциал Танга-Тённиса. [9] Хотя некоторые из них могут больше подходить для моделирования реальных жидкостей , [10] простота потенциала Леннарда-Джонса, а также его часто удивительная способность точно фиксировать поведение реальной жидкости, исторически сделали его парным потенциалом наибольшего общего значения. [11]
История
В 1924 году, в год, когда Леннард-Джонс получил докторскую степень в Кембриджском университете , он опубликовал [6] [12] ряд знаменательных статей о парных потенциалах, которые в конечном итоге были названы в его честь. [2] [3] [13] [1] В этих статьях он скорректировал параметры потенциала, а затем использовал результат в модели вязкости газа, ища набор значений, согласующихся с экспериментом. Его первоначальные результаты предполагали отталкивающий и притягивающий .
До Леннарда-Джонса, еще в 1903 году, Густав Ми работал над эффективными теориями поля; Эдуард Грюнайзен развил работу Ми для твердых тел, показав, что и требуется для твердых тел. В результате этой работы потенциал Леннарда-Джонса иногда называют потенциалом Ми-Грюнайзена в физике твердого тела . [3]
В 1930 году, после открытия квантовой механики , Фриц Лондон показал, что теория предсказывает, что сила притяжения на больших расстояниях должна иметь . В 1931 году Леннард-Джонс применил эту форму потенциала для описания многих свойств жидкостей, заложив основу для многих последующих исследований. [1]
Безразмерные (приведенные единицы)
Безразмерные приведенные единицы могут быть определены на основе параметров потенциала Леннарда-Джонса, что удобно для молекулярного моделирования. С числовой точки зрения преимущества этой системы единиц включают вычисление значений, которые ближе к единице, использование упрощенных уравнений и возможность легкого масштабирования результатов. [14] [15] Эта система приведенных единиц требует указания параметра размера и параметра энергии потенциала Леннарда-Джонса и массы частицы . Все физические свойства могут быть преобразованы напрямую с учетом соответствующей размерности, см. таблицу. Приведенные единицы часто сокращаются и обозначаются звездочкой.
В общем случае редуцированные единицы могут быть построены и на других потенциалах молекулярного взаимодействия, которые состоят из параметра длины и параметра энергии.
Дальнодействующие взаимодействия
Потенциал Леннарда-Джонса, см. уравнение (1) и рисунок сверху, имеет бесконечный диапазон. Только при его рассмотрении исследуется «истинный» и «полный» потенциал Леннарда-Джонса. Для оценки наблюдаемой ансамбля частиц, взаимодействующих посредством потенциала Леннарда-Джонса с использованием молекулярного моделирования, взаимодействия могут быть оценены явно только до определенного расстояния — просто из-за того, что число частиц всегда будет конечным. Максимальное расстояние, применяемое в моделировании, обычно называют радиусом «обрезки» (потому что потенциал Леннарда-Джонса радиально симметричен). Чтобы получить теплофизические свойства (как макроскопические, так и микроскопические) «истинного» и «полного» потенциала Леннарда-Джонса (ЛД), необходимо учесть вклад потенциала за пределами радиуса обрезки.
Были разработаны различные схемы коррекции для учета влияния дальнодействующих взаимодействий в моделировании и для поддержания достаточно хорошего приближения «полного» потенциала. [16] [14] Они основаны на упрощающих предположениях относительно структуры жидкости. Для простых случаев, таких как исследования равновесия однородных жидкостей, простые поправочные члены дают превосходные результаты. В других случаях, таких как исследования неоднородных систем с различными фазами, учет дальнодействующих взаимодействий более утомителен. Эти поправки обычно называют «дальнодействующими поправками». Для большинства свойств известны и хорошо известны простые аналитические выражения. Для заданной наблюдаемой «скорректированный» результат моделирования затем просто вычисляется из фактически выбранного значения и значения дальнодействующей поправки , например, для внутренней энергии . [14] Гипотетическое истинное значение наблюдаемой потенциала Леннарда-Джонса на действительно бесконечном расстоянии отсечки (термодинамический предел) в общем случае может быть только оценено.
Кроме того, качество схемы коррекции на больших расстояниях зависит от радиуса отсечки. Предположения, сделанные в схемах коррекции, обычно не оправдываются при (очень) коротких радиусах отсечки. Это проиллюстрировано в примере, показанном на рисунке справа. Говорят, что схема коррекции на больших расстояниях сходится, если остаточная ошибка схемы коррекции достаточно мала на заданном расстоянии отсечки, см. рисунок.
Расширения и модификации
Потенциал Леннарда-Джонса — как архетип межмолекулярных потенциалов — многократно использовался в качестве отправной точки для разработки более сложных или более обобщенных межмолекулярных потенциалов. Различные расширения и модификации потенциала Леннарда-Джонса были предложены в литературе; более обширный список приведен в статье « межатомный потенциал ». Следующий список относится только к нескольким примерам потенциалов, которые напрямую связаны с потенциалом Леннарда-Джонса и имеют как историческое значение, так и по-прежнему актуальны для современных исследований.
Потенциал Ми Потенциал Ми является обобщенной версией потенциала Леннарда-Джонса, т. е. показатели 12 и 6 вводятся как параметрыи. Известно, что особенно термодинамические производные свойства, например, сжимаемость и скорость звука , очень чувствительны к крутизне отталкивающей части межмолекулярного потенциала, который, следовательно, может быть более сложно смоделирован потенциалом Ми. [17] Первая явная формулировка потенциала Ми приписывается Эдуарду Грюнайзену . [18] [19] Следовательно, потенциал Ми был фактически предложен до потенциала Леннарда-Джонса. Потенциал Ми назван в честь Густава Ми . [8]
Потенциал Букингема Потенциал Букингема был предложен Ричардом Букингемом . Отталкивающая часть потенциала Леннарда-Джонса в нем заменена экспоненциальной функцией и включает в себя дополнительный параметр.
Потенциал Штокмайера Потенциал Штокмайера назван в честь WH Stockmayer. [20] Потенциал Штокмайера представляет собой комбинацию потенциала Леннарда-Джонса, наложенного диполем. Следовательно, частицы Штокмайера не являются сферически симметричными, а скорее имеют важную ориентационную структуру.
Двухцентровой потенциал Леннарда-Джонса Двухцентровой потенциал Леннарда-Джонса состоит из двух идентичных сайтов взаимодействия Леннарда-Джонса (одинаковых , , ), которые связаны как твердое тело. Его часто сокращают до 2CLJ. Обычно удлинение (расстояние между сайтами Леннарда-Джонса) значительно меньше параметра размера . Следовательно, два сайта взаимодействия значительно слиты.
Усеченный и сплайновый потенциал Леннарда-Джонса Усеченный и сплайновый потенциал Леннарда-Джонса — это редко используемый, но полезный потенциал. Подобно более популярному потенциалу LJTS, он прочно усечен на определенном «конечном» расстоянии , и за его пределами не рассматриваются никакие дальнодействующие взаимодействия. В отличие от потенциала LJTS, который смещен таким образом, что потенциал непрерывен, усеченный и сплайновый потенциал Леннарда-Джонса становится непрерывным с помощью произвольной, но благоприятной сплайновой функции.
Усеченный и смещенный потенциал Леннарда-Джонса (LJTS)
Усеченный и смещенный потенциал Леннарда-Джонса (LJTS) часто используется в качестве альтернативы «полному» потенциалу Леннарда-Джонса (см. уравнение (1)). «Полный» и «усеченный и смещенный» потенциал Леннарда-Джонса должны строго разделяться. Это просто два разных межмолекулярных потенциала, дающих разные теплофизические свойства. Усеченный и смещенный потенциал Леннарда-Джонса определяется как
с
Следовательно, потенциал LJTS усекается при и смещается на соответствующее значение энергии . Последнее применяется, чтобы избежать скачка разрыва потенциала при . Для потенциала LJTS не требуются дальнодействующие взаимодействия за пределами — ни явно, ни неявно. Наиболее часто используемая версия усеченного и смещенного потенциала Леннарда-Джонса — это та, в которой . [ необходима цитата ] Тем не менее, в литературе использовались разные значения. [24] [25] [26] [27] Каждый потенциал LJTS с заданным радиусом усечения следует рассматривать как потенциал и, соответственно, как отдельную субстанцию.
Потенциал LJTS в вычислительном отношении значительно дешевле, чем «полный» потенциал Леннарда-Джонса, но все еще охватывает основные физические характеристики материи (наличие критической и тройной точек, мягкие отталкивающие и притягивающие взаимодействия, фазовые равновесия и т. д.). Поэтому потенциал LJTS используется для тестирования новых алгоритмов, методов моделирования и новых физических теорий. [28] [29] [30]
Интересно, что для гомогенных систем межмолекулярные силы, которые рассчитываются из LJ и LJTS потенциала на заданном расстоянии, одинаковы (поскольку то же самое), тогда как потенциальная энергия и давление зависят от смещения. Кроме того, свойства вещества LJTS могут быть дополнительно затронуты выбранным алгоритмом моделирования, т. е. выборкой MD или MC (это, как правило, не относится к «полному» потенциалу Леннарда-Джонса).
Для потенциала LJTS с , сдвиг потенциальной энергии составляет приблизительно 1/60 от энергии дисперсии в потенциальной яме: . На рисунке справа показано сравнение равновесия пар-жидкость «полного» потенциала Леннарда-Джонса и «усеченного и смещенного» потенциала Леннарда-Джонса. Результаты «полного» потенциала Леннарда-Джонса преобладают при значительно более высокой критической температуре и давлении по сравнению с результатами потенциала LJTS, но критическая плотность очень похожа. [31] [32] [26] Давление пара и энтальпия испарения сильнее зависят от дальнодействующих взаимодействий, чем от насыщенных плотностей. Это связано с тем, что потенциал в основном энергетически изменяется путем усечения и смещения.
Приложения
Потенциал Леннарда-Джонса имеет не только фундаментальное значение в вычислительной химии и физике мягкой материи , но и для моделирования реальных веществ. Потенциал Леннарда-Джонса используется для фундаментальных исследований поведения материи и для объяснения атомистических явлений. Он также часто используется для несколько специальных случаев использования, например, для изучения теплофизических свойств двумерных или четырехмерных веществ [33] [34] [35] (вместо классических трех пространственных направлений нашей Вселенной).
Существует два основных применения потенциалов Леннарда-Джонса: (i) для изучения гипотетической субстанции Леннарда-Джонса [13] и (ii) для моделирования взаимодействий в моделях реальных субстанций. [3] [2] Эти два применения обсуждаются далее.
вещество Леннарда-Джонса
Вещество Леннарда-Джонса или «Леннард-Джонсий» — это название, данное идеализированному веществу, которое было бы результатом взаимодействия атомов или молекул исключительно через потенциал Леннарда-Джонса. [13] Статистическая механика [36] и компьютерное моделирование [15] [16] могут быть использованы для изучения потенциала Леннарда-Джонса и получения теплофизических свойств «вещества Леннарда-Джонса». Вещество Леннарда-Джонса часто называют «Леннард-Джонсием», [13] предполагая, что оно рассматривается как (фиктивный) химический элемент . [21] Более того, его параметры энергии и длины можно настроить так, чтобы они соответствовали многим различным реальным веществам. Как потенциал Леннарда-Джонса, так и, соответственно, вещество Леннарда-Джонса являются упрощенными, но реалистичными моделями, например, они точно отражают основные физические принципы, такие как наличие критической и тройной точек , конденсация и замерзание . Отчасти из-за своей математической простоты потенциал Леннарда-Джонса широко использовался в исследованиях материи с первых дней компьютерного моделирования. [37] [38] [39] [40]
Теплофизические свойства вещества Леннарда-Джонса
Теплофизические свойства вещества Леннарда-Джонса [13] , т. е. частицы, взаимодействующие с потенциалом Леннарда-Джонса, могут быть получены с помощью статистической механики. Некоторые свойства могут быть вычислены аналитически, т. е. с машинной точностью, тогда как большинство свойств могут быть получены только путем выполнения молекулярного моделирования. [15] Последнее, как правило, будет наложено как статистическими, так и систематическими неопределенностями. [43] [21] [44] [45] Например, вириальные коэффициенты могут быть вычислены непосредственно из потенциала Леннарда с использованием алгебраических выражений [36] , и поэтому представленные данные не имеют неопределенности. Результаты молекулярного моделирования, например, давление при заданной температуре и плотности, имеют как статистические, так и систематические неопределенности. [43] [45] Молекулярное моделирование потенциала Леннарда-Джонса, как правило, может быть выполнено с использованием либо моделирования молекулярной динамики (МД), либо моделирования Монте-Карло (МК). Для моделирования МК напрямую используется потенциал Леннарда-Джонса , тогда как моделирование МД всегда основано на производной потенциала, т.е. силе . Эти различия в сочетании с различиями в обработке дальнодействующих взаимодействий (см. ниже) могут влиять на вычисляемые теплофизические свойства. [46] [32]
Поскольку Леннард-Джонсиум является архетипом для моделирования простых, но реалистичных межмолекулярных взаимодействий, большое количество теплофизических свойств было изучено и описано в литературе. [21] Данные компьютерных экспериментов по потенциалу Леннарда-Джонса в настоящее время считаются наиболее точными известными данными в классической механике вычислительной химии. Следовательно, такие данные также в основном используются в качестве эталона для проверки и тестирования новых алгоритмов и теорий. Потенциал Леннарда-Джонса постоянно использовался с первых дней молекулярного моделирования. Первые результаты компьютерных экспериментов для потенциала Леннарда-Джонса были представлены Розенблютом и Розенблютом [38] и Вудом и Паркером [37] после того, как молекулярное моделирование на « быстрых вычислительных машинах » стало доступно в 1953 году. [47] С тех пор во многих исследованиях сообщалось о данных вещества Леннарда-Джонса; [21] приблизительно 50 000 точек данных находятся в открытом доступе. Текущее состояние исследований теплофизических свойств вещества Леннарда-Джонса обобщено Стефаном и др. [21] (которые не охватывают свойства переноса и смешивания). Национальный институт стандартов и технологий США (NIST) приводит примеры молекулярной динамики и кодов Монте-Карло вместе с результатами, полученными с их помощью. [48] Данные о свойствах переноса жидкостей Леннарда-Джонса были собраны Беллом и др. [49] и Лаутеншегером и Хассе. [50]
На рисунке справа показана фазовая диаграмма жидкости Леннарда-Джонса. Фазовые равновесия потенциала Леннарда-Джонса изучались много раз и, соответственно, известны сегодня с хорошей точностью. [41] [21] [51] На рисунке показаны результаты корреляций, полученные из результатов компьютерного эксперимента (поэтому показаны линии вместо точек данных).
Среднее межмолекулярное взаимодействие частицы Леннарда-Джонса сильно зависит от термодинамического состояния, т. е. температуры и давления (или плотности). Для твердых состояний притягивающее взаимодействие Леннарда-Джонса играет доминирующую роль, особенно при низких температурах. Для жидких состояний упорядоченная структура отсутствует по сравнению с твердыми состояниями. Средняя потенциальная энергия на частицу отрицательна. Для газообразных состояний притягивающие взаимодействия потенциала Леннарда-Джонса играют второстепенную роль, поскольку они находятся на большом расстоянии. Основная часть внутренней энергии запасается в виде кинетической энергии для газообразных состояний. В сверхкритических состояниях притягивающее взаимодействие Леннарда-Джонса играет второстепенную роль. С ростом температуры средняя кинетическая энергия частиц увеличивается и превышает энергетическую яму потенциала Леннарда-Джонса. Следовательно, частицы в основном взаимодействуют посредством мягких отталкивающих взаимодействий потенциалов, и средняя потенциальная энергия на частицу соответственно положительна.
В целом, из-за большого временного интервала, в течение которого изучался потенциал Леннарда-Джонса, и в литературе были представлены данные о теплофизических свойствах, а вычислительные ресурсы были недостаточны для точного моделирования (по современным стандартам), известно, что заметное количество данных является сомнительным. [21] Тем не менее, во многих исследованиях такие данные используются в качестве справочных. Отсутствие хранилищ данных и оценки данных является важнейшим элементом для будущей работы в долговременной области исследования потенциала Леннарда-Джонса.
Характерные точки и кривые
Наиболее важными характерными точками потенциала Леннарда-Джонса являются критическая точка и тройная точка пар-жидкость-твердое тело . Они были изучены много раз в литературе и собраны в [21] . Таким образом, критическая точка была оценена как расположенная в
Указанные неопределенности были рассчитаны на основе стандартного отклонения критических параметров, полученных из наиболее надежных доступных наборов данных о равновесии пар-жидкость . [21] Эти неопределенности можно принять как нижний предел точности, с которой критическая точка жидкости может быть получена из результатов молекулярного моделирования.
В настоящее время предполагается, что тройная точка находится в
Неопределенности представляют собой разброс данных разных авторов. [41] Критическая точка вещества Леннарда-Джонса изучалась гораздо чаще, чем тройная точка. Как для критической точки, так и для тройной точки пар-жидкость-твердое тело в нескольких исследованиях сообщалось о результатах, выходящих за пределы указанных выше диапазонов. Указанные выше данные в настоящее время считаются правильными и надежными. Тем не менее, определенность критической температуры и температуры тройной точки все еще неудовлетворительна.
Очевидно, что кривые сосуществования фаз (ср. рисунки) имеют фундаментальное значение для характеристики потенциала Леннарда-Джонса. Кроме того, характеристические кривые Брауна [55] дают наглядное описание существенных особенностей потенциала Леннарда-Джонса. Характеристические кривые Брауна определяются как кривые, на которых определенное термодинамическое свойство вещества соответствует свойству идеального газа . Для реальной жидкости и ее производных могут соответствовать значениям идеального газа для специальных комбинаций только в результате правила фаз Гиббса. Полученные точки в совокупности составляют характеристическую кривую. Определены четыре основные характеристические кривые: одна 0-го порядка (названная кривой Зенона ) и три кривые 1-го порядка (названные кривыми Амагата , Бойля и Чарльза ). Характеристическая кривая должна иметь отрицательную или нулевую кривизну на всем протяжении и один максимум на двойной логарифмической диаграмме давление-температура. Более того, характеристические кривые Брауна и вириальные коэффициенты напрямую связаны в пределе идеального газа и поэтому точно известны при . Результаты компьютерного моделирования и уравнения состояния были опубликованы в литературе для потенциала Леннарда-Джонса. [53] [21] [52] [56] [57]
Точки на кривой Зенона Z имеют коэффициент сжимаемости , равный единице . Кривая Зенона начинается при температуре Бойля , окружает критическую точку и имеет наклон, равный единице в пределе низких температур. [52] Точки на кривой Бойля B имеют . Кривая Бойля начинается с кривой Зенона при температуре Бойля, слабо окружает критическую точку и заканчивается на кривой давления пара. Точки на кривой Чарльза (также известной как кривая инверсии Джоуля-Томсона ) имеют и, что более важно , , т. е. никакого изменения температуры при изоэнтальпийном дросселировании. Она начинается при в пределе идеального газа, пересекает кривую Зенона и заканчивается на кривой давления пара. Точки на кривой Амагата A имеют . Она также начинается в пределе идеального газа при , окружает критическую точку и три другие характеристические кривые и переходит в область твердой фазы. Подробное обсуждение характеристических кривых потенциала Леннарда-Джонса дано Стефаном и Дейтерсом. [52]
Свойства жидкости Леннарда-Джонса
Свойства жидкости Леннарда-Джонса были подробно изучены в литературе из-за выдающейся важности потенциала Леннарда-Джонса в физике мягкой материи и смежных областях. [13] На сегодняшний день опубликовано около 50 наборов данных компьютерных экспериментов для равновесия пар-жидкость . [21] Кроме того, за эти годы было опубликовано более 35 000 точек данных в однородных состояниях жидкости, которые недавно были скомпилированы и оценены на предмет выбросов в базе данных с открытым доступом. [21]
Равновесие пар-жидкость вещества Леннарда-Джонса в настоящее время известно с точностью, т.е. взаимное согласие термодинамически согласованных данных для давления пара, плотности насыщенной жидкости, плотности насыщенного пара, энтальпии испарения и поверхностного натяжения. [21] Этот статус-кво нельзя считать удовлетворительным, учитывая тот факт, что статистические неопределенности, обычно сообщаемые для отдельных наборов данных, значительно ниже вышеуказанных значений (даже для гораздо более сложных молекулярных силовых полей).
Как свойства фазового равновесия, так и свойства однородного состояния при произвольной плотности в общем случае могут быть получены только из молекулярных симуляций, тогда как вириальные коэффициенты могут быть вычислены непосредственно из потенциала Леннарда-Джонса. [36] Численные данные для второго и третьего вириальных коэффициентов доступны в широком диапазоне температур. [58] [52] [21] Для более высоких вириальных коэффициентов (до шестнадцатого) количество доступных точек данных уменьшается с увеличением номера вириального коэффициента. [59] [60] Также были изучены транспортные свойства (вязкость, теплопроводность и коэффициент самодиффузии) жидкости Леннарда-Джонса, [61] [62] но база данных значительно менее плотная, чем для однородных равновесных свойств, таких как – или данные о внутренней энергии. Более того, было разработано большое количество аналитических моделей ( уравнений состояния ) для описания жидкости Леннарда-Джонса (подробности см. ниже).
Свойства тела Леннарда-Джонса
База данных и знания для твердого тела Леннарда-Джонса значительно хуже, чем для жидких фаз. Было рано осознано, что взаимодействия в твердых фазах не должны быть аппроксимированы как попарно аддитивные – особенно для металлов. [63] [64]
Тем не менее, потенциал Леннарда-Джонса используется в физике твердого тела из-за его простоты и вычислительной эффективности. Следовательно, основные свойства твердых фаз и равновесия твердое тело-жидкость были исследованы несколько раз, например, [51] [41] [42] [65] [66] [54]
Вещество Леннарда-Джонса образует fcc (гранецентрированную кубическую), hcp (гексагональную плотноупакованную) и другие плотноупакованные политипные решетки — в зависимости от температуры и давления, см. рисунок выше с фазовой диаграммой. При низкой температуре и вплоть до умеренного давления энергетически выгодна решетка hcp и, следовательно, равновесная структура. Структура решетки fcc энергетически выгодна как при высокой температуре, так и при высоком давлении и, следовательно, в целом равновесная структура в более широком диапазоне состояний. Линия сосуществования между фазами fcc и hcp начинается примерно при , проходит через температурный максимум примерно при , а затем заканчивается на границе фаз пар–твердое тело примерно при , что тем самым образует тройную точку. [65] [41] Следовательно, только твердая фаза fcc демонстрирует фазовые равновесия с жидкой и сверхкритической фазой, см. рисунок выше с фазовой диаграммой.
Сообщается, что тройная точка двух твердых фаз (ГЦК и ГПУ) и паровой фазы находится в: [65] [41]
пока не сообщается
Обратите внимание, что в литературе также сообщалось о других и существенно отличающихся значениях. Следовательно, база данных для тройной точки fcc-hcp–vapor должна быть дополнительно укреплена в будущем.
Смеси веществ Леннарда-Джонса
Смеси частиц Леннарда-Джонса в основном используются в качестве прототипа для разработки теорий и методов растворов, а также для изучения свойств растворов в целом. Это восходит к фундаментальной работе теории конформных растворов Лонге-Хиггинса [67] и Леланда и Роулинсона и их коллег. [68] [69] Сегодня они являются основой большинства теорий смесей. [70] [71]
Смеси двух или более компонентов Леннарда-Джонса устанавливаются путем изменения по крайней мере одного параметра потенциального взаимодействия ( или ) одного из компонентов по отношению к другому. Для бинарной смеси это дает три типа парных взаимодействий, которые все моделируются потенциалом Леннарда-Джонса: взаимодействия 1-1, 2-2 и 1-2. Для перекрестных взаимодействий 1–2 требуются дополнительные предположения для спецификации параметров или из , и , . Различные варианты (все более или менее эмпирические и не основанные строго на физических аргументах) могут быть использованы для этих так называемых правил комбинирования. [72] Наиболее широко используемым [72] правилом комбинирования является правило Лоренца и Бертло [73]
Параметр является дополнительным параметром взаимодействия, независимым от состояния, для смеси. Параметр обычно устанавливается равным единице, поскольку среднее арифметическое можно считать физически правдоподобным для параметра размера перекрестного взаимодействия. С другой стороны, параметр часто используется для корректировки среднего геометрического , чтобы воспроизвести фазовое поведение модельной смеси. Для аналитических моделей, например, уравнений состояния , параметр отклонения обычно записывается как . Для энергия дисперсии перекрестного взаимодействия и, соответственно, сила притяжения между разнородными частицами усиливаются, а силы притяжения между разнородными частицами уменьшаются для .
Для смесей Леннарда-Джонса можно изучать как равновесия жидкой, так и твердой фаз , то есть пар–жидкость , жидкость–жидкость , газ–газ, твердое тело–пар, твердое тело–жидкость и твердое тело–твердое тело. Соответственно, могут существовать различные типы тройных точек (трехфазные равновесия) и критических точек , а также различные эвтектические и азеотропные точки . [74] [71] Бинарные смеси Леннарда-Джонса в области жидкости (различные типы равновесий жидких и газовых фаз) [31] [75 ] [76] [77] [78] были изучены более всесторонне, чем фазовые равновесия, включающие твердые фазы. [79] [80] [81] [82] [83] В литературе было изучено большое количество различных смесей Леннарда-Джонса. На сегодняшний день не установлен стандарт для них. Обычно параметры бинарного взаимодействия и параметры двух компонентов выбираются таким образом, чтобы получить смесь со свойствами, подходящими для данной задачи. Однако это часто затрудняет сравнение.
Для поведения жидкой фазы смеси демонстрируют практически идеальное поведение (в смысле закона Рауля ) для . Для преобладают притягивающие взаимодействия, и смеси имеют тенденцию образовывать высококипящие азеотропы, т. е. для стабилизации равновесия пар-жидкость требуется более низкое давление, чем давление паров чистых компонентов. Для преобладают отталкивающие взаимодействия, и смеси имеют тенденцию образовывать низкокипящие азеотропы, т. е. для стабилизации равновесия пар-жидкость требуется более высокое давление, чем давление паров чистых компонентов, поскольку средние дисперсионные силы уменьшаются. Особенно низкие значения, кроме того, приведут к зазорам смешиваемости жидкость-жидкость. Также в литературе изучались различные типы фазовых равновесий, включающих твердые фазы, например, Кэрол и ее коллегами. [81] [83] [80] [79] Также существуют случаи, когда границы твердой фазы прерывают равновесия жидкой фазы. Однако для фазовых равновесий, включающих твердые фазы, объем опубликованных данных невелик.
Уравнения состояния
Большое количество уравнений состояния (УС) для потенциала/вещества Леннарда-Джонса было предложено с тех пор, как его характеристика и оценка стали доступны с первыми компьютерными симуляциями. [47] В связи с фундаментальной важностью потенциала Леннарда-Джонса, большинство имеющихся в настоящее время молекулярных УС построены вокруг жидкости Леннарда-Джонса. Они были всесторонне рассмотрены Стефаном и др. [11] [52]
Уравнения состояния для жидкости Леннарда-Джонса имеют особое значение в физике мягких веществ и физической химии , используются в качестве отправной точки для разработки УС для сложных жидкостей, например, полимеров и ассоциированных жидкостей. Мономерные единицы этих моделей обычно напрямую адаптируются из УС Леннарда-Джонса в качестве строительного блока, например, УС PHC, [84] УС BACKONE, [85] [86] и УС типа SAFT . [17] [87] [88] [89]
В литературе было предложено более 30 уравнений состояния Леннарда-Джонса. Всесторонняя оценка [11] [52] таких уравнений состояния показала, что несколько уравнений состояния [90] [91] [92] [93] описывают потенциал Леннарда-Джонса с хорошей и схожей точностью, но ни одно из них не является выдающимся. Три из этих уравнений состояния показывают неприемлемое нефизическое поведение в некоторой области жидкости, например, множественные петли Ван-дер-Ваальса, будучи в остальном достаточно точными. Только уравнение состояния Леннарда-Джонса Колафы и Незбеды [91] оказалось надежным и точным для большинства термодинамических свойств жидкости Леннарда-Джонса. [52] [11] Кроме того, уравнение состояния Леннарда-Джонса Джонсона и др. [94] оказалось менее точным для практически всех доступных справочных данных [21] [11], чем уравнения состояния Колафы и Незбеды. [91]
Потенциал Леннарда-Джонса как строительный блок для силовых полей
Потенциал Леннарда-Джонса широко используется для молекулярного моделирования реальных веществ. По сути, существует два способа использования потенциала Леннарда-Джонса для молекулярного моделирования: (1) Реальный атом или молекула вещества моделируются непосредственно потенциалом Леннарда-Джонса, что дает очень хорошие результаты для благородных газов и метана , т. е. дисперсионно взаимодействующих сферических частиц. В случае метана предполагается, что молекула сферически симметрична, а атомы водорода объединены с атомом углерода в общую единицу. Это упрощение в общем случае может быть применено и к более сложным молекулам, но обычно дает плохие результаты. (2) Реальная молекула вещества состоит из нескольких участков взаимодействия Леннарда-Джонса, которые могут быть связаны либо жесткими связями, либо гибкими дополнительными потенциалами (и в конечном итоге также состоит из других типов потенциалов, например, частичных зарядов). Молекулярные модели (часто называемые « силовыми полями ») практически для всех молекулярных и ионных частиц могут быть построены с использованием этой схемы, например, для алканов .
При использовании первого изложенного подхода молекулярная модель имеет только два параметра потенциала Леннарда-Джонса , которые можно использовать для подгонки, например , и можно использовать для аргона . После подгонки параметров модели ε и σ к реальным свойствам вещества потенциал Леннарда-Джонса можно использовать для описания простого вещества (например, благородных газов ) с хорошей точностью. Очевидно, что этот подход является лишь хорошим приближением для сферических и просто дисперсионно взаимодействующих молекул и атомов. Прямое использование потенциала Леннарда-Джонса имеет большое преимущество в том, что результаты моделирования и теории для потенциала Леннарда-Джонса можно использовать напрямую. Следовательно, доступные результаты для потенциала Леннарда-Джонса и вещества можно напрямую масштабировать с помощью соответствующих и (см. приведенные единицы). Параметры потенциала Леннарда-Джонса и в целом можно подогнать под любое желаемое свойство реального вещества. В физике мягкой материи обычно для параметризации используются экспериментальные данные для фазового равновесия пар-жидкость или критической точки; в физике твердого тела используются скорее сжимаемость, теплоемкость или постоянные решетки. [63] [64]
Второй изложенный подход использования потенциала Леннарда-Джонса в качестве строительного блока удлиненных и сложных молекул гораздо более сложен. Молекулярные модели , таким образом, изготавливаются на заказ в том смысле, что результаты моделирования применимы только для этой конкретной модели. Этот подход к разработке молекулярных силовых полей сегодня в основном применяется в физике мягкой материи и связанных с ней областях, таких как химическая инженерия , химия и вычислительная биология. Большое количество силовых полей основано на потенциале Леннарда-Джонса, например, силовое поле TraPPE [95] , силовое поле OPLS [96] и силовое поле MolMod [97] (обзор молекулярных силовых полей выходит за рамки настоящей статьи). Для современного моделирования твердотельных материалов используются более сложные многочастичные потенциалы (например, потенциалы EAM [98] ).
Потенциал Леннарда-Джонса дает хорошее приближение межмолекулярных взаимодействий для многих приложений: макроскопические свойства, вычисленные с использованием потенциала Леннарда-Джонса, хорошо согласуются с экспериментальными данными для простых веществ, таких как аргон, с одной стороны, а потенциальная функция неплохо согласуется с результатами квантовой химии , с другой стороны. Потенциал Леннарда-Джонса дает хорошее описание молекулярных взаимодействий в жидких фазах, тогда как молекулярные взаимодействия в твердых фазах описываются лишь приблизительно. Это в основном связано с тем, что многочастичные взаимодействия играют значительную роль в твердых фазах, которые не включены в потенциал Леннарда-Джонса. Поэтому потенциал Леннарда-Джонса широко используется в физике мягкой материи и связанных с ней областях, тогда как он реже используется в физике твердого тела . Благодаря своей простоте потенциал Леннарда-Джонса часто используется для описания свойств газов и простых жидкостей, а также для моделирования дисперсионных и отталкивательных взаимодействий в молекулярных моделях . Он особенно точен для атомов благородных газов и метана . Кроме того, он является хорошим приближением для молекулярных взаимодействий на больших и малых расстояниях для нейтральных атомов и молекул. Поэтому потенциал Леннарда-Джонса очень часто используется в качестве строительного блока молекулярных моделей сложных молекул, например, алканов или воды . [95] [99] [97] Потенциал Леннарда-Джонса также может использоваться для моделирования адсорбционных взаимодействий на границах раздела твердое тело-жидкость, т. е. физосорбции или хемосорбции .
Общепризнано, что основные ограничения потенциала Леннарда-Джонса заключаются в том, что потенциал является парным потенциалом (не охватывает многочастичные взаимодействия) и что показатель степени используется для отталкивания. Результаты квантовой химии показывают, что необходимо использовать показатель степени выше 12, т. е. более крутой потенциал. Кроме того, потенциал Леннарда-Джонса имеет ограниченную гибкость, т. е. только два параметра модели и может использоваться для подгонки для описания реального вещества.
Смотрите также
На Викискладе есть медиафайлы по теме « Потенциалы Леннарда-Джонса» .
^ abcd Фишер, Иоганн; Вендланд, Мартин (2023-10-01). «Об истории ключевых эмпирических межмолекулярных потенциалов». Fluid Phase Equilibria . 573 : 113876. Bibcode : 2023FlPEq.57313876F. doi : 10.1016/j.fluid.2023.113876 . ISSN 0378-3812.
^ abc Ленхард, Йоханнес; Стефан, Саймон; Хассе, Ганс (2024-05-16). «Об истории потенциала Леннарда-Джонса». Annalen der Physik . 536 (6). doi :10.1002/andp.202400115. ISSN 0003-3804.
^ abcde Швердтфегер, Питер; Уэльс, Дэвид Дж. (2024-04-26). «100 лет потенциалу Леннарда-Джонса». Журнал химической теории и вычислений . 20 (9): 3379–3405. doi :10.1021/acs.jctc.4c00135. ISSN 1549-9618. PMID 38669689.
^ Стоун, А. Дж. (2016). Теория межмолекулярных сил (Второе издание, опубликованное в мягкой обложке). Оксфорд: Oxford University Press. ISBN978-0-19-878915-4.
^ ab Lennard-Jones, JE (1931-09-01). «Сплоченность». Труды Физического общества . 43 (5): 461–482. Bibcode :1931PPS....43..461L. doi :10.1088/0959-5309/43/5/301. ISSN 0959-5309.
^ ab Jones, JE (1924). «Об определении молекулярных полей.—I. Из изменения вязкости газа с температурой». Труды Лондонского королевского общества. Серия A, содержащая статьи математического и физического характера . 106 (738): 441–462. Bibcode : 1924RSPSA.106..441J. doi : 10.1098/rspa.1924.0081 . ISSN 0950-1207.
^ Ван, Сипэн; Рамирес-Хинестроса, Симон; Добникар, Юре; Френкель, Даан (2020). «Потенциал Леннарда-Джонса: когда (не) его использовать». Физическая химия Химическая физика . 22 (19): 10624–10633. arXiv : 1910.05746 . Bibcode :2020PCCP...2210624W. doi :10.1039/C9CP05445F. ISSN 1463-9076. PMID 31681941.
^ аб Ми, Густав (1903). «Zur kinetischen Theorie der Einatomigen Körper». Аннален дер Физик (на немецком языке). 316 (8): 657–697. Бибкод : 1903АнП...316..657М. дои : 10.1002/andp.19033160802.
^ Tang, KT; Toennies, J. Peter (1984-04-15). «Улучшенная простая модель для потенциала Ван-дер-Ваальса на основе универсальных функций затухания для коэффициентов дисперсии». Журнал химической физики . 80 (8): 3726–3741. Bibcode : 1984JChPh..80.3726T. doi : 10.1063/1.447150. ISSN 0021-9606.
^ Лафитт, Томас; Апостолаку, Анастасия; Авенданьо, Карлос; Галиндо, Ампаро; Аджиман, Клэр С.; Мюллер, Эрих А.; Джексон, Джордж (2013-10-21). "Точная статистическая теория ассоциированной жидкости для цепных молекул, образованных из сегментов Ми". Журнал химической физики . 139 (15). Bibcode : 2013JChPh.139o4504L. doi : 10.1063/1.4819786 . hdl : 10044/1/12859 . ISSN 0021-9606. PMID 24160524.
^ abcdefghi Стефан, Саймон; Штаубах, Йенс; Хассе, Ганс (2020). «Обзор и сравнение уравнений состояния для жидкости Леннарда-Джонса». Fluid Phase Equilibria . 523 : 112772. Bibcode : 2020FlPEq.52312772S. doi : 10.1016/j.fluid.2020.112772. S2CID 224844789.
^ Джонс, Дж. Э. (1924). «Об определении молекулярных полей. —II. Из уравнения состояния газа». Труды Лондонского королевского общества. Серия A, содержащая статьи математического и физического характера . 106 (738): 463–477. Bibcode : 1924RSPSA.106..463J. doi : 10.1098/rspa.1924.0082 . ISSN 0950-1207.
^ abcdef Ленхард, Йоханнес; Стефан, Саймон; Хассе, Ганс (февраль 2024 г.). «Дитя предсказания. Об истории, онтологии и вычислениях Леннарда-Джонсиума». Исследования по истории и философии науки . 103 : 105–113. doi :10.1016/j.shpsa.2023.11.007. PMID 38128443. S2CID 266440296.
^ abc Allen, Michael P.; Tildesley, Dominic J. (2017-11-23). "Компьютерное моделирование жидкостей". Oxford Scholarship Online . doi :10.1093/oso/9780198803195.001.0001. ISBN9780198803195.
^ abc DC Rapaport (1 апреля 2004 г.). Искусство моделирования молекулярной динамики. Cambridge University Press. ISBN978-0-521-82568-9.
^ ab Френкель, Д.; Смит, Б. (2002), Понимание молекулярного моделирования (второе издание), Сан-Диего: Academic Press, ISBN0-12-267351-4
^ ab Лафитт, Томас; Апостолаку, Анастасия; Авенданьо, Карлос; Галиндо, Ампаро; Аджиман, Клэр С.; Мюллер, Эрих А.; Джексон, Джордж (2013-10-16). "Точная статистическая теория ассоциированной жидкости для цепных молекул, образованных из сегментов Ми". Журнал химической физики . 139 (15): 154504. Bibcode : 2013JChPh.139o4504L. doi : 10.1063/1.4819786. hdl : 10044/1/12859 . ISSN 0021-9606. PMID 24160524.
^ Грюнайзен, Эдвард (1911). «Das Verhältnis der thermischen Ausdehnung zur spezifischen Wärme Fester Elemente». Zeitschrift für Elektrochemie und angewandte Physikalische Chemie . 17 (17): 737–739. дои : 10.1002/bbpc.191100004. S2CID 178760389.
^ Грюнайзен, Э. (1912). «Теория праздников: элементы einatomiger Elemente». Аннален дер Физик (на немецком языке). 344 (12): 257–306. Бибкод : 1912АнП...344..257Г. дои : 10.1002/andp.19123441202.
^ Stockmayer, WH (1941-05-01). "Вторые вириальные коэффициенты полярных газов". Журнал химической физики . 9 (5): 398–402. Bibcode :1941JChPh...9..398S. doi :10.1063/1.1750922. ISSN 0021-9606.
^ abcdefghijklmnopqrst Стефан, Саймон; Тол, Моника; Врабец, Ядран; Хассе, Ханс (28.10.2019). «Теплофизические свойства жидкости Леннарда-Джонса: база данных и оценка данных». Журнал химической информации и моделирования . 59 (10): 4248–4265. doi :10.1021/acs.jcim.9b00620. ISSN 1549-9596. PMID 31609113. S2CID 204545481.
^ Vrabec, Jadran; Kedia, Gaurav Kumar; Fuchs, Guido; Hasse, Hans (2006-05-10). «Комплексное исследование сосуществования пара и жидкости усеченной и смещенной жидкости Леннарда–Джонса, включая свойства плоского и сферического интерфейса». Молекулярная физика . 104 (9): 1509–1527. Bibcode : 2006MolPh.104.1509V. doi : 10.1080/00268970600556774. ISSN 0026-8976. S2CID 96606562.
^ Хейер, Микаэла; Стефан, Саймон; Лю, Джинлу; Чепмен, Уолтер Г.; Хассе, Ханс; Лангенбах, Кай (18.08.2018). «Уравнение состояния для усеченной и смещенной по Леннард-Джонсу жидкости с радиусом отсечки 2,5 σ на основе теории возмущений и ее приложений к межфазной термодинамике». Молекулярная физика . 116 (15–16): 2083–2094. Bibcode : 2018MolPh.116.2083H. doi : 10.1080/00268976.2018.1447153. ISSN 0026-8976. S2CID 102956189.
^ Шауль, Кэтрин RS; Шульц, Эндрю Дж.; Кофке, Дэвид А. (2010). «Влияние усечения и сдвига на вириальные коэффициенты потенциалов Леннарда–Джонса». Сборник чехословацких химических сообщений . 75 (4): 447–462. doi :10.1135/cccc2009113. ISSN 1212-6950.
^ Ши, Вэй; Джонсон, Дж. Карл (2001). «Изучение перевеса гистограмм и конечно-размерного масштабирования жидкостей Леннарда–Джонса». Fluid Phase Equilibria . 187–188: 171–191. Bibcode : 2001FlPEq.187..171S. doi : 10.1016/S0378-3812(01)00534-9.
^ ab Дуников, ДО; Малышенко, СП; Жаховский, ВВ (2001-10-08). "Закон соответствующих состояний и моделирование молекулярной динамики жидкости Леннарда-Джонса". Журнал химической физики . 115 (14): 6623–6631. Bibcode :2001JChPh.115.6623D. doi :10.1063/1.1396674. ISSN 0021-9606.
^ Ливия Б. Партай, Кристоф Ортнер, Альберт П. Барток, Крис Дж. Пикард и Габор Чани «Политипизм в структуре основного состояния Леннарда-Джонсиума», Физическая химия Химическая физика 19 19369 (2017)
^ Стефан, Саймон; Лю, Джинлу; Лангенбах, Кай; Чепмен, Уолтер Г.; Хассе, Ханс (2018). «Интерфейс пар-жидкость усеченной и смещенной жидкости Леннарда-Джонса: сравнение молекулярного моделирования, теории градиента плотности и теории функционала плотности». Журнал физической химии C. 122 ( 43): 24705–24715. doi :10.1021/acs.jpcc.8b06332. ISSN 1932-7447. S2CID 105759822.
^ Коб, Вальтер; Андерсен, Ганс К. (1995-05-01). "Тестирование теории связи мод для переохлажденной бинарной смеси Леннарда-Джонса I: Корреляционная функция Ван Хоува". Physical Review E. 51 ( 5): 4626–4641. arXiv : cond-mat/9501102 . Bibcode : 1995PhRvE..51.4626K. doi : 10.1103/PhysRevE.51.4626. PMID 9963176. S2CID 17662741.
^ abcde Стефан, Саймон; Хассе, Ганс (2020-06-01). "Влияние дисперсионных дальнодействующих взаимодействий на свойства равновесий пар–жидкость и интерфейсов бинарных смесей Леннарда-Джонса". Молекулярная физика . 118 (9–10): e1699185. Bibcode : 2020MolPh.11899185S. doi : 10.1080/00268976.2019.1699185. ISSN 0026-8976. S2CID 214174102.
^ ab Smit, B. (1992), "Фазовые диаграммы жидкостей Леннарда-Джонса" (PDF) , Журнал химической физики , 96 (11): 8639–8640, Bibcode : 1992JChPh..96.8639S, doi : 10.1063/1.462271
^ Смит, Б.; Френкель, Д. (1991-04-15). "Парожидкостные равновесия двумерной жидкости Леннарда-Джонса". Журнал химической физики . 94 (8): 5663–5668. Bibcode :1991JChPh..94.5663S. doi :10.1063/1.460477. ISSN 0021-9606. S2CID 1580499.
^ Scalise, Osvaldo H (июнь 2001 г.). "Тип I газожидкостных равновесий двумерной бинарной смеси Леннарда–Джонса". Fluid Phase Equilibria . 182 (1–2): 59–64. Bibcode : 2001FlPEq.182...59S. doi : 10.1016/s0378-3812(01)00380-6. ISSN 0378-3812.
^ Hloucha, M.; Sandler, SI (ноябрь 1999). «Фазовая диаграмма четырехмерной жидкости Леннарда-Джонса». Журнал химической физики . 111 (17): 8043–8047. Bibcode : 1999JChPh.111.8043H. doi : 10.1063/1.480138. ISSN 0021-9606.
^ abc Хилл, Террелл Л. (1956). Статистическая механика: принципы и избранные приложения. Нью-Йорк: Dover Publications. ISBN0-486-65390-0. OCLC 15163657.
^ ab Wood, WW; Parker, FR (1957). «Уравнение состояния молекул, взаимодействующих с потенциалом Леннарда-Джонса, методом Монте-Карло. I. Сверхкритическая изотерма при температуре, примерно в два раза превышающей критическую». Журнал химической физики . 27 (3): 720–733. Bibcode : 1957JChPh..27..720W. doi : 10.1063/1.1743822. ISSN 0021-9606.
^ ab Rosenbluth, Marshall N.; Rosenbluth, Arianna W. (1954). «Дальнейшие результаты по уравнениям состояния Монте-Карло». Журнал химической физики . 22 (5): 881–884. Bibcode : 1954JChPh..22..881R. doi : 10.1063/1.1740207. ISSN 0021-9606.
^ Alder, BJ; Wainwright, TE (1959). «Исследования по молекулярной динамике. I. Общий метод». Журнал химической физики . 31 (2): 459–466. Bibcode : 1959JChPh..31..459A. doi : 10.1063/1.1730376. ISSN 0021-9606.
^ Рахман, А. (1964-10-19). "Корреляции в движении атомов в жидком аргоне". Physical Review . 136 (2A): A405–A411. Bibcode :1964PhRv..136..405R. doi :10.1103/PhysRev.136.A405. ISSN 0031-899X.
^ abcdefgh Шульц, Эндрю Дж.; Кофке, Дэвид А. (2018-11-28). «Комплексное высокоточное высокоточное уравнение состояния и свойства сосуществования для классических кристаллов Леннарда-Джонса и низкотемпературных жидких фаз». Журнал химической физики . 149 (20): 204508. Bibcode : 2018JChPh.149t4508S. doi : 10.1063/1.5053714. ISSN 0021-9606. PMID 30501268. S2CID 54629914.
^ ab Шульц, Эндрю Дж.; Кофке, Дэвид А. (2020-08-07). "Erratum: "Комплексное высокоточное высокоточное уравнение состояния и свойства сосуществования для классических кристаллов Леннарда-Джонса и низкотемпературных жидких фаз" [J. Chem. Phys. 149, 204508 (2018)]". Журнал химической физики . 153 (5): 059901. doi : 10.1063/5.0021283 . ISSN 0021-9606. PMID 32770918.
^ ab Schappals, Michael; Mecklenfeld, Andreas; Kröger, Leif; Botan, Vitalie; Köster, Andreas; Stephan, Simon; García, Edder J.; Rutkai, Gabor; Raabe, Gabriele; Klein, Peter; Leonhard, Kai (2017-09-12). "Round Robin Study: Molecular Simulation of Thermodynamic Properties from Models with Internal Degrees of Freedom". Journal of Chemical Theory and Computation . 13 (9): 4270–4280. doi :10.1021/acs.jctc.7b00489. ISSN 1549-9618. PMID 28738147.
^ Лёффлер, Ханнес Х.; Босисио, Стефано; Дуарте Рамос Матос, Гильерме; Су, Донхёк; Ру, Бенуа; Мобли, Дэвид Л.; Мишель, Жюльен (13.11.2018). «Воспроизводимость расчетов свободной энергии в различных пакетах программного обеспечения для молекулярного моделирования». Журнал химической теории и вычислений . 14 (11): 5567–5582. doi :10.1021/acs.jctc.8b00544. hdl : 20.500.11820/52d85d71-d3df-468b-8f88-9c52e83da1f1 . ISSN 1549-9618. PMID 30289712. S2CID 52923832.
^ ab Ленхард, Йоханнес; Кюстер, Уве (2019). «Воспроизводимость и концепция численного решения». Minds and Machines . 29 (1): 19–36. doi :10.1007/s11023-019-09492-9. ISSN 0924-6495. S2CID 59159685.
^ Ши, Вэй; Джонсон, Дж. Карл (2001-09-15). "Изучение перевеса гистограмм и конечно-размерного масштабирования жидкостей Леннарда–Джонса". Fluid Phase Equilibria . 187–188: 171–191. Bibcode :2001FlPEq.187..171S. doi :10.1016/S0378-3812(01)00534-9. ISSN 0378-3812.
^ ab Метрополис, Николас; Розенблут, Арианна В.; Розенблут, Маршалл Н.; Теллер, Августа Х.; Теллер, Эдвард (1953). «Вычисления уравнений состояния с помощью быстрых вычислительных машин». Журнал химической физики . 21 (6): 1087–1092. Bibcode : 1953JChPh..21.1087M. doi : 10.1063/1.1699114. ISSN 0021-9606. OSTI 4390578. S2CID 1046577.
^ Дэниел Сидериус (01.09.2017). Сайт эталонного моделирования стандартов NIST — SRD 173. Национальный институт стандартов и технологий. doi : 10.18434/T4M88Q.
^ "Модифицированное масштабирование энтропии транспортных свойств жидкости Леннарда-Джонса". dx.doi.org . doi :10.1021/acs.jpcb.9b05808.s001 . Получено 20 мая 2024 г. .
^ Лаутеншлегер, Мартин П.; Хассе, Ганс (29.04.2021). «Транспортные свойства усеченной и смещенной жидкости Леннарда-Джонса из неравновесного молекулярно-динамического моделирования». dx.doi.org . doi :10.31219/osf.io/67w4u . Получено 20.05.2024 .
^ ab Köster, Andreas; Mausbach, Peter; Vrabec, Jadran (2017-10-10). «Предплавление, равновесие твердое тело–жидкость и термодинамические свойства в области высокой плотности на основе потенциала Леннарда-Джонса». Журнал химической физики . 147 (14): 144502. Bibcode : 2017JChPh.147n4502K. doi : 10.1063/1.4990667. ISSN 0021-9606. PMID 29031254.
^ abcdefghi Стефан, Саймон; Дейтерс, Ульрих К. (2020-08-20). "Характерные кривые жидкости Леннарда-Джонса". International Journal of Thermophysics . 41 (10): 147. Bibcode : 2020IJT....41..147S. doi : 10.1007/s10765-020-02721-9. ISSN 1572-9567. PMC 7441092. PMID 32863513 .
^ ab Deiters, Ulrich K.; Neumaier, Arnold (2016-08-11). «Компьютерное моделирование характеристических кривых чистых жидкостей». Journal of Chemical & Engineering Data . 61 (8): 2720–2728. doi :10.1021/acs.jced.6b00133. ISSN 0021-9568.
^ ab Агравал, Рупал; Кофке, Дэвид А. (1995). «Термодинамические и структурные свойства модельных систем при сосуществовании твердого тела и жидкости: II. Плавление и сублимация системы Леннарда-Джонса». Молекулярная физика . 85 (1): 43–59. doi :10.1080/00268979500100921. ISSN 0026-8976.
^ Браун, Э. Х. (1960). «О термодинамических свойствах жидкостей». Бюллетень Международного института холода . Приложение 1960-1: 169–178.
^ Апфельбаум, Э.М.; Воробьев, В.С. (2020-06-18). «Линия единичного фактора сжимаемости (Zeno-Line) для кристаллических состояний». Журнал физической химии B. 124 ( 24): 5021–5027. doi :10.1021/acs.jpcb.0c02749. ISSN 1520-6106. PMID 32437611. S2CID 218835048.
^ Апфельбаум, EM; Воробьев, VS; Мартынов, GA (2008). «О теории линии Зенона». Журнал физической химии A. 112 ( 26): 6042–6044. Bibcode :2008JPCA..112.6042A. doi :10.1021/jp802999z. ISSN 1089-5639. PMID 18543889.
^ Николас, Дж. Дж.; Габбинс, К. Э.; Стритт, В. Б.; Тилдесли, Д. Д. (1979). «Уравнение состояния для жидкости Леннарда-Джонса». Молекулярная физика . 37 (5): 1429–1454. Bibcode : 1979MolPh..37.1429N. doi : 10.1080/00268977900101051. ISSN 0026-8976.
^ Фэн, Чао; Шульц, Эндрю Дж.; Чаудхари, Випин; Кофке, Дэвид А. (2015-07-28). "Восьмой-шестнадцатый вириальные коэффициенты модели Леннарда-Джонса". Журнал химической физики . 143 (4): 044504. Bibcode : 2015JChPh.143d4504F. doi : 10.1063/1.4927339 . ISSN 0021-9606. PMID 26233142.
^ Шульц, Эндрю Дж.; Кофке, Дэвид А. (2009-11-10). «Шестой, седьмой и восьмой вириальные коэффициенты модели Леннарда-Джонса». Молекулярная физика . 107 (21): 2309–2318. Bibcode : 2009MolPh.107.2309S. doi : 10.1080/00268970903267053. ISSN 0026-8976. S2CID 94811614.
^ Белл, Ян Х.; Мессерли, Ричард; Тол, Моника; Костильола, Лоренцо; Дайр, Йеппе К. (2019-07-25). «Модифицированное масштабирование энтропии транспортных свойств жидкости Леннарда-Джонса». Журнал физической химии B. 123 ( 29): 6345–6363. doi :10.1021/acs.jpcb.9b05808. ISSN 1520-6106. PMC 7147083. PMID 31241958 .
^ Лаутеншлегер, Мартин П.; Хассе, Ганс (2019). «Транспортные свойства усеченной и смещенной жидкости Леннарда-Джонса из неравновесных симуляций молекулярной динамики». Fluid Phase Equilibria . 482 : 38–47. Bibcode : 2019FlPEq.482...38L. doi : 10.1016/j.fluid.2018.10.019. S2CID 106113718.
^ ab Zhen, Shu; Davies, GJ (1983-08-16). «Расчет параметров потенциальной энергии Леннарда-Джонса–m для металлов». Physica Status Solidi A (на немецком языке). 78 (2): 595–605. Bibcode :1983PSSAR..78..595Z. doi :10.1002/pssa.2210780226.
^ ab Halicioglu, T.; Pound, GM (1975-08-16). «Вычисление параметров потенциальной энергии из свойств кристаллического состояния». Physica Status Solidi A. 30 ( 2): 619–623. Bibcode :1975PSSAR..30..619H. doi :10.1002/pssa.2210300223.
^ abc Travesset, Alex (2014-10-28). "Фазовая диаграмма степенного закона и систем Леннарда-Джонса: кристаллические фазы". Журнал химической физики . 141 (16): 164501. Bibcode : 2014JChPh.141p4501T. doi : 10.1063/1.4898371. ISSN 0021-9606. PMID 25362319.
^ Лонге-Хиггинс, ХК (1951-02-07). "Статистическая термодинамика многокомпонентных систем". Труды Лондонского королевского общества. Серия A. Математические и физические науки . 205 (1081): 247–269. Bibcode :1951RSPSA.205..247L. doi :10.1098/rspa.1951.0028. ISSN 0080-4630. S2CID 202575459.
^ Leland, TW; Rowlinson, JS; Sather, GA (1968). "Статистическая термодинамика смесей молекул разных размеров". Transactions of the Faraday Society . 64 : 1447. doi :10.1039/tf9686401447. ISSN 0014-7672.
^ Мансури, Г. Али; Леланд, Томас В. (1972). «Статистическая термодинамика смесей. Новая версия теории конформного решения». Журнал химического общества, Faraday Transactions 2. 68 : 320. doi :10.1039/f29726800320. ISSN 0300-9238.
^ Роулинсон, Дж. С.; Суинтон, Ф. Л. (1982). Жидкости и жидкие смеси (Третье изд.). Лондон: Баттерворт.
^ ab Deiters, Ulrich K.; Kraska, Thomas (2012). Высоконапорные фазовые равновесия жидкости: феноменология и вычисления (1-е изд.). Амстердам: Elsevier. ISBN978-0-444-56354-5. OCLC 787847134.
^ аб Шнабель, Торстен; Врабец, Ядран; Хассе, Ганс (2007). «В отличие от параметров Леннарда – Джонса для равновесия пар-жидкость». Журнал молекулярных жидкостей . 135 (1–3): 170–178. arXiv : 0904.4436 . doi :10.1016/j.molliq.2006.12.024. S2CID 16111477.
^ Лоренц, HA (1881). «Ueber die Anwendung des Satzes vom Virial in der kinetischen Theorie der Gase». Аннален дер Физик (на немецком языке). 248 (1): 127–136. Бибкод : 1881АнП...248..127Л. дои : 10.1002/andp.18812480110.
^ van Konynenburg, PH; Scott, RL (1980-12-18). «Критические линии и фазовые равновесия в бинарных смесях Ван-дер-Ваальса». Philosophical Transactions of the Royal Society of London. Серия A, Математические и физические науки . 298 (1442): 495–540. Bibcode :1980RSPTA.298..495K. doi :10.1098/rsta.1980.0266. ISSN 0080-4614. S2CID 122538015.
^ Потофф, Джеффри Дж.; Панагиотопулос, Афанассиос З. (1998-12-22). «Критическая точка и фазовое поведение чистой жидкости и смеси Леннарда-Джонса». Журнал химической физики . 109 (24): 10914–10920. Bibcode : 1998JChPh.10910914P. doi : 10.1063/1.477787. ISSN 0021-9606.
^ Проценко, Сергей П.; Байдаков, Владимир Г. (2016). «Бинарные смеси Леннарда-Джонса с высокоасимметричными взаимодействиями компонентов. 1. Влияние энергетических параметров на фазовые равновесия и свойства интерфейсов жидкость–газ». Fluid Phase Equilibria . 429 : 242–253. Bibcode : 2016FlPEq.429..242P. doi : 10.1016/j.fluid.2016.09.009.
^ Проценко, Сергей П.; Байдаков, Владимир Г.; Брюханов, Василий М. (2016). «Бинарные смеси Леннарда-Джонса с высокоасимметричными взаимодействиями компонентов. 2. Влияние размера частиц на фазовые равновесия и свойства интерфейсов жидкость–газ». Fluid Phase Equilibria . 430 : 67–74. Bibcode :2016FlPEq.430...67P. doi :10.1016/j.fluid.2016.09.022.
^ Стефан, Саймон; Хассе, Ганс (2020-01-23). «Молекулярные взаимодействия на границах раздела пар–жидкость: бинарные смеси простых жидкостей». Physical Review E. 101 ( 1): 012802. Bibcode : 2020PhRvE.101a2802S. doi : 10.1103/PhysRevE.101.012802. ISSN 2470-0045. PMID 32069593. S2CID 211192904.
^ ab Lamm, Monica H.; Hall, Carol K. (2002). "Равновесия между твердой, жидкой и паровой фазами в бинарных смесях Леннарда–Джонса". Fluid Phase Equilibria . 194–197: 197–206. Bibcode :2002FlPEq.194..197L. doi :10.1016/S0378-3812(01)00650-1.
^ ab Lamm, Monica H.; Hall, Carol K. (2001). "Моделирование методом Монте-Карло полных фазовых диаграмм для бинарных смесей Леннарда–Джонса". Fluid Phase Equilibria . 182 (1–2): 37–46. Bibcode : 2001FlPEq.182...37L. doi : 10.1016/S0378-3812(01)00378-8.
^ ab Hitchcock, Monica R.; Hall, Carol K. (1999-06-15). "Равновесие твердой и жидкой фаз для бинарных смесей Леннарда-Джонса". Журнал химической физики . 110 (23): 11433–11444. Bibcode : 1999JChPh.11011433H. doi : 10.1063/1.479084. ISSN 0021-9606.
^ Юнгблут, Светлана; Деллаго, Кристоф (2011-03-14). «Кристаллизация бинарной смеси Леннарда-Джонса». Журнал химической физики . 134 (10): 104501. Bibcode : 2011JChPh.134j4501J. doi : 10.1063/1.3556664. ISSN 0021-9606. PMID 21405169.
^ ab Lamm, Monica H.; Hall, Carol K. (2004). «Влияние давления на полное фазовое поведение бинарных смесей». AIChE Journal . 50 (1): 215–225. Bibcode : 2004AIChE..50..215L. doi : 10.1002/aic.10020 . ISSN 0001-1541.
^ Коттерман, Р. Л.; Праусниц, Дж. М. (1986). «Молекулярная термодинамика жидкостей при низких и высоких плотностях. Часть II: Фазовые равновесия для смесей, содержащих компоненты с большими различиями в размере молекул или потенциальной энергии». Журнал AIChE . 32 (11): 1799–1812. Bibcode : 1986AIChE..32.1799C. doi : 10.1002/aic.690321105. ISSN 0001-1541. S2CID 96417239.
^ Мюллер, Андреас; Винкельманн, Йохен; Фишер, Иоганн (1996). «Семейство уравнений состояния Backone: 1. Неполярные и полярные чистые жидкости». Журнал AIChE . 42 (4): 1116–1126. Bibcode : 1996AIChE..42.1116M. doi : 10.1002/aic.690420423. ISSN 0001-1541.
^ Blas, FJ; Vega, LF (1997). "Термодинамическое поведение гомоядерных и гетероядерных цепей Леннарда-Джонса с сайтами ассоциации из моделирования и теории". Молекулярная физика . 92 (1): 135–150. Bibcode : 1997MolPh..92..135F. doi : 10.1080/002689797170707. ISSN 0026-8976.
^ Краска, Томас; Габбинс, Кит Э. (1996). «Расчеты фазовых равновесий с модифицированным уравнением состояния SAFT. 1. Чистые алканы, алканолы и вода». Industrial & Engineering Chemistry Research . 35 (12): 4727–4737. doi :10.1021/ie9602320. ISSN 0888-5885.
^ Гонасги, Д.; Чепмен, Уолтер Г. (1994). «Прогнозирование свойств модельных полимерных растворов и смесей». Журнал AIChE . 40 (5): 878–887. Bibcode : 1994AIChE..40..878G. doi : 10.1002/aic.690400514. ISSN 0001-1541.
^ Мекке, М.; Мюллер, А.; Винкельманн, Й.; Врабец, Й.; Фишер, Й.; Спан, Р.; Вагнер, В. (1996-03-01). «Точное уравнение состояния типа Ван-дер-Ваальса для жидкости Леннарда-Джонса». International Journal of Thermophysics . 17 (2): 391–404. Bibcode : 1996IJT....17..391M. doi : 10.1007/BF01443399. ISSN 1572-9567. S2CID 123304062.
^ Тол, Моника; Руткай, Габор; Кёстер, Андреас; Люстиг, Рольф; Спан, Роланд; Врабец, Ядран (2016). «Уравнение состояния для жидкости Леннарда-Джонса». Журнал физических и химических справочных данных . 45 (2): 023101. Bibcode : 2016JPCRD..45b3101T. doi : 10.1063/1.4945000. ISSN 0047-2689.
^ Готтшалк, Маттиас (2019-12-01). "Уравнение состояния для жидкости Леннарда-Джонса: подход вириального расширения". AIP Advances . 9 (12): 125206. Bibcode : 2019AIPA....9l5206G. doi : 10.1063/1.5119761 . ISSN 2158-3226.
^ Джонсон, Дж. Карл; Цоллвег, Джон А.; Габбинс, Кит Э. (1993-02-20). «Повторное рассмотрение уравнения состояния Леннарда-Джонса». Молекулярная физика . 78 (3): 591–618. Bibcode : 1993MolPh..78..591J. doi : 10.1080/00268979300100411. ISSN 0026-8976.
^ ab Eggimann, Becky L.; Sunnarborg, Amara J.; Stern, Hudson D.; Bliss, Andrew P.; Siepmann, J. Ilja (2014-01-02). "Онлайн-база данных параметров и свойств для силового поля TraPPE". Molecular Simulation . 40 (1–3): 101–105. doi :10.1080/08927022.2013.842994. ISSN 0892-7022. S2CID 95716947.
^ Йоргенсен, Уильям Л.; Максвелл, Дэвид С.; Тирадо-Райвс, Джулиан (январь 1996 г.). «Разработка и тестирование силового поля OPLS All-Atom Force Field на конформационной энергетике и свойствах органических жидкостей». Журнал Американского химического общества . 118 (45): 11225–11236. doi :10.1021/ja9621760. ISSN 0002-7863.
^ ab Стефан, Саймон; Хорш, Мартин Т.; Врабец, Ядран; Хассе, Ханс (2019-07-03). "MolMod – открытая база данных силовых полей для молекулярного моделирования жидкостей". Молекулярное моделирование . 45 (10): 806–814. arXiv : 1904.05206 . doi : 10.1080/08927022.2019.1601191. ISSN 0892-7022. S2CID 119199372.
^ Менделев, MI; Хан, S.; Сроловиц, DJ; Экленд, GJ; Сан, DY; Аста, M. (2003). «Разработка новых межатомных потенциалов, подходящих для кристаллического и жидкого железа». Philosophical Magazine . 83 (35): 3977–3994. Bibcode :2003PMag...83.3977A. doi :10.1080/14786430310001613264. ISSN 1478-6435. S2CID 4119718.
^ Абаскал, Дж. Л. Ф.; Вега, К. (15.12.2005). "Универсальная модель для конденсированных фаз воды: TIP4P/2005". Журнал химической физики . 123 (23): 234505. Bibcode : 2005JChPh.123w4505A. doi : 10.1063/1.2121687. ISSN 0021-9606. PMID 16392929.