Search

The long-term dynamics of particles involved in an incompressible flow with a small viscosity ($\epsilon>0$) and slow chemical reactions, is depicted by a class of stochastic reaction-diffusion-advection (RDA) equations with a fast advection term of magnitude $1/\epsilon$. It has been shown in [7] the fast advection asymptotics of stochastic RDA equation in $\mathbb{R}^2$ can be characterized through a stochastic partial differential equation (SPDE) on the graph associated with certain Hamiltonian. To simulate such fast advection asymptotics, we introduce and study an asymptotic-preserving (AP) exponential Euler approximation for the multiscale stochastic RDA equation. There are three key ingredients in proving asymptotic-preserving property of the proposed approximation. First, a strong error estimate, which depends on $1/\epsilon$ linearly, is obtained via a variational argument. Second, we prove the consistency of exponential Euler approximations on the fast advection asymptotics between the original problem and the SPDE on graph. Last, a graph weighted space is introduced to quantify the approximation error for SPDE on graph, which avoids the possible singularity near the vertices. Numerical experiments are carried out to support the theoretical results.

We propose a supervised learning scheme for the first order Hamilton-Jacobi PDEs in high dimensions. The scheme is designed by using the geometric structure of Wasserstein Hamiltonian flows via a density coupling strategy. It is equivalently posed as a regression problem using the Bregman divergence, which provides the loss function in learning while the data is generated through the particle formulation of Wasserstein Hamiltonian flow. We prove a posterior estimate on $L^1$ residual of the proposed scheme based on the coupling density. Furthermore, the proposed scheme can be used to describe the behaviors of Hamilton-Jacobi PDEs beyond the singularity formations on the support of coupling density.Several numerical examples with different Hamiltonians are provided to support our findings., Comment: 29 pages

To capture and simulate geometric surface evolutions, one effective approach is based on the phase field methods. Among them, it is important to design and analyze numerical approximations whose error bound depends on the inverse of the diffuse interface thickness (denoted by $\frac 1\epsilon$) polynomially. However, it has been a long-standing problem whether such numerical error bound exists for stochastic phase field equations. In this paper, we utilize the regularization effect of noise to show that near sharp interface limit, there always exists the weak error bound of numerical approximations, which depends on $\frac 1\epsilon$ at most polynomially. To illustrate our strategy, we propose a polynomial taming fully discrete scheme and present novel numerical error bounds under various metrics. Our method of proof could be also extended to a number of other fully numerical approximations for semilinear stochastic partial differential equations (SPDEs)., Comment: We need a lot of time to modify it and then update in the future

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.

Understanding the properties of the stochastic phase field models is crucial to model processes in several practical applications, such as soft matters and phase separation in random environments. To describe such random evolution, this work proposes and studies two mathematical models and their numerical approximations for parabolic stochastic partial differential equation (SPDE) with a logarithmic Flory--Huggins energy potential. These multiscale models are built based on a regularized energy technique and thus avoid possible singularities of coefficients. According to the large deviation principle, we show that the limit of the proposed models with small noise naturally recovers the classical dynamics in deterministic case. Moreover, when the driving noise is multiplicative, the Stampacchia maximum principle holds which indicates the robustness of the proposed model. One of the main advantages of the proposed models is that they can admit the energy evolution law and asymptotically preserve the Stampacchia maximum bound of the original problem. To numerically capture these asymptotic behaviors, we investigate the semi-implicit discretizations for regularized logrithmic SPDEs. Several numerical results are presented to verify our theoretical findings., Comment: 26 pages, 5 figures

We study the optimal control formulation for stochastic nonlinear Schrodinger equation (SNLSE) on a finite graph. By viewing the SNLSE as a stochastic Wasserstein Hamiltonian flow on density manifold, we show the global existence of a unique strong solution for SNLSE with a linear drift control or a linear diffusion control on graph. Furthermore, we provide the gradient formula, the existence of the optimal control and a description on the optimal condition via the forward and backward stochastic differential equations.

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.

We study a class of fully-discrete schemes for the numerical approximation of solutions of stochastic Cahn--Hilliard equations with cubic nonlinearity and driven by additive noise. The spatial (resp. temporal) discretization is performed with a spectral Galerkin method (resp. a tamed exponential Euler method). We consider two situations: space-time white noise in dimension $d=1$ and trace-class noise in dimensions $d=1,2,3$. In both situations, we prove weak error estimates, where the weak order of convergence is twice the strong order of convergence with respect to the spatial and temporal discretization parameters. To prove these results, we show appropriate regularity estimates for solutions of the Kolmogorov equation associated with the stochastic Cahn--Hilliard equation, which have not been established previously and may be of interest in other contexts.

We study the Wasserstein Hamiltonian flow with a common noise on the density manifold of a finite graph. Under the framework of stochastic variational principle, we first develop the formulation of stochastic Wasserstein Hamiltonian flow and show the local existence of a unique solution. We also establish a sufficient condition for the global existence of the solution. Consequently, we obtain the global well-posedness for the nonlinear Schr\"odinger equations with common noise on graph. In addition, using Wong-Zakai approximation of common noise, we prove the existence of the minimizer for an optimal control problem with common noise. We show that its minimizer satisfies the stochastic Wasserstein Hamiltonian flow on graph as well.

This paper focuses on the construction and analysis of explicit numerical methods of high dimensional stochastic nonlinear Schrodinger equations (SNLSEs). We first prove that the classical explicit numerical methods are unstable and suffer from the numerical divergence phenomenon. Then we propose a kind of explicit splitting numerical methods and prove that the structure-preserving splitting strategy is able to enhance the numerical stability. Furthermore, we establish the regularity analysis and strong convergence analysis of the proposed schemes for SNLSEs based on two key ingredients. One ingredient is proving new regularity estimates of SNLSEs by constructing a logarithmic auxiliary functional and exploiting the Bourgain space. Another one is providing a dedicated error decomposition formula and a novel truncated stochastic Gronwall's lemma, which relies on the tail estimates of underlying stochastic processes. In particular, our result answers the strong convergence problem of numerical methods for 2D SNLSEs emerged from [C. Chen, J. Hong and A. Prohl, Stoch. Partial Differ. Equ. Anal. Comput. 4 (2016), no. 2, 274-318] and [J. Cui and J. Hong, SIAM J. Numer. Anal. 56 (2018), no. 4, 2045-2069]., Comment: 30 pages

In this paper, we study the stochastic Hamiltonian flow in Wasserstein manifold, the probability density space equipped with $L^2$-Wasserstein metric tensor, via the Wong--Zakai approximation. We begin our investigation by showing that the stochastic Euler-Lagrange equation, regardless it is deduced from either variational principle or particle dynamics, can be interpreted as the stochastic kinetic Hamiltonian flows in Wasserstein manifold. We further propose a novel variational formulation to derive more general stochastic Wassersetin Hamiltonian flows, and demonstrate that this new formulation is applicable to various systems including the stochastic Schr\"odinger equation, Schr\"odinger equation with random dispersion, and Schr\"odinger bridge problem with common noise., Comment: 34 pages

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 propose a numerical method to solve the classic $L^2$-optimal transport problem. Our algorithm is based on use of multiple shooting, in combination with a continuation procedure, to solve the boundary value problem associated to the transport problem. We exploit the viewpoint of Wasserstein Hamiltonian flow with initial and target densities, and our method is designed to retain the underlying Hamiltonian structure. Several numerical examples are presented to illustrate the performance of the method., Comment: 25 pages

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

We present a definition of stochastic Hamiltonian process on finite graph via its corresponding density dynamics in Wasserstein manifold. We demonstrate the existence of stochastic Hamiltonian process in many classical discrete problems, such as the optimal transport problem, Schr\"odinger equation and Schr\"odinger bridge problem (SBP). The stationary and periodic properties of Hamiltonian processes are also investigated in the framework of SBP., Comment: 29 pages

This paper investigates the convergence of density approximations for stochastic heat equation in both uniform convergence topology and total variation distance. The convergence order of the densities in uniform convergence topology is shown to be exactly $1/2$ in the nonlinear case and nearly $1$ in the linear case. This result implies that the distributions of the approximations always converge to the distribution of the origin equation in total variation distance. As far as we know, this is the first result on the convergence of density approximations to the stochastic partial differential equation.

We study discretizations of Hamiltonian systems on the probability density manifold equipped with the $L^2$-Wasserstein metric. Based on discrete optimal transport theory, several Hamiltonian systems on graph (lattice) with different weights are derived, which can be viewed as spatial discretizations to the original Hamiltonian systems. We prove the consistency and provide the approximate orders for those discretizations. By regularizing the system using Fisher information, we deduce an explicit lower bound for the density function, which guarantees that symplectic schemes can be used to discretize in time. Moreover, we show desirable long time behavior of these schemes, and demonstrate their performance on several numerical examples., Comment: 34 pages

In this paper, we present an energy-preserving exponentially integrable numerical method for stochastic wave equation with cubic nonlinearity and additive noise. We first apply the spectral Galerkin method to discretizing the original equation and show that this spatial discretization possesses an energy evolution law and certain exponential integrability property. Then the exponential integrability property of the exact solution is deduced by proving the strong convergence of the semi-discretization. To propose a full discrete numerical method which could inherit both the energy evolution law and the exponential integrability, we use the splitting technique and averaged vector field method in the temporal direction. Combining these structure-preserving properties with regularity estimates of the exact and the numerical solutions, we obtain the strong convergence rate of the numerical method. Numerical experiments coincide with these theoretical results.

In this article, we develop and analyze a full discretization, based on the spatial spectral Galerkin method and the temporal drift implicit Euler scheme, for the stochastic Cahn--Hilliard equation driven by multiplicative space-time white noise. By introducing an appropriate decomposition of the numerical approximation, we first use the factorization method to deduce the a priori estimate and regularity estimate of the proposed full discretization. With the help of the variation approach, we then obtain the sharp spatial and temporal convergence rate in negative Sobolev space in mean square sense. Furthermore, the sharp mean square convergence rates in both time and space are derived via the Sobolev interpolation inequality and semigroup theory. To the best of our knowledge, this is the first result on the convergence rate of temporally and fully discrete numerical methods for the stochastic Cahn--Hilliard equation driven by multiplicative space-time white noise.

In this article, we consider the stochastic Cahn--Hilliard equation driven by multiplicative space-time white noise with diffusion coefficient of sublinear growth. By introducing the spectral Galerkin method, we first obtain the well-posedness of the approximated equation in finite dimension. Then with the help of the semigroup theory and the factorization method, the approximation processes is shown to possess many desirable properties. Further, we show that the approximation process is strongly convergent in certain Banach space via the interpolation inequality and variational approach. Finally, the global existence and regularity estimate of the unique solution process are proven by means of the strong convergence of the approximation process., Comment: 23 pages

In this article, we study the density function of the numerical solution of the splitting averaged vector field (AVF) scheme for the stochastic Langevin equation. To deal with the non-globally monotone coefficient in the considered equation, we first present the exponential integrability properties of the exact and numerical solutions. Then we show the existence and smoothness of the density function of the numerical solution by proving its uniform non-degeneracy in Malliavin sense. In order to analyze the approximate error between the density function of the exact solution and that of the numerical solution, we derive the optimal strong convergence rate in every Malliavin--Sobolev norm of the numerical scheme via Malliavin calculus. Combining the approximation result of Donsker's delta function and the smoothness of the density functions, we prove that the convergence rate in density coincides with the optimal strong convergence rate of the numerical scheme., Comment: 39 pages

This article presents explicit exponential integrators for stochastic Maxwell's equations driven by both multiplicative and additive noises. By utilizing the regularity estimate of the mild solution, we first prove that the strong order of the numerical approximation is $\frac 12$ for general multiplicative noise. Combing a proper decomposition with the stochastic Fubini's theorem, the strong order of the proposed scheme is shown to be $1$ for additive noise. Moreover, for linear stochastic Maxwell's equation with additive noise, the proposed time integrator is shown to preserve exactly the symplectic structure, the evolution of the energy as well as the evolution of the divergence in the sense of expectation. Several numerical experiments are presented in order to verify our theoretical findings., Comment: 21 Pages

In this article, we consider the stochastic Cahn--Hilliard equation driven by space-time white noise. We discretize this equation by using a spatial spectral Galerkin method and a temporal accelerated implicit Euler method. The optimal regularity properties and uniform moment bounds of the exact and numerical solutions are shown. Then we prove that the proposed numerical method is strongly convergent with the sharp convergence rate in a negative Sobolev space. By using an interpolation approach, we deduce the spatial optimal convergence rate and the temporal super-convergence rate of the proposed numerical method in strong convergence sense. To the best of our knowledge, this is the first result on the strong convergence rates of numerical methods for the stochastic Cahn--Hilliard equation driven by space-time white noise. This interpolation approach is also applied to the general noise and high dimension cases, and strong convergence rate results of the proposed scheme are given., Comment: 31 pages

Approximating the invariant measure and the expectation of the functionals for parabolic stochastic partial differential equations (SPDEs) with non-globally Lipschitz coefficients is an active research area and is far from being well understood. In this article, we study such problem in terms of a full discretization based on the spectral Galerkin method and the temporal implicit Euler scheme. By deriving the a priori estimates and regularity estimates of the numerical solution via a variational approach and Malliavin calculus, we establish the sharp weak convergence rate of the full discretization. When the SPDE admits a unique $V$-uniformly ergodic invariant measure, we prove that the invariant measure can be approximated by the full discretization. The key ingredients lie on the time-independent weak convergence analysis and time-independent regularity estimates of the corresponding Kolmogorov equation. Finally, numerical experiments confirm the theoretical findings., Comment: 47 pages

Strong and weak approximation errors of a spatial finite element method are analyzed for stochastic partial differential equations(SPDEs) with one-sided Lipschitz coefficients, including the stochastic Allen--Cahn equation, driven by additive noise. In order to give the strong convergence rate of the finite element method, we present an appropriate decomposition and some a priori estimates of the discrete stochastic convolution. To the best of our knowledge, there has been no essentially sharp weak convergence rate of spatial approximation for parabolic SPDEs with non-globally Lipschitz coefficients. To investigate the weak error, we first regularize the original equation by the splitting technique and derive the regularity estimates of the corresponding regularized Kolmogorov equation. Meanwhile, we present the regularity estimate in Malliavin sense and the refined estimate of the finite element method. Combining the regularity estimates of regularized Kolmogorov equation with Malliavin integration by parts formula, the weak convergence rate is shown to be twice the strong convergence rate., Comment: 26pages

This article analyzes an explicit temporal splitting numerical scheme for the stochastic Allen-Cahn equation driven by additive noise, in a bounded spatial domain with smooth boundary in dimension $d\le 3$. The splitting strategy is combined with an exponential Euler scheme of an auxiliary problem. When $d=1$ and the driving noise is a space-time white noise, we first show some a priori estimates of this splitting scheme. Using the monotonicity of the drift nonlinearity, we then prove that under very mild assumptions on the initial data, this scheme achieves the optimal strong convergence rate $\OO(\delta t^{\frac 14})$. When $d\le 3$ and the driving noise possesses some regularity in space, we study exponential integrability properties of the exact and numerical solutions. Finally, in dimension $d=1$, these properties are used to prove that the splitting scheme has a strong convergence rate $\OO(\delta t)$., Comment: 33 pages

In this paper, we consider the well-posedness of the weakly damped stochastic nonlinear Schr\"odinger(NLS) equation driven by multiplicative noise. First, we show the global existence of the unique solution for the damped stochastic NLS equation in critical case. Meanwhile, the exponential integrability of the solution is proved, which implies the continuous dependence on the initial data. Then, we analyze the effect of the damped term and noise on the blow-up phenomenon. By modifying the associated energy, momentum and variance identity, we deduce a sharp blow-up condition for damped stochastic NLS equation in supercritical case. Moreover, we show that when the damped effect is large enough, the damped effect can prevent the blow-up of the solution with high probability., Comment: 19 pages

In this paper, we investigate the damped stochastic nonlinear Schr\"odinger(NLS) equation with multiplicative noise and its splitting-based approximation. When the damped effect is large enough, we prove that the solutions of the damped stochastic NLS equation and the splitting scheme are exponential stable and possess some exponential integrability. These properties lead that the strong order of the scheme is $\frac 12$ and independent of time. Meanwhile, we analyze the regularity of the Kolmogorov equation with respect to the equation. As a consequence, the weak order of the scheme is shown to be twice the strong order and independent of time., Comment: 24 pages

We prove the optimal strong convergence rate of a fully discrete scheme, based on a splitting approach, for a stochastic nonlinear Schr\"odinger (NLS) equation. The main novelty of our method lies on the uniform a priori estimate and exponential integrability of a sequence of splitting processes which are used to approximate the solution of the stochastic NLS equation. We show that the splitting processes converge to the solution with strong order $1/2$. Then we use the Crank--Nicolson scheme to temporally discretize the splitting process and get the temporal splitting scheme which also possesses strong order $1/2$. To obtain a full discretization, we apply this splitting Crank--Nicolson scheme to the spatially discrete equation which is achieved through the spectral Galerkin approximation. Furthermore, we establish the convergence of this fully discrete scheme with optimal strong convergence rate $\mathcal{O}(N^{-2}+\tau^\frac12)$, where $N$ denotes the dimension of the approximate space and $\tau$ denotes the time step size. To the best of our knowledge, this is the first result about strong convergence rates of temporally numerical approximations and fully discrete schemes for stochastic NLS equations, or even for stochastic partial differential equations (SPDEs) with non-monotone coefficients. Numerical experiments verify our theoretical result., Comment: 39 pages

In this paper, we derive a strong convergence rate of spatial finite difference approximations for both focusing and defocusing stochastic cubic Schr\"odinger equations driven by a multiplicative $Q$-Wiener process. Beyond the uniform boundedness of moments for high order derivatives of the exact solution, the key requirement of our approach is the exponential integrability of both the exact and numerical solutions. By constructing and analyzing a Lyapunov functional and its discrete correspondence, we derive the uniform boundedness of moments for high order derivatives of the exact solution and the first order derivative of the numerical solution, which immediately yields the well-posedness of both the continuous and discrete problems. The latter exponential integrability is obtained through a variant of a criterion given by [Cox, Hutzenthaler and Jentzen, arXiv:1309.5595]. As a by-product of this exponential integrability, we prove that the exact and numerical solutions depend continuously on the initial data and obtain a large deviation-type result on the dependence of the noise with first order strong convergence rate., Comment: 26 pages

We indicate that the nonlinear Schr\"odinger equation with white noise dispersion possesses stochastic symplectic and multi-symplectic structures. Based on these structures, we propose the stochastic symplectic and multi-symplectic methods, which preserve the continuous and discrete charge conservation laws, respectively. Moreover, we show that the proposed methods are convergent with temporal order one in probability. Numerical experiments are presented to verify our theoretical results.