[论文解读] Crystallization of Simple Fluids: Relative Stability of f.c.c. and b.c.c Structures

Swarn Lata Singh, Yashwant Singh|arXiv (Cornell University)|Feb 27, 2009

Material Dynamics and Properties被引用 11

一句话总结

该论文提出了一种用于结晶流体的自由能泛函，其包含直接对相关函数（DPCF）的对称性保持部分与对称性破缺部分，利用带有Roger-Young闭合的Ornstein-Zernike方程及对称性破缺部分的微扰展开。该理论正确预测了n=12的反幂律流体凝固为面心立方（f.c.c.），而较软的势场（n≤6）则倾向于形成体心立方（b.c.c.），解决了密度泛函理论中长期存在的挑战。

ABSTRACT

A free-energy functional for a crystal that contains both the symmetry conserved and symmetry broken parts of the direct pair correlation function is developed. The free-energy functional is used to investigate the crystallization of fluids interacting via the inverse power potential ; $u(r)=ε{(σ/r)}^n$. In agreement with simulation results we find that for $n=12$ the freezing is into close packed f.c.c structure while for soft repulsions $(n\leq 6)$ b.c.c phase is more stable.

研究动机与目标

• 开发一种基于第一性原理的统计力学理论，准确捕捉流体结晶过程中面心立方（f.c.c.）与体心立方（b.c.c.）相的相对稳定性。
• 解决传统密度泛函理论在描述流体-固相变时的失败，原因在于忽略了直接对相关函数中对称性破缺的贡献。
• 引入在晶体相中出现、在熔化时消失的、定性上全新的对称性破缺DPCF成分。
• 通过反幂律势（n=12, 6, 4）的蒙特卡罗模拟数据，检验该改进泛函的预测能力。

提出的方法

• 构建自由能泛函 A[ρ] = A_id[ρ] + A_ex[ρ]，其中 A_ex 包含DPCF的对称性保持部分 c^(0) 和对称性破缺部分 c^(b)。
• 利用带有Roger-Young闭合的Ornstein-Zernike方程计算 c^(0)(r)，确保在不同密度下热力学一致性。
• 通过泛函泰勒展开建模 c^(b)(r₁,r₂;[ρ])，其主导项表示为各向同性流体三体DPCF与密度涨落 ρ_b(r) 的积分。
• 在倒空间中使用倒格矢 G 表达 c^(b)，其中系数 c^(G)(r) 依赖于差向量 r 和晶体结构。
• 通过双参数化（λ 和 ξ）对密度场与序参量场进行泛函积分，计算晶体的过剩自由能 A_ex^(b)。
• 利用巨热力学势 ΔW = W_l - W 确定流体-固相共存点，其中在相变点 ΔW/N = 0。

实验结果

研究问题

• RQ1 是否包含对称性保持与对称性破缺DPCF成分的自由能泛函能正确预测反幂律流体中f.c.c.与b.c.c.相的相对稳定性？
• RQ2 对称性破缺DPCF对流体与固相之间自由能差的贡献有多大，特别是在软排斥势情况下？
• RQ3 通过各向同性流体的三体DPCF对 c^(b) 进行微扰处理，能否可靠估计固相的稳定性？
• RQ4 为何传统DFT方法无法正确预测软排斥势下的结晶路径？
• RQ5 在n=12、6和4的情况下，引入 c^(b) 后与蒙特卡罗模拟结果的符合程度提高了多少？

主要发现

• 对于n=12，f.c.c.结构被预测为更稳定，与蒙特卡罗模拟结果一致，这得益于对称性保持与对称性破缺贡献的协同作用。
• 对于n=6和n=4，b.c.c.结构比f.c.c.更稳定，正确再现了模拟结果，其中对称性破缺DPCF起到了关键的稳定作用。
• 对称性破缺部分对自由能差 ΔW/N 的贡献不可忽视——在n=12时约占对称性保持部分的25%，在n=4时接近50%。
• 若不包含对称性破缺项，理论会高估流体的稳定性，尤其在软势场下，导致错误预测。
• 通过各向同性流体三体DPCF对 c^(b) 的微扰表达式，对较软势场（n=6,4）高估了其大小，导致预测的相变发生在比模拟更低的 γ 值处。
• f.c.c.与b.c.c.结构在 c^(G)(r) 的函数形式上存在显著差异，其大小与 r 依赖关系不同，从而解释了二者稳定性的本质区别。