Анализ роли межмолекулярных сил в фазовых переходах жидкость–газ и жидкость–твердое тело: как эти силы влияют на критические параметры и структуру конденсата, и какие модели описывают изменения при переходах.
Краткий ответ с ключевыми выводами и нужными формулами. 1) Роль межмолекулярных сил — общая картина - Притяжения (дисперсионные, водородные, ионные) определяют стабильность конденсатов, величину температур и давлений фазовых переходов; более сильное среднее притяжение → выше критическая температура TcT_cTc, выше энергия конденсации и больше плотность жидкости. - Отталкивания (короткий радиус отталкивания, «объём занятых частиц») задают уплотняемость, структуру ближнего окружения и влияют на PcP_cPc и объём конденсата. - Дальность действия сил (короткий vs длинный диапазон) влияет на универсальность поведения при критической точке и на кинетику переходов (например, длиннодальние кулоновские взаимодействия меняют класс универсальности). 2) Жидкость–газ: влияние на критические параметры и структуру - Простейшая модель: уравнение Ван‑дер‑Ваальса (мольные параметры a,ba,ba,b) даёт явную зависимость критических параметров: Tc=8a27Rb,Pc=a27b2,Vc=3b.T_c=\frac{8a}{27R b},\qquad P_c=\frac{a}{27 b^2},\qquad V_c=3b.Tc=27Rb8a,Pc=27b2a,Vc=3b. Здесь aaa отражает среднее притяжение, bbb — объём исключённый отталкиванием; рост aaa повышает TcT_cTc и PcP_cPc. - Микропотенциалы: например Леннард‑Джонc U(r)=4ε[(σr)12−(σr)6],U(r)=4\varepsilon\Big[\Big(\frac{\sigma}{r}\Big)^{12}-\Big(\frac{\sigma}{r}\Big)^6\Big],U(r)=4ε[(rσ)12−(rσ)6], где глубина ямы ε\varepsilonε задаёт шкалу TcT_cTc (Tc∝ε/kBT_c\propto\varepsilon/k_BTc∝ε/kB). - Структура описывается радиальной корреляцией g(r)g(r)g(r) и структурным фактором S(k)S(k)S(k): S(k)=1+ρ∫e−ik⋅r(g(r)−1) dr.S(k)=1+\rho\int e^{-i\mathbf{k}\cdot\mathbf{r}}\big(g(r)-1\big)\,d\mathbf{r}.S(k)=1+ρ∫e−ik⋅r(g(r)−1)dr. Ближе к TcT_cTc растёт корреляционная длина ξ\xiξ, флуктуации плотности усиливаются. - Критические явления: порядок (плотность) ведёт себя как Δρ=ρℓ−ρv∼(Tc−T)β,\Delta\rho=\rho_\ell-\rho_v\sim (T_c-T)^\beta,Δρ=ρℓ−ρv∼(Tc−T)β, с универсальным показателем β≈0.325\beta\approx0.325β≈0.325 для 3D‑Ising‑класса; корреляционная длина ξ∼∣T−Tc∣−ν,ν≈0.63.\xi\sim|T-T_c|^{-\nu},\qquad \nu\approx0.63.ξ∼∣T−Tc∣−ν,ν≈0.63. В области далеко от TcT_cTc применимы mean‑field/вандерваальсовы приближения; вблизи TcT_cTc важны флуктуации (Ginzburg‑критерий). 3) Жидкость–твердое тело (кристаллизация): влияние сил и структура конденсата - Структура кристалла определяется форме потенциальной ямы и симметрией взаимодействий: жёсткие сферические отталкивания → плотноупакованные fcc/hcp (жёсткие сферы); потенциалы с направленными связями (водородные, ковалентные) → открытые структуры (лед, алмаз). - Баланс коротко‑диапазонного сильного отталкивания и умеренных притяжений определяет плотность и тип упорядоченности; сила притяжения влияет на температуру плавления TmT_mTm и величину скрытой теплоты плавления. В общем: сильнее притяжение → выше TmT_mTm и большая энергия кристаллообразования. - Модели: классические DFT (density functional theory) и фазефилд‑модели/Phase‑Field Crystal описывают переход жидкость→кристалл и предсказывают пространственную плотность ρ(r)\rho(\mathbf r)ρ(r) с образованием периодической решётки. 4) Модели и методы, описывающие изменения при переходах - Макроскопические/средние модели: уравнение Ван‑дер‑Ваальса, теория законов соответствия (law of corresponding states). - Латтис‑газ ↔ Изинг: дискретная модель, полезна для изучения универсальности жидкость–газ. - Статистические методы: виральное разложение, интегральные уравнения Орнштейна‑Цернике с замыканиями (HNC, Percus‑Yevick) для корреляций. - Пертурбационные подходы: разделение потенциала на жёсткую часть + слабые притяжения (WCA, Barker‑Henderson). - Молекулярная динамика / Монте‑Карло — численные симуляции с реальным потенциалом (LJ, полярные, ковариантные). - Классическая DFT — строит функционал свободной энергии F[ρ(r)]F[\rho(\mathbf r)]F[ρ(r)] и позволяет получать фазовые диаграммы, межфазное натяжение и структуру конденсата. - Критические явления — теория ренормгруппы (RG) для определения показателей и масштабирования. 5) Кинетика переходов: нуклеация и спинодальная декомпозиция - Классическая нуклеационная теория (сфера): критический радиус и барьер: r∗=2γΔp,ΔG∗=16πγ33(Δp)2,r^*=\frac{2\gamma}{\Delta p},\qquad \Delta G^*=\frac{16\pi\gamma^3}{3(\Delta p)^2},r∗=Δp2γ,ΔG∗=3(Δp)216πγ3, где γ\gammaγ — межфазная поверхность, Δp\Delta pΔp — переохлаждение/супернасыщение (разность термодинамических потенциалов). - При сильном пересыщении/переохлаждении вместо нуклеации возможна спинодальная декомпозиция: условие спинодали — вторая производная свободной энергии по плотности = 0: (∂2f∂ρ2)T=0.\left(\frac{\partial^2 f}{\partial\rho^2}\right)_{T}=0.(∂ρ2∂2f)T=0. 6) Особые случаи и влияние анизотропных/долгих сил - Водородные связки, ассоциативные взаимодействия и полярность дают аномалии (например, вода: расширение при замерзании, двузначная кристаллическая структура). - Длинно‑дальние дипольные/кулоновские силы меняют масштабирование и могут приводить к иным классам универсальности. Краткое резюме: межмолекулярные силы контролируют и термодинамику (критические параметры Tc,PcT_c,P_cTc,Pc), и микроструктуру конденсатов (g(r),S(k)g(r),S(k)g(r),S(k), тип кристаллической решётки), а также кинетику переходов (нуклеация vs спинодал). Для описания используют и простые уравнения состояния (ван‑дер‑ваальс, латтис‑газ), и микроскопические потенциалы (LJ), и более точные методы (классическая DFT, интегральные уравнения, симуляции, RG для критических явлений).
1) Роль межмолекулярных сил — общая картина
- Притяжения (дисперсионные, водородные, ионные) определяют стабильность конденсатов, величину температур и давлений фазовых переходов; более сильное среднее притяжение → выше критическая температура TcT_cTc , выше энергия конденсации и больше плотность жидкости.
- Отталкивания (короткий радиус отталкивания, «объём занятых частиц») задают уплотняемость, структуру ближнего окружения и влияют на PcP_cPc и объём конденсата.
- Дальность действия сил (короткий vs длинный диапазон) влияет на универсальность поведения при критической точке и на кинетику переходов (например, длиннодальние кулоновские взаимодействия меняют класс универсальности).
2) Жидкость–газ: влияние на критические параметры и структуру
- Простейшая модель: уравнение Ван‑дер‑Ваальса (мольные параметры a,ba,ba,b) даёт явную зависимость критических параметров:
Tc=8a27Rb,Pc=a27b2,Vc=3b.T_c=\frac{8a}{27R b},\qquad P_c=\frac{a}{27 b^2},\qquad V_c=3b.Tc =27Rb8a ,Pc =27b2a ,Vc =3b.
Здесь aaa отражает среднее притяжение, bbb — объём исключённый отталкиванием; рост aaa повышает TcT_cTc и PcP_cPc .
- Микропотенциалы: например Леннард‑Джонc
U(r)=4ε[(σr)12−(σr)6],U(r)=4\varepsilon\Big[\Big(\frac{\sigma}{r}\Big)^{12}-\Big(\frac{\sigma}{r}\Big)^6\Big],U(r)=4ε[(rσ )12−(rσ )6],
где глубина ямы ε\varepsilonε задаёт шкалу TcT_cTc (Tc∝ε/kBT_c\propto\varepsilon/k_BTc ∝ε/kB ).
- Структура описывается радиальной корреляцией g(r)g(r)g(r) и структурным фактором S(k)S(k)S(k):
S(k)=1+ρ∫e−ik⋅r(g(r)−1) dr.S(k)=1+\rho\int e^{-i\mathbf{k}\cdot\mathbf{r}}\big(g(r)-1\big)\,d\mathbf{r}.S(k)=1+ρ∫e−ik⋅r(g(r)−1)dr.
Ближе к TcT_cTc растёт корреляционная длина ξ\xiξ, флуктуации плотности усиливаются.
- Критические явления: порядок (плотность) ведёт себя как
Δρ=ρℓ−ρv∼(Tc−T)β,\Delta\rho=\rho_\ell-\rho_v\sim (T_c-T)^\beta,Δρ=ρℓ −ρv ∼(Tc −T)β,
с универсальным показателем β≈0.325\beta\approx0.325β≈0.325 для 3D‑Ising‑класса; корреляционная длина
ξ∼∣T−Tc∣−ν,ν≈0.63.\xi\sim|T-T_c|^{-\nu},\qquad \nu\approx0.63.ξ∼∣T−Tc ∣−ν,ν≈0.63.
В области далеко от TcT_cTc применимы mean‑field/вандерваальсовы приближения; вблизи TcT_cTc важны флуктуации (Ginzburg‑критерий).
3) Жидкость–твердое тело (кристаллизация): влияние сил и структура конденсата
- Структура кристалла определяется форме потенциальной ямы и симметрией взаимодействий: жёсткие сферические отталкивания → плотноупакованные fcc/hcp (жёсткие сферы); потенциалы с направленными связями (водородные, ковалентные) → открытые структуры (лед, алмаз).
- Баланс коротко‑диапазонного сильного отталкивания и умеренных притяжений определяет плотность и тип упорядоченности; сила притяжения влияет на температуру плавления TmT_mTm и величину скрытой теплоты плавления. В общем: сильнее притяжение → выше TmT_mTm и большая энергия кристаллообразования.
- Модели: классические DFT (density functional theory) и фазефилд‑модели/Phase‑Field Crystal описывают переход жидкость→кристалл и предсказывают пространственную плотность ρ(r)\rho(\mathbf r)ρ(r) с образованием периодической решётки.
4) Модели и методы, описывающие изменения при переходах
- Макроскопические/средние модели: уравнение Ван‑дер‑Ваальса, теория законов соответствия (law of corresponding states).
- Латтис‑газ ↔ Изинг: дискретная модель, полезна для изучения универсальности жидкость–газ.
- Статистические методы: виральное разложение, интегральные уравнения Орнштейна‑Цернике с замыканиями (HNC, Percus‑Yevick) для корреляций.
- Пертурбационные подходы: разделение потенциала на жёсткую часть + слабые притяжения (WCA, Barker‑Henderson).
- Молекулярная динамика / Монте‑Карло — численные симуляции с реальным потенциалом (LJ, полярные, ковариантные).
- Классическая DFT — строит функционал свободной энергии F[ρ(r)]F[\rho(\mathbf r)]F[ρ(r)] и позволяет получать фазовые диаграммы, межфазное натяжение и структуру конденсата.
- Критические явления — теория ренормгруппы (RG) для определения показателей и масштабирования.
5) Кинетика переходов: нуклеация и спинодальная декомпозиция
- Классическая нуклеационная теория (сфера): критический радиус и барьер:
r∗=2γΔp,ΔG∗=16πγ33(Δp)2,r^*=\frac{2\gamma}{\Delta p},\qquad \Delta G^*=\frac{16\pi\gamma^3}{3(\Delta p)^2},r∗=Δp2γ ,ΔG∗=3(Δp)216πγ3 ,
где γ\gammaγ — межфазная поверхность, Δp\Delta pΔp — переохлаждение/супернасыщение (разность термодинамических потенциалов).
- При сильном пересыщении/переохлаждении вместо нуклеации возможна спинодальная декомпозиция: условие спинодали — вторая производная свободной энергии по плотности = 0:
(∂2f∂ρ2)T=0.\left(\frac{\partial^2 f}{\partial\rho^2}\right)_{T}=0.(∂ρ2∂2f )T =0.
6) Особые случаи и влияние анизотропных/долгих сил
- Водородные связки, ассоциативные взаимодействия и полярность дают аномалии (например, вода: расширение при замерзании, двузначная кристаллическая структура).
- Длинно‑дальние дипольные/кулоновские силы меняют масштабирование и могут приводить к иным классам универсальности.
Краткое резюме: межмолекулярные силы контролируют и термодинамику (критические параметры Tc,PcT_c,P_cTc ,Pc ), и микроструктуру конденсатов (g(r),S(k)g(r),S(k)g(r),S(k), тип кристаллической решётки), а также кинетику переходов (нуклеация vs спинодал). Для описания используют и простые уравнения состояния (ван‑дер‑ваальс, латтис‑газ), и микроскопические потенциалы (LJ), и более точные методы (классическая DFT, интегральные уравнения, симуляции, RG для критических явлений).