Read--Shockley formula for a general Bravais lattice in two dimensions

以下是基于论文《READ–SHOCKLEY FORMULA FOR A GENERAL BRAVAIS LATTICE IN TWO DIMENSIONS》（二维一般布拉维晶格的 Read-Shockley 公式）的详细技术总结：

1. 研究背景与问题 (Problem Statement)

• 核心问题：研究二维晶体中两个具有微小取向差（misorientation angle, $\theta$）的晶粒之间形成的**小角度晶界（grain boundary）**的能量。
• 物理模型：采用由 Lauteri 和 Luckhaus 提出的非线性半离散位错能量模型（semi-discrete dislocation energy）。该模型结合了连续介质力学（弹性应变能）和离散结构（位错核心能）。
• 目标：
  1. 为**任意布拉维晶格（General Bravais Lattice）**构建一个简单且物理意义明确的晶界构造。
  2. 证明该构造的能量上界符合著名的 Read-Shockley 公式，即能量随取向差 $\theta$ 呈 $\theta |\log \theta|$ 的标度关系。
  3. 显式计算公式中的常数，并验证其与正方形晶格特例及经典文献结果的一致性。

2. 方法论 (Methodology)

2.1 能量泛函定义

定义在区域 $\Omega$ 上的能量泛函 $F_\varepsilon(\beta, S)$ 包含两项：

1. 弹性应变能：$\int_{\Omega \setminus B_{\lambda\varepsilon}(S)} \text{dist}^2(\beta, SO(2)) dx$，衡量应变场 $\beta$ 偏离刚体旋转的程度。
2. 核心能（Core Energy）：$\frac{\tau^2}{\lambda^2} \mathcal{L}^2(B_{\lambda\varepsilon}(S))$，惩罚位错核心区域的长度，其中 $S$ 是位错集合，$\varepsilon$ 是微观尺度，$\tau$ 是晶格间距，$\lambda$ 是核心半径参数。

2.2 构造策略 (The Construction)

作者提出了一种比现有文献（如 [4, 9]）更简单的构造方法，用于连接两个旋转角度分别为 $R_{-\theta}$ 和 $R_{\theta}$ 的晶粒：

• 域分割：将晶界区域分割为四个垂直条带。中间两个条带（$\Sigma_1, \Sigma_2$）用于容纳位错阵列，两侧条带（$\Sigma_{-\theta}, \Sigma_{\theta}$）匹配边界条件。
• 分形/层级结构：
  • 在中间条带内，将区域划分为一系列正方形 $Q_i$。
  • 每个正方形进一步细分为二进环形（dyadic square annuli） $D_n$。
  • 利用凸多边形（三角形 $\Delta$ 和梯形 $T$）对正方形进行细分，形成**分片仿射（piecewise affine）**变形场。
• 位错模拟：
  • 通过引入Burgers 矢量（晶格生成元 $b_1$ 或 $b_2$）的累积平移来模拟位错。
  • 随着向下移动（沿晶界方向），累积的平移量增加，模拟了晶格为了容纳多余半原子面而发生的“张开”。
  • 应变场 $\beta$ 在大部分区域是常数（分片常数），仅在位错核心附近（$D_0$ 的边界）发生跳变。
• Rank-one 连接：构造的关键在于确保相邻区域的应变值在公共界面上满足 Rank-one 连接（即 $\beta_1 - \beta_2 = a \otimes n$），从而最小化界面能。作者通过调整条带宽度参数 $r_{\bar{n}_i}$ 来满足这一条件。

2.3 可行性验证

• 验证了构造满足半离散能量模型的所有假设（H1-H3），特别是绕闭合曲线的应变积分（Burgers 矢量）必须属于晶格 $\tau\varepsilon B$。
• 证明了构造中的应变场是容许的（Admissible），且 Curl $\beta$ 仅集中在位错核心边界上。

3. 主要结果 (Key Results)

3.1 能量上界估计

通过计算弹性能和核心能，作者证明了总能量满足以下上界：
$$F_\varepsilon(\beta, S) \leq C \varepsilon L E_0 \theta (A - \log \theta)$$
其中：

• $C$ 是常数。
• $L$ 是晶界长度。
• $E_0$ 和 $A$ 是依赖于晶格方向（由角度 $\phi, \eta$ 定义）和核心半径的常数。
• 该结果严格符合 Read-Shockley 公式 的对数标度行为。

3.2 常数显式表达

作者显式推导了常数 $E_0$ 和 $A$ 的表达式（公式 20-21, 25）：

• $E_0$：与晶格生成元 $b_1, b_2$ 相对于晶界法线的角度有关。
• $A$：包含对数项，依赖于核心半径 $\lambda$ 和晶格参数。

3.3 特例验证：正方形晶格

在正方形晶格（$B=\mathbb{Z}^2$）的特殊情况下：

• 推导出的位错间距 $2r_{\bar{n}_i}$ 与 Read 和 Shockley 在 1950 年经典论文中的结果完全一致。
• 常数 $E_0$ 与 Read-Shockley 公式中的系数完全吻合。
• 常数 $A$ 与经典结果存在细微差异（核心区域形状为正方形而非圆形导致的几何因子差异，以及核心能模型的简化），但在主导阶（leading order）上是一致的。

4. 主要贡献 (Key Contributions)

1. 通用性：首次为任意二维布拉维晶格（不仅仅是正方形晶格）提供了显式的晶界构造和能量上界证明。
2. 简化构造：相比之前的文献（如 [4] 需要三个条带），本文提出了仅需两个条带的更简单构造。
3. 分片常数应变：构造中的应变场是**分片常数（piecewise constant）**的，这使得问题可以转化为矩阵连接问题，且满足更严格的变分问题（公式 4），即 $\beta$ 为分片常数时的能量极小化。
4. 物理直观性：构造清晰地展示了 Burgers 矢量的累积效应和晶格的几何适应过程，具有更强的物理直观性。

5. 意义与影响 (Significance)

• 理论验证：从严格的数学分析角度，验证了 Read-Shockley 公式在一般晶格和半离散能量模型下的有效性，填补了从离散位错模型到宏观晶界能量公式之间的理论空白。
• 模型扩展性：这种简化的构造方法被认为易于推广，例如用于处理多势阱能量（multi-well energy）、孪晶（twins）和层状结构（laminates）等问题。
• 计算材料学基础：为理解小角度晶界的形成机制、能量标度律以及位错阵列的稳定性提供了坚实的数学基础，有助于指导多晶体材料的微观结构设计。

总结：该论文通过构建一个巧妙的、基于分形细分和分片仿射变形的几何构造，成功证明了在一般二维布拉维晶格下，小角度晶界的能量遵循 Read-Shockley 公式，并给出了精确的常数表达式，是位错理论和晶界能量研究中的重要进展。