Search

In this manuscript, we propose efficient stochastic semi-explicit symplectic schemes tailored for nonseparable stochastic Hamiltonian systems (SHSs). These semi-explicit symplectic schemes are constructed by introducing augmented Hamiltonians and using symmetric projection. In the case of the artificial restraint in augmented Hamiltonians being zero, the proposed schemes also preserve quadratic invariants, making them suitable for developing semi-explicit charge-preserved multi-symplectic schemes for stochastic cubic Schr\"odinger equations with multiplicative noise. Through numerical experiments that validate theoretical results, we demonstrate that the proposed stochastic semi-explicit symplectic scheme, which features a straightforward Newton iteration solver, outperforms the traditional stochastic midpoint scheme in terms of effectiveness and accuracy.

It is known that when the diffuse interface thickness $\epsilon$ vanishes, the sharp interface limit of the stochastic reaction-diffusion equation is formally a stochastic geometric flow. To capture and simulate such geometric flow, it is crucial to develop numerical approximations whose error bounds depends on $\frac 1\epsilon$ polynomially. However, due to loss of spectral estimate of the linearized stochastic reaction-diffusion equation, how to get such error bound of numerical approximation has been an open problem. In this paper, we solve this weak error bound problem for stochastic reaction-diffusion equations near sharp interface limit. We first introduce a regularized problem which enjoys the exponential ergodicity. Then we present the regularity analysis of the regularized Kolmogorov and Poisson equations which only depends on $\frac 1{\epsilon}$ polynomially. Furthermore, we establish such weak error bound. This phenomenon could be viewed as a kind of the regularization effect of noise on the numerical approximation of stochastic partial differential equation (SPDE). As a by-product, a central limit theorem of the weak approximation is shown near sharp interface limit. Our method of proof could be extended to a number of other spatial and temporal numerical approximations for semilinear SPDEs., Comment: 58 pages

This paper is concerned with the random effect of the noise dispersion for stochastic logarithmic Schr\"odinger equation emerged from the optical fibre with dispersion management. The well-posedness of the logarithmic Schr\"odinger equation with white noise dispersion is established via the regularization energy approximation and a spatial scaling property. For the small noise case, the effect of the noise dispersion is quantified by the proven large deviation principle under additional regularity assumptions on the initial datum. As an application, we show that for the regularized model, the exit from a neighborhood of the attractor of deterministic equation occurs on a sufficiently large time scale. Furthermore, the exit time and exit point in the small noise case, as well as the effect of large noise dispersion, is also discussed for the stochastic logarithmic Schr\"odinger equation.

Stochastic Klein--Gordon--Schr\"odinger (KGS) equations are important mathematical models and describe the interaction between scalar nucleons and neutral scalar mesons in the stochastic environment. In this paper, we propose novel structure-preserving schemes to numerically solve stochastic KGS equations with additive noise, which preserve averaged charge evolution law, averaged energy evolution law, symplecticity, and multi-symplecticity. By applying central difference, sine pseudo-spectral method, or finite element method in space and modifying finite difference in time, we present some charge and energy preserved fully-discrete scheme for the original system. In addition, combining the symplectic Runge-Kutta method in time and finite difference in space, we propose a class of multi-symplectic discretizations preserving the geometric structure of the stochastic KGS equation. Finally, numerical experiments confirm theoretical findings., Comment: This is a duplicate submission and has been withdrawn. Please see arXiv:2305.05117

In this paper, we propose and analyze semi-implicit numerical schemes for the stochastic wave equation (SWE) with general nonlinearity and multiplicative noise. These numerical schemes, called stochastic scalar auxiliary variable (SAV) schemes, are constructed by transforming the considered SWE into a higher dimensional stochastic system with a stochastic SAV. We prove that they can be solved explicitly and preserve the modified energy evolution law and the regularity structure of the original system. These structure-preserving properties are the keys to overcoming the mutual effect of the noise and nonlinearity. By proving new regularity estimates of the introduced SAV, we establish the strong convergence rate of stochastic SAV schemes and the further fully-discrete schemes with the finite element method in spatial direction. To the best of our knowledge, this is the first result on the construction and strong convergence of semi-implicit energy-preserving schemes for nonlinear SWE.

Stochastic Hamiltonian partial differential equations, which possess the multi-symplectic conservation law, are an important and fairly large class of systems. The multi-symplectic methods inheriting the geometric features of stochastic Hamiltonian partial differential equations provide numerical approximations with better numerical stability, and are of vital significance for obtaining correct numerical results. In this paper, we propose three novel multi-symplectic methods for stochastic Hamiltonian partial differential equations based on the local radial basis function collocation method, the splitting technique, and the partitioned Runge-Kutta method. Concrete numerical methods are presented for nonlinear stochastic wave equations, stochastic nonlinear Schr\"odinger equations, stochastic Korteweg-de Vries equations and stochastic Maxwell equations. We take stochastic wave equations as examples to perform numerical experiments, which indicate the validity of the proposed methods.

In this paper, we study two kinds of structure-preserving splitting methods, including the Lie--Trotter type splitting method and the finite difference type method, for the stochasticlogarithmic Schr\"odinger equation (SlogS equation) via a regularized energy approximation. We first introduce a regularized SlogS equation with a small parameter $0<\epsilon\ll1$ which approximates the SlogS equation and avoids the singularity near zero density. Then we present a priori estimates, the regularized entropy and energy, and the stochastic symplectic structure of the proposed numerical methods. Furthermore, we derive both the strong convergence rates and the convergence rates of the regularized entropy and energy. To the best of our knowledge, this is the first result concerning the construction and analysis of numerical methods for stochastic Schr\"odinger equations with logarithmic nonlinearities., Comment: 25

In this paper, we focus on constructing numerical schemes preserving the averaged energy evolution law for nonlinear stochastic wave equations driven by multiplicative noise. We first apply the compact finite difference method and the interior penalty discontinuous Galerkin finite element method to discretize space variable and present two semi-discrete schemes, respectively. Then we make use of the discrete gradient method and the Pad\'e approximation to propose efficient fully-discrete schemes. These semi-discrete and fully-discrete schemes are proved to preserve the discrete averaged energy evolution law. In particular, we also prove that the proposed fully-discrete schemes exactly inherit the averaged energy evolution law almost surely if the considered model is driven by additive noise. Numerical experiments are given to confirm theoretical findings.

In this paper, we prove the global existence and uniqueness of the solution of the stochastic logarithmic Schr\"odinger (SlogS) equation driven by additive noise or multiplicative noise. The key ingredient lies on the regularized stochastic logarithmic Schr\"odinger (RSlogS) equation with regularized energy and the strong convergence analysis of the solutions of (RSlogS) equations. In addition, temporal H\"older regularity estimates and uniform estimates in energy space $\mathbb H^1(\mathcal O)$ and weighted Sobolev space $L^2_{\alpha}(\mathcal O)$ of the solutions for both SlogS equation and RSlogS equation are also obtained., Comment: 32 pages

In this paper, we focus on two kinds of large deviations principles (LDPs) of the invariant measures of Langevin equations and their numerical methods, as the noise intensity $\epsilon\to 0$ and the dissipation intensity $\nu\to\infty$ respectively. First, by proving the weak LDP and the exponential tightness, we conclude that the invariant measure $\{\mu_{\nu,\epsilon}\}$ of the exact solution satisfies the LDPs as $\epsilon\to0$ and $\nu\to\infty$ respectively. Then, we study whether there exist numerical methods asymptotically preserving these two LDPs of $\{\mu_{\nu,\epsilon}\}$ in the sense that the rate functions of invariant measures of numerical methods converge pointwise to the rate function of $\{\mu_{\nu,\epsilon}\}$ as the step-size tends to zero. The answer is positive for the linear Langevin equation. For the small noise case, we show that a large class of numerical methods can asymptotically preserve the LDP of $\{\mu_{\nu,\epsilon}\}_{\epsilon>0}$ as $\epsilon\to0$. For the strong dissipation case, we study the stochastic $\theta$-method ($\theta\in[1/2,1]$) and show that only the midpoint scheme ($\theta=1/2$) can asymptotically preserve the LDP of $\{\mu_{\nu,\epsilon}\}_{\nu>0}$ as $\nu\to\infty$. These results indicate that in the linear case, the LDP as $\epsilon\to0$ and the LDP as $\nu\to\infty$ for the invariant measures of numerical methods have intrinsic differences: the common numerical methods can asymptotically preserve the LDP of $\{\mu_{\nu,\epsilon}\}_{\epsilon>0}$ as $\epsilon\to0$ while the asymptotical preservation of numerical methods for the LDP of $\{\mu_{\nu,\epsilon}\}_{\nu>0}$ as $\nu\to\infty$ depends on the choice of numerical methods. To the best of our knowledge, this is the first result of investigating the relationship between the LDPs of invariant measures of stochastic differential equations and those of their numerical methods.

We propose a class of numerical integration methods for stochastic Poisson systems (SPSs) of arbitrary dimensions. Based on the Darboux-Lie theorem, we transform the SPSs to their canonical form, the generalized stochastic Hamiltonian systems (SHSs), via canonical coordinate transformations found by solving certain PDEs defined by the Poisson brackets of the SPSs. An a-generating function approach with \alpha\in [0,1] is then used to create symplectic discretizations of the SHSs, which are then transformed back by the inverse coordinate transformation to numerical integrators for the SPSs. These integrators are proved to preserve both the Poisson structure and the Casimir functions of the SPSs. Applications to a three-dimensional stochastic rigid body system and a three-dimensional stochastic Lotka-Volterra system show efficiency of the proposed methods.

In this paper, we consider the large deviations principles (LDPs) for the stochastic linear Schr\"odinger equation and its symplectic discretizations. These numerical discretizations are the spatial semi-discretization based on spectral Galerkin method, and the further full discretizations with symplectic schemes in temporal direction. First, by means of the abstract G\"artner--Ellis theorem, we prove that the observable $B_T=\frac{u(T)}{T}$, $T>0$ of the exact solution $u$ is exponentially tight and satisfies an LDP on $L^2(0, \pi; \mathbb C)$. Then, we present the LDPs for both $\{B^M_T\}_{T>0}$ of the spatial discretization $\{u^M\}_{M\in\mathbb N}$ and $\{B^M_N\}_{N\in \mathbb N}$ of the full discretization $\{u^M_N\}_{M,N\in\mathbb N}$, where $B^M_T=\frac{u^M(T)}{T}$ and $B^M_N=\frac{u^M_N}{N\tau}$ are the discrete approximations of $B_T$. Further, we show that both the semi-discretization $\{u^M\}_{M\in \mathbb N}$ and the full discretization $\{u^M_N\}_{M,N\in \mathbb N}$ based on temporal symplectic schemes can weakly asymptotically preserve the LDP of $\{B_T\}_{T>0}$. These results show the ability of symplectic discretizations to preserve the LDP of the stochastic linear \xde equation, and first provide an effective approach to approximating the LDP rate function in infinite dimensional space based on the numerical discretizations.