2. 微结构演化物理基础

相场法之所以能模拟微结构演化，根本原因在于它把“系统往哪里变、变得多快”这件事，写成了自由能随成分场与序参量场变化的数学问题。理解这一节，需要抓住一条主线：自由能决定方向，化学势决定扩散，界面能决定形貌，曲率与热力学驱动力共同决定长大或粗化。下面按物理量之间的因果关系展开。

一、自由能：演化的“地形图”

在恒温恒压下，系统自发演化的判据是总自由能降低。相场模型通常把总自由能写成体自由能与梯度能之和：

F = ∫ [ f(c, η) + (κ_c/2)|∇c|² + (κ_η/2)|∇η|² ] dV

其中 c 是成分场，η 是序参量（如有序度、晶体取向或相标识），f 是均匀体自由能密度，κ 是梯度能系数。体自由能 f 常写成双阱形式，例如 f = W·c²(1−c)²，两个极小值分别对应两相平衡成分；W 是势垒高度，控制相分离驱动力的大小。梯度能项则惩罚成分或序参量的空间突变，是界面具有有限厚度的来源。

需要强调：自由能不是“能量数值”本身，而是演化方向的地形图。相场方程本质上是让系统沿自由能下降方向运动，因此自由能函数的形状直接决定最终微结构是层状、滴状还是连通网络。

二、化学势：扩散的真正驱动力

很多人误以为浓度梯度驱动扩散，严格来说，化学势梯度才是扩散的驱动力。对非均匀系统，化学势由变分导数给出：

μ = δF/δc = ∂f/∂c − κ_c ∇²c

扩散通量遵循 J = −M ∇μ，其中 M 是原子迁移率。代入连续性方程 ∂c/∂t = −∇·J，得到 Cahn–Hilliard 方程：

∂c/∂t = ∇·[ M ∇( ∂f/∂c − κ_c ∇²c ) ]

这个方程说明：即使浓度梯度存在，若化学势梯度为零，净扩散仍为零；反之，在自由能曲率变号（∂²f/∂c² < 0）的失稳区，微小浓度起伏会自发放大，这就是调幅分解的物理根源。迁移率 M 只控制“多快”，不改变“往哪走”。

三、界面能：形貌选择的“成本项”

界面能来自梯度能项，可近似写为 γ ≈ ∫ √(2κ_c Δf) dc，其中 Δf 是势垒高度。界面能的存在使系统倾向于减少界面面积，这是粗化的根本原因。但界面能并非只起“收缩”作用：在弹性、各向异性或外场耦合下，界面能各向异性会诱导枝晶、片层或立方状析出。相场中常用 κ(θ) 或 W(θ) 随界面法向角变化来引入各向异性。

四、曲率驱动与 Gibbs–Thomson 效应

弯曲界面的局部平衡浓度高于平界面，这就是 Gibbs–Thomson 效应：

c_eq(R) ≈ c_eq(∞) · exp( 2γΩ / (R T R_curv) )

其中 Ω 是原子体积，R_curv 是曲率半径。曲率越大（半径越小），平衡浓度越高，因此小颗粒周围化学势更高，原子从小颗粒向大颗粒扩散，表现为Ostwald 熟化。在相场模拟中，这一效应不需要额外写规则，只要正确设置 κ 与 M，曲率驱动会自然从 −κ∇²c 项中涌现。

五、热力学驱动力：从“能不能”到“多快”

热力学驱动力通常定义为 ΔG = G_母相 − G_新相，它决定相变能否发生以及新相长大速度。在相场中，驱动力通过体自由能差 Δf 进入演化方程。若 Δf 为正，新相长大；若为零，界面处于平衡；若为负，新相回溶。实际模拟中，常把驱动力写成温度、成分或外场的函数，例如 Δf = L·(T − T₀)/T₀，其中 L 是潜热，T₀ 是平衡温度。

六、各物理量的耦合关系

把上述概念串起来，可以得到一条清晰的因果链：自由能地形决定化学势分布；化学势梯度驱动扩散通量；扩散改变成分场与序参量场；界面能通过梯度项约束形貌；曲率通过 Gibbs–Thomson 效应调节局部平衡浓度；热力学驱动力决定相变方向与速率。相场方程只是把这六个环节写成一个偏微分方程组，并用数值方法迭代求解。

七、工程实现中的关键参数与检查

在具体写代码时，建议先固定以下量纲与参数，再调演化行为：

一个实用的稳定性判据是：Δt ≤ Δx⁴ / (4 M κ_c)。若模拟中出现成分场振荡或自由能上升，优先减小 Δt 或增大 κ_c。此外，界面宽度应至少覆盖 4～6 个网格点，否则曲率驱动会被数值噪声淹没。

最后，把本节内容压缩成一句话：相场法用自由能定义“地形”，用化学势定义“坡度”，用扩散定义“流速”，用界面能和曲率定义“河道形状”，用热力学驱动力定义“水流方向”。理解这五者的分工与耦合，后续编写 Cahn–Hilliard 与 Allen–Cahn 方程时，才能清楚每一个参数在物理上到底控制什么。