A Nonlocal Orientation Field Phase-Field Model for Misorientation- and Inclination- Dependent Grain Boundaries

本文提出了一种非局部取向场相场模型，该模型利用单一取向场结合了取决于误取向和倾斜角的晶界各向异性，从而在简化拟合过程的同时，实现了对晶界能量的精确调控，并能准确重现线性晶粒生长和三叉点平衡等关键微观结构行为。

要点

想象一块金属块或一块陶瓷砖。在显微镜下，你看到的并不是单一、均匀的材料。相反，你看到的是由许多微小晶体（称为晶粒）组成的拼布被面。当两个晶粒相遇时，那里会有一个被称为晶界的边界。

可以将这些晶粒想象成拥挤房间里的行人。每个人的朝向都略有不同。晶界就是两个朝向不同的个体并排而立的那条线。

问题所在：“地图”缺失了

科学家使用计算机模拟（称为相场模型）来预测材料随时间的变化方式——例如金属如何变强或晶体如何生长。为了实现这一点，他们需要一张数学“地图”，告诉计算机形成一个晶界需要多少能量。

问题在于，边界的能量取决于两个棘手的因素：

1. 取向差（Misorientation）： 两个邻居彼此错开的角度（比如两个人的朝向相差10度与相差90度的区别）。
2. 倾角（Inclination）： 边界线本身切割材料的角度（比如一条南北走向的直线围栏与一条斜穿田野的围栏之间的区别）。

之前的计算机模型就像是在试图用一张只显示街道却没显示建筑物的地图来导航城市。它们可以处理简单的情况，但在处理晶粒旋转方式复杂或边界发生倾斜的情况时，难以准确预测能量。它们要么需要过高的计算能力，要么做出了过多的简化假设。

解决方案：一个“非局部”望远镜

本文作者提出了一种构建这张地图的新方法。他们称之为非局部取向场相场模型。

这里有一个类比：
想象你正站在两个社区交界处的边界上（即晶界）。在旧模型中，你只能看到你脚下的那条街道。你不知道另一侧的社区是什么样子的。

在这个新模型中，计算机给了你一个望远型望远镜。尽管你站在边界线上，望远镜能瞬间“向左看”一段距离，并“向右看”一段距离。它会立即告诉你：

• “好的，左边的晶粒朝向北。”
• “右边的晶粒朝向东。”

因为计算机现在能同时感知两侧的取向，所以无论这个特定的边界如何扭曲或倾斜，它都能计算出该特定边界的确切能量成本。

它是如何工作的（“智能围栏”）

该模型使用一条单一且平滑的线来表示晶粒之间的边界。

• 核心内部： 在边界的正中间，模型使用了一个特殊的“能量函数”，它了解倾斜和扭转的信息。这就像一个智能围栏，清楚地知道将两个特定的人固定在一起需要多少努力。
• 外部边缘： 当你远离边界进入坚实的晶粒内部时，模型会切换到一种更简单的规则，以确保晶粒保持固体状态而不变得“模糊”。

作者用几种场景测试了这种“望远镜”方法：

1. 稳定性： 他们检查了边界是否能稳定在正确的形状上。结果是肯定的。
2. 能量准确性： 他们测试了当旋转晶粒或倾斜边界时，能量是否会发生正确变化。结果与数学计算完全吻合。
3. 生长： 他们模拟了一个小晶粒在大晶粒内部缩小的过程（就像气泡破裂一样）。该模型正确预测了缩小的速度。
4. 复杂形状： 他们展示了该模型如何预测晶粒在试图最小化其能量时所呈现出的奇特、非圆形形状（称为 Wulff 形态），这取决于能量的各向异性（方向依赖性）。

为什么这很重要

其主要成就在于简洁与精准。

• 旧方法： 要模拟一个拥有100个不同晶粒的材料，你可能需要同时运行100个不同的数学方程，这既慢又笨重。
• 新方法： 该模型无论有多少个晶粒，都只使用一个方程来处理整个系统。它捕捉到了每个晶界的复杂“个性”，而无需为每一个晶界准备单独的方程。

简而言之，作者为计算机构建了一种更聪明、更高效的方式，让计算机能够“看见”将晶体连接在一起的无形力量，从而在不需要超级计算机进行繁重计算的情况下，实现对材料行为更准确的预测。

技术摘要

问题陈述
相场建模是研究多晶材料微观结构演化的标准工具，然而准确描述晶界（GB）能量仍然是一个重大挑战。现有方法面临着截然不同的局限性：多相场（MPF）方法为每个晶粒分配一个单独的标量场，其计算规模随系统尺寸增加而迅速恶化。经典的 KWC 模型通过将单一取向场与一个序参数耦合，本质上限制了晶界能量必须是误取向度的单调递增函数，从而无法实现通用的误取向度依赖性。其他试图引入误取向度依赖性的方法通常需要复杂的反向拟合程序，或者放弃平滑的时间演化方程而转向非平滑优化，亦或是缺乏能够参数化任意各向异性的数学保证。目前需要一种紧凑、高效的相场框架，能够在不引入大量辅助场的情况下，实现任意误取向度及倾角依赖的界面过剩自由能描述。

方法论
作者提出了一种非局部取向场相场模型，该模型扩展了标准的偏微分方程（PDE）表述形式。其核心创新在于使用单一标量取向场 $\theta(\mathbf{x}, t)$，并结合了一个显式包含界面两侧晶粒取向变化的非局部泛函。

• 非局部表述： 为了确定晶界内相邻晶粒的取向（$\theta^+$ 和 $\theta^-$），该模型采用了“非局部”外推法。模型并非仅仅依赖局部场值，而是通过在晶界附近（具体为沿梯度方向 $\nabla\theta/|\nabla\theta|$ 的距离 $d$ 处）优化的位置对取向场进行采样。这使得模型能够直接从场构型中计算出误取向度（$\Delta\theta = |\theta^+ - \theta^-|$）和倾角。
• 自由能泛函： 总自由能 $F$ 由两项组成，并由一个权重函数 $w(|\nabla\theta|)$ 进行加权，该函数用于区分晶界内部区域与外部过渡/体相区域。
  • 内部项（$f_1$）在晶界内占主导地位，利用一个各向异性的晶界函数 $B_{aniso}(\theta^+, \theta^-, \mathbf{v})$，该函数取决于误取向度以及通过单位法向量 $\mathbf{v}$ 表示的倾角。
  • 外部项（$f_2$）控制晶界宽度，并使用各向同性的晶界函数 $B_{iso}(\theta^+, \theta^-)$ 来强制执行体相取向约束。
• 晶界函数： 模型定义了特定的 $B_{iso}$ 和 $B_{aniso}$ 函数形式，以捕捉晶格对称性和倾角各向异性。这些函数经过设计具有灵活性，允许在独立于特定材料参数的情况下系统地研究模型行为。
• 时间演化： 系统通过 Allen-Cahn 型方程 $\partial\theta/\partial t = -M \delta F/\delta \theta$ 进行演化。泛函导数通过解析方式推导，非局部采样通过梯度搜索算法处理。对于三叉点（triple junctions），由于搜索方向具有歧义，模型设定了最大搜索距离以确保数值稳定性。
• 实现： 使用有限差分法（FDM）结合二阶中心差分和显式前向欧拉时间积分来求解偏微分方程。

主要贡献

1. 通用各向异性： 该模型利用单一取向场成功引入了任意误取向度和倾角依赖的晶界能量，克服了 MPF 的扩展性问题以及 KWC 模型的单调性限制。
2. 非局部采样： 通过引入一个在特定距离处采样取向场的非局部泛函，模型无需辅助场即可获取必要的晶粒取向信息，从而计算界面能量。
3. 精确调控： 该形式化方法通过显式的晶界函数实现了对晶界能量各向异性的简单且精确的调控，减少了与反向问题相关的广泛拟合需求。
4. 框架不变性： 作者证明了该模型具有内在的框架不变性，确保能量泛函在参考系发生刚性旋转时保持不变。

结果
模型通过一系列解析和数值测试得到了验证：

• 平衡剖面： 对准一维系统的模拟表明，晶界剖面能快速收敛至解析稳态解。晶界宽度可以通过材料参数 $\alpha$ 和 \beta} 进行调控。
• 晶界能量依赖性： 模型重现了预期的误取向度相关能量密度，包括基于晶格对称性（$n$）的周期性以及 Read-Shockley 型行为。至关重要的是，它成功模拟了倾角依赖的能量，展示了能量密度随晶界倾角变化的正弦波动，这与预设的各向异性函数一致。
• 迁移率与生长： 对收缩圆形晶粒的模拟表明，晶界迁移率（$M_{GB}$）与取向场迁移率（$M$）呈线性相关，并遵循经典的晶粒生长关系（$R_0^2 - R^2 = Kt$）。
• Wulff 形貌： 模型成功重现了具有不同误取向度和各向异性系数的晶粒 Wulff 形貌。模拟的晶粒形状与解析 Wulff 构建相匹配，证实了模型捕捉各向异性界面能量的能力。
• 多晶系统： 多晶系统的模拟展示了真实的晶粒生长、合并及三叉点行为，并有明确证据表明各向异性影响了晶粒形貌。

意义
本文声称，这项工作提供了一套紧凑且高效的相场方程组，能够描述任意误取向度和倾角依赖的界面过剩自由能。通过避免引入大量辅助相场，该模型弥合了原子尺度晶界能量计算与介观尺度微观结构演化描述之间的鸿沟。作者断言，该模型提供了一个透明且通用的理论框架，在简化晶界能量各向异性调控的同时，保持了传统的平滑时间演化方程，从而解决了多晶材料相场建模中的长期局限性。