- Open Access
Identification of mechanisms of magnetic transitions using an efficient method for converging on first-order saddle points
Phys. Rev. B 112, 104433 – Published 24 September, 2025
DOI: https://doi.org/10.1103/z673-hhnp
Abstract
An efficient and scalable implementation of a method for locating first-order saddle points on the energy surface of a magnetic system is presented, along with several applications in which the mechanisms of various magnetic transitions are identified. The starting point for the iterative search algorithm can be anywhere, even close to a local energy minimum representing an initial state of the system, and the final state need not be specified. Convergence on a saddle point is obtained by inverting the component of the gradient along the minimum mode, thereby effectively transforming the neighborhood of the saddle point to that of a local minimum. The method requires only the lowest two eigenvalues and corresponding eigenvectors of the Hessian of the system's energy and they are found using a quasi-Newton limited-memory Broyden-Fletcher-Goldfarb-Shanno solver for the minimization of the Rayleigh quotient without explicit evaluation of the Hessian. The method is applicable to large systems, as it does not introduce additional scaling overhead to the computational complexity determined by the interactions present in the system. Applications are presented to transitions in systems that reveal significant complexity of coexisting magnetic states, such as skyrmions, skyrmion bags, skyrmion tubes, chiral bobbers, and globules. The identification of new metastable three-dimensional (3D) textures, such as magnetic bobbers with extended equilibrium distance between the base and the terminating Bloch point, and magnetic globules appearing as isolated states in 3D due to magnetostatic interactions, demonstrates the usefulness of the method for the characterization of complex energy surfaces of magnetic systems. When combined with rate theory within the harmonic approximation, the method can be used for simulations of the long timescale dynamics of complex magnetic systems characterized by multiple metastable states.
Physics Subject Headings (PhySH)
Article Text
Supplemental Material
References (77)
- E. Wigner, The transition state method, Trans. Faraday Soc. 34, 29 (1938).
- G. H. Vineyard, Frequency factors and isotope effects in solid state rate processes, J. Phys. Chem. Solids 3, 121 (1957).
- H. A. Kramers, Brownian motion in a field of force and the diffusion model of chemical reactions, Physica (Amsterdam) 7, 284 (1940).
- J. Langer, Statistical theory of the decay of metastable states, Ann. Phys. 54, 258 (1969).
- D. V. Berkov, Magnetization dynamics including thermal fluctuations: Basic phenomenology, fast remagnetization processes and transitions over high-energy barriers, in Handbook of Magnetism and Advanced Magnetic Materials (Wiley, Hoboken, NJ, 2007), Vol. 2, pp. 795–823.
- H. Jónsson, G. Mills, and K. W. Jacobsen, Nudged elastic band method for finding minimum energy paths of transitions, in Classical and Quantum Dynamics in Condensed Phase Simulations, edited by B. J. Berne, G. Ciccotti, and D. F. Coker (World Scientific, Singapore, 1998), pp. 385–404.
- P. F. Bessarab, V. M. Uzdin, and H. Jónsson, Method for finding mechanism and activation energy of magnetic transitions, applied to skyrmion and antivortex annihilation, Comput. Phys. Commun. 196, 335 (2015).
- G. Henkelman and H. Jónsson, Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table, J. Chem. Phys. 115, 9657 (2001).
- H. Jónsson, Simulation of surface processes, Proc. Natl. Acad. Sci. USA 108, 944 (2011).
- M. Plasencia, A. Pedersen, A. Arnaldsson, J.-C. Berthet, and H. Jónsson, Geothermal model calibration using a global minimization algorithm based on finding saddle points and minima of the objective function, Comput. Geosci. 65, 110 (2014).
- D. M. Einarsdóttir, A. Arnaldsson, F. Óskarsson, and H. Jónsson, Path optimization with application to tunneling, in Applied Parallel and Scientific Computing, edited by K. Jónasson (Springer, Berlin, 2012), pp. 45–55.
- V. Ásgeirsson, A. Arnaldsson, and H. Jónsson, Efficient evaluation of atom tunneling combined with electronic structure calculations, J. Chem. Phys. 148, 102334 (2018).
- S. M. Vlasov, P. F. Bessarab, I. S. Lobanov, M. N. Potkina, V. M. Uzdin, and H. Jónsson, Magnetic skyrmion annihilation by quantum mechanical tunneling, New J. Phys. 22, 083013 (2020).
- I. A. Nosikov, M. V. Klimenko, G. A. Zhbankov, A. V. Podlesnyi, V. A. Ivanova, and P. F. Bessarab, Generalized force approach to point-to-point ionospheric ray tracing and systematic identification of high and low rays, IEEE Trans. Antennas Propag. 68, 455 (2020).
- B. Peters, Reaction Rate Theory and Rare Events (Elsevier, Amsterdam, 2017).
- G. Henkelman and H. Jónsson, A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives, J. Chem. Phys. 111, 7010 (1999).
- R. Malek and N. Mousseau, Dynamics of Lennard-Jones clusters: A characterization of the activation-relaxation technique, Phys. Rev. E 62, 7723 (2000).
- A. Pedersen, S. F. Hafstein, and H. Jónsson, Efficient sampling of saddle points with the minimum-mode following method, SIAM J. Sci. Comput. 33, 633 (2011).
- M. P. Gutiérrez, C. Argáez, and H. Jónsson, Improved minimum mode following method for finding first order saddle points, J. Chem. Theory Comput. 13, 125 (2017).
- C. Lanczos, An iteration method for the solution of the eigenvalue problem of linear differential and integral operators, J. Res. Natl. Bureau Stand. 45, 255 (1950).
- E. R. Davidson, The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices, J. Comput. Phys. 17, 87 (1975).
- W. T. Coffey, D. A. Garanin, and D. J. Mccarthy, Crossover formulas in the kramers theory of thermally activated escape rates–application to spin systems, in Advances in Chemical Physics (Wiley, Hoboken, NJ, 2001), pp. 483–765.
- P. F. Bessarab, V. M. Uzdin, and H. Jónsson, Harmonic transition-state theory of thermal spin transitions, Phys. Rev. B 85, 184409 (2012).
- P. F. Bessarab, V. M. Uzdin, and H. Jónsson, Potential energy surfaces and rates of spin transitions, Z. Phys. Chem. 227, 1543 (2013).
- I. S. Lobanov, H. Jónsson, and V. M. Uzdin, Mechanism and activation energy of magnetic skyrmion annihilation obtained from minimum energy path calculations, Phys. Rev. B 94, 174418 (2016).
- D. Cortés-Ortuño, W. Wang, M. Beg, R. A. Pepper, M.-A. Bisotti, R. Carey, M. Vousden, T. Kluyver, O. Hovorka, and H. Fangohr, Thermal stability and topological protection of skyrmions in nanotracks, Sci. Rep. 7, 4060 (2017).
- S. von Malottki, B. Dupé, P. F. Bessarab, A. Delin, and S. Heinze, Enhanced skyrmion stability due to exchange frustration, Sci. Rep. 7, 12299 (2017).
- P. F. Bessarab, G. P. Müller, I. S. Lobanov, F. N. Rybakov, N. S. Kiselev, H. Jónsson, V. M. Uzdin, S. Blügel, L. Bergqvist, and A. Delin, Lifetime of racetrack skyrmions, Sci. Rep. 8, 3433 (2018).
- B. Heil, A. Rosch, and J. Masell, Universality of annihilation barriers of large magnetic skyrmions in chiral and frustrated magnets, Phys. Rev. B 100, 134424 (2019).
- D. Cortés-Ortuño, N. Romming, M. Beg, K. von Bergmann, A. Kubetzka, O. Hovorka, H. Fangohr, and R. Wiesendanger, Nanoscale magnetic skyrmions and target states in confined geometries, Phys. Rev. B 99, 214408 (2019).
- H. Schrautzer, S. von Malottki, P. F. Bessarab, and S. Heinze, Effects of interlayer exchange on collapse mechanisms and stability of magnetic skyrmions, Phys. Rev. B 105, 014414 (2022).
- M. A. Goerzen, S. von Malottki, S. Meyer, P. F. Bessarab, and S. Heinze, Lifetime of coexisting sub-10 nm zero-field skyrmions and antiskyrmions, npj Quantum Mater. 8, 54 (2023).
- J. Hagemeister, A. Siemens, L. Rózsa, E. Y. Vedmedenko, and R. Wiesendanger, Controlled creation and stability of skyrmions on a discrete lattice, Phys. Rev. B 97, 174436 (2018).
- L. Desplat, J.-V. Kim, and R. L. Stamps, Paths to annihilation of first- and second-order (anti)skyrmions via (anti)meron nucleation on the frustrated square lattice, Phys. Rev. B 99, 174409 (2019).
- V. M. Kuchkin, B. Barton-Singer, F. N. Rybakov, S. Blügel, B. J. Schroers, and N. S. Kiselev, Magnetic skyrmions, chiral kinks, and holomorphic functions, Phys. Rev. B 102, 144422 (2020).
- V. M. Kuchkin, N. S. Kiselev, F. N. Rybakov, and P. F. Bessarab, Tailed skyrmions—An obscure branch of magnetic solitons, Front. Phys. 11, 1171079 (2023).
- F. N. Rybakov, A. B. Borisov, S. Blügel, and N. S. Kiselev, New type of stable particlelike states in chiral magnets, Phys. Rev. Lett. 115, 117201 (2015).
- G. P. Müller, F. N. Rybakov, H. Jónsson, S. Blügel, and N. S. Kiselev, Coupled quasimonopoles in chiral magnets, Phys. Rev. B 101, 184405 (2020).
- V. M. Kuchkin and N. S. Kiselev, Homotopy transitions and 3D magnetic solitons, APL Mater. 10, 071102 (2022).
- I. S. Lobanov and V. M. Uzdin, Lifetime, collapse, and escape paths for hopfions in bulk magnets with competing exchange interactions, Phys. Rev. B 107, 104405 (2023).
- M. Sallermann, H. Jónsson, and S. Blügel, Stability of hopfions in bulk magnets with competing exchange interactions, Phys. Rev. B 107, 104404 (2023).
- Y. Li, Y. Zang, R. Chen, and C. Moutafis, Tailoring energy barriers of bloch-point-mediated transitions between topological spin textures, Phys. Rev. B 109, 024418 (2024).
- Y. L. A. Schmerwitz, G. Levi, and H. Jónsson, Calculations of excited electronic states by converging on saddle points using generalized mode following, J. Chem. Theory Comput. 19, 3634 (2023).
- G. P. Müller, P. F. Bessarab, S. M. Vlasov, F. Lux, N. S. Kiselev, S. Blügel, V. M. Uzdin, and H. Jónsson, Duplication, collapse, and escape of magnetic skyrmions revealed using a systematic saddle point search method, Phys. Rev. Lett. 121, 197202 (2018).
- P. F. Bessarab, V. M. Uzdin, and H. Jónsson, Calculations of magnetic states and minimum energy paths of transitions using a noncollinear extension of the alexander-anderson model and a magnetic force theorem, Phys. Rev. B 89, 214424 (2014).
- A. V. Ivanov, D. Dagbartsson, J. Tranchida, V. M. Uzdin, and H. Jónsson, Efficient optimization method for finding minimum energy paths of magnetic transitions, J. Phys.: Condens. Matter 32, 345901 (2020).
- H. Bocquet and P. M. Derlet, Searching for activated transitions in complex magnetic systems, Phys. Rev. B 108, 174419 (2023).
- Intel math kernel library (2024.2.2) (2024).
- V. P. Antropov, M. I. Katsnelson, B. N. Harmon, M. van Schilfgaarde, and D. Kusnezov, Spin dynamics in magnets: Equation of motion and finite temperature effects, Phys. Rev. B 54, 1019 (1996).
- A. Edelman, T. A. Arias, and S. T. Smith, The geometry of algorithms with orthogonality constraints, SIAM J. Matrix Anal. Appl. 20, 303 (1998).
- A. S. Varentcova, S. von Malottki, M. N. Potkina, G. Kwiatkowski, S. Heinze, and P. F. Bessarab, Toward room-temperature nanoscale skyrmions in ultrathin films, npj Comput. Mater. 6, 193 (2020).
- K. Ohno and S. Maeda, A scaled hypersphere search method for the topography of reaction pathways on the potential energy surface, Chem. Phys. Lett. 384, 277 (2004).
- S. Maeda, Y. Watanabe, and K. Ohno, A scaled hypersphere interpolation technique for efficient construction of multidimensional potential energy surfaces, Chem. Phys. Lett. 414, 265 (2005).
- P.-A. Absil, R. Mahony, and R. Sepulchre, Optimization Algorithms on Matrix Manifolds (Princeton University Press, Princeton, NJ, 2008).
- D. C. Liu and J. Nocedal, On the limited memory BFGS method for large scale optimization, Math. Program. 45, 503 (1989).
- A. V. Ivanov, V. M. Uzdin, and H. Jónsson, Fast and robust algorithm for energy minimization of spin systems applied in an analysis of high temperature spin configurations in terms of skyrmion density, Comput. Phys. Commun. 260, 107749 (2021).
- G. L. G. Sleijpen and H. A. Van der Vorst, A Jacobi–Davidson iteration method for linear eigenvalue problems, SIAM Rev. 42, 267 (2000).
- G. W. Stewart, A Krylov–Schur algorithm for large eigenproblems, SIAM J. Matrix Anal. Appl. 23, 601 (2002).
- C. Brezinski, Some pioneers of extrapolation methods, in The Birth of Numerical Analysis (World Scientific, Singapoire, 2010), pp. 1–22.
- L. F. Richardson, IX. The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam, Philos. Trans. R. Soc. London, Ser. A 210, 307 (1911).
- L. F. Richardson and J. A. Gaunt, VIII. The deferred approach to the limit, Philos. Trans. R. Soc. London, Ser. A 226, 299 (1927).
- A. Pedersen, J.-C. Berthet, and H. Jónsson, Simulated annealing with coarse graining and distributed computing, in Applied Parallel and Scientific Computing, edited by K. Jónasson (Springer, Berlin, 2012), pp. 34–44.
- S. T. Chill, M. Welborn, R. Terrell, L. Zhang, J.-C. Berthet, A. Pedersen, H. Jónsson, and G. Henkelman, Eon: Software for long time simulations of atomic scale systems, Modell. Simul. Mater. Sci. Eng. 22, 055002 (2014).
- O. Eriksson, A. Bergman, L. Bergqvist, and J. Hellsvik, Atomistic Spin Dynamics: Foundations and Applications (Oxford University Press, Oxford, 2017).
- Q. Xu, I. P. Miranda, M. Pereiro, F. N. Rybakov, D. Thonig, E. Sjöqvist, P. F. Bessarab, A. Bergman, O. Eriksson, P. Herman, and A. Delin, Metaheuristic conditional neural network for harvesting skyrmionic metastable states, Phys. Rev. Res. 5, 043199 (2023).
- F. Zheng, F. N. Rybakov, A. B. Borisov, D. Song, S. Wang, Z.-A. Li, H. Du, N. S. Kiselev, J. Caron, A. Kovács, M. Tian, Y. Zhang, S. Blügel, and R. E. Dunin-Borkowski, Experimental observation of chiral magnetic bobbers in b20-type fege, Nat. Nanotechnol. 13, 451 (2018).
- J. W. Cooley and J. W. Tukey, An algorithm for the machine calculation of complex fourier series, Math. Comput. 19, 297 (1965).
- N. Hayashi, K. Saito, and Y. Nakatani, Calculation of demagnetizing field distribution based on fast fourier transform of convolution, Jpn. J. Appl. Phys. 35, 6065 (1996).
- H. J. G. Draaisma and W. J. M. de Jonge, Surface and volume anisotropy from dipole-dipole interactions in ultrathin ferromagnetic films, J. Appl. Phys. 64, 3610 (1988).
- V. L. Zhang, C. G. Hou, K. Di, H. S. Lim, S. C. Ng, S. D. Pollard, H. Yang, and M. H. Kuok, Eigenmodes of néel skyrmions in ultrathin magnetic films, AIP Adv. 7, 055212 (2017).
- L. Desplat, D. Suess, J.-V. Kim, and R. L. Stamps, Thermal stability of metastable magnetic skyrmions: Entropic narrowing and significance of internal eigenmodes, Phys. Rev. B 98, 134407 (2018).
- S. von Malottki, P. F. Bessarab, S. Haldar, A. Delin, and S. Heinze, Skyrmion lifetime in ultrathin films, Phys. Rev. B 99, 060409(R) (2019).
- See Supplemental Material at http://link.aps.org/supplemental/10.1103/z673-hhnp for details of the calculations illustrating the breakdown of the skyrmion escape mechanism.
- V. M. Kuchkin, N. S. Kiselev, A. Haller, I. C. V. Liščák, A. Michels, and T. L. Schmidt, Stability and nucleation of dipole strings in uniaxial chiral magnets, Phys. Rev. B 111, 174410 (2025).
- G. P. Müller, M. Hoffmann, C. Dißelkamp, D. Schürhoff, S. Mavros, M. Sallermann, N. S. Kiselev, H. Jónsson, and S. Blügel, Spirit: Multifunctional framework for atomistic spin simulations, Phys. Rev. B 99, 224414 (2019).
- V. Ásgeirsson, B. O. Birgisson, R. Bjornsson, U. Becker, F. Neese, C. Riplinger, and H. Jónsson, Nudged elastic band method for molecular reactions using energy-weighted springs combined with eigenvector following, J. Chem. Theory Comput. 17, 4929 (2021).
- G. Fiedler, J. Fidler, J. Lee, T. Schrefl, R. L. Stamps, H. B. Braun, and D. Suess, Direct calculation of the attempt frequency of magnetic structures using the finite element method, J. Appl. Phys. 111, 093917 (2012).