Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConvergence in Density of Splitting AVF Scheme for Stochastic Langevin Equation
134
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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 Donskers 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.
rate research
Read More
We consider a repeated quantum interaction model describing a small system $Hh_S$ in interaction with each one of the identical copies of the chain $bigotimes_{N^*}C^{n+1}$, modeling a heat bath, one after another during the same short time intervals $[0,h]$. We suppose that the repeated quantum interaction Hamiltonian is split in two parts: a free part and an interaction part with time scale of order $h$. After giving the GNS representation, we establish the relation between the time scale $h$ and the classical low density limit. We introduce a chemical potential $mu$ related to the time $h$ as follows: $h^2=e^{betamu}$. We further prove that the solution of the associated discrete evolution equation converges strongly, when $h$ tends to 0, to the unitary solution of a quantum Langevin equation directed by Poisson processes.
Strong convergence rates for (temporal, spatial, and noise) numerical approximations of semilinear stochastic evolution equations (SEEs) with smooth and regular nonlinearities are well understood in the scientific literature. Weak convergence rates for numerical approximations of such SEEs have been investigated since about 11 years and are far away from being well understood: roughly speaking, no essentially sharp weak convergence rates are known for parabolic SEEs with nonlinear diffusion coefficient functions; see Remark 2.3 in [A. Debussche, Weak approximation of stochastic partial differential equations: the nonlinear case, Math. Comp. 80 (2011), no. 273, 89-117] for details. In this article we solve the weak convergence problem emerged from Debussches article in the case of spectral Galerkin approximations and establish essentially sharp weak convergence rates for spatial spectral Galerkin approximations of semilinear SEEs with nonlinear diffusion coefficient functions. Our solution to the weak convergence problem does not use Malliavin calculus. Rather, key ingredients in our solution to the weak convergence problem emerged from Debussches article are the use of appropriately modifie
We consider a mono-dimensional two-velocities scheme used to approximate the solutions of a scalar hyperbolic conservative partial differential equation. We prove the convergence of the discrete solution toward the unique entropy solution by first estimating the supremum norm and the total variation of the discrete solution, and second by constructing a discrete kinetic entropy-entropy flux pair being given a continuous entropy-entropy flux pair of the hyperbolic system. We finally illustrate our results with numerical simulations of the advection equation and the Burgers equation.
We put forward and analyze an explicit finite difference scheme for the Camassa-Holm shallow water equation that can handle general $H^1$ initial data and thus peakon-antipeakon interactions. Assuming a specified condition restricting the time step in terms of the spatial discretization parameter, we prove that the difference scheme converges strongly in $H^1$ towards a dissipative weak solution of Camassa-Holm equation.
We propose an unbiased Monte-Carlo estimator for $mathbb{E}[g(X_{t_1}, cdots, X_{t_n})]$, where $X$ is a diffusion process defined by a multi-dimensional stochastic differential equation (SDE). The main idea is to start instead from a well-chosen simulatable SDE whose coefficients are updated at independent exponential times. Such a simulatable process can be viewed as a regime-switching SDE, or as a branching diffusion process with one single living particle at all times. In order to compensate for the change of the coefficients of the SDE, our main representation result relies on the automatic differentiation technique induced by Bismu-Elworthy-Li formula from Malliavin calculus, as exploited by Fournie et al.(1999) for the simulation of the Greeks in financial applications. In particular, this algorithm can be considered as a variation of the (infinite variance) estimator obtained in Bally and Kohatsu-Higa [Section 6.1](2014) as an application of the parametrix method.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
