Jianyuan Yin

Jianyuan Yin, Lei Zhang, Zhiyi Zhang

We propose a subspace curvature-scaling high-index saddle dynamics (SCS-HiSD) method to accelerate high-index saddle dynamics (HiSD) for locating ill-conditioned saddle points. The key observation is that HiSD already computes approximations of the unstable Hessian eigenvectors during iteration, which can be used at negligible additional cost to construct an inverse-Hessian approximation on the unstable subspace. This subspace curvature information is incorporated to adaptively scale the dynamics along each unstable direction, eliminating the dependence of the convergence rate on the smallest-magnitude negative eigenvalues and thereby substantially accelerating the convergence for ill-conditioned saddle points. We establish the linear stability of the continuous SCS-HiSD system and provide a local convergence analysis for the discrete iterative scheme. This method extends naturally to address slow convergence caused by small positive eigenvalues. Numerical experiments on benchmark problems and a liquid-crystal model demonstrate that SCS-HiSD substantially accelerates the computation of ill-conditioned saddle points, particularly in severely ill-conditioned cases.

Yuchen Xie, Jianyuan Yin, Lei Zhang

The discovery of ordered structures in pattern-forming systems, such as the Landau-Brazovskii (LB) model, is often limited by the sensitivity of numerical solvers to the prescribed computational domain size. Incompatible domains induce artificial stress, frequently trapping the system in high-energy metastable configurations. To resolve this issue, we propose a Geometry-Adaptive Deep Variational Framework (GeoDVF) that jointly optimizes the infinite-dimensional order parameter, which is parameterized by a neural network, and the finite-dimensional geometric parameters of the computational domain. By explicitly treating the domain size as trainable variables within the variational formulation, GeoDVF naturally eliminates artificial stress during training. To escape the attraction basin of the disordered phase under small initializations, we introduce a warmup penalty mechanism, which effectively destabilizes the disordered phase, enabling the spontaneous nucleation of complex three-dimensional ordered phases from random initializations. Furthermore, we design a guided initialization protocol to resolve topologically intricate phases associated with narrow basins of attraction. Extensive numerical experiments show that GeoDVF provides a robust and geometry-consistent variational solver capable of identifying both stable and metastable states without prior knowledge.