Improved Energy Stable Symmetric Gauss-Seidel Projection Method for Micromagnetics Simulations

技术摘要：用于微磁学模拟的改进型能量稳定对称高斯-赛德尔投影法

问题陈述
由兰道-利夫希茨-吉尔伯特（LLG）方程控制的微磁学模拟，对于研究铁磁材料中的磁化动力学至关重要。这些模拟面临着显著的数值挑战，主要源于方程的刚性、非线性和磁化矢量严格的单位长度约束。虽然高斯-赛德尔投影法（GSPM）已成为一种鲁棒的半隐式框架，具有无条件稳定性和计算效率，但它也存在固有的局限性。具体而言，标准的 GSPM 依赖于单向（单侧）高斯-赛德尔迭代。这种方法会引入方向性的数值偏差，并可能导致能量稳定性下降，特别是在具有超低吉尔伯特阻尼参数（$\alpha \approx 10^{-3} \sim 10^{-4}$）的系统中，这类参数在钇铁石榴石（YIG）和外延 FeCo 合金等高质量磁性材料中十分常见。在无阻尼（$\alpha=0$）或弱阻尼机制下，标准 GSPM 可能无法保证离散能量稳定性，从而导致非物理性的能量增长或振荡。

方法论
作者提出了一种对称高斯-赛德尔投影法（SGSPM），以解决标准 GSPM 的稳定性缺陷。该方法保留了 GSPM 的核心结构——包括分数步分解、流形投影以及用于求解线性系统的快速傅里叶变换（FFT）技术——但修改了其迭代求解器。

1. 对称迭代： 与执行单向前向高斯-赛德尔扫描（$1 \to n$）的标准 GSPM 不同，SGSPM 采用双向对称迭代。这包括一个前向扫描紧接着一个后向扫描（$n \to 1$）。
2. 算法步骤：
   • 步骤 1（隐式对称高斯-赛德尔）： 使用对称迭代过程更新磁化场，该过程在正向和反向方向上充分利用热扩散子问题中最新的更新状态变量。
   • 步骤 2（无约束热流）： 执行一个无约束的热流步骤，以处理阻尼项和有效场项。
   • 步骤 3（投影）： 通过流形投影步骤强制执行磁化矢量的单位长度约束（$|m|=1$）。
3. 理论基础： 本文证明，对于严格对角占优矩阵，对称迭代矩阵 $G_{SGS}$ 满足 $\rho(G_{SGS}) \leq \rho(G_{GS}) < 1$，从理论上确保了其相比于单向方案具有更优越的收敛性和稳定性特性。

核心贡献

• 开发 SGSPM： 主要贡献在于 SGSPM 的构建，这是 GSPM 的一种变体，它在保持原 GSPM 一阶时间精度和二阶空间精度的同时，显著增强了稳定性。
• 离散能量稳定性： 该方法严格保证了离散能量稳定性。作者证明，对称迭代消除了标准高斯-赛德尔扫描中固有的方向性偏差，而这种偏差正是弱阻尼系统中能量不稳定的根源。
• 低阻尼机制下的鲁棒性： 该方案专门设计用于处理超低阻尼参数（$\alpha \to 0$）的情况，在这一机制下，标准显式或半隐式方法往往会失效，或者需要极小的步长。

数值结果
通过涉及 1D 和 3D 设置下的磁化动力学演化及畴壁传播的广泛数值实验，验证了 SGSPM 的性能。

• 能量稳定性： 在零阻尼（$\alpha=0$）模拟中，标准 GSPM 表现出非物理性的能量增长和振荡。相比之下，SGSPM 展示了严格的单调能量衰减，符合物理能量耗散定律。
• 阻尼敏感性： 对于正阻尼值（$\alpha > 0$），两种方法都产生了耗散能量曲线。然而，与 GSPM 相比，SGSPM 始终实现了更快的能量衰减速率和更低的残余能量水平。随着阻尼增加，两者的性能差距缩小，但 SGSPM 依然保持领先。
• 精度与效率： 数值测试确认，SGSPM 保留了一阶时间精度和二阶空间精度。虽然由于增加了后向扫描，SGSPM 的计算成本略高于 GSPM，但通过使用 FFT 求解底层线性系统，保持了整体的高效率。
• 畴壁动力学： 在外部磁场下的 Néel 畴壁运动模拟中，SGSPM 产生的磁化解与 GSPM 相比在数值一致性和视觉匹配度上均表现良好，证实了稳定性提升并未损害动力学的物理保真度。

意义与主张
本文声称，SGSPM 为微磁学模拟提供了一个更鲁棒的数值框架，特别是针对涉及超低阻尼材料的现代应用。通过采用对称高斯-赛德尔策略，该方法有效地抑制了残余振荡和方向性数值偏差，确保即使在缺乏吉尔伯特阻尼的情况下，离散能量演化仍保持物理有效。作者认为，这种改进使得 SGSPM 成为长时动态模拟和高频自旋波问题中更优的选择，因为在这些场景下能量稳定性至关重要。此外，作者指出，这种对称策略可以扩展到更高阶的 GSPM 变体和其他改进方案，以进一步增强在弱阻尼条件下的稳定性。