CC BY
34ПРИКЛАДНАЯ МАТЕМАТИКА И МЕХАНИКА
Вестник Сыктывкарского университета. Серия 1: Математика. Механика. Информатика. Выпуск 4 (25). 2017
УДК 548.52
ДИНАМИКА СЕТКИ МЕЖМОЛЕКУЛЯРНЫХ СВЯЗЕЙ И ФАЗОВЫЕ ПЕРЕХОДЫ В КОНДЕНСИРОВАННЫХ
СРЕДАХ
В. Н. Чередов, Л. А. Куратова
Предложен новый подход к исследованию молекулярной структуры жидкой и твердой фазы вещества — модель мерцающих связей. Данный подход основывается на развитии модели тепловых колебаний атомов (молекул) вещества и их влиянии на динамику молекулярной структуры и структуру сетки межмолекулярных связей твердой и жидкой фаз вещества.
Выявлена температурная зависимость динамики свойств сетки межмолекулярных связей твердой и жидкой фаз вещества, а также динамики свойств указанной сетки связей в фазовых переходах первого рода «твердое тело — жидкость» и «жидкость — газ». На основе построенной модели изучена динамика структуры Н20 и ее фазовых переходов.
Ключевые слова: межмолекулярные связи, фазовые переходы, кристаллизация, структура решетки.
1. Введение
Исследования организации вещества в молекулярные (атомные) структуры неизменно находятся в центре внимания физиков и химиков. В последнее время особое значение данные исследования получили в связи с повышенным вниманием науки к наноструктурам. В последние годы развитие получили оригинальные модели процессов организации вещества на наноуровне, которые используются для решения широкого круга физических проблем.
В данной работе в развитие теории самооорганизации вещества на наноуровне предлагается модель мерцающих связей для исследования
© Чередов В. Н., Куратова Л. А., 2017.
структуры межмолекулярных связей кристаллической и жидкой фаз вещества, фазовых переходов 1-го рода. Данная модель существенно расширяет наши возможности анализа фазовых переходов, понимания природы и структуры конденсированного состояния вещества.
2. Модель мерцающих связей
Рассмотрим однородный участок твердой или жидкой фазы в гипотетическом состоянии, когда все связи между молекулами (атомами) являются устойчивыми. Пусть данный участок вещества получает некоторый небольшой объем тепла АQ. Тогда внутренняя энергия молекул увеличивается на такую же величину АЕ = АQ. Вполне естественно, что увеличение внутренней энергии молекул может привести к разрыву некоторых межмолекулярных связей.
Природа разрыва связей между молекулами (атомами) зависит от типа этих связей. В общем виде разрыв связи можно представить как временное перераспределение электронной плотности у атомов, образующих молекулы и участвующих в связи, в состояние отсутствия связи. Например, для водородной и вандервальсовой связи разрыв связи между молекулами означает временное изменение электронной плотности, при котором одна из молекул, участвующих в связи, теряет на время дипольные свойства [1].
Может реализоваться ситуация, когда часть связей между молекулами в структуре вещества находится в устойчивом состоянии, а часть времени в разорванном состоянии. Если вещество находится в состоянии с постоянной температурой и давлением, т. е. с одинаковой суммарной внутренней энергией молекул, то часть межмолекулярных связей всегда находится в разорванном состоянии. Разорванные связи восстанавливаются, переходят в устойчивое состояние, а другие устойчивые связи в это же время рвутся. Происходит как бы «мерцание» или «осцилляция» разрывов связей внутри объема вещества.
Применяя принцип однородности пространства, приходим к понятию мерцающей связи. Распределение разорванных связей в пространстве и во времени должно быть равномерно. То есть каждая межмолекулярная связь во всем объеме вещества должна мерцать, время от времени разрываться и потом опять возникать.
При дальнейшем увеличении внутренней энергии молекул происходит увеличение количества разорванных связей молекулы Н и, соответственно, уменьшение количества устойчивых связей молекулы Г (Г + Н = Кт, где Кт — координационное число молекулы). То есть суммарная внутренняя энергия молекул и температура вещества определяются в данной модели количеством устойчивых межмолекулярных
связей молекулы вещества.
Разрыв связи, который произошел в одном месте решетки, прекращается, связь восстанавливается, и одновременно происходит разрыв связи в другом месте решетки, и т. д. Конкретная межмолекулярная связь находится в двух состояниях. Часть времени Tf межмолекулярная связь устойчива, а часть времени т9 неустойчива или разорвана, при этом время мерцания связи тт = Tf + т9. В данной формулировке все межмолекулярные связи мерцающие и все реально существующие, т. е. в идеале их количество в конкретный момент времени определяется эмпирическим путем. Тем не менее, внутренняя энергия вещества в конкретный момент времени определяется исключительно количеством устойчивых связей Г = KmTf /тт или в силу однородности вещества количеством устойчивых связей одной молекулы Н = Кттд/тт, где Нт — общее количество связей вещества, Кт — количество связей одной молекулы вещества.
Рассмотрим ряд условий и обстоятельств, характеризующих модель мерцающих связей. Во-первых, термодинамическое состояние вещества определяется в каждый момент времени исключительно соотношением количества устойчивых и разорванных связей между молекулами (атомами). Данное условие исключает рассмотрение случаев, при которых другие факторы могут изменять термодинамическое состояние вещества.
Рис. 1. Мгновенная структура твердого тела с простой кубической кристаллической решеткой с мерцающими связями (разорванные связи показаны
разорванной линией)
Во-вторых, в течение какого времени количество и, соответственно, доля разорванных связей в общей доле межмолекулярных связей в одном термодинамическом состоянии сохраняется? В данном представлении термодинамическое состояние твердого тела или внутренняя энергия тела определяется мгновенным «снимком» всех связей структуры: разорванных и устойчивых (рис. 1). В каждый последующий момент времени мгновенная структура вещества будет постоянно меняться.
В-третьих, связи между молекулами (атомами), как устойчивые, так и разорванные, распределены по всему объему вещества равномерно. Данное условие ограничивает исследуемые объекты модели мерцающих (осциллирующих) связей. Фактически в данной модели рассматриваются объекты вещества, весь объем которых находится в одном термодинамическом состоянии, исключая неоднородности как самого вещества, так и его состояния. Другими словами, если весь объект нагревается (или охлаждается), то нагреваются (или охлаждаются) одинаково все части объекта. Если объект переходит из одной фазы вещества в другую, то каждая его часть одновременно переходит из той же фазы в другую. Это ограничивает объекты рассмотрения до участков одной фазы при фазовых переходах или до размеров толщины границы фазового перехода (фронта кристаллизации), т. е. до наноразмеров.
В-четвертых, на поверхности твердой и жидкой фазы вещества количество разорванных связей молекул много больше, чем в объеме вещества. Соответственно, количество устойчивых связей молекул в поверхностном слое вещества меньше, чем у молекул объема. Это приводит к различию термодинамических свойств вещества на поверхности и в объеме. В частности, температура плавления тонких пленок должна быть меньше, чем в объеме вещества. Твердое тело должно начинать плавится на своей поверхности.
В-пятых, в твердом агрегатном состоянии вещества в среднем молекулы (атом) имеют не менее 3-х устойчивых межмолекулярных связей. В противном случае происходит фазовый переход из твердого агрегатного состояния вещества в жидкое и наоборот. В общем случае, обратное утверждение неверно. Учитывая, что равновесное (без учета тепловых колебаний) положение молекул в твердом теле (кристалле) устойчиво, следует предположить, что количество устойчивых межмолекулярных связей одной молекулы в твердом теле должно быть в среднем не менее 3-х. Так как в трехмерном пространстве существует три поступательные степени свободы, то 3-х неразорванных межмолекулярных связей для одной молекулы вполне достаточно для ее устойчивого положения в структуре твердого тела и сохранения структурой устойчивой решетки.
Таким образом, ограничение снизу на среднее количество устойчивых связей молекул твердой фазы вещества: Г > 3.
Для наглядности рассмотрим решетку твердого тела с однородной структурой в форме прямоугольного параллелепипеда с количеством связей у молекулы к, I, т по 3-м пространственным осям х, у, г. Общее количество межмолекулярных связей N = 3 * к * I * т. На рис. 2 представлен граф устойчивых связей, являющийся гипотетическим состоянием решетки твердого тела, при котором количество межмолекулярных связей у каждой молекулы равно 3-м. У молекул, количество связей которых меньше 3-х, отсутствующие связи не изображены, так как они соединяют данные молекулы с молекулами за пределами данного участка вещества. Подчеркнем, что данное состояние решетки нам необходимо исключительно для наглядности и представляет собой один из практически бесконечного множества вариантов пространственной решетки распределения устойчивых связей.
Рис. 2. Граф одномоментных устойчивых связей в простой кубической решетке твердого тела (размерами к * I * т), у которого каждая молекула имеет 3 межмолекулярные связи (изображены только устойчивые связи)
В-шестых, в жидком агрегатном состоянии вещества в среднем молекулы (атомы) имеют не менее 2-х устойчивых межмолекулярных связей: Г > 2. В случае снижения количества устойчивых межмолекулярных связей до меньшего числа происходит фазовый переход из жидкого агрегатного состояния вещества в газообразное и наоборот.
На рис. 3 представлено гипотетическое состояние решетки жидкой
фазы, при котором количество межмолекулярных связей у каждой молекулы равно 2-м. Фактически получилась замкнутая одномерная цепочка молекул. Подчеркнем, что данное состояние решетки нам необходимо исключительно для наглядности и оно реально невозможно в связи с необходимостью соблюдения принципа однородности пространства, а также в данном случае и в связи с невозможностью данного состояния из-за эффекта поверхностного натяжения жидкости. Тем не менее, при перераспределении устойчивых межмолекулярных связей однородно по всей трехмерной структуре общее их количество не изменится, а вышеперечисленные соображения удовлетворяются (см. рис.).
Рис. 3. Гипотетическое трансформированное состояние решетки жидкой фазы вещества, у которого каждая молекула имеет две устойчивые межмоле-кулярныв связи (изображены только устойчивые связи)
В зависимости от того, сколько связей молекулы в данный момент времени находится в устойчивом или в разорванном состоянии, молекула занимает более или менее возбужденное равновесное состояние (положение в структуре), характеризующееся, соответственно, большим или меньшим значением внутренней энергии молекулы. При этом под воздействием оставшихся связей молекул-соседей сама молекула немного смещается в пространстве. При восстановлении связи, т. е. при приобретении связью устойчивого характера, молекула опять занимает свое исходное пространственное положение и переходит в исходное менее возбужденное состояние с меньшей энергией.
То же самое происходит при «мерцании» всех остальных связей данной молекулы. При этом следует учитывать, что молекулы, с которыми выбранная молекула имеет связи, сами по себе также колеблются около положения равновесия при «мерцании» своих связей. Это придает колебаниям молекулы хаотический характер со степенями свободы, соответствующими степеням свободы тепловых колебаний.
Описанный механизм колебания молекул и является в модели мер-
цающих связей основным механизмом теплового колебания. При этом в данной модели первично не тепловое колебание молекул вещества, а «мерцание» связей молекул вещества. Модель мерцающих связей дает дополнительные возможности компьютерного моделирования динамики теплового колебания молекул вещества.
Осцилляция связей молекул жидкой и твердой фазы вещества имеет тепловую природу и определяется температурой вещества. В связи с этим мы предполагаем, что количество мерцающих связей в разорванном состоянии в общем количестве связей постоянно и определяется прежде всего температурой Т и, соответственно, внутренней энергией молекулы Ег. Чем больше Т, тем больше Н и тем меньше Г.
Молекулам вещества, находящегося при Т = 0 К, чтобы освободится от межмолекулярных связей, необходимо придать энергию, называемую энергией сублимации Е3 и равную:
г Е
Е3 = 2 кв Ть + + AQb = Кт^, (1)
где г — число степеней свободы сублимированной молекулы, к в — постоянная Больцмана, Ть — температура кипения, AQm — теплота плавления из расчета на одну молекулу, AQь — теплота парообразования из расчета на одну молекулу, Кт — координационное число молекулы в твердой фазе, Ет — энергия разрыва одной связи при Т = 0 К.
Формула (1) показывает, какую энергию вещество в твердой фазе при абсолютном нуле температуры должно получить, чтобы все молекулы стали свободными и молекула из твердой фазы перешла в газообразную, т. е. стала свободна от межмолекулярных связей. Данная энергия сублимации равняется энергии связи молекулы при нулевой температуре. Деление на 2 в правой части формулы (1) возникает вследствие того, что разрыв одной связи относится сразу к двум молекулам. По формуле (1) можно рассчитать энергию разрыва одной связи молекулы.
Энергия связи молекулы Еь = ГЕт = (Кт — Н)Ет будет определяться формулой:
ЕЕ Еь = Г= Е3 — Ег = Кт 2"" — Ег. (2)
Для твердой фазы внутренняя энергия молекул определяется энергией их тепловых колебаний (энергией образования и динамики дефектов кристаллической решетки [2] пренебрегаем, так как она значительно меньше последней).
Как уже отмечалось ранее, мерцание межмолекулярных связей приводит к тепловым колебаниям молекул, а последние описываются в
термодинамике моделями Дюлонга и Пти, моделью Энштейна и самой близкой к реальности моделью Дебая [3]. Поэтому для определения энергии связи молекул в условиях мерцания межмолекулярных связей, которая равна в модели мерцающих связей энергии тепловых колебаний молекул, используем модель Дебая [3].
Тогда энергия связи молекул твердой фазы определяется следующими формулами:
e E / T \ 3 г©D/T x3 dx
= F-'-f. = Km-f - T(_) I ex-i, <3>
где Fsp — количество устойчивых связей молекулы в твердой фазе при температуре T, OD — температура Дебая.
При температуре плавления T = Tm формула (3) позволяет рассчитать количество устойчивых связей молекулы в твердой фазе при температуре плавления Fm из выражения для энергии связи:
1 Em, т ^ Em, ~ * ^ ( Tm^\ f © / x dx
Eb = F^ = Km^- - 9ksTm[^) I (4)
b m 2 m 2 B m\ On L ex — 1' W
где Е^— количество устойчивых связей молекулы в твердой фазе при температуре плавления Т = Тм.
Подставляя в формулу (4) константы и известные данные, получаем количество устойчивых связей молекулы в твердой фазе при температуре плавления Е£!р. При температуре плавления в твердой фазе у молекулы имеется Е£Р устойчивых связей, а значит, Кт — Е^р связей разорвано.
В жидкой фазе при температуре плавления (кристаллизации) необходимо вычесть из выражения энергии связи в формуле (4) теплоту плавления, приведенную на одну молекулу. Соответственно, получаем выражение для количества устойчивых связей молекулы в жидкой фазе при температуре плавления Е^:
Е Е / Т \ 3 Г/Тт ]х
ЕЬ = Е1РЕт = КтЕт — 9квТт I еХ—1 — А®т> (5)
где Е}р — количество устойчивых связей молекулы в жидкой фазе при температуре плавления Т = Т .
Для газообразной фазы необходимо учесть теплоту парообразования, приведенную на одну молекулу. Для этого подставим формулу (1)
в формулу (2). Получаем формулу энергии связи молекулы газообразной фазы при температуре кипения Т = ТЬ:
Е Е ?
Еь = Е?^ = - 2квТь - - AQb = 0, (6)
где Е^р — количество устойчивых связей молекулы в газообразной фазе при температуре кипения Т = Ть (Е^р = 0). В жидкой фазе при температуре кипения:
Еь = Е1ьрЕт = КтЕт - 2квТь - AQm = AQь, (7)
где Е1ЬР — количество устойчивых связей молекулы в жидкой фазе при температуре кипения Т = Ть.
В жидкой фазе для Тт < Т < Ть количество устойчивых связей молекулы меняется в следующем диапазоне ЕьР < Е1р < Е^, причем при увеличении температуры количество связей падает. Так как в термодинамике на сегодняшний день отсутствует полноценная модель внутренней энергии молекул жидкости, поэтому аппроксимируем функцию внутренней энергии жидкой фазы от температуры линейной зависимостью, тем более что в реальных жидкостях, например воде, она близка к линейной. Тогда для количества устойчивых связей молекулы получим следующую формулу:
Е = Е1рЕт = Ет I е1Р - Т - Тт\ (8)
Еь = Е 2 = 2 ^ Ав'р ] ' (8)
где Е1р — количество устойчивых связей молекулы в жидкой фазе при температуре Тт < Т < Ть, Ав1рН — величина характеристической разности температур жидкой фазы вещества (Тт + Ав^) - Тт.
В формуле (8) вводится величина характеристической разности температур жидкой фазы вещества, которая определяется по следующему соотношению:
Ть - Тт Е1Р — Ед
Вт - Вт
Характеристическая разность температур жидкой фазы вещества Ав^ является разностью температур, при увеличении на которую температуры жидкости от температуры плавления, должна разорваться дополнительно ровно одна межмолекулярная связь молекулы. Данная характеристическая разность температур также вводится в настоящей работе и является уникальным параметром для каждого вещества.
Ав1р = —_— (9)
сН тЯ,р трдр' К '
В газообразной фазе количество среднее устойчивых связей молекулы можно считать равным нулю Едр = 0.
Сводя воедино все участки количества устойчивых связей молекулы Е воедино, т. е. формулы (3) и (8), получаем:
F
F sp
- 9Emквт (£)3 /с
©D /T x3dx
ex-l '
при 0 < T < Tm (твердая фаза); F1P = Flp — T -T
A©
ip
(10)
при Tm < T < Tb (жидкая фаза); Fgp = 0,
при T < Tb (газообразная фаза).
Для удобства и наглядности введем два показателя: долю устойчивых связей молекул f = F/Km и долю разорванных связей молекул h = H/Km. Координационное число соответствует координационному числу молекул в соответствующей фазе, т. е. для молекул твердой фазы Km, в жидкой фазе Klm (может быть как больше, так и меньше Km), в газообразной фазе K®p = 0. Область изменения показателей f и h определена от 0 до 1. При этом всегда f + h = 1.
3. Структура воды в модели мерцающих связей Отобразим полученные выше результаты на графике температурной зависимости показателей F из сводной формулы 10 для H2O (рис. 4).
Рис. 4. Фазовые переходы и температурная зависимость количества и доли устойчивых связей молекулы в решетке Н2О
Данные по воде брались из работ [4-6]. Расчеты по формулам 4, 5, 7 дают для воды следующие значения количества устойчивых связей для льда при температуре плавления Етр = 3, 64. При температуре плавления во льду у молекулы имеется 3,64 устойчивых связей, а значит, 0,36 связей разорвано. В жидкой фазе при температуре плавления (кристаллизации) для воды получаем Е]р = 3, 22. В жидкой фазе при температуре кипения для воды получается значение ЕьР = 2, 84.
4. Заключение
Впервые представлена модель мерцающих связей, в рамках которой построена модель структуры кристалла и жидкости как непрерывной сетки межмолекулярных связей однородной по количеству устойчивых межмолекулярных связей.
Выявлена температурная зависимость количества устойчивых связей молекул для твердой и жидкой фазы, его изменения при фазовых переходах первого рода «твердое тело — жидкость» и «жидкость — газ». На примере воды показана зависимость количества устойчивых межмолекулярных связей в зависимости от температуры и фазы вещества, ее изменение в процессе фазовых переходов.
Список литературы
1. Каплан И. Г. Межмолекулярные взаимодействия. Физическая интерпретация, компьютерные расчёты и модельный потенциал. М.: БИНОМ. Лаборатория знаний, 2012. 400 с.
2. Чередов В. Н. Статика и динамика дефектов в синтетических кристаллах флюорита. СПб.: Наука, 1993. 112 с.
3. Ландау Л. Д., Лифшиц Е. М. Статистическая физика. М.: Физ-матлит, 2010. Ч. 1. 616 с.
4. Енохович А. С. Справочник по физике и технике. М.: Просвещение, 1989. 224 с.
5. Зацепина Г. Н. Физические свойства и структура воды. М.: МГУ, 1998. 184 с.
6. Эйзенберг Д., Кауцман В. Структура и свойства воды. М.: Директ-медиа, 2012. 284 с.
Summary
Cheredov V. N., Kuratova L. A. Dynamics of a network of intermolecular bonds and phase transitions in condensed media
A new approach to the investigation of the molecular structure of the liquid and solid phases of matter — the model of flickering bonds — is proposed. This approach is based on the development of the model of thermal vibrations of atoms (molecules) of a matter and their effect on the dynamics of the molecular structure and the structure of the intermolecular bond network of the solid and liquid phases of matter.
The temperature dependence of the dynamics of the properties of the network of intermolecular bonds of the solid and liquid phases of matter, as well as the dynamics of the properties of this bond network in the first-order phase transitions «solid-liquid» and «liquid-gas» is revealed. On the basis of the constructed model, the dynamics of the structure of H2O and its phase transitions is studied.
Keywords: intermolecular bonds, phase transitions, crystallization, lattice structure.
References
1. Kaplan I. G. Mezhmolekuljarnye vzaimodejstvija. Fizicheskaja inter-pretacija, komp'juternye raschjoty i model'nye potencial (Intermolecular interactions. Physical interpretation, computer calculations and model potentials), Moscow: BINOM, Laboratory of Knowledge, 2012, 400 p.
2. Cheredov V. N. Statika i dinamika defektov v sinteticheskih kristal-lah fljuorita (Statics and dynamics of defects in synthetic fluorite crystals), Saint-Petersburg: Nauka, 1993, 112 p.
3. Landau L. D., Lifshitz E. M. Statisticheskaja fizika (Statistical physics), part 1, Moscow: Fizmatlit, 2010, 616 p.
4. Enochovich A. S. Spravochnik po fizike i tehnike (Reference book on physics and techniques), Moscow: Prosveshenie, 1989, 224 p.
5. Zatsepina G. N. Fizicheskie svojstva i struktura vody (Physical properties and structure of water), Moscow: Moscow State University, 1998, 184 p.
6. Eisenberg D., Kautsman V. Struktura i svojstva vody (Structure and properties of water), Moscow: Direct-Media, 2012, 284 p.
Для цитирования: Чередов В. Н., Куратова Л. А. Динамика сетки межмолекулярных связей и фазовые переходы в конденсированных средах // Вестник Сыктывкарского университета. Сер. 1: Математика. Механика. Информатика. 2017. Вып. 4 (25). C. 20-32.
For citation: Cheredov V. N., Kuratova L. A. Dynamics of a network of intermolecular bonds and phase transitions in condensed media, Bulletin of Syktyvkar University, Series 1: Mathematics. Mechanics. Informatics, 2017, №4 (25), pp. 20-32.
СГУ им. Питирима Сорокина
Поступила 19.12.2017
