{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,14]],"date-time":"2026-05-14T15:17:31Z","timestamp":1778771851608,"version":"3.51.4"},"reference-count":47,"publisher":"Springer Science and Business Media LLC","issue":"3","license":[{"start":{"date-parts":[[2026,5,14]],"date-time":"2026-05-14T00:00:00Z","timestamp":1778716800000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2026,5,14]],"date-time":"2026-05-14T00:00:00Z","timestamp":1778716800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"funder":[{"name":"Hong Kong RGC","award":["16302223"],"award-info":[{"award-number":["16302223"]}]},{"name":"Hong Kong RGC","award":["16300524"],"award-info":[{"award-number":["16300524"]}]}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["J Sci Comput"],"published-print":{"date-parts":[[2026,6]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>This paper introduces a class of high-order numerical schemes for solving differential equations on the unit sphere, called Spherical Runge-Kutta methods with Richardson Extrapolation (SRKRE). Traditionally, Richardson Extrapolation (RE) is applied in Cartesian spaces through simple linear combinations of numerical solutions. Our approach adapts this idea directly to spherical geometry by extending the Spherical Linear Interpolation (SLERP) operator to perform extrapolation on the sphere itself. This results in a geometry-preserving formulation of RE that guarantees all intermediate and final solutions lie exactly on the unit sphere, without the need for any additional projection. The SRKRE schemes are constructed by combining existing low-order spherical integrators with this intrinsic extrapolation to achieve higher-order accuracy. Furthermore, we introduce a generalized counterpart, SRKREg, which extends the framework by allowing different choices of reference points on the sphere. SRKREg not only underscores the geometric foundations of SRKRE but also provides valuable insight into its accuracy properties. We analyze the stability and computational complexity of these methods and demonstrate their superior performance in both accuracy and structure preservation. Numerical experiments confirm that SRKRE and SRKREg achieve the expected order of convergence while fully honoring the geometric constraints of the problem.<\/jats:p>","DOI":"10.1007\/s10915-026-03318-4","type":"journal-article","created":{"date-parts":[[2026,5,14]],"date-time":"2026-05-14T15:06:11Z","timestamp":1778771171000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["High-Order Spherical Runge-Kutta Schemes Based on Spherical Linear Extrapolation"],"prefix":"10.1007","volume":"107","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-2599-3892","authenticated-orcid":false,"given":"Young","family":"Kyu Lee","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Shingyu","family":"Leung","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"297","published-online":{"date-parts":[[2026,5,14]]},"reference":[{"key":"3318_CR1","doi-asserted-by":"publisher","first-page":"611","DOI":"10.1007\/BF01211069","volume":"104","author":"SL Adler","year":"1986","unstructured":"Adler, S.L.: Quaternionic Quantum Field Theory. Commun. Math. Phys. 104, 611\u2013656 (1986)","journal-title":"Commun. Math. Phys."},{"key":"3318_CR2","doi-asserted-by":"publisher","first-page":"1939","DOI":"10.1002\/mma.651","volume":"28","author":"L Banas","year":"2005","unstructured":"Banas, L.: A numerical method for the Landau-Lifshitz equation with magnetostriction. Math. Meth. Appl. Sci. 28, 1939\u20131954 (2005)","journal-title":"Math. Meth. Appl. Sci."},{"key":"3318_CR3","doi-asserted-by":"crossref","unstructured":"Bayleyegn, T., Farago, I., Havasi, A.: On the convergence of multiple Richardson extrapolation combined with explicit Runge\u2013Kutta methods. Periodica Mathematica Hungarica, 88(335-353), (2024)","DOI":"10.1007\/s10998-023-00557-y"},{"key":"3318_CR4","doi-asserted-by":"crossref","unstructured":"Cai, Y., Chen, J., Wang, C., Xie, C.; A second-order numerical method for Landau-Lifshitz-Gilbert equation with large damping parameters. J. Comput. Phys., 451(110831), (2022)","DOI":"10.1016\/j.jcp.2021.110831"},{"issue":"3","key":"3318_CR5","doi-asserted-by":"publisher","first-page":"868","DOI":"10.1137\/04061979X","volume":"28","author":"M Calvo","year":"2006","unstructured":"Calvo, M., Hernandez-Abreu, D., Montijano, J.I., Randez, L.: On the preservation of invariants by explicit Runge-Kutta methods. SIAM J. Sci. Comput. 28(3), 868\u2013885 (2006)","journal-title":"SIAM J. Sci. Comput."},{"key":"3318_CR6","doi-asserted-by":"publisher","first-page":"189","DOI":"10.1016\/S0168-9274(99)00127-0","volume":"34","author":"RPK Chan","year":"2000","unstructured":"Chan, R.P.K., Murua, A.: Extrapolation of symplectic methods for Hamiltonian problems. Appl. Num. Math. 34, 189\u2013205 (2000)","journal-title":"Appl. Num. Math."},{"key":"3318_CR7","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1007\/BF03024947","volume":"15","author":"I Cimr\u00e1k","year":"2007","unstructured":"Cimr\u00e1k, I.: A Survey on the numerics and computations for the Landau-Lifshitz equation of micromagnetism. Archives of Computational Methods in Engineering 15, 1\u201337 (2007)","journal-title":"Archives of Computational Methods in Engineering"},{"key":"3318_CR8","doi-asserted-by":"crossref","unstructured":"E, W., Wang, X.P.: Numerical Methods for the Landau-Lifshitz Equation. SIAM J. Numer. Anal., 38(5):1647\u20131665, (2000)","DOI":"10.1137\/S0036142999352199"},{"key":"3318_CR9","doi-asserted-by":"publisher","first-page":"160","DOI":"10.1006\/jcph.1997.5672","volume":"133","author":"J Frank","year":"1997","unstructured":"Frank, J., Huang, W., Leimkuhler, B.: Geometrical integrators for classical spin systems. J. Comput. Phys. 133, 160\u2013172 (1997)","journal-title":"J. Comput. Phys."},{"issue":"8","key":"3318_CR10","doi-asserted-by":"publisher","first-page":"1969","DOI":"10.1088\/0951-7715\/19\/8\/011","volume":"19","author":"JD Gibbon","year":"2006","unstructured":"Gibbon, J.D., Holm, D.D., Kerr, R.M., Roulstone, I.: Quaternions and particle dynamics in the Euler fluid equations. Nonlinearity 19(8), 1969\u20131983 (2006)","journal-title":"Nonlinearity"},{"issue":"1","key":"3318_CR11","doi-asserted-by":"publisher","first-page":"84","DOI":"10.1137\/080726926","volume":"2","author":"D Goldfarb","year":"2009","unstructured":"Goldfarb, D., Wen, Z., Yin, W.: A Curvilinear Search Method for $$p$$-Harmonic Flows on Spheres. SIAM J. Imaging Sciences 2(1), 84\u2013109 (2009)","journal-title":"SIAM J. Imaging Sciences"},{"key":"3318_CR12","doi-asserted-by":"publisher","first-page":"73","DOI":"10.1090\/S0025-5718-98-00913-2","volume":"67","author":"S Gottlieb","year":"1998","unstructured":"Gottlieb, S., Shu, C.-W.: Total variation diminishing Runge-Kutta schemes. Math. Comput. 67, 73\u201385 (1998)","journal-title":"Math. Comput."},{"key":"3318_CR13","doi-asserted-by":"publisher","first-page":"89","DOI":"10.1137\/S003614450036757X","volume":"43","author":"S Gottlieb","year":"2000","unstructured":"Gottlieb, S., Shu, C.W., Tadmor, E.: Strong stability preserving high order time discretization methods. SIAM Rev. 43, 89\u2013112 (2000)","journal-title":"SIAM Rev."},{"issue":"2","key":"3318_CR14","doi-asserted-by":"publisher","first-page":"174","DOI":"10.1093\/imamat\/4.2.174","volume":"4","author":"R Grimshaw","year":"1968","unstructured":"Grimshaw, R.: Propagation of surface waves at high frequencies. IMA J. Appl. Math. 4(2), 174\u2013193 (1968)","journal-title":"IMA J. Appl. Math."},{"issue":"4","key":"3318_CR15","doi-asserted-by":"publisher","first-page":"726","DOI":"10.1023\/A:1022344502818","volume":"40","author":"E Hairer","year":"2000","unstructured":"Hairer, E.: Symmetric projection methods for differential equations on manifolds. BIT Numer. Math. 40(4), 726\u2013734 (2000)","journal-title":"BIT Numer. Math."},{"issue":"5","key":"3318_CR16","doi-asserted-by":"publisher","first-page":"996","DOI":"10.1023\/A:1021989212020","volume":"41","author":"E Hairer","year":"2001","unstructured":"Hairer, E.: Geometric integration of ordinary differential equations on manifolds. BIT Numer. Math. 41(5), 996\u20131007 (2001)","journal-title":"BIT Numer. Math."},{"issue":"2","key":"3318_CR17","doi-asserted-by":"publisher","first-page":"164","DOI":"10.1109\/2945.468403","volume":"1","author":"AJ Hanson","year":"1995","unstructured":"Hanson, A.J., Ma, H.: Quaternion Frame Approach to Streamline Visualization. IEEE Trans. Visual Comput. Graphics 1(2), 164\u2013174 (1995)","journal-title":"IEEE Trans. Visual Comput. Graphics"},{"issue":"4","key":"3318_CR18","doi-asserted-by":"publisher","first-page":"5231","DOI":"10.1093\/mnras\/stz884","volume":"486","author":"DM Hernandez","year":"2019","unstructured":"Hernandez, D.M.: Should N-body integrators be symplectic everywhere in phase space? Mon. Not. R. Astron. Soc. 486(4), 5231\u20135238 (2019)","journal-title":"Mon. Not. R. Astron. Soc."},{"key":"3318_CR19","doi-asserted-by":"publisher","first-page":"613","DOI":"10.1016\/j.cam.2010.01.002","volume":"234","author":"D Jeong","year":"2010","unstructured":"Jeong, D., Kim, J.: A Crank-Nicolson scheme for the Landau-Lifshitz equation. J. Comput. and Appl. Math. 234, 613\u2013623 (2010)","journal-title":"J. Comput. and Appl. Math."},{"issue":"1","key":"3318_CR20","doi-asserted-by":"publisher","first-page":"3","DOI":"10.1111\/sapm.12211","volume":"141","author":"KI Kou","year":"2018","unstructured":"Kou, K.I., Xia, Y.-H.: Linear Quaternion Differential Equations: Basic Theory and Fundamental Results. Stud. Appl. Math. 141(1), 3\u201345 (2018)","journal-title":"Stud. Appl. Math."},{"key":"3318_CR21","doi-asserted-by":"publisher","first-page":"1013","DOI":"10.1016\/S0010-4655(98)00009-5","volume":"111","author":"M Krech","year":"1998","unstructured":"Krech, M., Bunker, A., Landau, D.: Fast spin dynamics algorithms for classical spin systems. Comput. Phys. Commun. 111, 1013 (1998)","journal-title":"Comput. Phys. Commun."},{"key":"3318_CR22","unstructured":"Leung, S.: A Spherical Crank-Nicolson Integrator Based on the Exponential Map and the Spherical Linear Interpolation. arXiv preprint arXiv:2503.17618, (2025)"},{"key":"3318_CR23","doi-asserted-by":"crossref","unstructured":"Leung, S., Chau, W.M., Lee, Y.K.: SLERP-TVDRK (STVDRK) Methods for Ordinary Differential Equations on Spheres. J. Sci. Comput. 101(63), (2024) (arXiv:2410.10420)","DOI":"10.1007\/s10915-024-02702-2"},{"key":"3318_CR24","doi-asserted-by":"publisher","first-page":"1345","DOI":"10.1109\/TIP.2004.834662","volume":"13","author":"M Lysaker","year":"2004","unstructured":"Lysaker, M., Osher, S., Tai, X.-C.: Noise removal using smoothed normals and surface fitting. IEEE Transactions on Image Procecssing 13, 1345\u20131357 (2004)","journal-title":"IEEE Transactions on Image Procecssing"},{"issue":"307","key":"3318_CR25","doi-asserted-by":"publisher","first-page":"2325","DOI":"10.1090\/mcom\/3153","volume":"86","author":"R McLachlan","year":"2017","unstructured":"McLachlan, R., Modin, K., Verdier, O.: A Minimal-variable Sympletic Integrator on Spheres. Math. Comp. 86(307), 2325\u20132344 (2017)","journal-title":"Math. Comp."},{"key":"3318_CR26","doi-asserted-by":"publisher","first-page":"L447","DOI":"10.1088\/0305-4470\/39\/27\/L01","volume":"39","author":"R McLachlan","year":"2006","unstructured":"McLachlan, R., O\u2019Neale, D.: Geometric integration for a two-spin system. J. Phys. A: Math. Gen. 39, L447\u2013L452 (2006)","journal-title":"J. Phys. A: Math. Gen."},{"issue":"3","key":"3318_CR27","doi-asserted-by":"publisher","first-page":"1105","DOI":"10.3934\/dcds.2014.34.1105","volume":"34","author":"PC Moan","year":"2014","unstructured":"Moan, P.C., Niesen, J.: On an asymtotic method for computing the modified energy for sympletic methods. Discrete Contin. Dynam. Systems 34(3), 1105\u20131120 (2014)","journal-title":"Discrete Contin. Dynam. Systems"},{"issue":"5","key":"3318_CR28","doi-asserted-by":"publisher","first-page":"1089","DOI":"10.1023\/A:1021957732493","volume":"41","author":"B Orel","year":"2001","unstructured":"Orel, B.: Extrapolated Magnus Methods. BIT 41(5), 1089\u20131100 (2001)","journal-title":"BIT"},{"key":"3318_CR29","doi-asserted-by":"crossref","unstructured":"Proskova, J.: Description of protein secondary structure using dual quaternions. Journal of Molecular Structure, 1076(89-93), (2014)","DOI":"10.1016\/j.molstruc.2014.07.031"},{"key":"3318_CR30","doi-asserted-by":"publisher","first-page":"306","DOI":"10.1016\/0021-9991(85)90009-9","volume":"60","author":"DC Rapaport","year":"1985","unstructured":"Rapaport, D.C.: Molecular Dynamics Simulation Using Quaternions. J. Comput. Phys. 60, 306\u2013314 (1985)","journal-title":"J. Comput. Phys."},{"issue":"459\u2013470","key":"3318_CR31","first-page":"307","volume":"210","author":"LF Richardson","year":"1911","unstructured":"Richardson, L.F.: The approximate arithmetical solution by finite differences of physical problems including differential equations, with an application to the stresses in a masonry dam. Phil. Trans. Roy. Soc. A 210(459\u2013470), 307\u2013357 (1911)","journal-title":"Phil. Trans. Roy. Soc. A"},{"key":"3318_CR32","doi-asserted-by":"crossref","unstructured":"Schoeller, S.F., Townsend, A.K., Westwood, T.A., Keaveny, E.E.: Methods for suspensions of passive and active filaments. J. Comput. Phys., 424(109846), (2021)","DOI":"10.1016\/j.jcp.2020.109846"},{"key":"3318_CR33","doi-asserted-by":"publisher","first-page":"1495","DOI":"10.1007\/s00214-014-1495-4","volume":"133","author":"A Sergi","year":"2014","unstructured":"Sergi, A.: Computer simulation of quantum dynamics in a classical spin enviornment. Theor. Chem. Acc. 133, 1495 (2014)","journal-title":"Theor. Chem. Acc."},{"issue":"3","key":"3318_CR34","doi-asserted-by":"publisher","first-page":"731","DOI":"10.1093\/imanum\/drn033","volume":"29","author":"T Shingel","year":"2009","unstructured":"Shingel, T.: Interpolation in special orthogonal groups. IMAJ Num. Analy. 29(3), 731\u2013745 (2009)","journal-title":"IMAJ Num. Analy."},{"key":"3318_CR35","doi-asserted-by":"crossref","unstructured":"Shoemake, K.: Animating rotation with quaternion curves. In Proceedings of the 12th annual conference on Computer graphics and interactive techniques, pages 245\u2013254, (1985)","DOI":"10.1145\/325334.325242"},{"key":"3318_CR36","doi-asserted-by":"publisher","first-page":"1073","DOI":"10.1137\/0909073","volume":"9","author":"CW Shu","year":"1988","unstructured":"Shu, C.W.: Total-Variation-Diminishing Time Discreatizations. SIAM. J. Sci. Stat. Compt. 9, 1073\u20131084 (1988)","journal-title":"SIAM. J. Sci. Stat. Compt."},{"key":"3318_CR37","doi-asserted-by":"publisher","first-page":"439","DOI":"10.1016\/0021-9991(88)90177-5","volume":"77","author":"CW Shu","year":"1988","unstructured":"Shu, C.W., Osher, S.J.: Efficient implementation of essentially non-oscillatory shock capturing schemes. J. Comput. Phys. 77, 439\u2013471 (1988)","journal-title":"J. Comput. Phys."},{"key":"3318_CR38","unstructured":"Sol\u00e0, J.: Quaternion kinematics for the error-state Kalman filter. arXiv:1711.02508 [CS.RO], (2017)"},{"key":"3318_CR39","doi-asserted-by":"crossref","unstructured":"Tang, B., Sapiro, G., Caselles: Diffusion of General Data on Non-Flat Manifolds via Harmonic Maps Theory: The Direction Diffusion Case. International Journal of Computer Vision, 36(2):149\u2013161, (2000)","DOI":"10.1023\/A:1008152115986"},{"key":"3318_CR40","doi-asserted-by":"crossref","unstructured":"Tang, B., Sapiro, G., Caselles: Color image enhancement via chromaticity diffusion. IEEE Trans. Image Process., 10:701\u2013707, (2001)","DOI":"10.1109\/83.918563"},{"key":"3318_CR41","doi-asserted-by":"crossref","unstructured":"Tschisgale, S., Frohlich, J.: An immersed boundary method for the fluid-structure interaction of slender flexible structures in viscous fluid. J. Comput. Phys., 423(109801), (2020)","DOI":"10.1016\/j.jcp.2020.109801"},{"key":"3318_CR42","doi-asserted-by":"crossref","unstructured":"Udwadia, F.E., Schutte, A.D.: An Alternative Derivation of the Quaternion Equations of Motion for Rigid-Body Rotational Dynamics. J. Applied Mechanics, 77(044505-1), (2010)","DOI":"10.1115\/1.4000917"},{"issue":"6","key":"3318_CR43","doi-asserted-by":"publisher","first-page":"2085","DOI":"10.1137\/S0036142901396715","volume":"40","author":"L Vese","year":"2002","unstructured":"Vese, L., Osher, J.: Numerical Methods for $$p$$-Harmonic Flows and Applications to Image Processing. SIAM J. Numer. Anal. 40(6), 2085\u20132104 (2002)","journal-title":"SIAM J. Numer. Anal."},{"issue":"3","key":"3318_CR44","doi-asserted-by":"publisher","first-page":"365","DOI":"10.1109\/TVCG.2006.48","volume":"12","author":"R Weinstein","year":"2006","unstructured":"Weinstein, R., Teran, J., Fedkiw, R.: Dynamic Simulation of Articulated Rigid Bodies with Contact and Collision. IEEE Trans. Visual Comput. Graphics 12(3), 365\u2013374 (2006)","journal-title":"IEEE Trans. Visual Comput. Graphics"},{"key":"3318_CR45","doi-asserted-by":"publisher","first-page":"591","DOI":"10.1007\/s211-001-8017-5","volume":"89","author":"J Wennsch","year":"2001","unstructured":"Wennsch, J.: Extrapolation Methods in Lie Groups. Numer. Math. 89, 591\u2013604 (2001)","journal-title":"Numer. Math."},{"key":"3318_CR46","doi-asserted-by":"publisher","first-page":"2163","DOI":"10.1016\/j.jde.2009.06.015","volume":"247","author":"P Wilczynski","year":"2009","unstructured":"Wilczynski, P.: Quaternionic-valued ordinary differential equations. The Riccati equation. J. Differential Equations 247, 2163\u20132187 (2009)","journal-title":"The Riccati equation. J. Differential Equations"},{"key":"3318_CR47","doi-asserted-by":"publisher","DOI":"10.1016\/j.jcp.2019.109104","volume":"404","author":"C Xie","year":"2020","unstructured":"Xie, C., Garcia-Cervera, C.J., Wang, C., Zhou, Z., Chen, J.: Second-order semi-implicit projection methods for micromagnetics simulations. J. Comput. Phys. 404, 109104 (2020)","journal-title":"J. Comput. Phys."}],"container-title":["Journal of Scientific Computing"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s10915-026-03318-4.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1007\/s10915-026-03318-4","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s10915-026-03318-4.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,5,14]],"date-time":"2026-05-14T15:06:19Z","timestamp":1778771179000},"score":1,"resource":{"primary":{"URL":"https:\/\/link.springer.com\/10.1007\/s10915-026-03318-4"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,5,14]]},"references-count":47,"journal-issue":{"issue":"3","published-print":{"date-parts":[[2026,6]]}},"alternative-id":["3318"],"URL":"https:\/\/doi.org\/10.1007\/s10915-026-03318-4","relation":{},"ISSN":["0885-7474","1573-7691"],"issn-type":[{"value":"0885-7474","type":"print"},{"value":"1573-7691","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,5,14]]},"assertion":[{"value":"1 July 2025","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"18 March 2026","order":2,"name":"revised","label":"Revised","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"29 April 2026","order":3,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"14 May 2026","order":4,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"Leung acknowledges the support of the Hong Kong RGC under grants 16302223 and 16300524.","order":1,"name":"Ethics","group":{"name":"EthicsHeading","label":"Declarations"}}],"article-number":"103"}}