Long-time behavior of exact and numerical solutions of stochastic evolution equations on the sphere
Pith reviewed 2026-05-10 18:59 UTC · model grok-4.3
The pith
The stochastic exponential integrator preserves long-time trace formulas for energy, mass, and momentum in linear stochastic PDEs on the sphere, while Euler-Maruyama schemes do not.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For linear stochastic evolution equations on the sphere, the exact solutions satisfy trace formulas for physically relevant quantities such as energy, mass, and momentum. The forward and backward Euler-Maruyama methods fail to reproduce these formulas in the long-time limit, whereas the stochastic exponential integrator preserves them exactly for the three models considered.
What carries the argument
The stochastic exponential integrator, which applies the exact semigroup operator to the linear deterministic part and adds the stochastic integral directly, is the mechanism that maintains the trace formulas.
Load-bearing premise
The equations must be linear so that explicit trace formulas for energy, mass, and momentum can be derived from the forcing terms.
What would settle it
A long-time numerical simulation of one of the three equations using the stochastic exponential integrator that shows growing deviation from the expected trace formula value would falsify the preservation result.
Figures
read the original abstract
We investigate the long-time behavior of exact solutions and numerical approximations of linear stochastic evolution equations defined on the sphere. We focus on three classical models arising in mathematical physics: the stochastic wave equation, the stochastic Schr\"odinger equation, and the stochastic Maxwell's equations. For these SPDEs, we analyze several widely used time integrators with respect to trace formulas describing the evolution of physically relevant quantities such as energy, mass, and momentum dependent on the forcing term. In particular, we prove that the forward and backward Euler-Maruyama schemes fail to reproduce the correct long-time behavior of the exact solutions. In addition, we prove that the stochastic exponential integrator preserves the correct long-time behavior of the physical quantities of interest. Finally, several numerical experiments are provided to illustrate our theoretical findings.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates the long-time behavior of exact solutions and numerical approximations for three linear stochastic evolution equations on the sphere: the stochastic wave equation, the stochastic Schrödinger equation, and stochastic Maxwell's equations. It derives trace formulas for physically relevant quantities such as energy, mass, and momentum that depend on the forcing term, proves that the forward and backward Euler-Maruyama schemes fail to reproduce the correct long-time asymptotics of the exact solutions, demonstrates that the stochastic exponential integrator preserves these trace formulas exactly, and illustrates the findings with numerical experiments.
Significance. If the central claims hold, the work provides a rigorous justification for preferring stochastic exponential integrators in long-time simulations of linear SPDEs on the sphere, where preservation of invariants is critical. The explicit derivations via mild solutions and Itô calculus, together with the numerical validation, constitute a clear strength; the restriction to linear problems with additive noise that admits closed-form traces further strengthens the results by enabling precise comparisons without additional modeling assumptions.
minor comments (3)
- Abstract: the statement that 'several numerical experiments are provided to illustrate our theoretical findings' would be more informative if it briefly indicated which of the three models and which trace formulas are tested numerically.
- The notation and precise definitions of the trace formulas (energy, mass, momentum) should be introduced with a dedicated preliminary subsection before their use in the proofs of preservation or non-preservation, to improve readability for readers outside the immediate subfield.
- The description of the spherical discretization (e.g., truncation level of spherical harmonics or finite-element spaces) in the numerical experiments section should include explicit parameter values and convergence checks to facilitate reproducibility.
Simulated Author's Rebuttal
We sincerely thank the referee for their positive and accurate summary of our manuscript on the long-time behavior of exact and numerical solutions for linear stochastic evolution equations on the sphere. We appreciate the recognition of the rigorous derivations via mild solutions and Itô calculus, as well as the numerical validation demonstrating the advantages of the stochastic exponential integrator in preserving trace formulas for energy, mass, and momentum. The recommendation for minor revision is noted, and we will incorporate improvements to enhance the manuscript accordingly.
Circularity Check
No significant circularity; proofs are direct verifications for linear SPDEs
full rationale
The paper derives exact trace formulas for energy, mass, and momentum from the mild solution and Itô calculus applied to the linear stochastic evolution equations. It then verifies that the stochastic exponential integrator reproduces the exact linear evolution operator and integrated stochastic convolution, thereby satisfying identical discrete trace relations. This is a self-contained mathematical argument relying on linearity and additive forcing to obtain closed-form expressions; no data fitting, self-definitional loops, or load-bearing self-citations are used to establish the preservation property. The failure of Euler-Maruyama schemes is shown by explicit counterexamples to the same trace relations. The derivation chain does not reduce any claimed result to its inputs by construction.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The stochastic evolution equations are linear.
- domain assumption Trace formulas exist that describe the evolution of energy, mass, and momentum depending on the forcing term.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
we prove that the stochastic exponential integrator preserves the correct long-time behavior of the physical quantities of interest
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
trace formula for the energy: E[E(t)] = E[E(0)] + t/2 Tr(Q)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
G. D. Akrivis. Finite difference discretization of the cu bic Schr¨ odinger equation.IMA Journal of Numerical Analysis, 13(1):115–124, January 1993
work page 1993
- [2]
-
[3]
V. V. Anh, P. Broadbridge, A. Olenko, and Y. G. Wang. On app roximation for fractional stochastic partial differential equations on the sphere. Stochastic Environmental Research and Risk Assessment , 32(9):2585–2603, Sep 2018
work page 2018
-
[4]
X. Antoine, W. Bao, and C. Besse. Computational methods f or the dynamics of the nonlinear Schr¨ odinger/Gross–Pitaevskii equations.Computer Physics Communications , 184(12):2621–2633, Decem- ber 2013
work page 2013
-
[5]
R. Anton and D. Cohen. Exponential integrators for stoch astic Schr¨ odinger equations driven by Itˆ o noise. J. Comput. Math. , 36(2):276–309, 2018
work page 2018
- [6]
- [7]
-
[8]
Vector spherical h armonics and their application to magneto- statics
R G Barrera, G A Estevez, and J Giraldo. Vector spherical h armonics and their application to magneto- statics. European Journal of Physics , 6(4):287, oct 1985. LONG-TIME BEHA VIOR OF STOCHASTIC EVOLUTION EQUATIONS ON TH E SPHERE 33
work page 1985
-
[9]
C. Besse. A Relaxation Scheme for the Nonlinear Schr¨ odi nger Equation. SIAM Journal on Numerical Analysis, 42(3):934–952, January 2004
work page 2004
- [10]
- [11]
-
[12]
Z. Brze´ zniak and S. Cerrai. Stochastic wave equations with constraints: well-posedness and Smoluchowski- Kramers diffusion approximation. Comm. Math. Phys. , 406(9):Paper No. 223, 59, 2025
work page 2025
-
[13]
Z. Brze´ zniak and J. Hussain. Global solution of nonlin ear stochastic heat equation with solutions in a Hilbert manifold. Stoch. Dyn. , 20(6):2040012, 29, 2020
work page 2020
-
[14]
Z. Brze´ zniak and M. Ondrej´ at. Strong solutions to sto chastic wave equations with values in Riemannian manifolds. J. Funct. Anal. , 253(2):449–481, 2007
work page 2007
-
[15]
K. Burrage and P. M. Burrage. Low rank Runge-Kutta metho ds, symplecticity and stochastic Hamiltonian problems with additive noise. J. Comput. Appl. Math. , 236(16):3920–3930, 2012
work page 2012
-
[16]
P. M. Burrage and K. Burrage. Structure-preserving Run ge-Kutta methods for stochastic Hamiltonian equations with additive noise. Numer. Algorithms , 65(3):519–532, 2014
work page 2014
-
[17]
W. E. Byerly. An Elementary Treatise on Fourier’s Series and Spherical, C ylindrical, and Ellipsoidal Harmonics: With Applications to Problems in Mathematical P hysics. Ginn, 1893
-
[18]
F. Cabral and F. S. N. Lobo. Electrodynamics and Spaceti me Geometry: Foundations. Foundations of Physics, 47(2):208–228, 2017
work page 2017
-
[19]
R. Carles. Semi-Classical Analysis For Nonlinear Schr¨ odinger Equations. World Scientific Publishing Com- pany, Singapore, 2008
work page 2008
-
[20]
C. Chen, D. Cohen, R. D’Ambrosio, and A. Lang. Drift-pre serving numerical integrators for stochastic Hamiltonian systems. Adv. Comput. Math. , 46(2):Paper No. 27, 22, 2020
work page 2020
-
[21]
C. Chen, J. Hong, and L. Zhang. Preservation of physical properties of stochastic Maxwell equations with additive noise via stochastic multi-symplectic methods. J. Comput. Phys. , 306:500–519, 2016
work page 2016
-
[22]
S. Chun. Method of moving frames to solve time-dependen t Maxwell’s equations on anisotropic curved surfaces: Applications to invisible cloak and ELF propagat ion. J. Comput. Phys. , 340:85–104, 2017
work page 2017
-
[23]
S. Chun, J. Marcon, J. Peir´ o, and S. J. Sherwin. Reducin g errors caused by geometrical inaccuracy to solve partial differential equations with moving frames on c urvilinear domain. Comput. Methods Appl. Mech. Engrg., 398:Paper No. 115261, 24, 2022
work page 2022
-
[24]
D. Cohen. On the numerical discretisation of stochasti c oscillators. Math. Comput. Simulation , 82(8):1478– 1495, 2012
work page 2012
- [25]
-
[26]
D. Cohen and L. Gauckler. One-stage exponential integr ators for nonlinear Schr¨ odinger equations over long times. BIT Numerical Mathematics , 52(4):877–903, December 2012
work page 2012
- [27]
-
[28]
D. Cohen and A. Lang. Numerical approximation and simul ation of the stochastic wave equation on the sphere. Calcolo, 59(3):Paper No. 32, 32, 2022
work page 2022
- [29]
-
[30]
D. Cohen and M. Sigg. Convergence analysis of trigonome tric methods for stiff second-order stochastic differential equations. Numer. Math. , 121(1):1–29, 2012
work page 2012
-
[31]
D. Cohen and G. Vilmart. Drift-preserving numerical in tegrators for stochastic Poisson systems. Int. J. Comput. Math. , 99(1):4–20, 2022
work page 2022
-
[32]
R. D’Ambrosio, M. Moccaldi, and B. Paternoster. Numeri cal preservation of long-term dynamics by stochastic two-step methods. Discrete Contin. Dyn. Syst. Ser. B , 23(7):2763–2773, 2018
work page 2018
-
[33]
H. de la Cruz, J. C. Jimenez, and J. P. Zubelli. Locally li nearized methods for the simulation of stochastic oscillators driven by random forces. BIT, 57(1):123–151, 2017
work page 2017
-
[34]
K. D. Elworthy. Stochastic Differential Equations on Manifolds . London Mathematical Society Lecture Note Series. Cambridge University Press, 1982
work page 1982
-
[35]
L. C. Evans. Partial differential equations , volume 19 of Graduate Studies in Mathematics . American Mathematical Society, Providence, RI, second edition, 201 0. 34 DA VID COHEN, BJ ¨ORN M ¨ULLER, AND ANDREA PAPINI
-
[36]
Y. Fang and A. Waters. Dispersive estimates for Maxwell ’s equations in the exterior of a sphere. Journal of Differential Equations , 415:855–885, 2025
work page 2025
- [37]
- [38]
- [39]
-
[40]
J. Hong, R. Scherer, and L. Wang. Predictor-corrector m ethods for a linear stochastic oscillator with additive noise. Math. Comput. Modelling , 46(5-6):738–764, 2007
work page 2007
-
[41]
E. Jansson, A. Lang, and M. Pereira. Non-stationary Gau ssian random fields on hypersurfaces: Sampling and strong error analysis, 2024
work page 2024
- [42]
-
[43]
Y. Kazashi and Q. T. Le Gia. A non-uniform discretizatio n of stochastic heat equations with multiplicative noise on the unit sphere. J. Complexity , 50:43–65, 2019
work page 2019
-
[44]
P. E. Kloeden and E. Platen. Numerical solution of stochastic differential equations , volume 23 of Appli- cations of Mathematics (New York) . Springer-Verlag, Berlin, 1992
work page 1992
-
[45]
A. Lang and I. Motschan-Armen. Euler-Maruyama approxi mations of the stochastic heat equation on the sphere. J. Comput. Dyn. , 11(1):23–42, 2024
work page 2024
-
[46]
A. Lang and B. M¨ uller. Isotropic Q-fractional Brownian motion on the sphere: regularity and f ast simu- lation. Philos. Trans. Roy. Soc. A , 383(2298):Paper No. 20240238, 13, 2025
work page 2025
-
[47]
A. Lang, A. Papini, and V. Schwarz. Approximation of the L´ evy-driven stochastic heat equation on the sphere, 2025
work page 2025
-
[48]
A. Lang and C. Schwab. Isotropic Gaussian random fields o n the sphere: regularity, fast simulation and stochastic partial differential equations. Ann. Appl. Probab. , 25(6):3047–3094, 2015
work page 2015
-
[49]
Q. T. Le Gia and J. Peach. A spectral method to the stochas tic Stokes equations on the sphere. In B. Lamichhane, T. Tran, and J. Bunder, editors, Proceedings of the 18th Biennial Computational Tech- niques and Applications Conference , CTAC-2018 , volume 60 of ANZIAM J. , pages C52–C64, June 2019
work page 2018
-
[50]
J. M. Lee. Introduction to Riemannian Manifolds , volume 176 of Graduate Texts in Mathematics . Springer International Publishing, Cham, 2018
work page 2018
-
[51]
G. J. Lord, C. E. Powell, and T. Shardlow. An introduction to computational stochastic PDEs . Cambridge Texts in Applied Mathematics. Cambridge University Press, New York, 2014
work page 2014
-
[52]
C. Lubich. From Quantum to Classical Molecular Dynamics: Reduced Mode ls and Numerical Analysis , volume 12 of Zurich Lectures in Advanced Mathematics . EMS Press, 1 edition, September 2008
work page 2008
-
[53]
M. Malzoumati-Khiaban, A. Foroush Bastani, and M. R. Ya ghouti. Long-term adaptive symplectic nu- merical integration of linear stochastic oscillators driv en by additive white noise. Numer. Algorithms , 80(3):1059–1095, 2019
work page 2019
-
[54]
D. Marinucci and G. Peccati. Random Fields on the Sphere: Representation, Limit Theorem s and Cos- mological Applications. London Mathematical Society Lecture Note Series. Cambrid ge University Press, Cambridge, 2011
work page 2011
-
[55]
A. Milani and R. Picard. Decomposition theorems and their application to non-linea r electro- and magneto- static boundary value problems , pages 317–340. Springer Berlin Heidelberg, Berlin, Heide lberg, 1988
work page 1988
-
[56]
G. N. Milstein and M. V. Tretyakov. Stochastic numerics for mathematical physics . Scientific Computation. Springer, Cham, [2021] ©2021. Second edition [of 2069903]
work page 2021
-
[57]
Finite Element Methods for Maxwell’s Equations
Peter Monk. Finite Element Methods for Maxwell’s Equations . Oxford University Press, 04 2003
work page 2003
-
[58]
S. Munshi and R. Yang. Complex solutions to Maxwell’s eq uations. Complex Analysis and its Synergies , 8(1), 2022
work page 2022
-
[59]
S. Ohkuro. Differential Forms and Maxwell’s Field: An Ap plication of Harmonic Integrals. Journal of Mathematical Physics , 11(6):2005–2012, 06 1970
work page 2005
-
[60]
S. Peszat and J. Zabczyk. Stochastic Partial Differential Equations with L´ evy noise: An Evolution Equa- tion Approach. Cambridge University Press, 2007
work page 2007
-
[61]
J. M. Sanz-Serna and J. G. Verwer. Conerservative and No nconservative Schemes for the Solution of the Nonlinear Schr¨ odinger Equation.IMA Journal of Numerical Analysis , 6(1):25–42, January 1986
work page 1986
-
[62]
H. Schurz. A numerical method for nonlinear stochastic wave equations in R1. Dyn. Contin. Discrete Impuls. Syst. Ser. A Math. Anal. , 14(Advances in Dynamical Systems, suppl. S2):74–78, 2007 . LONG-TIME BEHA VIOR OF STOCHASTIC EVOLUTION EQUATIONS ON TH E SPHERE 35
work page 2007
-
[63]
H. Schurz. Analysis and discretization of semi-linear stochastic wave equations with cubic nonlinearity and additive space-time noise. Discrete Contin. Dyn. Syst. Ser. S , 1(2):353–363, 2008
work page 2008
-
[64]
H. Schurz. Stochastic wave equations with cubic nonlin earity and Q-regular additive noise in R2. Discrete Contin. Dyn. Syst. , (Dynamical systems, differential equations and applicati ons. 8th AIMS Conference. Suppl. Vol. II):1299–1308, 2011
work page 2011
-
[65]
M. J. Senosiain and A. Tocino. A review on numerical sche mes for solving a linear stochastic oscillator. BIT, 55(2):515–529, 2015
work page 2015
-
[66]
M. J. Senosiain and A. Tocino. On the numerical integrat ion of the undamped harmonic oscillator driven by independent additive gaussian white noises. Appl. Numer. Math. , 137:49–61, 2019
work page 2019
-
[67]
M. H. Stone. On One-Parameter Unitary Groups in Hilbert Space. Annals of Mathematics , 33(3):643–648, 1932
work page 1932
-
[68]
R. S. Strichartz. Analysis of the Laplacian on the compl ete Riemannian manifold. Journal of Functional Analysis, 1983
work page 1983
-
[69]
A. H. Strømmen Melbø and D. J. Higham. Numerical simulat ion of a linear stochastic oscillator with additive noise. Appl. Numer. Math. , 51(1):89–99, 2004
work page 2004
-
[70]
W. Su. On the peaks of a stochastic heat equation on a sphe re with a large radius. Electron. J. Probab. , 25:Paper No. 5, 38, 2020
work page 2020
-
[71]
J. Sun, C. Shu, and Y. Xing. Multi-symplectic discontin uous Galerkin methods for the stochastic Maxwell equations with additive noise. J. Comput. Phys. , 461:Paper No. 111199, 30, 2022
work page 2022
-
[72]
G. Szeg. Orthogonal polynomials, volume 23. American Mathematical Soc., 1939
work page 1939
-
[73]
A. Tocino. On preserving long-time features of a linear stochastic oscillator. BIT, 47(1):189–196, 2007
work page 2007
-
[74]
P. Wang, J. Hong, and D. Xu. Construction of symplectic R unge-Kutta methods for stochastic Hamilton- ian systems. Commun. Comput. Phys. , 21(1):237–270, 2017
work page 2017
-
[75]
Y. Wang, Z. Lang, X. Yin, and Z. Zhao. Exponentially fitte d midpoint scheme for a stochastic oscillator. Mathematics, 14(1), 2026
work page 2026
-
[76]
N. Weck. Maxwell’s boundary value problem on riemannia n manifolds with nonsmooth boundaries. Jour- nal of Mathematical Analysis and Applications , 46(2):410–437, 1974
work page 1974
- [77]
-
[78]
K. Yee. Numerical solution of initial boundary value pr oblems involving Maxwell’s equations in isotropic media. IEEE Transactions on Antennas and Propagation , 14(3):302–307, 1966
work page 1966
-
[79]
X. Yin, C. Zhang, and Y. Liu. Exponentially fitted trapez oidal scheme for a stochastic oscillator. J. Comput. Math. , 35(6):801–813, 2017
work page 2017
-
[80]
E ” ∥∇ S2uptq∥2 L2pS2;C3q ı “ E
Y. Zhou and D. Liang. Modeling and FDTD discretization o f stochastic Maxwell’s equations with Drude dispersion. J. Comput. Phys. , 509:Paper No. 113033, 29, 2024. Department of Mathematical Sciences, Chalmers University of Technology and University of Gothenburg, 41296 Gothenburg, Sweden Email address : david.cohen@chalmers.se Department of Mathematical ...
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.