Optimal-order error estimates for Euler discretization of high-index saddle dynamics

. High-index saddle dynamics provides an effective means to compute the any-index saddle points and construct the solution landscape. In this paper we prove error estimates for Euler discretization of high-index saddle dynamics with respect to the time step size, which remains untreated in the literature. We overcome the main difficulties that lie in the strong nonlinearity of the saddle dynamics and the orthonormalization procedure in the numerical scheme that is uncommon in standard discretization of differential equations. The derived methods are further extended to study the generalized high-index saddle dynamics for non-gradient systems and provide theoretical support for the accuracy of numerical implementations.

[1]  Zhen Huang,et al.  Constrained High-Index Saddle Dynamics for the Solution Landscape with Equality Constraints , 2020, J. Sci. Comput..

[2]  Yucheng Hu,et al.  Solution landscapes of the simplified Ericksen–Leslie model and its comparisonwith the reduced Landau–deGennes model , 2021, Proceedings of the Royal Society A.

[3]  Wei Wang,et al.  Modelling and computation of liquid crystals , 2021, Acta Numerica.

[4]  Y. Nishiura,et al.  Solution landscapes of the diblock copolymer-homopolymer model under two-dimensional confinement. , 2021, Physical review. E.

[5]  Pingwen Zhang,et al.  Transition pathways connecting crystals and quasicrystals , 2020, Proceedings of the National Academy of Sciences.

[6]  Pingwen Zhang,et al.  Solution landscape of a reduced Landau–de Gennes model on a hexagon , 2020, Nonlinearity.

[7]  Jianyuan Yin,et al.  Solution landscape of Onsager functional identifies non-axisymmetric critical points , 2021, ArXiv.

[8]  Lei Zhang,et al.  Searching the solution landscape by generalized high-index saddle dynamics , 2020, Science China Mathematics.

[9]  Pingwen Zhang,et al.  Construction of a Pathway Map on a Complicated Energy Landscape. , 2019, Physical review letters.

[10]  Pingwen Zhang,et al.  Transition pathways between defect patterns in confined nematic liquid crystals , 2019, J. Comput. Phys..

[11]  A. Shi,et al.  Pathways connecting two opposed bilayers with a fusion pore: a molecularly-informed phase field approach. , 2019, Soft matter.

[12]  Pingwen Zhang,et al.  High-Index Optimization-Based Shrinking Dimer Method for Finding High-Index Saddle Points , 2019, SIAM J. Sci. Comput..

[13]  Tong Zhang,et al.  Near-optimal stochastic approximation for online principal component estimation , 2016, Math. Program..

[14]  Christoph Ortner,et al.  Convergence and Cycling in Walker-type Saddle Search Algorithms , 2016, SIAM J. Numer. Anal..

[15]  Lei Zhang,et al.  Recent developments in computational modelling of nucleation in phase transformations , 2016 .

[16]  Nicholas I. M. Gould,et al.  A dimer-type saddle search algorithm with preconditioning and linesearch , 2016, Math. Comput..

[17]  Qiang Du,et al.  Optimization-based Shrinking Dimer Method for Finding Transition States , 2016, SIAM J. Sci. Comput..

[18]  Guigang Zhang,et al.  Deep Learning , 2016, Int. J. Semantic Comput..

[19]  Simon W. Funke,et al.  Deflation Techniques for Finding Distinct Solutions of Nonlinear Partial Differential Equations , 2014, SIAM J. Sci. Comput..

[20]  Xiang Zhou,et al.  An Iterative Minimization Formulation for Saddle Point Search , 2014, SIAM J. Numer. Anal..

[21]  Christina L. Ting,et al.  Nucleation in polymers and soft matter. , 2014, Annual review of physical chemistry.

[22]  Qiang Du,et al.  Shrinking Dimer Dynamics and Its Applications to Saddle Point Search , 2012, SIAM J. Numer. Anal..

[23]  Dhagash Mehta,et al.  Finding all the stationary points of a potential-energy landscape via numerical polynomial-homotopy-continuation method. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[24]  E Weinan,et al.  The gentlest ascent dynamics , 2010, 1011.0042.

[25]  Eric Vanden-Eijnden,et al.  Transition-path theory and path-finding algorithms for the study of rare events. , 2010, Annual review of physical chemistry.

[26]  Qiang Du,et al.  Simultaneous prediction of morphologies of a critical nucleus and an equilibrium precipitate in solids , 2009 .

[27]  Qiang Du,et al.  Morphology of critical nuclei in solid-state phase transformations. , 2007, Physical review letters.

[28]  Yurii Nesterov,et al.  Introductory Lectures on Convex Optimization - A Basic Course , 2014, Applied Optimization.

[29]  J. Doye,et al.  Saddle Points and Dynamics of Lennard-Jones Clusters, Solids and Supercooled Liquids , 2001, cond-mat/0108310.

[30]  Yongxin Li,et al.  A Minimax Method for Finding Multiple Critical Points and Its Applications to Semilinear PDEs , 2001, SIAM J. Sci. Comput..

[31]  G. Henkelman,et al.  A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives , 1999 .

[32]  Jesse Freeman,et al.  in Morse theory, , 1999 .

[33]  P. Schleyer,et al.  Internal conrotation and disrotation in H2BCH2BH2 and diborylmethane 1,3 H exchange , 1997 .

[34]  B. Eckhardt.  Irregular scattering , 1988 .

[35]  J. Baker.  An algorithm for the location of transition states , 1986 .