An electromechanically coupled multiphase-field model with generalized kinetic relations for ferroelectrics

技术摘要：一种具有广义动力学关系的铁电体电机械耦合多相场模型

问题陈述
经典的铁电材料相场模型通常依赖于 Allen–Cahn 演化律，该定律将界面运动描述为自由能的梯度下降。这种表述本质上意味着畴壁速度与其热力学驱动力之间存在线性动力学关系。然而，实验证据和分子动力学模拟表明，铁电畴壁通常表现出非线性动力学行为，例如 Merz–Stadler 定律所描述的混合指数-幂律关系，或阈值型响应。此外，现有模型往往无法在单一多相框架内解释不同类型畴壁（例如 180° 壁与 90° 壁）所具有的截然不同的动力学特性和界面能。虽然水平集（level-set）方法允许预设动力学关系，但它们在处理多相演化时面临挑战，且需要频繁进行重新初始化。因此，需要一种能够直接引入任意、具有物理意义的动力学关系，同时保持电机械耦合和多相能力的相场框架。

方法论
作者提出了一种受水平集运动学启发并由热力学原理推导出的通用动力学多相场框架。该方法包括：

1. 热力学推导： 演化律是从局域在界面上的 Clausius-Duhem 不等式中推导出来的。该模型并非通过梯度下降来最小化自由能，而是定义了一个基于相之间体自由能、势能和梯度能差异的正则化热力学驱动力 ($f_{\beta\alpha}$)。
2. 广义演化律： 序参数 ($\phi_\alpha$) 的时间演化遵循类似于水平集方法的运动学方程：$\partial_t \phi_\alpha = \frac{1}{m} \sum_{\beta} V_{\beta\alpha} |\nabla(\phi_\alpha - \phi_\beta)|$。至关重要的是，法向速度 $V_{\beta\alpha}$ 并非由驱动力导出，而是被预设为一个显式函数 $V_{\beta\alpha} = \check{V}_{\beta\alpha}(f_{\beta\alpha})$。这使得任何具有物理意义的动力学关系都能直接集成到演化方程中。
3. 应力驱动表述： 该模型采用了吉布斯自由能密度的应力驱动表述。这种方法有助于直接应用平均应力，并利用快速傅里叶变换 (FFT) 谱解法实现电机械场的均匀化，从而提高效率。
4. 多相扩展： 该框架采用序参数的线性变换来处理多个畴变体（$M$ 个相）之间的相互作用。这种表述允许为特定的序参数对分配不同的动力学关系和梯度系数（代表界面能），从而区分不同的畴壁类型（如 180° 和 90° 壁）。
5. 正则化： 使用多障碍势（multi-obstacle potential）以确保有限的界面厚度，并防止在交界处出现伪相。正则化参数 $z$ 控制界面能和厚度，并具有一个导出的下界，以确保数值稳定性和物理有效性。

核心贡献

• 通用动力学框架： 本文引入了一种相场模型，它将演化律与特定形式的动力学关系解耦，从而可以直接指定任意非线性动力学（如 Merz–Stadler 定律、阈值定律）。
• 差异化畴壁特性： 多相表述使得在同一模拟中为不同类型的畴壁（180° 对比 90°）分配不同的动力学关系和界面能成为可能，反映了这些壁具有不同结构和能量特性的物理事实。
• 电机械耦合： 该模型在应力驱动框架内充分整合了电机械耦合，捕捉了电场和机械应力对畴演化的影响。
• 数值实现： 作者使用基于 FFT 的谱解法结合交错方案进行场均匀化，并使用显式时间积分进行畴演化，展示了该模型的效率和稳定性。

结果
针对单晶钛酸钡 (BaTiO$_3$) 的数值模拟验证了该模型：

• 动力学恢复： 模型准确地重现了预设的非线性动力学关系，包括从指数到幂律行为（Merz–Stadler）以及阈值型平方根关系的转变，证实了畴壁速度遵循分配给它的驱动力函数。
• 选择性动力学控制： 在一个秩-2 层状结构中，通过将特定类型的畴壁动力学关系设为零，模型成功抑制了特定畴壁的运动（例如冻结 180° 壁而允许 90° 壁运动），展示了对不同界面类型的独立控制。
• 界面能效应： 通过改变梯度系数比率 ($g = \gamma_{90}/\gamma_{180}$)，模型显示出不同的界面能会稳定四重接点，而相同的能量则会导致接点分裂，这突显了各向异性界面能在微观结构演化中的重要性。
• 各向异性诱导的反转： 模拟捕捉到了随着施加电场强度的增加，90° 畴壁驱动力发生符号改变的现象，这一现象是由 BaTiO$_3$ 中极强的介电常数各向异性驱动的。这导致在高电场下出现了优选畴取向的反转。
• 正则化分析： 研究确定了维持有限界面厚度并防止体相区域偏离平衡态的正则化参数 $z$ 的下界，将数值稳定性与物理长度尺度联系起来。

意义与主张
作者声称，这项工作提供了一个灵活且热力学一致的框架，用于模拟铁电体中的畴壁动力学，克服了经典 Allen–Cahn 模型的局限性。通过允许直接引入实验测定或理论推导的非线性动力学关系，该模型为畴壁动力学提供了更准确的表示。为不同类型的界面分配不同属性的能力，使其特别适用于复杂的多元相铁电系统。作者指出，虽然该模型已在铁电体中得到演示，但该框架广泛适用于涉及非线性动力学的其他界面运动材料系统，如马氏体相变。论文结论较为谦逊，承认诸如显式成核方案和降低正则化要求等进一步的改进仍是开放性的挑战。