Related papers: An explicit splitting SAV scheme for the kinetic L…
There has been considerable interest in designing Markov chain Monte Carlo algorithms by exploiting numerical methods for Langevin dynamics, which includes Hamiltonian dynamics as a deterministic case. A prominent approach is Hamiltonian…
We carry out a rigorous error analysis of the first-order semi-discrete (in time) consistent splitting scheme coupled with a generalized scalar auxiliary variable (GSAV) approach for the Navier-Stokes equations with no-slip boundary…
In this paper we propose and analyze a second order accurate (in time) numerical scheme for the square phase field crystal (SPFC) equation, a gradient flow modeling crystal dynamics at the atomic scale in space but on diffusive scales in…
In this paper, two efficient and magnetization norm preserving numerical schemes based on the scalar auxiliary variable (SAV) method are developed for calculating the ground state in micromagnetic structures. The first SAV scheme is based…
For a class of fourth order gradient flow problems, integration of the scalar auxiliary variable (SAV) time discretization with the penalty-free discontinuous Galerkin (DG) spatial discretization leads to SAV-DG schemes. These schemes are…
Consistent splitting schemes are among the most accurate pressure segregation methods, incurring no splitting errors or spurious boundary conditions. Nevertheless, their theoretical properties are not yet fully understood, especially when…
In this paper, we consider numerical approximations for solving the inductionless magnetohydrodynamic (MHD) equations. By utilizing the scalar auxiliary variable (SAV) approach for dealing with the convective and coupling terms, we propose…
We propose a new algorithm---Stochastic Proximal Langevin Algorithm (SPLA)---for sampling from a log concave distribution. Our method is a generalization of the Langevin algorithm to potentials expressed as the sum of one stochastic smooth…
Langevin dynamics has found a large number of applications in sampling, optimization and estimation. Preconditioning the gradient in the dynamics with the covariance - an idea that originated in literature related to solving estimation and…
The Langevin dynamics is a diffusion process extensively used, in particular in molecular dynamics simulations, to sample Gibbs measures. Some alternatives based on (piecewise deterministic) kinetic velocity jump processes have gained…
In this paper, we consider a novel auxiliary variable method to obtain energy stable schemes for gradient flows. The auxiliary variable based on energy bounded above does not limited to the hypothetical conditions adopted in previous…
We study the design and implementation of numerical methods to solve the generalized Langevin equation (GLE) focusing on canonical sampling properties of numerical integrators. For this purpose, we cast the GLE in an extended phase space…
In this paper, we propose a new class of splitting methods to solve the stochastic Langevin equation, which can simultaneously preserve the ergodicity and exponential integrability of the original equation. The central idea is to extract a…
We present an implicit Split-Step explicit Euler type Method (dubbed SSM) for the simulation of McKean-Vlasov Stochastic Differential Equations (MV-SDEs) with drifts of superlinear growth in space, Lipschitz in measure and non-constant…
This article is concerned with sampling from Gibbs distributions $\pi(x)\propto e^{-U(x)}$ using Markov chain Monte Carlo methods. In particular, we investigate Langevin dynamics in the continuous- and the discrete-time setting for such…
We establish a general framework for developing, efficient energy stable numerical schemes for gradient flows and develop three classes of generalized scalar auxiliary variable approaches (G-SAV). Numerical schemes based on the G-SAV…
We carry out the convergence analysis of the Scalar Auxiliary Variable (SAV) method applied to the nonlinear Schr\"odinger equation which preserves a modified Hamiltonian on the discrete level. We derive a weak and strong convergence…
We develop several efficient numerical schemes which preserve exactly the global constraints for constrained gradient flows. Our schemes are based on the SAV approach combined with the Lagrangian multiplier approach. They are as efficient…
Monte Carlo simulations are widely used to simulate complex molecular systems, but standard approaches suffer from metastability. Lately, the use of non-local proposal updates in a collective-variable (CV) space has been proposed in several…
This paper proposes a finite element scheme, based on the Scalar Auxiliary Variable (SAV) approach, for the Cahn-Hilliard equation--a model that possesses significant physical relevance and a rich mathematical structure. A convergence…