[论文解读] Convex Splitting Method for the Calculation of Transition States of Energy Functional

Shuting Gu, Xiang Zhou|arXiv (Cornell University)|Oct 23, 2016

Block Copolymer Self-Assembly参考文献 26被引用 3

一句话总结

该论文在迭代最小化公式（IMF）中引入了一种凸分裂方法，以高效计算能量泛函的指标-1鞍点（过渡态），例如Allen-Cahn和Cahn-Hilliard系统中的情况。通过在每个IMF循环中对辅助泛函应用凸分裂，该方法确保了无条件能量稳定性，并支持大时间步长，显著提高了1D和2D模型（包括Landau-Brazovskii泛函）的计算效率。

ABSTRACT

Among numerical methods for partial differential equations arising from steepest descent dynamics of energy functionals (e.g., Allen-Cahn and Cahn-Hilliard equations), the convex splitting method is well-known to maintain unconditional energy stability for a large time step size. In this work, we show how to use the convex splitting idea to find transition states, i.e., index-1 saddle points of the same energy functionals. Based on the iterative minimization formulation (IMF) for saddle points (SIAM J. Numer. Anal., vol. 53, p1786, 2015), we introduce the convex splitting method to minimize the auxiliary functional at each cycle of the IMF. We present a general principle of constructing convex splitting forms for these auxiliary functionals and show how to avoid solving nonlinear equations. The new numerical scheme based on the convex splitting method allows for large time step sizes. The new methods are tested for the one dimensional Ginzburg-Landau energy functional in the search of the Allen-Cahn or Cahn-Hilliard types of transition states. We provide the numerical results of transition states for the two dimensional Landau-Brazovskii energy functional for diblock copolymers.

研究动机与目标

• 开发一种高效的数值方法，用于计算空间扩展系统中能量泛函的指标-1鞍点（过渡态）。
• 将已知在梯度流模拟中具有能量稳定性的凸分裂方法，适配到迭代最小化公式（IMF）中，用于鞍点计算。
• 在IMF中求解辅助最小化子问题时，实现大时间步长，克服标准方案的稳定性限制。
• 在$H^{-1}$和$L^2$度量下，针对1D Ginzburg-Landau和2D Landau-Brazovskii能量泛函，展示该方法的鲁棒性和高效性。
• 通过强制最小模式与状态梯度的正交性，解决由于平移模式导致的虚假收敛等实际挑战。

提出的方法

• 该方法将凸分裂方法整合到迭代最小化公式（IMF）中，其中每个IMF循环通过辅助泛函的凸分裂求解子问题。
• 提出了一般性原则，将辅助泛函$L$分解为凸部分和凹部分，确保大时间步长下的无条件能量稳定性。
• 凸分裂格式通过使用半隐式、线性化格式，避免在每个时间步长求解非线性方程，从而提高计算效率。
• 该方法确保用于辅助泛函的最小模式与空间梯度$\partial_x \phi$正交，防止收敛到平凡的平移模式。
• 在子问题中应用自适应停止准则，以平衡精度与计算成本，同时保持IMF的二次收敛速率。
• 该方法在1D Ginzburg-Landau和2D Landau-Brazovskii模型上得到验证，数值结果表明即使在大时间步长下也能实现稳定收敛。

实验结果

研究问题

• RQ1 凸分裂方法能否有效适配到迭代最小化公式（IMF）中，用于计算指标-1鞍点？
• RQ2 在IMF中对辅助泛函采用凸分裂方法是否能确保无条件能量稳定性并支持大时间步长？
• RQ3 该方法在避免因系统平移不变性导致的虚假解收敛方面表现如何？
• RQ4 与标准方案相比，IMF中使用凸分裂方法在过渡态计算中的计算效率提升如何？
• RQ5 该方法能否推广到其他能量泛函（如二嵌段共聚物的Landau-Brazovskii模型）？

主要发现

• 在IMF框架内应用凸分裂方法，实现了无条件能量稳定的时间积分，支持大时间步长而不损失稳定性。
• 对于1D Ginzburg-Landau能量，该方法分别在$L^2$和$H^{-1}$度量下成功计算了Allen-Cahn和Cahn-Hilliard类型的过渡态。
• 在2D Landau-Brazovskii模型中，即使采用大时间步长，该方法仍能收敛到稳定的过渡态，显示出在高维情况下的鲁棒性。
• 通过强制最小模式与状态空间梯度的正交性，该方法避免了收敛到平凡的平移模式，确保了物理上有意义的鞍点搜索。
• 大量数值实验表明，该方法保持了二次收敛性，并且与标准方案相比显著提高了计算效率。
• 该方法提供了一个通用框架，可将现有的梯度流能量稳定方案重新用于高效鞍点计算。