基于梯度能量各向异性微扰模型的铁电/反铁电复合材料相场模拟新方法

《Materials》：A Perturbation Model of Gradient Energy Anisotropy for Phase-Field Simulation of Ferroelectrics Xiaoming Shi, Jiecheng Liu, Ke Xu, Haoyu Wang, Zheng Wang, Nan Wang, Houbing Huang and Zhuhong Liu

铁电/反铁电复合材料体系的相场模拟面临梯度能量各向异性难以高效精确描述的挑战。为解决此问题，本文发展了一种基于傅里叶谱方法的梯度能量各向异性微扰模型。该模型通过傅里叶空间微扰方案，实现了对具有空间变化的全各向异性梯度能量张量的处理，克服了以往常数系数或各向同性近似的局限。应用于铁电/反铁电复合材料体系，该方法展现出卓越的鲁棒性和收敛效率，可精准复现反铁电相图及AFE-FE相变路径，并具有优越的可扩展性，便于揭示复杂纳米畴结构及三维极化织构的演化机制。

　　

随着高功率脉冲功率技术的快速发展，对介电材料提出了日益苛刻的要求，需要它们具备更高的能量密度、充放电效率及响应速度。在这一背景下，弛豫铁电体与反铁电体因其独特的极化响应机制和优异的能量存储能力，成为了前沿的候选材料。弛豫铁电体以其扩散相变、频散介电响应及极性纳米区的存在，展现出细长的电滞回线、高能量密度和低损耗。反铁电体则可在电场诱导下发生反铁电态与铁电态之间的相变，实现快速的极化翻转，从而获得高可恢复能量密度和出色的循环稳定性。这些特性使它们在诸如高能量密度电容器和脉冲功率系统等领域极具应用前景。

然而，决定这些材料性能的根本物理机制在于其内部复杂的局域异质纳米结构。纳米尺度上的特征，例如极性纳米区、化学有序-无序分布以及晶格畸变，直接调制了极化行为和相变动力学。因此，深刻理解并精确描述这些局域结构的形成与演化，是揭示其优越性能来源的关键。当前的理论建模，特别是基于化学掺杂引起的涨落进行的相场模拟，通常使用无序场来模拟掺杂对相行为的影响。但现有模型主要聚焦于相涨落、弹性能涨落和电场能涨落，而普遍忽略了梯度能量涨落的关键作用。这种简化在铁电与反铁电相共存或竞争的复杂体系中，构成了显著的理论缺陷。由于氧八面体的协同倾斜和晶格应变，极化的空间变化具有很强的方向依赖性，导致梯度能量呈现出显著的各向异性。忽略这种效应会阻碍对非均匀极化演化的精确捕获，从而影响对相变路径、畴结构演化及储能机制的精确预测。

为了准确描述此类各向异性效应，发展一种能够纳入梯度能量各向异性的相场理论至关重要。在数值求解策略方面，傅里叶谱方法因其在周期性边界条件下的高效率，已成为求解相场方程的主流方法。它通过在频域直接计算空间微分算子，展现出谱收敛性，与实空间有限差分法相比，允许更大的时间步长，并显著提高了计算效率。尤其对于那些在弹性、介电、热导或电导率方面表现出各向异性的材料体系，傅里叶空间微扰迭代方案通常仅需数次迭代即可收敛到稳定解，表现出良好的计算鲁棒性。尽管如此，主流的傅里叶空间算法框架要么忽略了梯度能量系数的空间非均匀性，要么依赖于实空间离散化，未能充分利用傅里叶谱方法的效率。具体而言，现有的谱方法实现均无法同时处理一个同样在空间上变化的、完全各向异性的梯度系数张量。这一算法空白严重限制了在复杂的铁电-反铁电共存体系中，对极化织构演化、多相竞争及场响应机制进行高保真模拟的能力。

为应对这一关键科学挑战，本项研究提出了一种用于求解铁电材料中梯度能量各向异性的傅里叶空间微扰模型。通过在倒空间构建各向异性梯度能量算子，该模型充分考虑了极化梯度在不同晶向上的非对称响应，从而实现了对梯度能量各向异性的精确离散化和高效求解。该模型的发展不仅弥补了现有相场模拟的理论短板，也为深入理解弛豫铁电/反铁电材料中纳米结构与性能的关系提供了强大的计算工具，为合理设计下一代高性能储能材料铺平了道路。

本研究主要应用了以下几个关键技术方法：1. 相场模型，以全局极化_P和位移场_u作为序参量，采用时间相关的Ginzburg-Landau方程和应力/电场平衡方程描述极化演化。2. 傅里叶谱方法，用于高效求解相场方程，特别适用于周期性边界条件。3. 傅里叶空间微扰方案，核心创新在于构造了倒空间中的各向异性梯度能量算子，以处理梯度能量系数的全各向异性张量及其空间变化，克服了已有方法的局限。4. 数值模拟与分析，将所提出的微扰模型应用于铁电/反铁电复合材料体系，评估其鲁棒性、收敛效率和可扩展性，并进行三维模拟以可视化复杂纳米畴结构和极化织构的演化。

研究结果

梯度能量各向异性微扰模型的构建与验证

研究人员构建了一个傅里叶空间微扰模型，用于求解铁电材料中的各向异性梯度能量。通过在倒易空间构建各向异性梯度能量算子，该模型能够精确离散化并高效处理梯度能量系数的全各向异性张量及其空间变化。数值验证表明，所提出的微扰方案在求解具有空间变化梯度能量系数的相场方程时，表现出良好的收敛性和鲁棒性。

应用于铁电/反铁电复合材料的模拟结果

将该模型应用于铁电/反铁电复合材料体系，成功再现了反铁电相图以及在外加电场或梯度能量参数变化下的AFE-FE相变路径。模拟结果精确捕捉了不同梯度能量各向异性程度对相变行为和畴结构的影响，证明了该模型在描述复杂多相系统方面的有效性。

算法的可扩展性与三维模拟能力

所提出的算法展现出优越的可扩展性，能够无缝扩展到三维模拟区域。这为可视化复杂的纳米畴结构（如极性纳米区PNRs的形貌和排列）提供了可能，并揭示了极化织构在三维空间中的演化机制，这是以往基于简化梯度能量模型的模拟所难以实现的。

研究结论与讨论

本研究成功发展了一种基于傅里叶空间微扰的梯度能量各向异性模型，有效解决了铁电/反铁电复合材料相场模拟中长期存在的理论和技术瓶颈。该模型的核心优势在于能够同时处理具有空间变化的、完全各向异性的梯度能量张量，从而显著提高了对极化非均匀演化、相变路径和纳米畴结构预测的准确性。在铁电/反铁电复合材料体系中的应用表明，该模型不仅计算稳健、收敛高效，而且具备出色的三维扩展能力，能够揭示极化织构复杂的三维演化过程。

此项工作的意义重大。在理论层面，它弥补了现有相场模拟忽略梯度能量各向异性及空间变化的理论缺陷，为更精确地描述铁电/反铁电材料中的复杂物理现象提供了坚实的理论基础。在技术层面，所发展的傅里叶空间微扰算法为高性能、高保真的相场模拟提供了强大的计算工具，极大地提升了模拟的效率和尺度。最终，该模型为深入理解弛豫铁电体和反铁电体中纳米结构（如PNRs）与宏观性能（如能量存储）之间的构效关系提供了新途径，将有力推动下一代高性能储能介电材料的理性设计与开发。论文发表于《Materials》期刊。