Error Estimate for Semi-implicit Method of Sphere-Constrained High-Index Saddle Dynamics

Abstract

The authors prove error estimates for the semi-implicit numerical scheme of sphere-constrained high-index saddle dynamics, which serves as a powerful instrument in finding saddle points and constructing the solution landscapes of constrained systems on the high-dimensional sphere. Due to the semi-implicit treatment and the novel computational procedure, the orthonormality of numerical solutions at each time step could not be fully employed to simplify the derivations, and the computations of the state variable and directional vectors are coupled with the retraction, the vector transport and the orthonormalization procedure, which significantly complicates the analysis. They address these issues to prove error estimates for the proposed semi-implicit scheme and then carry out numerical experiments to substantiate the theoretical findings.

Keywords

Saddle point / Constrained saddle dynamics / Solution landscape / Semi-implicit / Numerical analysis

References

[1] Baker J. An algorithm for the location of transition states. J. Comput. Chew., 1986, 7: 385-395

[2] Bao W, Cai Y. Mathematical theory and numerical methods for Bose-Einstein condensation. Kinet. Relat. Models, 2013, 6: 1-135

[3] Doye J, Wales D. Saddle points and dynamics of Lennard-Jones clusters, solids, and supercooled liquids. J. Chem. Phys., 2002, 116: 3777-3788

[4] E W, Vanden-Eijnden E. Transition-path theory and path-finding algorithms for the study of rare events. Annu. Rev. Phys. Chem., 2010, 61: 391-420

[5] E W, Zhou X. The gentlest ascent dynamics. Nonlinearity, 2011, 24: 1831-1842

[6] Farrell P, Birkisson Á, Funke S. Deflation techniques for finding distinct solutions of nonlinear partial differential equations. SIAM J. Sci. Comput., 2015, 37: A2026-A2045

[7] Gao W, Leng J, Zhou X. An iterative minimization formulation for saddle point search. SIAM J. Numer. Anal, 2015, 53: 1786-1805

[8] Gould N, Ortner C, Packwood D. A dimer-type saddle search algorithm with preconditioning and linesearch. Math. Comp., 2016, 85: 2939-2966

[9] Grantham W. Gradient transformation trajectory following algorithms for determining stationary min-max saddle points. Advances in Dynamic Game Theory, 2007, Boston, MA: Birkhauser Boston 639-657 9

[10] Han Y, Hu Y, Zhang P Transition pathways between defect patterns in confined nematic liquid crystals. J. Comput. Phys., 2019, 396: 1-11

[11] Han Y, Xu Z, Shi A, Zhang L. Pathways connecting two opposed bilayers with a fusion pore: A molecularly-informed phase field approach. Soft Matter, 2020, 16: 366-374

[12] Henkelman G, Jonsson H. A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives. J. Chew. Phys., 1999, 111: 7010-7022

[13] Levitt A, Ortner C. Convergence and cycling in walker-type saddle search algorithms. SIAM J. Numer. Anal, 2017, 55: 2204-2227

[14] Li Y, Zhou J. A minimax method for finding multiple critical points and its applications to semilinear PDEs. SIAM J. Sci. Comput., 2001, 23: 840-865

[15] Luo, Y., Zhang, L., Zhang, P., et al., Semi-implicit method of high-index saddle dynamics and application to construct solution landscape, Submitted.

[16] Mehta D. Finding all the stationary points of a potential-energy landscape via numerical polynomial-homotopy-continuation method. Phys. Rev. E, 2011, 84: 025702

[17] Milnor J. Morse Theory. Ann. of Math. Stud., 1963, NJ: Princeton University Press 51

[18] Thomson J. XXIV, On the structure of the atom: An investigation of the stability and periods of oscillation of a number of corpuscles arranged at equal intervals around the circumference of a circle; with application of the results to the theory of atomic structure. London, Edinburgh, Dublin Phil. Mag. J. Sci., 1904, 7: 237-265

[19] Wang W, Zhang L, Zhang P. Modelling and computation of liquid crystals. Acta Numerica, 2021, 30: 765-851

[20] Xie Z, Yuan Y, Zhou J. On solving semilinear singularly perturbed Neumann problems for multiple solutions. SIAM J. Sci. Comput., 2022, 44: A501-A523

[21] Yin J, Huang Z, Zhang L. Constrained high-index saddle dynamics for the solution landscape with equality constraints. J. Sci. Comput., 2022, 91: 62

[22] Yin J, Jiang K, Shi A Transition pathways connecting crystals and quasicrystals. Proc. Natl. Acad. Sci. U.S.A., 2021, 118: e2106230118

[23] Yin J, Wang Y, Chen J Construction of a pathway map on a complicated energy landscape. Phys. Rev. Lett., 2020, 124: 090601

[24] Yin J, Yu B, Zhang L. Searching the solution landscape by generalized high-index saddle dynamics. Sci. China Math., 2021, 64: 1801

[25] Yin J, Zhang L, Zhang P. High-index optimization-based shrinking dimer method for finding high-index saddle points. SIAM J. Sci. Comput., 2019, 41: A3576-A3595

[26] Yu B, Zheng X, Zhang P, Zhang L. Computing solution landscape of nonlinear space-fractional problems via fast approximation algorithm. J. Comput. Phys., 2022, 468: 111513

[27] Zhang L, Chen L, Du Q. Morphology of critical nuclei in solid-state phase transformations. Phys. Rev. Lett., 2007, 98: 265703

[28] Zhang L, Chen L, Du Q. Simultaneous prediction of morphologies of a critical nucleus and an equilibrium precipitate in solids. Commun. Comput. Phys., 2010, 7: 674-682

[29] Zhang J, Du Q. Shrinking dimer dynamics and its applications to saddle point search. SIAM J. Numer. Anal, 2012, 50: 1899-1921

[30] Zhang L, Ren W, Samanta A, Du Q. Recent developments in computational modelling of nucleation in phase transformations. npj Comput. Mater., 2016, 2: 16003

[31] Zhang L, Zhang P, Zheng X. Error estimates of Euler discretization to high-index saddle dynamics. SIAM J. Numer. Anal., 2022, 60: 2925-2944

[32] Zhang L, Zhang P, Zheng X. Mathematical and numerical analysis to shrinking-dimer saddle dynamics with local Lipschitz conditions. CSIAM Trans. Appl. Math., 2023, 4: 157-176

[33] Zhang, L., Zhang, P. and Zheng, X., Discretization and index-robust error analysis for constrained high-index saddle dynamics on high-dimensional sphere, Sci. China Math., 2023, DOI: https://doi.org/10.1007/s11425-022-2149-2.