Searching the solution landscape by generalized high-index saddle dynamics

Abstract

We introduce a generalized numerical algorithm to construct the solution landscape, which is a pathway map consisting of all stationary points and their connections. Based on the high-index optimization-based shrinking dimer (HiOSD) method for gradient systems, a generalized high-index saddle dynamics (GHiSD) is proposed to compute any-index saddles of dynamical systems. Linear stability of the index-$k$ saddle point can be proved for the GHiSD system. A combination of the downward search algorithm and the upward search algorithm is applied to systematically construct the solution landscape, which not only provides a powerful and efficient way to compute multiple solutions without tuning initial guesses, but also reveals the relationships between different solutions. Numerical examples, including a three-dimensional example and the phase field model, demonstrate the novel concept of the solution landscape by showing the connected pathway maps.

References in corpus (9)

• Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms
• String Method for the Study of Rare Events
• Essentially No Barriers in Neural Network Energy Landscape
• The Gentlest Ascent Dynamics
• Construction of a pathway map on a complicated energy landscape
• A Comparative Analysis of the Optimization and Generalization Property of Two-layer Neural Network and Random Feature Models Under Gradient Descent Dynamics
• Finding All the Stationary Points of a Potential Energy Landscape via Numerical Polynomial Homotopy Continuation Method
• Pathways connecting two opposed bilayers with a fusion pore: A molecularly-informed phase field approach
• An Inversion-Relaxation Approach for Sampling Stationary Points of Spin Model Hamiltonians

Cited by in corpus (13)

• Constrained high-index saddle dynamics for the solution landscape with equality constraints
• Solution landscape of the Onsager model identifies non-axisymmetric critical points
• Computing solution landscape of nonlinear space-fractional problems via fast approximation algorithm
• Solution landscapes of the diblock copolymer-homopolymer model under two-dimensional confinement
• Solution Landscapes of the Simplified Ericksen--Leslie Model and its Comparison with the Reduced Landau--de Gennes Model
• Hierarchies of Critical Points of a Landau-de Gennes Free Energy on Three-Dimensional Cuboids
• Normalized Wolfe-Powell-type local minimax method for finding multiple unstable solutions of nonlinear elliptic PDEs
• Transition pathways in Cylinder-Gyroid interface
• Spectral Deferred Correction Method for Landau-Brazovskii Model with Convex Splitting Technique
• Companion-Based Multi-Level Finite Element Method for Computing Multiple Solutions of Nonlinear Differential Equations
• Modeling and Computation of Liquid Crystals
• Numerical Discretization of Variational Phase Field Model for Phase Transitions in Ferroelectric Thin Films
• Construction of a minimum energy path for the VT flash model by an exponential time differencing scheme with the string method