FR39 · 晶体缺陷

曲率驱动的晶粒长大：小晶粒为什么会消失

为什么边数少于六的晶粒会缩小，而边数多于六的会长大？ 每条晶界都朝自己的曲率中心移动，三条晶界相交处的夹角趋于 120 度。于是边数少于六的晶粒，晶界向外凸，晶粒缩小；边数多于六的，晶界向内凹，晶粒长大。在二维中，这个速率只取决于边数，与晶粒大小无关。

把相场多晶退火，检验冯·诺依曼–马林斯定律：晶粒面积的变化速率只由它的边数决定。

方盒中的晶粒数 $N=70$

平均边数 $\langle n\rangle=5.94$

测得的冯·诺依曼–马林斯斜率 $\frac{\mathrm{slope}}{Lk}=1.016\pm0.017$

平均晶粒面积 $\langle A\rangle=234\ \text{格}$

把多晶退火，看它粗化

晶粒着色依据

起始晶粒数

70

N_0

放入的圆形晶粒半径

24

R/\Delta x

晶界迁移率

1

L

梯度能系数

1

k

把圆环拖到任意晶粒上跟踪它；或者把加号拖进方盒，放下一个新的圆形晶粒。

播放速度

每帧时间步数

4

n_{\mathrm{step}}

模型、方法与来源

$\frac{\partial\eta_i}{\partial t}=-L\Bigl(-\eta_i+\eta_i^3+2\eta_i\sum_{j\neq i}\eta_j^2-k\nabla^2\eta_i\Bigr),\qquad \frac{dA}{dt}=\frac{\pi}{3}Lk\,(n-6)$

每个晶粒都有自己的序参量，它们一起沿 Kazaryan、Wang、Dregia 与 Patton（1999）多序参量模型的自由能梯度弛豫，晶界能和迁移率各向同性。这样，晶界以正比于曲率的速度移动：边数少于六的晶粒缩小，多于六的长大，速率只由边数决定。方盒是周期性的；两个晶粒至少共享两个格面才算相邻。

Kazaryan, Wang, Dregia & Patton (1999), arXiv:cond-mat/9912044v1 · Mason, Lazar, MacPherson & Srolovitz (2015), arXiv:1507.03379v1

推导冯·诺依曼–马林斯定律

模型假设

每个晶粒都有自己的序参量，整个场按 Kazaryan、Wang、Dregia 与 Patton（1999）的多序参量模型弛豫，晶界能和迁移率各向同性。晶界以正比于曲率的速度移动；三条晶界相交处，夹角趋于 120 度。散点图收集了在短时间窗内边数不变的每个晶粒的面积变化率。虚线是各向同性二维晶粒长大的冯·诺依曼–马林斯定律，是参照，不是拟合。方盒是周期性的，宽 128 格；组织是模拟出来的，不是显微照片。

面积变化率对边数

所选晶粒的面积

更多数值与计算方法

所选晶粒的边数: $0$
所选晶粒的面积变化率: $0$
冯·诺依曼–马林斯斜率: $0$
已收集的晶粒速率: $0$
模型时间: $0$
场的自由能: $0$
运行次数与种子: $0$

晶粒面积是它在每个格点所占序参量份额的总和，所以变化是平滑的。速率是在边数保持不变的时间窗内做直线拟合得到的；斜率读数是把所有这样的速率对“边数减六”做加权拟合。导出文件列出种子、设置和收集到的每个速率。

物理教程

晶界在变直，代价由小晶粒来付

背景金属和陶瓷通常是多晶：由许多小晶体，即晶粒，拼成，每个晶粒的晶格取向不同。晶粒之间的晶界要花能量，所以材料热到原子能跳过晶界时，晶界就会移动，让总长度变短。

为什么值得理解晶粒尺寸决定强度、韧性和抗蠕变能力，所以退火工艺要围绕晶粒粗化的快慢来设计。在二维中，这种长大服从材料科学里少有的精确定律之一。

先抓重点

核心问题: 为什么边数少于六的晶粒会缩小，而边数多于六的会长大？
一句话直觉: 每条晶界朝自己的曲率中心移动，速度正比于曲率。三条晶界相交处，夹角弛豫到 120 度。把晶粒周界转过的总角度加起来，面积变化率就只依赖边数： $dA/dt=(\pi/3)\,M\gamma\,(n-6)$ 。

核心数学模型

曲率驱动的晶界运动

v_n=M\gamma\,\kappa,\qquad \frac{dA}{dt}=-2\pi M\gamma\ \ \text{(isolated round grain)}

半径为 $R$ 的圆形晶粒有 $\kappa=1/R$ ，所以无论多大，它的面积都以恒定速率减少。在相场模型里，乘积 $M\gamma$ 对应 $Lk$ 。

冯·诺依曼–马林斯定律

\frac{dA_n}{dt}=-M\gamma\Bigl(2\pi-n\,\frac{\pi}{3}\Bigr)=\frac{\pi}{3}\,M\gamma\,(n-6)

闭合周界的切线一共转过 2π。每个三叉点贡献 π/3 的跳变，剩下的由弯曲的边来补，正是这一部分驱动面积变化。

常见难点

决定命运的不是大小

容易误解大晶粒长大、小晶粒缩小，是因为它们的大小。

正确理解在二维中，速率取决于边数，而不是面积。小晶粒往往边数少，所以通常会缩小；但一个小的七边形晶粒照样会长大。

动手实验

01
看多晶粗化
从“多晶粗化”开始，看面积变化率对边数的散点图。
你会看到： 这些点排在虚线的冯·诺依曼–马林斯直线上；晶粒数在减少，平均面积在增大。
02
让一个孤立的圆形晶粒缩小
选择“圆形晶粒匀速缩小”，看所选晶粒的曲线。
你会看到： 它的面积沿一条直线下降直到消失，正如圆形晶粒速率恒定的预测。
03
跟踪蜂窝里的缺陷
选择“蜂窝里的一个缺陷”，把圆环拖到中间的小晶粒上。
你会看到： 六边形晶粒几乎不变；边数少的小晶粒缩小直至消失，它的邻居各少掉一条边。

曲率驱动的晶粒长大：小晶粒为什么会消失

晶界在变直，代价由小晶粒来付

核心数学模型

曲率驱动的晶界运动

冯·诺依曼–马林斯定律

常见难点

决定命运的不是大小

动手实验

看多晶粗化

让一个孤立的圆形晶粒缩小

跟踪蜂窝里的缺陷

相关实验与文章

从伯格斯矢量到合金枝晶