Description

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

Одним из современных методов исследования сильно неоднородных систем является метод функционала молекулярной плотности, который в разных вариантах в последнее время применялся для описания поверхностных слоев как на границе твердое тело–жидкость, так и на границе жидкость–пар. Другими широко используемыми методами является молекулярно-динамическое моделирование и метод Монте-Карло. Кроме того, в последнее время эффективно развивается вычислительный метод подталкивания упругой ленты (Nudged Elastic Band, NEB) в применении к задачам, где требуется анализ энергетической гиперповерхности в пространстве высокой или бесконечной размерности. Этот численный метод позволяет найти путь минимального перепада энергии (ПМПЭ или Minimum Energy Path, MEP в англоязычной литературе), который «проходит» система от исходного нестабильного или метастабильного состояния к конечному, стабильному. Этот метод применим к системам с различной геометрией.

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

Key findings for the stage (in detail)

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

Рассмотрены различные уравнения состояния системы жестких сфер для описания однокомпонентного флюида: уравнения Карнахана–Старлинга, Перкуса–Йевика, Колафы, уравнения Русанова 6-го и 7-го порядков. В качестве «эталонного» использовано модифицированное уравнение Колафы–Лабика–Малиевского, описывающее с высокой точностью, но в ограниченном диапазоне плотностей результаты молекулярного моделирования системы жестких сфер. Проведено сравнение того, насколько точно описываются зависимости фактора сжимаемости, химического потенциала и плотности свободной энергии с помощью этих уравнений. Для уравнений Карнахана–Старлинга, Колафы, усеченного уравнения Русанова 6-го порядка, а также модифицированного уравнения Колафы–Лабика–Малиевского, в рамках интегрального метода функционала плотности в приближении случайной фазы получены профили молекулярной плотности в радиально неоднородных сферических малых капельках и пузырьках аргоноподобного вещества, являющихся критическими каплями/пузырьками при нуклеации. Для описания флюида все использованные уравнения состояния дополнены слагаемым, описывающим притяжение частиц в среднеполевом приближении. При сравнении уравнений состояния системы жестких сфер для описания как однородного флюида, так и критических капель/пузырьков уравнение Колафы оказалось лучшим из рассмотренных уравнений состояния системы жестких сфер, с высокой точностью повторяющим результаты для «эталонного» уравнения Колафы–Лабика–Малиевского. При этом уравнение Колафы имеет намного более простой вид и напоминает широко используемое уравнение Карнахана–Старлинга, заметно превосходя его по точности при описании как однородных, так и неоднородных флюидов жестких сфер.

Проведено молекулярно-динамическое моделирование образования сидячей капели на плоской поверхности с вычислением ее краевого угла. Реализовано два способа вычисления краевого угла, отличающихся процедурой кластеризации. Рассмотрены капли аргона и капли воды в полноатомной модели TIP3P. В качестве плоской поверхности была использована непроницаемая стенка, взаимодействующая с каплей в рамках потенциала Стила (потенциал «3-9»), а также полноатомные модели графита и кристобалита (SiO₂). Для сидячих капель, состоящих из молекул аргона, была определена зависимость краевого угла от параметра взаимодействия потенциала Стила, при помощи которого регулировалась смачиваемость поверхности. Наблюдалось уменьшение краевого угла капли при увеличении параметра взаимодействия. Была определена зависимость краевого угла сидячей капли аргона от температуры при постоянном параметре взаимодействия молекул аргона с поверхностью. Наблюдалось уменьшение краевого угла смачивания при возрастании температуры от 65 К до 100 К, сопровождающееся увеличением числа молекул аргона, находящихся в газовой фазе. Проведено моделирование сидячих капель воды на поверхности графита. Для моделирования различной степени смачиваемости поверхности варьировался параметр ε потенциала Леннард-Джонса для атомов, составляющих поверхность. Наблюдалось уменьшение краевого угла капли воды при увеличении силы взаимодействия молекул воды с поверхностью. Было проверено, что изменение числа слоёв графита в диапазоне от 4 до 10 не оказывает существенного влияния на краевой угол капли. Для определения зависимости угла смачивания от температуры была проведена серия расчётов с каплей воды из 1000 молекул на графитовой поверхности при температурах от 280 до 360 К. Краевой угол несколько уменьшился при увеличении температуры в рассмотренном диапазоне, при этом его изменение осталось в пределах статистической погрешности. Для исследования зависимости краевого угла сидячей капли воды от ее размера были рассмотрены системы, содержащие 500, 1000, 2000 и 8000 молекул воды на графитовой поверхности, при этом большинство молекул воды принадлежали сидячей капле и лишь небольшое их число оказалось в газовой фазе. Результаты расчётов показали, что краевой угол оказался практически неизменным в рассмотренном диапазоне размеров капель и составил, приблизительно, 52 градуса. Также было проведено моделирование сидячей капли воды на поверхности кристобалита (SiO₂) в полноатомной модели. Смачиваемость поверхности варьировалась путем присоединения к атомам кислорода, находящимся на поверхности, атомов водорода (с образованием OH групп), при этом частичные заряды всех атомов поверхности пересчитывались для сохранения электронейтральности системы. При отсутствии OH групп на поверхности краевой угол образовавшейся капли составил, приблизительно, 88 градусов. В системе с максимально возможной плотностью OH групп наблюдалась полная смачиваемость поверхности водой.

Рассмотрена теория эффекта Казимира для двух дифракционных решеток, разделенных вакуумной щелью, обсуждается нарушение симметрии при геометрических переходах. Следствием нарушения трансляционной инвариантности при вращении бесконечной дифракционной решетки на конечный угол является нарушение симметрии для векторов обратного векторного пространства, что приводит к несохранению компоненты волнового вектора обратного векторного пространства и, как следствие, к фундаментальной перестройке структуры обратного векторного пространства. В результате перестройки обратного векторного пространства в системе двух дифракционных решеток с совпадающими периодами имеет место геометрический переход – разрыв энергии Казимира. В системе из двух конечных дифракционных решеток крутящий момент принимает сколь угодно большие значения при одновременном увеличении размеров дифракционных решеток. Предсказанный эффект может представлять значительный интерес в связи с новым механизмом нарушения симметрии и возможностью проявления данного механизма в различных физических системах. Разработан новый, калибровочно-инвариантный по построению, метод вывода электрических и магнитных функций Грина в теории эффекта Казимира. Функции Грина находятся решением системы уравнений, полученной из граничных условий для тангенциальных компонент электрического и магнитного полей на двух границах диэлектрик-вакуум и условий поперечности полей. Подчеркнем, приведенная впервые в данной работе система из 12 уравнений ранее не использовалась в теории эффекта Казимира. Разработанный метод эффективен для решения различных задач теории эффекта Казимира, особенно при наличии геометрий с плоскопараллельными поверхностями и заданными граничными условиями на границах раздела сред. С использованием найденных электрической и магнитной функций Грина в работе также получена классическая формула Лифшица для давления в системе из двух диэлектрических полупространств, разделенных вакуумной щелью.

Key findings for the stage (summarized)

Recently, using the gradient and integral molecular density functional method, we have shown that not only stable droplets are formed from vapor around solid lyophilic particles, but also thin stable concentric vapor shells-bubbles appear around lyophobic particles in a stretched liquid. Such droplets and bubbles correspond to local minima of the grand thermodynamic potential of the system with the new phase nucleus, the initial phase, and the nucleation seed. One of the objectives of our work within the framework of the RSF grant is to study stable droplets and bubbles on lyophilic and lyophobic nucleation seeds using the modern molecular density functional method with the most complete allowance for hard-sphere correlations in terms of fundamental measure theory. Within this approach, equilibrium three-dimensional density profiles of droplets around solid lyophilic spherical particles in undersaturated and supersaturated vapor and thin concentric vapor shells on lyophobic particles in a stretched and stable liquid were calculated in 2022. The existence of enveloping stable droplets and bubbles on nanosized lyophilic and lyophobic nucleation seeds has been confirmed, the structure of droplets and bubbles has been described, dependences of the chemical potential of molecules in a droplet and in a bubble on the thickness of a liquid or vapor film have been plotted, and the threshold values of the chemical potential of vapor and liquid for barrier-free heterogeneous nucleation have been found. With the help of the expression for a grand thermodynamic potential as a molecular density functional, the disjoining pressures in thin liquid films around nanosized wettable spherical particles and in thin vapor layers around nonwettable particles are calculated depending on the energy parameters of molecular interactions, film thickness, and particle size. Although the results show a qualitative agreement between the calculated dependences of the disjoining pressure and those obtained in the framework of a simpler gradient method of the molecular density functional, the new results differ significantly quantitatively. It has been confirmed that the disjoining pressure in a liquid film around a nanosized lyophilic particle increases with an increase in particle size and lyophilicity.

When studying inhomogeneous fluids using the density functional theory in many of its variants, one employs expressions for the local free energy, as well as for the chemical potential in a system of hard spheres without attraction as functions of the concentration of the hard spheres. Exact expressions for these functions are unknown, but there are various approximations that demonstrate different accuracy. In 2022, we have considered and compare the results given by the Carnahan–Starling, the Percus–Yevik, the Kolafa equations of state, and the 6th and the 7th-order Rusanov equations of state for a system of hard spheres without attraction. The dependences of the compressibility factor, chemical potential, and free energy density on the particle number density have been considered and compared with the modified Kolafa–Labík–Malijevský equation describing results of molecular dynamics simulations of the compressibility factor with high accuracy, though in a limited density range. We have calculated, within the integral density functional theory in random-phase approximation, the molecular density profiles in radially inhomogeneous small spherical droplets and bubbles of an argon-like substance representing critical droplets/bubbles in nucleation. We have such calculations for the Carnahan–Starling, Kolafa and truncated 6th-order Rusanov equations of state for hard spheres, as well as for the modified Kolafa–Labík–Malijevský equation. To describe the fluid, all the employed equations of state were supplemented with a term describing the attraction of particles in a mean-field approximation. When comparing the equations of state for a system of hard spheres to describe both a homogeneous fluid and critical droplets/bubbles, the Kolafa equation turned out to be the best of the considered equations of state for a system of hard spheres, repeating the results for the “reference” Kolafa–Labík–Malijevský equation with high accuracy. At the same time, the Kolafa equation is much simpler and resembles the widely used Carnahan–Starling equation, significantly exceeding it in accuracy in describing both homogeneous and inhomogeneous fluids of hard spheres.

The increase in performance of computing equipment as well as its accessibility to scientific community inevitably leads to more and more frequent use of molecular modeling methods, such as molecular dynamics and Monte Carlo methods, for tackling the fundamental and applied scientific problems. These problems certainly include the studies of the properties of sessile droplets on the surfaces of different structure and chemical composition. One of the most important properties of the surface is its wettability, because its modification opens up wide opportunities for practical applications. In 2022, using the molecular dynamics method, we simulated the formation of a sessile droplet on a flat surface and calculated its contact angle which is the macroscopic characteristic of the surface wettability. Two methods for calculating the contact angle differing in the clustering procedure (i.e., separation of the molecules constituting the droplet from the molecules of the gas phase) were used. The droplets consisting of argon or water molecules were considered. An impenetrable wall interacting with a drop according to the Steele potential (“3-9” potential) was used as a flat surface, as well as the all-atom models of graphite and cristobalite (SiO₂). For all the considered models of the surface we varied the degree of its wettability. The dependences of the contact angle of a droplet on its size and temperature were studied.

As is known, the question of rigorous calculation of the stress tensor of a fluctuating electromagnetic field, free energy, and the Casimir–Polder potential in inhomogeneous systems in the presence of thin films, drops, bubbles, vacuum and air slots is of fundamental importance. Simple molecular field integrations for solids can introduce uncontrollable errors. An important task of this project is to study the dependence of free energy and the Casimir–Polder potential on the geometric characteristics of systems and on the physical and chemical properties of materials. In 2022, the theory of the Casimir effect was considered for two diffraction gratings separated by a vacuum slit. Discontinuity in the Casimir energy during rotation of the grating is shown to be the consequence of the fundamental transformation of structure of reciprocal vector space, which is a result of translational symmetry breaking.

A novel explicitly gauge-invariant method for derivation of electric and magnetic Green functions in the Casimir effect theory is developed. The method is based on application of Weyl formula and makes it possible to derive the Casimir pressure and the Casimir–Polder potential of an anisotropic atom by a solution of an explicit gauge-invariant system of equations following from transversality of fields and boundary conditions.
AcronymRSF_RG_2022 - 1
StatusFinished
Effective start/end date18/05/2231/12/22

ID: 95229566