Surface Phase-Field-Crystal-Helfrich model for out-of-plane deformations in thin crystalline sheets with lattice mismatch

技术摘要：用于面外变形的表面相场晶体-赫尔菲希（sPFC-Helfrich）模型

问题陈述
薄且具有柔性的晶体片表现出由面内弹性与面外变形耦合驱动的复杂力学行为。描述此类系统需要一个综合框架，既能体现晶体结构及其缺陷（如位错），也能体现异质结构中的晶格失配。虽然以往的方法利用了原子论描述、连续介质模型或多尺度方法，但仍需要一种既能保留原子细节，又能捕捉宏观弯曲效应以及空间变化平衡晶格参数的框架。具体而言，理解局域化晶格本征应变（模拟晶格失配）与缺陷如何相互作用从而驱动面外变形，仍然是一个关键挑战。

方法论
作者通过将表面相场晶体（sPFC）模型与宏观 Helfrich 能量描述相结合来处理弯曲效应，从而扩展了 sPFC 模型。该方法的核心组成部分包括：

1. 模型构建： 系统由自由能泛函 $F = F_L + F_P + F_b$ 定义，其中 $F_L$ 和 $F_P$ 代表 sPFC 能量（表征晶格结构和相空间），而 $F_b$ 代表 Helfrich 弯曲能量。
   • 一个关键的扩展是允许使用空间变化的平衡晶格参数 $q(x)$，从而能够表示局域化的晶格本征应变，以模拟异质结构中的晶格失配。
   • 该模型采用高度观测者参数化方案 $X(t, x, y) = [x, y, h(t, x, y)]^T$ 来描述表面变形。
2. 数学框架： 系统演化受 $(L^2, H^{-1})$ 表面梯度流控制。为了确保在演化表面上的能量耗散和质量守恒（$\frac{d}{dt}\int_S \psi dS = 0$），作者使用了 Truesdell 时间导数和 Truesdell 表面独立性规范。这种方法解决了以往推导中关于表面局部扩张或收缩的局限性。
3. 数值实现： 使用傅里叶伪谱法结合半隐式时间步进方案来求解密度场 $\psi$ 和高度场 $h$ 的耦合方程。非线性项在物理空间进行计算，并通过快速傅里叶变换（FFT）进行转换。

主要贡献

• 扩展模型： 开发了一种多尺度 sPFC-Helfrich 模型，该模型整合了空间变化的平衡晶格间距，从而能够在晶体片内模拟晶格失配和本征应变。
• 理论精炼： 利用 Truesdell 规范对表面梯度流进行了严谨的推导，这引入了以往简化模型中所没有的附加项（例如 $\psi \Delta_S h$）。
• 验证： 该模型通过两个经典的解析基准进行了验证：
  • Föppl–von Kármán (FvK) 方程： 在单轴压缩下，模型的屈曲振幅预测与解析解一致，包括高阶修正。
  • Eshelby 夹杂问题： 对于具有本征应变的球形夹杂物，模型重现了夹杂内部的均匀应力场以及外部的衰减过程，当弯曲刚度较高时，这与经典弹性理论一致。

结果
经过验证的框架被应用于研究局域化晶格失配、缺陷与面外变形之间的相互作用：

• 弯曲刚度效应： 降低弯曲刚度 ($\beta$) 会引入面外变形的尺寸依赖性。与经典的 Eshelby 问题（或高 $\beta$ 情况）不同，较低的刚度允许薄片发生弯曲，从而提供了一种额外的松弛机制，部分释放了夹杂物内的应力。
• 缺陷驱动变形： 相对于基体旋转的夹杂物（产生了带有孤立位错的小角度晶界）会诱发面外凸起。这些凸起形成于静水压力为压缩应力的位错处。凸起的方向取决于初始扰动，但相邻位错之间的相互作用会影响最终的高度轮廓。
• 综合效应： 在结合了本征应变（晶格失配）和旋转诱导位错的模拟中，来自夹杂物的应力释放与位错引起的凸起发生叠加。结果表明，缺陷作为应变释放的额外来源，降低了与具有相同本征应变的无缺陷情况相比的整体屈曲振幅。因此，夹杂物平均承受的压缩应力较小，而其与周围基体的界面则承受较大的拉应力。

意义与主张
本文声称，开发的 sPCL-Helfrich 模型为预测和控制二维晶体片的形貌及力学响应提供了一个稳健的框架。通过成功地将原子细节与宏观弯曲耦合在一起，该模型阐明了缺陷形成（塑性变形）与面外变形作为内部应力松弛机制之间的竞争关系。作者指出，虽然目前的工作侧重于特定的验证案例和简单的异质结构，但该方法适用于更大的系统和更复杂的结构，例如多晶金属箔。论文谦虚地总结道，若要重现存在大角度晶界下的高度调制以及动态条件下的表现，仍需进一步评估，这将是未来工作的课题。