Из Википедии, бесплатной энциклопедии
Перейти к навигации Перейти к поиску
Рисунок 1 . График потенциальной функции Леннарда-Джонса: межмолекулярная потенциальная энергия как функция расстояния пары частиц. Потенциальный минимум составляет .

Потенциал Леннард-Джонс (также называется потенциалом ЛДА или 12-6 потенциала ) представляет собой межмолекулярный парный потенциал. Среди межмолекулярных потенциалов потенциал Леннарда-Джонса играет центральную роль, как вода среди реальных жидкостей: это потенциал, который изучен наиболее широко и тщательно. Он считается прототипом простых, но реалистичных межмолекулярных взаимодействий.

Потенциал Леннарда-Джонса моделирует мягкие отталкивающие и притягивающие взаимодействия. Следовательно, потенциал Леннарда-Джонса описывает электронно-нейтральные атомы или молекулы. Он назван в честь Джона Леннард-Джонса . [1] [2] [3] Обычно используется выражение для потенциала Леннарда-Джонса:

где - расстояние между двумя взаимодействующими частицами, - глубина потенциальной ямы (обычно называемая «энергией дисперсии»), и - расстояние, на котором потенциальная энергия частица-частица равна нулю (часто называемая «размером частица '). Потенциал Леннарда-Джонса имеет минимум на расстоянии , где потенциальная энергия имеет значение .

Потенциал Леннарда-Джонса - это упрощенная модель, которая все же описывает существенные особенности взаимодействия между простыми атомами и молекулами: две взаимодействующие частицы отталкиваются друг от друга на очень близком расстоянии, притягиваются друг к другу на умеренном расстоянии и не взаимодействуют на бесконечном расстоянии, см. Рисунок 1. Потенциал Леннарда-Джонса - это парный потенциал, т.е. никакие трех- или многочастичные взаимодействия не покрываются потенциалом.

Статистическая механика [4] и компьютерное моделирование [5] [6] могут быть использованы для изучения потенциала Леннарда-Джонса и получения теплофизических свойств «вещества Леннарда-Джонса». И потенциал Леннарда-Джонса, и, соответственно, субстанция Леннарда-Джонса являются упрощенными, но реалистичными моделями, например, они точно отражают основные физические принципы, такие как наличие критической и тройной точки , конденсация и замерзание и т. Д. Потенциал математически прост и поэтому широко используется в исследованиях материи с первых дней компьютерного моделирования. [7] [8] [9] [10] Благодаря своей математической простоте и универсальным возможностям моделирования потенциал Леннарда-Джонса, вероятно, до сих пор остается наиболее изучаемым модельным потенциалом. [11] [12] Вещество Леннарда-Джонса часто даже называют «Леннард-Джонсиум», предполагая, что оно рассматривается как химический элемент . Потенциал Леннарда-Джонса обычно является стандартным выбором для разработки теорий материи (особенно мягкой), а также для разработки и тестирования вычислительных методов и алгоритмов. После настройки параметров модели и реальных свойств вещества потенциал Леннарда-Джонса можно использовать для описания простых веществ (например, благородных газов.) с хорошей точностью. Кроме того, потенциал Леннарда-Джонса часто используется в качестве строительного блока в молекулярных моделях (так называемых силовых полях ) для более сложных веществ. [13] [14] [15] [16] [17]

Физические основы и математические детали [ править ]

Потенциал Леннарда-Джонса моделирует два наиболее важных и фундаментальных молекулярных взаимодействия: член отталкивания ( член) описывает отталкивание Паули на коротких расстояниях взаимодействующих частиц из-за перекрытия электронных орбиталей, а член притяжения ( член) описывает притяжение при дальнодействующих взаимодействиях. ( дисперсионная сила ), которые исчезают на бесконечном расстоянии между двумя частицами. Крутые отталкивающие взаимодействия на малых расстояниях приводят к низкой сжимаемости твердой и жидкой фаз; притягивающие дисперсионные взаимодействия стабилизируют конденсированную фазу, особенно парожидкостное равновесие .

Функциональная форма члена притяжения, то есть показатель степени «6», имеет физическое обоснование, которое не выполняется так строго для члена отталкивания с показателем степени «12». Притягивающие дисперсионные взаимодействия между простыми атомами и молекулами являются результатом флуктуирующих частичных зарядов. Квантово-химические расчеты показали, что этот дисперсионный вклад должен затухать с . [18]

Этот термин используется в основном потому, что он может быть реализован с вычислительной точки зрения очень эффективно как квадрат , который не выполняется в той же степени для значений, отличных от «12». Кроме того, достаточно хорошо аппроксимирует отталкивание Паули . Потенциал Леннарда-Джонса можно обобщить, используя произвольные показатели вместо 12 и 6. Полученный потенциал называется потенциалом Ми. В данной статье рассматривается исключительно классический (12-6) потенциал Леннарда-Джонса.

Потенциал Леннарда-Джонса имеет полюс при , т. Е. Потенциальная энергия расходится , что может вызвать нестабильность в молекулярном моделировании, например, при измерении химического потенциала. Потенциал Леннарда-Джонса сходится к для . Следовательно, с математической точки зрения, притягивающие взаимодействия присутствуют для бесконечно удаленных частиц. Эти дисперсионные «дальнодействующие» взаимодействия имеют важное влияние на несколько свойств вещества Леннард-Джонса, например, давление или теплоемкость в окрестности критической точки и самой критической точки. Важность дальнодействующих взаимодействий была замечена уже на ранних этапах статистической механики . [19]Для компьютерного моделирования можно использовать только конечное число частиц, что приводит к тому, что потенциал может быть оценен только до конечного радиуса , что является так называемым эффектом конечного размера. Существуют хорошо зарекомендовавшие себя методы для неявного учета игнорируемого таким образом вклада дальнего действия для данной наблюдаемой (подробности приведены ниже).

Часто утверждают, что существуют множественные потенциалы Леннарда-Джонса и соответствующие вещества в зависимости от того, как обращаться с дальнодействующими взаимодействиями. Это заблуждение. Существует только один «потенциал Леннарда-Джонса», который точно определяется формулой. (1). Потенциал Леннарда-Джонса требует рассмотрения и оценки дальнодействующих взаимодействий вплоть до очень больших (фактически бесконечных) расстояний - по крайней мере, так, чтобы влияние усечения не влияло на наблюдаемую, представляющую интерес для указанных десятичных знаков.

Потенциал Леннарда-Джонса подразумевает, что частицы являются точечными массами с массой . Хотя этот параметр часто называют «размером частицы», частицы, взаимодействующие с потенциалом Леннарда-Джонса, не имеют однозначно определенного «размера» - в отличие от потенциала твердой сферы . Частицы, взаимодействующие с потенциалом Леннарда-Джонса, скорее имеют мягкие ядра отталкивания.

Модель Леннарда-Джонса описывает потенциальную межмолекулярную энергию между двумя частицами на основе изложенных принципов. После механики Ньютона , фактическая сила между двумя взаимодействующими частицами просто получается дифференцирование потенциала Леннард-Джонс относительно , то есть . В зависимости от расстояния между двумя частицами результирующая сила может быть либо притягивающей, либо отталкивающей.

Потенциал Леннарда-Джонса дает хорошее приближение межмолекулярных взаимодействий для многих приложений: макроскопические свойства, рассчитанные с использованием потенциала Леннарда-Джонса, с одной стороны, хорошо согласуются с экспериментальными данными для простых веществ, таких как аргон, с одной стороны, а потенциальная функция удовлетворительно согласуется с результаты квантовой химии с другой стороны. Потенциал Леннарда-Джонса дает хорошее описание молекулярных взаимодействий в жидких фазах , тогда как молекулярные взаимодействия в твердых фазах описываются лишь приблизительно. В основном это связано с тем, что многочастичные взаимодействия играют важную роль в твердых фазах, которые не входят в потенциал Леннарда-Джонса. Поэтому потенциал Леннарда-Джонса широко используется вфизика мягкой материи и связанные с ней области, тогда как в физике твердого тела она используется реже . Из-за своей простоты потенциал Леннарда-Джонса часто используется для описания свойств газов и простых жидкостей, а также для моделирования дисперсионных и отталкивающих взаимодействий в молекулярных моделях . Это особенно верно для атомов благородных газов и метана . Кроме того, это хорошее приближение для молекулярных взаимодействий на больших и малых расстояниях для нейтральных атомов и молекул. Поэтому потенциал Леннарда-Джонса очень часто используется в качестве строительного блока молекулярных моделей сложных молекул, например алканов или воды . [16][20] [15] Потенциал Леннарда-Джонса также можно использовать для моделирования адсорбционных взаимодействий на границах раздела твердое тело-жидкость, то есть физадсорбции или хемосорбции .

Хорошо известно, что основные ограничения потенциала Леннарда-Джонса заключаются в том, что потенциал является парным (не охватывает взаимодействия нескольких тел), и что показатель степени используется для отталкивания. Результаты квантовой химии предполагают, что следует использовать более высокий показатель степени, чем 12, т.е. более крутой потенциал. Кроме того, потенциал Леннарда-Джонса имеет ограниченную гибкость, т.е. только два параметра модели, и может использоваться для подгонки для описания реального вещества.

Многочисленные межмолекулярные потенциалы были предложены в прошлом для моделирования простых мягких отталкивающих и привлекательных взаимодействий между сферический симметричными частицами, то есть общая форма показана на рисунке 1. Примеров других потенциалов являются потенциалом Морзе , потенциал Ми, [21] Потенциал Букингема и потенциал Тан-Тонни. [22] Тем не менее, ни один из них не имеет такого общего значения, как потенциал Леннарда-Джонса.

Применение потенциала Леннарда-Джонса в молекулярном моделировании [ править ]

Потенциал Леннарда-Джонса имеет фундаментальное значение не только в вычислительной химии и физике мягкой материи , но и для моделирования реальных веществ. По существу, есть два способа использования потенциала Леннарда-Джонса для этой цели: (1) реальный атом или молекула вещества моделируется непосредственно потенциалом Леннарда-Джонса, который дает очень хорошие результаты для благородных газов и метана., т.е. дисперсно взаимодействующие сферические частицы. В случае метана предполагается, что молекула сферически симметрична, а атомы водорода слиты с атомом углерода в общую единицу. Это упрощение, как правило, может применяться и к более сложным молекулам, но обычно дает плохие результаты. (2) Реальная молекула вещества построена из множества сайтов леннард-джонсовского взаимодействия, которые могут быть связаны либо жесткими связями, либо гибкими дополнительными потенциалами (и в конечном итоге также состоит из других типов потенциалов, например, частичных зарядов). Молекулярные модели (часто называемые « силовым полем ») практически для всех молекулярных и ионных частиц могут быть построены с использованием этой схемы, например, для алканов .

При использовании первого подхода , описанной, молекулярная модель имеет только два параметра потенциала Леннард-Джонса , и которые могут быть использованы для присоединения, например , и часто используются для аргона . Очевидно, что этот подход является лишь хорошим приближением для сферических и просто дисперсионно взаимодействующих молекул и атомов. Прямое использование потенциала Леннарда-Джонса имеет большое преимущество в том, что результаты моделирования и теории для потенциала Леннарда-Джонса могут использоваться напрямую. Следовательно, доступные результаты для потенциала Леннарда-Джонса и вещества могут быть напрямую масштабированы с использованием соответствующих и (см. Сокращенные единицы). Параметры потенциала Леннарда-Джонса ив целом может быть приспособлено к любому желаемому объекту недвижимости. В физике мягкой материи для параметризации обычно используются экспериментальные данные для фазового равновесия пар-жидкость или критической точки; в физике твердого тела используются, скорее, сжимаемость, теплоемкость или постоянные решетки. [23] [24]

Второй изложенный подход к использованию потенциала Леннарда-Джонса в качестве строительного блока вытянутых и сложных молекул является гораздо более сложным. Таким образом, молекулярные модели адаптированы в том смысле, что результаты моделирования применимы только для этой конкретной модели. Этот подход к развитию молекулярных силовых полей сегодня в основном применяется в физике мягкой материи и связанных областях, таких как химическая инженерия . Большое количество силовых полей основано на потенциале Леннарда-Джонса, например, силовое поле TraPPE , [16] силовое поле OPLS, [25] и силовое поле MolMod [15] (обзормолекулярные силовые поля выходят за рамки данной статьи). Для современного моделирования твердотельных материалов используются более сложные многочастичные потенциалы (например, потенциалы EAM [26] ).

Альтернативные обозначения потенциала Леннарда-Джонса [ править ]

Есть несколько различных способов сформулировать потенциал Леннарда-Джонса, помимо уравнения. (1). Альтернативы:

Форма AB

Форма AB часто используется в реализациях программного обеспечения для моделирования, поскольку она удобна с вычислительной точки зрения. Потенциал Леннарда-Джонса можно записать как

где, и . Наоборот, и . Это форма, в которой Леннард-Джонс написал потенциал, названный его именем. [27]

форма n-exp

Форма n-exp является математически более общей формой и может быть записана как

где и - энергия связи молекулы (энергия, необходимая для разделения атомов). Применяя гармоническое приближение в минимуме потенциала (at ), показатель степени и параметр энергии можно связать с жесткостью пружины:

откуда можно рассчитать, если известно. Обычно известны гармонические состояния,, где . также может быть связано с групповой скоростью в кристалле,

где - расстояние решетки, а - масса частицы.

Безразмерные (уменьшенные единицы) [ править ]

Безразмерные приведенные единицы могут быть определены на основе параметров потенциала Леннарда-Джонса, что удобно для молекулярного моделирования. С числовой точки зрения преимущества этой системы единиц включают вычисление значений, которые ближе к единице, использование упрощенных уравнений и возможность легко масштабировать результаты. [28] [5] Эта сокращенная система единиц требует указания параметра размера и параметра энергии потенциала Леннарда-Джонса и массы частицы . Все физические свойства могут быть легко преобразованы с учетом соответствующих размеров, см. Таблицу. Сокращенные единицы часто сокращаются и обозначаются звездочкой.

В общем, восстановленные единицы также могут быть построены на других потенциалах взаимодействия молекул, которые состоят из параметра длины и параметра энергии.

Теплофизические свойства вещества Леннарда-Джонса [ править ]

Рисунок 2 . Фазовая диаграмма вещества Леннард-Джонса. Корреляции и числовые значения для критической точки и тройной точки (точек) взяты из работ. [12] [29] [11] Звездочка указывает на критическую точку. [12] Кружком обозначена тройная точка пар-жидкость-твердое тело, а треугольником - тройная точка пар-твердое тело (ГЦК) -твердое (ГПУ). [29] [30] Сплошными линиями обозначены линии сосуществования двух фаз. [12] [29] Пунктирными линиями обозначена спинодаль пар-жидкость. [11]

Теплофизические свойства вещества Леннарда-Джонса, т.е. частицы, взаимодействующие с потенциалом Леннарда-Джонса, могут быть получены с помощью статистической механики. Некоторые свойства можно вычислить аналитически, то есть с машинной точностью, тогда как большинство свойств можно получить только путем молекулярного моделирования. [5] На последнее обычно накладываются как статистические, так и систематические неопределенности. [31] [12] [32] [33] Вириальные коэффициенты можно, например, вычислить непосредственно из потенциала Леннарда, используя алгебраические выражения [4]поэтому представленные данные не содержат неопределенности. Результаты молекулярного моделирования, например давление при заданной температуре и плотности, имеют как статистические, так и систематические погрешности. [31] [33] Молекулярное моделирование потенциала Леннарда-Джонса, как правило, может быть выполнено с использованием моделирования молекулярной динамики (МД) или моделирования методом Монте-Карло (МК). Для моделирования MC напрямую используется потенциал Леннарда-Джонса , в то время как модели MD всегда основаны на производной потенциала, то есть силе . Эти различия в сочетании с различиями в трактовке дальнодействующих взаимодействий (см. Ниже) могут влиять на расчетные теплофизические свойства. [34] [35]

Поскольку Леннард-Джонсиум является прототипом для моделирования простых, но реалистичных межмолекулярных взаимодействий, большое количество теплофизических свойств было изучено и опубликовано в литературе. [12] Данные компьютерных экспериментов с потенциалом Леннарда-Джонса в настоящее время считаются наиболее точно известными данными в вычислительной химии классической механики. Следовательно, такие данные также в основном используются в качестве эталона для проверки и тестирования новых алгоритмов и теорий. Потенциал Леннарда-Джонса постоянно использовался с первых дней молекулярного моделирования. Первые результаты компьютерных экспериментов для потенциала Леннарда-Джонса были представлены Розенблютом и Розенблютом [8] и Вудом и Паркером [7].после того, как в 1953 г. стало доступно молекулярное моделирование на « быстрых вычислительных машинах » [36] С тех пор во многих исследованиях были представлены данные о веществе Леннарда-Джонса; [12] публично доступно около 50 000 точек данных. Текущее состояние исследований теплофизических свойств вещества Леннарда-Джонса можно обобщить следующим образом. Наиболее полное резюме и цифровая база данных были предоставлены Stephan et al. [12] В настоящее время ни один репозиторий данных не покрывает и не поддерживает эту базу данных (или любую другую потенциальную модель) - даже к данным и результатам, заявленным на веб-сайте NIST, следует относиться с осторожностью (не воспроизводимые и вводящие в заблуждение ссылки [12] ).

На рисунке 2 показана фазовая диаграмма флюида Леннард-Джонса. Фазовые равновесия потенциала Леннарда-Джонса изучались много раз и, соответственно, известны сегодня с хорошей точностью. [29] [12] [37] На рисунке 2 показаны корреляции результатов, полученные из результатов компьютерных экспериментов (следовательно, вместо точек данных показаны линии).

Среднее межмолекулярное взаимодействие леннард-джонсовской частицы сильно зависит от термодинамического состояния, то есть температуры и давления (или плотности). Для твердых тел притягивающее леннард-джонсовское взаимодействие играет доминирующую роль, особенно при низких температурах. Для жидких состояний упорядоченная структура отсутствует по сравнению с твердыми состояниями. Средняя потенциальная энергия на частицу отрицательна. Для газообразных состояний притягивающие взаимодействия потенциала Леннарда-Джонса играют второстепенную роль, поскольку они далеки друг от друга. Основная часть внутренней энергии сохраняется в виде кинетической энергии для газообразных состояний. В сверхкритических состояниях притягивающее взаимодействие Леннарда-Джонса играет второстепенную роль. С повышением температуры средняя кинетическая энергия частиц увеличивается и превышает энергетическую яму потенциала Леннарда-Джонса. Следовательно,частицы в основном взаимодействуют посредством мягкого отталкивающего взаимодействия потенциалов, и средняя потенциальная энергия на частицу соответственно положительна.

В целом, из-за большого промежутка времени, когда потенциал Леннарда-Джонса был изучен, данные о теплофизических свойствах были представлены в литературе, а вычислительные ресурсы были недостаточны для точного моделирования (по современным стандартам), заметный объем данных, как известно, является сомнительным. [12] Тем не менее, во многих исследованиях данные уклонения используются в качестве справочных. Отсутствие репозиториев данных и оценки данных - важный элемент для будущей работы в долгосрочной области потенциальных исследований Леннарда-Джонса.

Характерные точки и кривые [ править ]

Наиболее важными характеристиками потенциала Леннарда-Джонса являются критическая точка и тройная точка пар-жидкость-твердое тело . Они неоднократно изучались в литературе и собраны в [2]. [12] Таким образом, критическая точка находилась в

Указанные неопределенности были рассчитаны на основе стандартного отклонения критических параметров, полученных из наиболее надежных имеющихся наборов данных о равновесии пар-жидкость . [12] Эти неопределенности можно принять как нижний предел точности, с которой критическая точка жидкости может быть получена из результатов молекулярного моделирования.

Рисунок 3 . Характеристические кривые вещества Леннард-Джонса. Толстая черная линия указывает на парожидкостное равновесие; звездочка указывает на критическую точку. Коричневая линия указывает на равновесие твердое тело-жидкость. Другие черные сплошные линии и символы обозначают характеристические кривые Брауна (подробности см. В тексте) вещества Леннарда-Джонса: линии представляют собой результаты уравнения состояния, символы - модели молекулярного моделирования, а треугольники - точные данные в пределе идеального газа, полученные из вириальных коэффициентов. . Данные взяты из. [38] [39] [40]

В настоящее время предполагается, что тройная точка расположена в

Неопределенности представляют собой разброс данных от разных авторов. [29] Критическая точка вещества Леннарда-Джонса изучалась гораздо чаще, чем тройная точка. Как для критической точки, так и для тройной точки пар-жидкость-твердое тело в нескольких исследованиях сообщалось о результатах, выходящих за пределы указанных выше диапазонов. Приведенные выше данные являются в настоящее время правильными и надежными. Тем не менее детерминированность критической температуры и температуры тройной точки все еще неудовлетворительна.

Очевидно, что кривые сосуществования фаз (см. Рис. 2) имеют фундаментальное значение для характеристики потенциала Леннарда-Джонса. Кроме того, характеристические кривые Брауна [41] дают иллюстративное описание существенных особенностей потенциала Леннарда-Джонса. Характеристические кривые Брауна определяются как кривые, на которых определенные термодинамические свойства вещества соответствуют характеристикам идеального газа . Для реальной жидкости, и ее производные могут соответствовать значениям идеального газа для специального , комбинации только в результате фазового правило Гиббса. Полученные точки вместе составляют характеристическую кривую. Определены четыре основные характеристические кривые: одна 0-го порядка (называемая кривой Зенона).) и три кривые 1-го порядка (названные кривой Амагата , Бойля и Чарльза ). Требуется, чтобы характеристическая кривая имела отрицательную или нулевую кривизну на всем протяжении и единственный максимум на диаграмме двойного логарифма давление-температура. Кроме того, характеристические кривые Брауна и вириальные коэффициенты непосредственно связаны в пределе идеального газа и поэтому известны точно при . Как результаты компьютерного моделирования, так и результаты уравнения состояния описаны в литературе для потенциала Леннарда-Джонса. [39] [12] [38] [42] [43]

Точки на кривой Зенона Z имеют коэффициент сжимаемости, равный единице . Кривая Зенона начинается при температуре Бойля , окружает критическую точку и имеет наклон, равный единице в нижнем пределе температуры. [38] Точки на кривой Бойля B имеют . Кривая Бойля начинается с кривой Зенона при температуре Бойля, слабо окружает критическую точку и заканчивается кривой давления пара. Точки на кривой Чарльза (также известной как кривая инверсии Джоуля-Томсона ) имеют и, что более важно , отсутствие изменения температуры при изэнтальпическом дросселировании. Он берет свое начало в в пределе идеального газа пересекает кривую Зенона и заканчивается кривой давления пара. Точки на кривой Амагат имеют . Он также начинается в пределе идеального газа при , окружает критическую точку и другие три характеристические кривые и переходит в область твердой фазы. Подробное обсуждение характеристических кривых потенциала Леннарда-Джонса дано Стефаном и Дейтерсом. [38]

Рисунок 4 . Вириальные коэффициенты потенциала Леннарда-Джонса как функция температуры: второй вириальный коэффициент (вверху) и третий вириальный коэффициент (внизу). Круг показывает температуру Бойля . Результаты взяты из. [38]

Свойства жидкости Леннарда-Джонса [ править ]

Рисунок 5 . Парожидкостное равновесие вещества Леннарда-Джонса: давление пара (вверху), плотности насыщения (в центре) и межфазное натяжение (внизу). Символы указывают на результаты молекулярного моделирования. [44] [12] Линии показывают результаты из уравнения состояния (и теории квадратного градиента для межфазного натяжения). [44] [11]

Свойства жидкости Леннарда-Джонса широко изучались в литературе из-за выдающейся важности потенциала Леннарда-Джонса в физике мягкой материи и связанных областях. К настоящему времени опубликовано около 50 наборов данных компьютерных экспериментов для парожидкостного равновесия . [12] Более того, более 35 000 точек данных по однородным жидкостным состояниям были опубликованы за эти годы и недавно были скомпилированы и оценены на предмет выбросов в базе данных с открытым доступом. [12]

Парожидкостное равновесие вещества Леннард-Джонса в настоящее время известно с точностью, т. Е. Взаимным согласием термодинамически согласованных данных, для давления пара, для плотности насыщенной жидкости, для плотности насыщенного пара, для энтальпии парообразования, и для поверхностного натяжения. [12] Это статус-кво нельзя считать удовлетворительным, учитывая тот факт, что статистические неопределенности, обычно сообщаемые для отдельных наборов данных, значительно ниже указанных выше значений (даже для гораздо более сложных полей молекулярных сил).

Как свойства фазового равновесия, так и свойства однородного состояния при произвольной плотности можно в общем случае получить только из молекулярного моделирования, тогда как вириальные коэффициенты можно вычислить непосредственно из потенциала Леннарда-Джонса. [4] Численные данные для второго и третьего вириальных коэффициентов доступны в широком диапазоне температур. [45] [38] [12] Для более высоких вириальных коэффициентов (до шестнадцатого) количество доступных точек данных уменьшается с увеличением номера вириального коэффициента. [46] [47] Также часто изучаются транспортные свойства (вязкость, теплопроводность и коэффициент самодиффузии) леннард-джонсовской жидкости, [48] [49]но база данных значительно менее плотна, чем для однородных свойств равновесия, таких как данные о внутренней энергии. Кроме того, было разработано большое количество аналитических моделей ( уравнений состояния ) для описания жидкости Леннарда-Джонса (подробности см. Ниже).

Свойства твердого тела Леннарда-Джонса [ править ]

База данных и знания для твердого тела Леннарда-Джонса значительно хуже, чем для жидких фаз, что в основном связано с тем, что потенциал Леннарда-Джонса менее часто используется в приложениях для моделирования твердых веществ. Было рано осознавать, что взаимодействия в твердых фазах не следует приближать к попарно аддитивным, особенно для металлов. [23] [24]

Тем не менее, потенциал Леннарда-Джонса по-прежнему часто используется в физике твердого тела из-за его простоты и вычислительной эффективности. Таким образом, основные свойства твердых фаз и фазовых равновесий твердое тело-жидкость исследовались несколько раз, например, в работах [5,11]. [37] [29] [30] [50] [51] [40]

Вещество Леннарда-Джонса образует как ГЦК (гранецентрированная кубическая), так и ГПУ (гексагональная плотноупакованная) решетки - в зависимости от температуры и давления, ср. Рисунок 2. При низкой температуре и умеренном давлении решетка ГПУ энергетически благоприятствует и, следовательно, является равновесной структурой. Структура ГЦК-решетки энергетически выгодна как при высокой температуре, так и при высоком давлении и, следовательно, в целом является равновесной структурой в более широком диапазоне состояний. Линия сосуществования между ГЦК- и ГПУ-фазой начинается примерно в точке , проходит через максимум температуры примерно в точке и затем заканчивается на границе пар-твердая фаза примерно в точке , что, таким образом, образует тройную точку. [50] [29]Следовательно, только твердая фаза с ГЦК имеет фазовое равновесие с жидкой и сверхкритической фазой, ср. фигура 2.

Сообщается, что тройная точка двух твердых фаз (ГЦК и ГПУ) и паровой фазы расположена по адресу: [50] [29]

  • еще не сообщается

Обратите внимание, что в литературе также приводятся другие и значительно отличающиеся значения. Следовательно, база данных для тройной точки ГЦК-ГПУ-пар должна быть дополнительно укреплена в будущем.

Рисунок 6 . Парожидкостные равновесия бинарных смесей Леннарда-Джонса. Во всех показанных случаях компонент 2 является более летучим компонентом (обогащается паровой фазой). Единицы приведены в и компонента 1, который является одинаковым во всех четырех показанных смесей. Температура такая . Символы - это результаты молекулярного моделирования, а линии - результаты уравнения состояния. Данные взяты из Ref. [44]

Смеси веществ Леннард-Джонса [ править ]

Смеси частиц Леннарда-Джонса в основном используются в качестве прототипа для разработки теорий и методов решений, но также и для изучения свойств растворов в целом. Это восходит к фундаментальной работе Лонге-Хиггинса [52], Леланда и Роулинсона с сотрудниками по теории конформных решений . [53] [54] Сегодня это основа большинства теорий смесей. [55] [56]

Смеси двух или более компонентов Леннарда-Джонса устанавливают путем изменения, по меньшей мере, одного потенциального параметра взаимодействия ( или ) одного из компонентов по отношению к другому. Для бинарной смеси это дает три типа парных взаимодействий, которые все моделируются потенциалом Леннарда-Джонса: взаимодействия 1-1, 2-2 и 1-2. Для перекрестных взаимодействий 1-2, дополнительные предположения необходимы для определения параметров или из , и , . Для этих так называемых правил комбинирования можно использовать различные варианты (все более или менее эмпирические и не основанные строго на физических аргументах). [57] Наиболее часто используемым правилом комбинации является правило Лоренца и Бертло.[58]

Параметр является дополнительным параметром взаимодействия смеси, не зависящим от состояния. Параметр обычно устанавливается равным единице, так как среднее арифметическое можно считать физически правдоподобным для параметра размера перекрестного взаимодействия. С другой стороны, параметр часто используется для регулировки фазового поведения модельной смеси. Для аналитических моделей, например, уравнений состояния , параметр отклонения обычно записывается как . Для дисперсионной энергии перекрестного взаимодействия и, соответственно, силы притяжения между разнородными частицами усиливаются. Наоборот, силы притяжения между разнородными частицами уменьшаются при .

Для смесей Леннарда-Джонса могут быть изучены как жидкие, так и твердые фазовые равновесия , т.е. пар-жидкость , жидкость-жидкость , газ-газ, твердое тело-пар, твердое тело-жидкость и твердое тело-твердое тело. Соответственно, могут существовать разные типы тройных точек (трехфазное равновесие) и критические точки, а также разные эвтектические и азеотропные точки . [59] [56] Бинарные смеси Леннарда-Джонса в жидкой области (различные типы равновесий жидкой и газовой фаз) [44] [60] [61] [62] [63]изучены более полно, чем фазовые равновесия, содержащие твердые фазы. [64] [65] [66] [67] [68]

Что касается поведения жидкой фазы, смеси демонстрируют практически идеальное поведение (в смысле закона Рауля ) для . Поскольку преобладают взаимодействия притяжения, и смеси имеют тенденцию образовывать высококипящие азеотропы, то есть для стабилизации парожидкостного равновесия требуется более низкое давление, чем давление паров чистых компонентов. Поскольку преобладают отталкивающие взаимодействия, и смеси имеют тенденцию образовывать низкокипящие азеотропы, то есть для стабилизации парожидкостного равновесия требуется более высокое давление, чем давление пара чистых компонентов, поскольку средние дисперсионные силы уменьшаются. Кроме того, особенно низкие значения могут привести к разрывам смешиваемости жидкость-жидкость. Также различные типы фазовых равновесий, содержащих твердые фазы, были изучены в литературе, например,Кэрол и коллеги. [66] [68] [65] [64] Также существуют случаи, когда границы твердой фазы нарушают равновесие жидкой фазы. Однако для фазовых равновесий, которые включают твердые фазы, количество опубликованных данных скудно.

Уравнения состояния для потенциала Леннарда-Джонса [ править ]

Большое количество уравнений состояния (EOS) для потенциала / вещества Леннарда-Джонса было предложено с тех пор, как его характеристики стали доступны после первого компьютерного моделирования. [36] Из-за фундаментальной важности потенциала Леннарда-Джонса большинство доступных в настоящее время EOS описывают жидкость Леннарда-Джонса. Они были всесторонне рассмотрены Stephan et al. [11] [38]

Уравнения состояния для жидкости Леннарда-Джонса имеют особое значение в физике мягкой материи и физической химии, поскольку они часто используются в качестве отправной точки для разработки EOS для сложных жидкостей, например полимеров и связанных жидкостей. Мономерные звенья этих моделей обычно напрямую адаптированы из EOS Леннарда-Джонса в качестве строительного блока, например, PHC EOS, [69] BACKONE EOS, [70] [71] и EOS типа SAFT. [72] [73] [74] [75]

В литературе предложено более 30 EOS Леннарда-Джонса. Всесторонняя оценка [11] [38] таких УС показала, что несколько УС [76] [77] [78] [79] описывают потенциал Леннарда-Джонса с хорошей и аналогичной точностью, но ни одно из них не является выдающимся. Три из этих EOS демонстрируют неприемлемое нефизическое поведение в некоторой области жидкости, например, множественные петли Ван-дер-Ваальса, хотя в остальном они являются достаточно точными. Только уравнение состояния Леннард-Джонса Колафы и Незбеды [77] оказалось надежным и точным для большинства термодинамических свойств леннард-джонсовской жидкости. [38] [11] Следовательно, уравнение состояния Леннарда-Джонса Колафы и Незбеды [77]в настоящее время считается наиболее полезным выбором, поскольку он надежен и точен. Кроме того, EOS Леннарда-Джонса Johnson et al. [80] оказался менее точным практически для всех имеющихся справочных данных [12] [11], чем EOS Колафы и Незбеды. [77] Интересно отметить, что LJ EOS Johnson et al. [80] все же используется гораздо чаще, чем у Колафы и Незбеды. [77]

Дальнодействующие взаимодействия потенциала Леннарда-Джонса [ править ]

Рисунок 7 . Наглядный пример сходимости схемы коррекции для учета дальнодействующих взаимодействий потенциала Леннарда-Джонса. В нем указывается примерная наблюдаемая и применяемый радиус отсечения. Скорректированное значение на большом расстоянии обозначается как (символы и линия в качестве ориентира для глаза); гипотетическое «истинное» значение как (пунктирная линия).

Потенциал Леннарда-Джонса, ср. Уравнение (1) и рисунок 1 имеет бесконечный диапазон. Только при его рассмотрении исследуется «истинный» и «полный» потенциал Леннарда-Джонса. Для оценки наблюдаемого ансамбля частиц, взаимодействующих с помощью потенциала Леннарда-Джонса, с использованием молекулярного моделирования, взаимодействия могут быть оценены явно только до определенного расстояния - просто из-за того, что количество частиц всегда будет конечным. Максимальное расстояние, применяемое при моделировании, обычно называется радиусом отсечения.(поскольку потенциал Леннарда-Джонса радиально симметричен). Чтобы получить теплофизические свойства (как макроскопические, так и микроскопические) «истинного» и «полного» потенциала Леннарда-Джонса (ЛД), необходимо учитывать вклад потенциала за пределы радиуса отсечки.

Различные схемы коррекции были разработаны для учета влияния дальнодействующих взаимодействий при моделировании и для поддержания достаточно хорошего приближения к «полному» потенциалу. [6] [28] Они основаны на упрощающих предположениях относительно структуры жидкости. Для простых случаев, таких как исследования равновесия однородных жидкостей, простые поправочные члены дают отличные результаты. В других случаях, например, при исследовании неоднородных систем с разными фазами, учет дальнодействующих взаимодействий более утомителен. Эти поправки обычно называют «дальнодействующими поправками». Для большинства свойств известны и хорошо установлены простые аналитические выражения. Для данной наблюдаемой «исправленный» результат моделированиязатем просто вычисляется из фактически выбранного значения и значения дальнодействующей поправки , например, для внутренней энергии . [28] Гипотетическое истинное значение наблюдаемого потенциала Леннарда-Джонса на действительно бесконечном расстоянии отсечки (термодинамический предел) можно в целом только оценить.

Кроме того, качество схемы дальней коррекции зависит от радиуса отсечки. Допущения, сделанные в схемах коррекции, обычно не оправдываются при (очень) коротких радиусах отсечки. Это проиллюстрировано в примере, показанном на рисунке 7. Схема коррекции с большим радиусом действия называется сходящейся, если остающаяся ошибка схемы коррекции достаточно мала на данном расстоянии отсечки, см. рисунок 7.

Усеченный и сдвинутый потенциал Леннарда-Джонса (LJTS) [ править ]

Рисунок 8 . Сравнение парожидкостного равновесия «полного» потенциала Леннарда-Джонса (черный) и «усеченного и смещенного» потенциала Леннарда-Джонса (синий). Символы указывают на результаты молекулярного моделирования; [12] [81] линии показывают результаты из уравнений состояния [11] [82] .

Усеченный и сдвинутый потенциал Леннарда-Джонса (LJTS) - часто используемая альтернатива «полному» потенциалу Леннарда-Джонса (см. Уравнение (1)). «Полный» и «усеченный и сдвинутый» потенциал Леннарда-Джонса следует строго разделять. Это просто два разных потенциала, дающих разные теплофизические свойства. Усеченный и сдвинутый потенциал Леннарда-Джонса определяется как

с

Следовательно, потенциал LJTS устойчиво усечен и сдвинут на соответствующее значение энергии . Последний применяется, чтобы избежать скачка скачка потенциала при . Для потенциала LJTS не рассматриваются дальнейшие взаимодействия - ни явно, ни неявно. Потенциал просто резко заканчивается на . Наиболее часто используемая версия усеченного и сдвинутого потенциала Леннарда-Джонса - это вариант с . Тем не менее в литературе использовались разные значения. [83] [84] [85] Каждый потенциал LJTS с заданным радиусом усечения следует рассматривать как потенциал и, соответственно, как отдельную субстанцию.

Потенциал LJTS в вычислительном отношении значительно дешевле, чем «полный» потенциал Леннарда-Джонса, но все же охватывает основные физические характеристики материи (наличие критической и тройной точек, мягкое отталкивающее и притягивающее взаимодействия, фазовые равновесия и т. Д.). Поэтому потенциал LJTS очень часто используется для тестирования новых алгоритмов, методов моделирования и новых физических теорий. [86] [87]

Интересно, что для однородных систем межмолекулярные силы, которые рассчитываются на основе LJ и потенциала LJTS на заданном расстоянии, одинаковы (так как это одно и то же), тогда как смещение влияет на потенциальную энергию и давление. Кроме того, свойства LJTS-вещества могут быть затронуты выбранным алгоритмом моделирования, т.е. выборкой MD или MC (это, как правило, не относится к «полному» потенциалу Леннарда-Джонса).

Для потенциала LJTS с , потенциальный сдвиг энергии составляет приблизительно 1/60 от энергии дисперсии при потенциальной яме: . На рисунке 8 показано сравнение парожидкостного равновесия «полного» потенциала Леннарда-Джонса и «усеченного и смещенного» потенциала Леннарда-Джонса. Результаты «полного» потенциала Леннарда-Джонса имеют значительно более высокую критическую температуру и давление по сравнению с результатами потенциала LJTS, но критическая плотность очень похожа. [44] [35] [85]Давление пара и энтальпия парообразования находятся под более сильным влиянием дальнодействующих взаимодействий, чем насыщенные плотности. Это связано с тем, что потенциалом манипулируют в основном энергетически путем усечения и сдвига.

Расширения и модификации потенциала Леннарда-Джонса [ править ]

Потенциал Леннарда-Джонса - как архетип межмолекулярных потенциалов - много раз использовался в качестве отправной точки для разработки более сложных межмолекулярных потенциалов. В литературе предлагались различные расширения и модификации потенциала Леннарда-Джонса. Можно утверждать, что все силовые поля (существуют сотни) можно проследить до потенциала Леннарда-Джонса. Более обширный список приведен в статье о функциях «межатомный потенциал». Следующий список относится только к потенциалам, которые напрямую связаны с потенциалом Леннарда-Джонса и имеют как историческое значение, так и все еще актуальны для настоящего исследования.

  • Потенциал Ми Потенциал Ми - это обобщенная версия потенциала Леннарда-Джонса, т.е. показатели 12 и 6 вводятся как параметрыи.Известно, чтоособенно термодинамические производные свойства, например сжимаемость и скорость звука , очень чувствительны к крутизне отталкивающей части межмолекулярного потенциала, которая, следовательно, может быть смоделирована более сложным образом с помощью потенциала Ми. [72] Первая явная формулировка потенциала Ми приписывается Эдуарду Грюнайзену . [88] [89] Следовательно, потенциал Ми был фактически предложен до потенциала Леннарда-Джонса. Потенциал Ми назван в честьГустав Мие . [21]
  • Букингемский потенциал . Букингемский потенциал был предложен Ричардом Букингемом . В нем отталкивающая часть потенциала Леннарда-Джонса заменена экспоненциальной функцией и включает дополнительный параметр.
  • Потенциал Стокмайера Потенциал Стокмайера назван в честь У. С. Стокмайера. [90] Потенциал Стокмайера представляет собой комбинацию потенциала Леннарда-Джонса, наложенного диполем. Следовательно, частицы Стокмайера не сферически симметричны, а имеют важную ориентационную структуру.
  • Два центра ЛДА-Джонс Два центра ЛДА-Джонс состоит из двух идентичных участков взаимодействия ЛДОВ-( то же самых , , ), которые связаны как твердое тело. Часто его обозначают сокращенно 2CLJ. Обычно удлинение (расстояние между леннард-джонсовскими площадками) значительно меньше параметра размера . Следовательно, два сайта взаимодействия существенно слиты.
  • Усеченный и шлицевый потенциал Леннарда-Джонса Усеченный и шлицевый потенциал Леннарда-Джонса - редко используемый, но полезный потенциал. Подобно более популярному потенциалу LJTS, он надежно усечен на определенном «конечном» расстоянии, и за его пределами дальнодействующие взаимодействия не рассматриваются. В отличие от потенциала LJTS, который смещен таким образом, что потенциал является непрерывным, усеченный и шлицевый потенциал Леннарда-Джонса становится непрерывным за счет использования произвольной, но подходящей функции сплайна.

См. Также [ править ]

  • Молекулярная механика
  • Встроенная модель атома
  • Силовое поле (химия)
  • Сравнение реализаций силового поля
  • Потенциал Морзе и Morse / Long диапазон потенциала
  • Вириальное расширение

Ссылки [ править ]

  1. ^ Джонс, JE (1924). «Об определении молекулярных полей. - I. Из изменения вязкости газа с температурой» . Труды Лондонского королевского общества. Серия A, содержащая статьи математического и физического характера . 106 (738): 441–462. Bibcode : 1924RSPSA.106..441J . DOI : 10.1098 / RSPA.1924.0081 . ISSN  0950-1207 .
  2. ^ Джонс, JE (1924). «Об определении молекулярных полей. —II. Из уравнения состояния газа» . Труды Лондонского королевского общества. Серия A, содержащая статьи математического и физического характера . 106 (738): 463–477. Bibcode : 1924RSPSA.106..463J . DOI : 10.1098 / RSPA.1924.0082 . ISSN 0950-1207 . 
  3. ^ Леннард-Джонс, JE (1931-09-01). «Сплоченность» . Труды физического общества . 43 (5): 461–482. Bibcode : 1931PPS .... 43..461L . DOI : 10.1088 / 0959-5309 / 43/5/301 . ISSN 0959-5309 . 
  4. ^ a b c Хилл, Террелл Л. (1956). Статистическая механика: принципы и избранные приложения . Нью-Йорк: Dover Publications. ISBN 0-486-65390-0. OCLC  15163657 .
  5. ^ a b c Д. К. Рапапорт (1 апреля 2004 г.). Искусство моделирования молекулярной динамики . Издательство Кембриджского университета. ISBN 978-0-521-82568-9.
  6. ^ а б Френкель, Д .; Смит Б. (2002), Understanding Molecular Simulation (Second ed.), San Diego: Academic Press, ISBN 0-12-267351-4
  7. ^ a b Дерево, WW; Паркер, FR (1957). "Уравнение Монте-Карло состояния молекул, взаимодействующих с потенциалом Леннарда-Джонса. I. Сверхкритическая изотерма примерно при вдвое большей критической температуре" . Журнал химической физики . 27 (3): 720–733. Bibcode : 1957JChPh..27..720W . DOI : 10.1063 / 1.1743822 . ISSN 0021-9606 . 
  8. ^ a b Розенблют, Маршалл Н .; Розенблут, Арианна В. (1954). «Дальнейшие результаты по уравнениям состояния Монте-Карло» . Журнал химической физики . 22 (5): 881–884. Полномочный код : 1954JChPh..22..881R . DOI : 10.1063 / 1.1740207 . ISSN 0021-9606 . 
  9. Ольга, Би Джей; Уэйнрайт, Т. Е. (1959). «Исследования молекулярной динамики. I. Общий метод» . Журнал химической физики . 31 (2): 459–466. Bibcode : 1959JChPh..31..459A . DOI : 10.1063 / 1.1730376 . ISSN 0021-9606 . 
  10. ^ Рахман, А. (1964-10-19). «Корреляции в движении атомов в жидком аргоне» . Физический обзор . 136 (2A): A405 – A411. Полномочный код : 1964PhRv..136..405R . DOI : 10.1103 / PhysRev.136.A405 . ISSN 0031-899X . 
  11. ^ a b c d e f g h я Стефан, Саймон; Штаубах, Йенс; Хассе, Ганс (2020). «Обзор и сравнение уравнений состояния жидкости Леннарда-Джонса» . Равновесия жидкой фазы . 523 : 112772. дои : 10.1016 / j.fluid.2020.112772 .
  12. ^ Б с д е е г ч я J к л м п о р д т ы т у Стефана, Саймон; Тол, Моника; Врабец, Ядран; Хассе, Ханс (2019-10-28). «Теплофизические свойства жидкости Леннард-Джонса: база данных и оценка данных» . Журнал химической информации и моделирования . 59 (10): 4248–4265. DOI : 10.1021 / acs.jcim.9b00620 . ISSN 1549-9596 . PMID 31609113 .  
  13. ^ Йоргенсен, Уильям Л .; Максвелл, Дэвид С .; Тирадо-Ривес, Джулиан (январь 1996 г.). «Разработка и испытание полностью атомного силового поля OPLS по конформационной энергии и свойствам органических жидкостей». Журнал Американского химического общества . 118 (45): 11225–11236. DOI : 10.1021 / ja9621760 . ISSN 0002-7863 . 
  14. ^ Ван, Цзюньмэй; Вольф, Ромен М .; Колдуэлл, Джеймс У .; Коллман, Питер А .; Дело, Дэвид А. (2004-07-15). «Разработка и тестирование общего янтарного силового поля». Журнал вычислительной химии . 25 (9): 1157–1174. DOI : 10.1002 / jcc.20035 . ISSN 0192-8651 . PMID 15116359 . S2CID 18734898 .   
  15. ^ a b c Стефан, Саймон; Хорш, Мартин Т .; Врабец, Ядран; Хассе, Ганс (2019-07-03). «MolMod - открытая база данных силовых полей для молекулярного моделирования жидкостей» . Молекулярное моделирование . 45 (10): 806–814. arXiv : 1904.05206 . DOI : 10.1080 / 08927022.2019.1601191 . ISSN 0892-7022 . S2CID 119199372 .  
  16. ^ a b c Eggimann, Becky L .; Sunnarborg, Amara J .; Stern, Hudson D .; Блисс, Эндрю П .; Зипманн, Й. Илья (02.01.2014). «Онлайн-база данных параметров и свойств силового поля TraPPE» . Молекулярное моделирование . 40 (1–3): 101–105. DOI : 10.1080 / 08927022.2013.842994 . ISSN 0892-7022 . S2CID 95716947 .  
  17. ^ Чжэнь, Шу; Дэвис, GJ (16 августа 1983 г.). «Расчет n - m параметров потенциальной энергии Леннард-Джонса для металлов». Physica Status Solidi . 78 (2): 595–605. Bibcode : 1983PSSAR..78..595Z . DOI : 10.1002 / pssa.2210780226 .
  18. ^ Eisenschitz, R .; Лондон, Ф. (1930-07-01). "Über das Verhältnis der van der Waalsschen Kräfte zu den homöopolaren Bindungskräften" . Zeitschrift für Physik (на немецком языке). 60 (7): 491–527. Bibcode : 1930ZPhy ... 60..491E . DOI : 10.1007 / BF01341258 . ISSN 0044-3328 . S2CID 125644826 .  
  19. ^ Роулинсон, JS (2006-11-20). «Эволюция некоторых статистических механических идей» . Молекулярная физика . 104 (22–24): 3399–3410. Bibcode : 2006MolPh.104.3399R . DOI : 10.1080 / 00268970600965835 . ISSN 0026-8976 . S2CID 119942778 .  
  20. ^ Abascal, JLF; Вега, К. (15 декабря 2005 г.). «Универсальная модель конденсированных фаз воды: TIP4P / 2005» . Журнал химической физики . 123 (23): 234505. Bibcode : 2005JChPh.123w4505A . DOI : 10.1063 / 1.2121687 . ISSN 0021-9606 . PMID 16392929 .  
  21. ^ a b Ми, Густав (1903). "Zur kinetischen Theorie der einatomigen Körper" . Annalen der Physik (на немецком языке). 316 (8): 657–697. Bibcode : 1903AnP ... 316..657M . DOI : 10.1002 / andp.19033160802 .
  22. ^ Тан, KT; Тоеннис, Дж. Питер (1984-04-15). «Улучшенная простая модель для потенциала Ван-дер-Ваальса, основанная на универсальных функциях демпфирования для коэффициентов дисперсии» . Журнал химической физики . 80 (8): 3726–3741. Bibcode : 1984JChPh..80.3726T . DOI : 10.1063 / 1.447150 . ISSN 0021-9606 . 
  23. ^ а б Чжэнь, Шу; Дэвис, GJ (1983-08-16). «Расчет параметров потенциальной энергии Леннарда-Джонса – м для металлов» . Physica Status Solidi (A) (на немецком языке). 78 (2): 595–605. Bibcode : 1983PSSAR..78..595Z . DOI : 10.1002 / pssa.2210780226 .
  24. ^ a b Halicioglu, T .; Паунд, GM (1975-08-16). «Расчет параметров потенциальной энергии по свойствам кристаллического состояния» . Physica Status Solidi (А) . 30 (2): 619–623. Bibcode : 1975PSSAR..30..619H . DOI : 10.1002 / pssa.2210300223 .
  25. ^ Йоргенсен, Уильям Л .; Максвелл, Дэвид С .; Тирадо-Ривес, Джулиан (январь 1996 г.). «Разработка и испытание полностью атомного силового поля OPLS по конформационной энергии и свойствам органических жидкостей» . Журнал Американского химического общества . 118 (45): 11225–11236. DOI : 10.1021 / ja9621760 . ISSN 0002-7863 . 
  26. ^ Менделев, М.И.; Han, S .; Сроловиц, диджей; Ackland, GJ; Sun, DY; Аста, М. (2003). «Разработка новых межатомных потенциалов, подходящих для кристаллического и жидкого железа» . Философский журнал . 83 (35): 3977–3994. Bibcode : 2003PMag ... 83.3977A . DOI : 10.1080 / 14786430310001613264 . ISSN 1478-6435 . S2CID 4119718 .  
  27. ^ Леннард-Джонс, JE (1931). «Сплоченность». Труды физического общества . 43 (5): 461–482. Bibcode : 1931PPS .... 43..461L . DOI : 10.1088 / 0959-5309 / 43/5/301 .
  28. ^ a b c Аллен, Майкл П .; Тилдесли, Доминик Дж. (23.11.2017). «Компьютерное моделирование жидкостей». Оксфордская стипендия онлайн . DOI : 10.1093 / oso / 9780198803195.001.0001 . ISBN 9780198803195.
  29. ^ a b c d e f g h Шульц, Эндрю Дж .; Кофке, Дэвид А. (28.11.2018). «Комплексное высокоточное высокоточное уравнение состояния и свойств сосуществования для классических кристаллов Леннард-Джонса и низкотемпературных жидких фаз» . Журнал химической физики . 149 (20): 204508. DOI : 10,1063 / 1,5053714 . ISSN 0021-9606 . PMID 30501268 .  
  30. ^ a b Шульц, Эндрю Дж .; Кофке, Дэвид А. (07.08.2020). "Erratum:" Комплексное высокоточное уравнение состояния и сосуществования для классических кристаллов Леннард-Джонса и низкотемпературных жидких фаз "[J. Chem. Phys. 149, 204508 (2018)]" . Журнал химической физики . 153 (5): 059901. Bibcode : 2020JChPh.153e9901S . DOI : 10.1063 / 5.0021283 . ISSN 0021-9606 . PMID 32770918 .  
  31. ^ a b Шаппалс, Майкл; Мекленфельд, Андреас; Крегер, Лейф; Ботан, Виталий; Кёстер, Андреас; Стефан, Саймон; Гарсия, Эддер Дж .; Руткаи, Габор; Раабе, Габриэле; Кляйн, Питер; Леонхард, Кай (2017-09-12). «Круговое исследование: молекулярное моделирование термодинамических свойств на основе моделей с внутренними степенями свободы» . Журнал химической теории и вычислений . 13 (9): 4270–4280. DOI : 10.1021 / acs.jctc.7b00489 . ISSN 1549-9618 . PMID 28738147 .  
  32. ^ Loeffler, Hannes H .; Бозисио, Стефано; Дуарте Рамос Матос, Гильерме; Сух, Донхёк; Ру, Бенуа; Мобли, Дэвид Л .; Мишель, Жюльен (13.11.2018). «Воспроизводимость расчетов свободной энергии в различных пакетах программного обеспечения для молекулярного моделирования» . Журнал химической теории и вычислений . 14 (11): 5567–5582. DOI : 10.1021 / acs.jctc.8b00544 . ISSN 1549-9618 . PMID 30289712 .  
  33. ^ а б Ленхард, Йоханнес; Кюстер, Уве (2019). «Воспроизводимость и концепция численного решения» . Умы и машины . 29 (1): 19–36. DOI : 10.1007 / s11023-019-09492-9 . ISSN 0924-6495 . S2CID 59159685 .  
  34. ^ Ши, Вэй; Джонсон, Дж. Карл (15 сентября 2001 г.). «Изучение гистограммы повторного взвешивания и масштабирования конечного размера флюидов Леннарда – Джонса» . Равновесия жидкой фазы . 187–188: 171–191. DOI : 10.1016 / S0378-3812 (01) 00534-9 . ISSN 0378-3812 . 
  35. ^ a b Смит Б. (1992), «Фазовые диаграммы жидкостей Леннарда-Джонса» (PDF) , Журнал химической физики , 96 (11): 8639–8640, Bibcode : 1992JChPh..96.8639S , doi : 10.1063 / 1,462271
  36. ^ a b Метрополис, Николай; Розенблут, Арианна В .; Rosenbluth, Marshall N .; Teller, Augusta H .; Теллер, Эдвард (1953). «Уравнение состояний на быстрых вычислительных машинах» . Журнал химической физики . 21 (6): 1087–1092. Bibcode : 1953JChPh..21.1087M . DOI : 10.1063 / 1.1699114 . ISSN 0021-9606 . 
  37. ^ а б Кёстер, Андреас; Маусбах, Питер; Врабец, Ядран (10.10.2017). «Предплавление, равновесие твердой и жидкой фаз и термодинамические свойства в области высокой плотности на основе потенциала Леннарда-Джонса» . Журнал химической физики . 147 (14): 144502. Bibcode : 2017JChPh.147n4502K . DOI : 10.1063 / 1.4990667 . ISSN 0021-9606 . PMID 29031254 .  
  38. ^ a b c d e f g h я Стефан, Саймон; Дейтерс, Ульрих К. (20.08.2020). «Характеристические кривые жидкости Леннарда-Джонса» . Международный журнал теплофизики . 41 (10): 147. Bibcode : 2020IJT .... 41..147S . DOI : 10.1007 / s10765-020-02721-9 . ISSN 1572-9567 . PMC 7441092 . PMID 32863513 .   
  39. ^ a b Deiters, Ulrich K .; Ноймайер, Арнольд (2016-08-11). «Компьютерное моделирование характеристических кривых чистых жидкостей» . Журнал химических и технических данных . 61 (8): 2720–2728. DOI : 10.1021 / acs.jced.6b00133 . ISSN 0021-9568 . 
  40. ^ а б Агравал, Рупал; Кофке, Дэвид А. (1995). «Термодинамические и структурные свойства модельных систем при сосуществовании твердой и жидкой фаз: II. Плавление и сублимация системы Леннарда-Джонса» . Молекулярная физика . 85 (1): 43–59. DOI : 10.1080 / 00268979500100921 . ISSN 0026-8976 . 
  41. ^ Браун, EH (1960). «О термодинамических свойствах жидкостей». Bulletin de l'Institut International du Froid . Приложение 1960-1: 169–178.
  42. ^ Апфельбаум, EM; Воробьев, В.С. (18.06.2020). «Линия единичного коэффициента сжимаемости (зенолинии) для кристаллических состояний» . Журнал физической химии B . 124 (24): 5021–5027. DOI : 10.1021 / acs.jpcb.0c02749 . ISSN 1520-6106 . PMID 32437611 .  
  43. ^ Апфельбаум, EM; Воробьев, ВС; Мартынов, Г.А. (2008). «По поводу теории линии Зенона» . Журнал физической химии . 112 (26): 6042–6044. Bibcode : 2008JPCA..112.6042A . DOI : 10.1021 / jp802999z . ISSN 1089-5639 . PMID 18543889 .  
  44. ^ a b c d e Стефан, Саймон; Хассе, Ганс (2020-06-01). «Влияние дисперсионных дальнодействующих взаимодействий на свойства парожидкостных равновесий и границ раздела бинарных леннард-джонсовских смесей» . Молекулярная физика . 118 (9–10): e1699185. Bibcode : 2020MolPh.11899185S . DOI : 10.1080 / 00268976.2019.1699185 . ISSN 0026-8976 . S2CID 214174102 .  
  45. ^ Николас, JJ; Gubbins, KE; Streett, WB; Тилдесли, DJ (1979). «Уравнение состояния жидкости Леннарда-Джонса» . Молекулярная физика . 37 (5): 1429–1454. Bibcode : 1979MolPh..37.1429N . DOI : 10.1080 / 00268977900101051 . ISSN 0026-8976 . 
  46. ^ Фэн, Чао; Шульц, Эндрю Дж .; Чаудхари, Випин; Кофке, Дэвид А. (2015-07-28). «От восьмого до шестнадцатого вириальных коэффициентов модели Леннарда-Джонса» . Журнал химической физики . 143 (4): 044504. Bibcode : 2015JChPh.143d4504F . DOI : 10.1063 / 1.4927339 . ISSN 0021-9606 . PMID 26233142 .  
  47. ^ Шульц, Эндрю Дж .; Кофке, Дэвид А. (10 ноября 2009 г.). «Шестой, седьмой и восьмой вириальные коэффициенты модели Леннарда-Джонса» . Молекулярная физика . 107 (21): 2309–2318. Bibcode : 2009MolPh.107.2309S . DOI : 10.1080 / 00268970903267053 . ISSN 0026-8976 . S2CID 94811614 .  
  48. ^ Белл, Ян Х .; Мессерли, Ричард; Тол, Моника; Костильола, Лоренцо; Дайр, Джепп К. (25.07.2019). "Модифицированное масштабирование энтропии транспортных свойств жидкости Леннард-Джонса" . Журнал физической химии B . 123 (29): 6345–6363. DOI : 10.1021 / acs.jpcb.9b05808 . ISSN 1520-6106 . PMC 7147083 . PMID 31241958 .   
  49. ^ Lautenschlaeger, Мартин П .; Хассе, Ганс (2019). «Транспортные свойства леннард-джонсовской усеченной и сдвинутой жидкости из неравновесного молекулярно-динамического моделирования». Равновесия жидкой фазы . 482 : 38–47. DOI : 10.1016 / j.fluid.2018.10.019 .
  50. ^ a b c Травессет, Алекс (2014-10-28). «Фазовая диаграмма степенного закона и системы Леннарда-Джонса: кристаллические фазы» . Журнал химической физики . 141 (16): 164501. Bibcode : 2014JChPh.141p4501T . DOI : 10.1063 / 1.4898371 . ISSN 0021-9606 . PMID 25362319 .  
  51. ^ Хансен, Жан-Пьер; Верле, Лу (1969-08-05). «Фазовые переходы системы Леннарда-Джонса» . Физический обзор . 184 (1): 151–161. Bibcode : 1969PhRv..184..151H . DOI : 10.1103 / PhysRev.184.151 . ISSN 0031-899X . 
  52. ^ Лонг-Хиггинс, HC (1951-02-07). «Статистическая термодинамика многокомпонентных систем» . Труды Лондонского королевского общества. Серия А. Математические и физические науки . 205 (1081): 247–269. Bibcode : 1951RSPSA.205..247L . DOI : 10.1098 / rspa.1951.0028 . ISSN 0080-4630 . S2CID 202575459 .  
  53. ^ Леланд, TW; Роулинсон, Дж. С.; Sather, GA (1968). «Статистическая термодинамика смесей молекул разных размеров» . Труды общества Фарадея . 64 : 1447. DOI : 10.1039 / tf9686401447 . ISSN 0014-7672 . 
  54. ^ Мансури, Г. Али; Лиланд, Томас В. (1972). «Статистическая термодинамика смесей. Новая версия теории конформного раствора» . Журнал химического общества, Труды Фарадея 2 . 68 : 320. DOI : 10.1039 / f29726800320 . ISSN 0300-9238 . 
  55. ^ Rowlinson, JS; Суинтон, Флорида (1982). Жидкости и жидкие смеси (Третье изд.). Лондон: Баттерворт.
  56. ^ a b Deiters, Ulrich K .; Краска, Томас (2012). Фазовые равновесия жидкости высокого давления: феноменология и вычисления (1-е изд.). Амстердам: Эльзевир. ISBN 978-0-444-56354-5. OCLC  787847134 .
  57. ^ Шнабель, Торстен; Врабец, Ядран; Хассе, Ганс (2007). «В отличие от параметров Леннарда – Джонса для парожидкостного равновесия». Журнал молекулярных жидкостей . 135 (1–3): 170–178. arXiv : 0904.4436 . DOI : 10.1016 / j.molliq.2006.12.024 . S2CID 16111477 . 
  58. ^ Лоренц, HA (1881). "Ueber die Anwendung des Satzes vom Virial in der kinetischen Theorie der Gase" . Annalen der Physik (на немецком языке). 248 (1): 127–136. Bibcode : 1881AnP ... 248..127L . DOI : 10.1002 / andp.18812480110 .
  59. ^ ван Кониненбург, PH; Скотт, Р.Л. (1980-12-18). «Критические линии и фазовые равновесия в бинарных ван-дер-ваальсовых смесях» . Философские труды Лондонского королевского общества. Серия А, Математические и физические науки . 298 (1442): 495–540. Bibcode : 1980RSPTA.298..495K . DOI : 10,1098 / rsta.1980.0266 . ISSN 0080-4614 . S2CID 122538015 .  
  60. ^ Potoff, Джеффри Дж .; Панагиотопулос, Афанассиос З. (1998-12-22). «Критическая точка и фазовое поведение чистой жидкости и смеси Леннарда-Джонса» . Журнал химической физики . 109 (24): 10914–10920. Bibcode : 1998JChPh.10910914P . DOI : 10.1063 / 1.477787 . ISSN 0021-9606 . 
  61. ^ Проценко, Сергей П .; Байдаков, Владимир Г. (2016). «Бинарные смеси Леннарда-Джонса с сильно асимметричным взаимодействием компонентов. 1. Влияние энергетических параметров на фазовые равновесия и свойства границ раздела жидкость-газ». Равновесия жидкой фазы . 429 : 242–253. DOI : 10.1016 / j.fluid.2016.09.009 .
  62. ^ Проценко, Сергей П .; Байдаков, Владимир Г .; Брюханов, Василий М. (2016). «Бинарные смеси Леннарда-Джонса с сильно асимметричным взаимодействием компонентов. 2. Влияние размера частиц на фазовые равновесия и свойства границ раздела жидкость-газ». Равновесия жидкой фазы . 430 : 67–74. DOI : 10.1016 / j.fluid.2016.09.022 .
  63. ^ Стефан, Саймон; Хассе, Ганс (2020-01-23). «Молекулярные взаимодействия на границе раздела пар-жидкость: бинарные смеси простых жидкостей» . Physical Review E . 101 (1): 012802. Bibcode : 2020PhRvE.101a2802S . DOI : 10.1103 / PhysRevE.101.012802 . ISSN 2470-0045 . PMID 32069593 .  
  64. ^ а б Ламм, Моника Х .; Холл, Кэрол К. (2002). «Равновесия между твердой, жидкой и паровой фазами в бинарных смесях Леннарда – Джонса». Равновесия жидкой фазы . 194–197: 197–206. DOI : 10.1016 / S0378-3812 (01) 00650-1 .
  65. ^ а б Ламм, Моника Х .; Холл, Кэрол К. (2001). «Моделирование методом Монте-Карло полных фазовых диаграмм для бинарных смесей Леннарда – Джонса». Равновесия жидкой фазы . 182 (1–2): 37–46. DOI : 10.1016 / S0378-3812 (01) 00378-8 .
  66. ^ a b Хичкок, Моника Р .; Холл, Кэрол К. (1999-06-15). «Равновесие твердой и жидкой фаз для бинарных леннард-джонсовских смесей» . Журнал химической физики . 110 (23): 11433–11444. Bibcode : 1999JChPh.11011433H . DOI : 10.1063 / 1.479084 . ISSN 0021-9606 . 
  67. ^ Юнгблут, Swetlana; Деллаго, Кристоф (14 марта 2011 г.). «Кристаллизация бинарной смеси Леннард-Джонса» . Журнал химической физики . 134 (10): 104501. Bibcode : 2011JChPh.134j4501J . DOI : 10.1063 / 1.3556664 . ISSN 0021-9606 . PMID 21405169 .  
  68. ^ а б Ламм, Моника Х .; Холл, Кэрол К. (2004). «Влияние давления на полное фазовое поведение бинарных смесей» . Журнал Айше . 50 (1): 215–225. DOI : 10.1002 / aic.10020 . ISSN 0001-1541 . 
  69. ^ Коттерман, RL; Праусниц, JM (1986). «Молекулярная термодинамика жидкостей при низкой и высокой плотности. Часть II: Фазовые равновесия для смесей, содержащих компоненты с большими различиями в размерах молекул или потенциальной энергии» . Журнал Айше . 32 (11): 1799–1812. DOI : 10.1002 / aic.690321105 . ISSN 0001-1541 . 
  70. ^ Мюллер, Андреас; Винкельманн, Йохен; Фишер, Иоганн (1996). «Основа семейства уравнений состояния: 1. Неполярные и полярные чистые жидкости» . Журнал Айше . 42 (4): 1116–1126. DOI : 10.1002 / aic.690420423 . ISSN 0001-1541 . 
  71. ^ Weingerl, Ulrike; Вендланд, Мартин; Фишер, Иоганн; Мюллер, Андреас; Винкельманн, Йохен (2001). «Базовое семейство уравнений состояния: 2. Неполярные и полярные смеси жидкостей» . Журнал Айше . 47 (3): 705–717. DOI : 10.1002 / aic.690470317 .
  72. ^ a b Лафит, Томас; Апостолакова Анастасия; Авенданьо, Карлос; Галиндо, Ампаро; Adjiman, Claire S .; Мюллер, Эрих А .; Джексон, Джордж (2013-10-16). «Точная статистическая ассоциативная теория жидкости для цепных молекул, образованных из сегментов Ми» . Журнал химической физики . 139 (15): 154504. Bibcode : 2013JChPh.139o4504L . DOI : 10.1063 / 1.4819786 . hdl : 10044/1/12859 . ISSN 0021-9606 . PMID 24160524 .  
  73. ^ Blas, FJ; Вега, LF (1997). «Термодинамическое поведение гомоядерных и гетероядерных цепей Леннарда-Джонса с сайтами ассоциации из моделирования и теории» . Молекулярная физика . 92 (1): 135–150. Bibcode : 1997MolPh..92..135F . DOI : 10.1080 / 002689797170707 . ISSN 0026-8976 . 
  74. ^ Краска, Томас; Габбинс, Кейт Э. (1996). «Расчеты фазового равновесия с модифицированным уравнением состояния SAFT. 1. Чистые алканы, алканолы и вода» . Промышленные и инженерные химические исследования . 35 (12): 4727–4737. DOI : 10.1021 / ie9602320 . ISSN 0888-5885 . 
  75. ^ Ghonasgi, D .; Чепмен, Уолтер Г. (1994). «Прогнозирование свойств модельных полимерных растворов и смесей» . Журнал Айше . 40 (5): 878–887. DOI : 10.1002 / aic.690400514 . ISSN 0001-1541 . 
  76. ^ Mecke, M .; Мюллер, А .; Winkelmann, J .; Vrabec, J .; Фишер, Дж .; Span, R .; Вагнер, В. (1996-03-01). «Точное уравнение состояния типа Ван-дер-Ваальса для жидкости Леннарда-Джонса» . Международный журнал теплофизики . 17 (2): 391–404. Bibcode : 1996IJT .... 17..391M . DOI : 10.1007 / BF01443399 . ISSN 1572-9567 . S2CID 123304062 .  
  77. ^ a b c d e Колафа, Иржи; Незбеда, Иво (1994). «Жидкость Леннарда-Джонса: точное аналитическое и теоретически обоснованное уравнение состояния». Равновесия жидкой фазы . 100 : 1–34. DOI : 10.1016 / 0378-3812 (94) 80001-4 .
  78. ^ Тол, Моника; Руткаи, Габор; Кёстер, Андреас; Люстиг, Рольф; Спан, Роланд; Врабец, Ядран (2016). "Уравнение состояния жидкости Леннарда-Джонса" . Журнал физических и химических справочных данных . 45 (2): 023101. Bibcode : 2016JPCRD..45b3101T . DOI : 10.1063 / 1.4945000 . ISSN 0047-2689 . 
  79. ^ Готшальк, Маттиас (2019-12-01). «УСЛ для жидкости Леннарда-Джонса: подход вириального расширения» . Продвижение AIP . 9 (12): 125206. Bibcode : 2019AIPA .... 9l5206G . DOI : 10.1063 / 1.5119761 . ISSN 2158-3226 . 
  80. ^ a b Джонсон, Дж. Карл; Zollweg, John A .; Губбинс, Кейт Э. (1993-02-20). "Уравнение состояния Леннарда-Джонса повторно" . Молекулярная физика . 78 (3): 591–618. Bibcode : 1993MolPh..78..591J . DOI : 10.1080 / 00268979300100411 . ISSN 0026-8976 . 
  81. ^ Врабец, Ядран; Кедия, Гаурав Кумар; Fuchs, Guido; Хассе, Ганс (10 мая 2006 г.). «Комплексное исследование сосуществования пара и жидкости в усеченной и смещенной жидкости Леннарда-Джонса, включая свойства плоской и сферической границы раздела» . Молекулярная физика . 104 (9): 1509–1527. Bibcode : 2006MolPh.104.1509V . DOI : 10.1080 / 00268970600556774 . ISSN 0026-8976 . S2CID 96606562 .  
  82. ^ Хейер, Микаэла; Стефан, Саймон; Лю, Цзиньлу; Чепмен, Уолтер Дж .; Хассе, Ганс; Лангенбах, Кай (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 .  
  83. ^ Шауль, Кэтрин RS; Шульц, Эндрю Дж .; Кофке, Дэвид А. (2010). «Влияние усечения и сдвига на вириальные коэффициенты потенциалов Леннарда – Джонса» . Сборник чехословацких химических сообщений . 75 (4): 447–462. DOI : 10.1135 / cccc2009113 . ISSN 1212-6950 . 
  84. ^ Ши, Вэй; Джонсон, Дж Карл (2001). «Изучение гистограммы повторного взвешивания и масштабирования конечного размера флюидов Леннарда – Джонса». Равновесия жидкой фазы . 187–188: 171–191. DOI : 10.1016 / S0378-3812 (01) 00534-9 .
  85. ^ а б Дуников Д.О .; Малышенко, ИП; Жаховский, В.В. (2001-10-08). "Соответствующие законы состояний и моделирование молекулярной динамики жидкости Леннарда-Джонса" . Журнал химической физики . 115 (14): 6623–6631. Bibcode : 2001JChPh.115.6623D . DOI : 10.1063 / 1.1396674 . ISSN 0021-9606 . 
  86. ^ Чипев, Никола; Секлер, Штеффен; Хайнен, Маттиас; Врабец, Ядран; Гратль, Фабио; Хорш, Мартин; Бернройтер, Мартин; Гласс, Колин В; Нитхаммер, Кристоф; Молоток, Николай; Крищок, Бернд (2019). «TweTriS: Моделирование двадцати триллионов атомов» . Международный журнал приложений высокопроизводительных вычислений . 33 (5): 838–854. DOI : 10.1177 / 1094342018819741 . ISSN 1094-3420 . S2CID 59345875 .  
  87. ^ Стефан, Саймон; Лю, Цзиньлу; Лангенбах, Кай; Чепмен, Уолтер Дж .; Хассе, Ганс (2018). "Граница раздела пар-жидкость усеченной и сдвинутой по Леннарду-Джонсу жидкости: сравнение молекулярного моделирования, теории градиента плотности и функциональной теории плотности" . Журнал физической химии C . 122 (43): 24705–24715. DOI : 10.1021 / acs.jpcc.8b06332 . ISSN 1932-7447 . 
  88. ^ Grüneisen, Эдвард (1911). "Das Verhältnis der thermischen Ausdehnung zur spezifischen Wärme fester Elemente" . Zeitschrift für Elektrochemie und angewandte Physikalische Chemie . 17 : 737–739. doi : 10.1002 / bbpc.191100004 (неактивный 2021-01-10).CS1 maint: DOI inactive as of January 2021 (link)
  89. ^ Грюнайзена, Е. (1912). "Theorie des festen Zustandes einatomiger Elemente" . Annalen der Physik (на немецком языке). 344 (12): 257–306. Bibcode : 1912AnP ... 344..257G . DOI : 10.1002 / andp.19123441202 .
  90. ^ Штокмайер, WH (1941-05-01). «Вторые вириальные коэффициенты полярных газов» . Журнал химической физики . 9 (5): 398–402. Bibcode : 1941JChPh ... 9..398S . DOI : 10.1063 / 1.1750922 . ISSN 0021-9606 . 

Внешние ссылки [ править ]

  • Модель Леннарда-Джонса на SklogWiki .