Related papers: An explicit splitting SAV scheme for the kinetic L…
The Smectic-A (SmA) phase is modeled by a modified Landau-de Gennes (mLdG) model proposed by Xia et al. [Phys. Rev. Lett., 126 (2021), 177801], in which a tensor order parameter $\mathbf{Q}$ for the orientational order is coupled with a…
The paper studies a scalar auxiliary variable (SAV) method to solve the Cahn-Hilliard equation with degenerate mobility posed on a smooth closed surface {\Gamma}. The SAV formulation is combined with adaptive time stepping and a…
We develop an efficient sampling method by simulating Langevin dynamics with an artificial force rather than a natural force by using the gradient of the potential energy. The standard technique for sampling following the predetermined…
In this paper, a novel high order semi-Lagrangian (SL) spectral volume (SV) method is proposed and studied for nonlinear Vlasov-Poisson (VP) simulations via operator splitting. The proposed algorithm combines both advantages of…
In this paper, a new scheme of arbitrary high order accuracy in both space and time is proposed to solve hyperbolic conservative laws. Based on the idea of flux vector splitting(FVS) scheme, we split all the space and time derivatives in…
This paper proposes a new sampling scheme based on Langevin dynamics that is applicable within pseudo-marginal and particle Markov chain Monte Carlo algorithms. We investigate this algorithm's theoretical properties under standard…
In this paper, we develop a symmetric accelerated stochastic Alternating Direction Method of Multipliers (SAS-ADMM) for solving separable convex optimization problems with linear constraints. The objective function is the sum of a possibly…
We present an algorithm for the simulation of the exact real-time dynamics of classical many-body systems with discrete energy levels. In the same spirit of kinetic Monte Carlo methods, a stochastic solution of the master equation is found,…
We show how the Hamiltonian Monte Carlo algorithm can sometimes be speeded up by "splitting" the Hamiltonian in a way that allows much of the movement around the state space to be done at low computational cost. One context where this is…
This paper applies several well-known tricks from the numerical treatment of deterministic differential equations to improve the efficiency of the Multilevel Monte Carlo (MLMC) method for stochastic differential equations (SDEs) and…
In this paper, we propose a regularized auxiliary variable (RAV) approach and construct accurate and robust time-discrete schemes for a large class of gradient flows. By introducing an auxiliary variable $r=0$ and constructing an auxiliary…
We propose a stochastic method to generate exactly the overdamped Langevin dynamics of semi-flexible Gaussian chains, conditioned to evolve between given initial and final conformations in a preassigned time. The initial and final…
Classical solvable stochastic volatility models (SVM) use a CEV process for instantaneous variance where the CEV parameter $\gamma$ takes just few values: 0 - the Ornstein-Uhlenbeck process, 1/2 - the Heston (or square root) process, 1-…
This paper studies a class of linear unconditionally energy stable schemes for the gradient flows. Such schemes are built on the SAV technique and the general linear time discretization (GLTD) as well as the linearization based on the…
We introduce Semi-Implicit Lagrangian Voronoi Approximation (SILVA), a novel numerical method for the solution of the incompressible Euler and Navier-Stokes equations, which combines the efficiency of semi-implicit time marching schemes…
We present an optimization algorithm that can identify a global minimum of a potentially nonconvex smooth function with high probability, assuming the Gibbs measure of the potential satisfies a logarithmic Sobolev inequality. Our…
We provide a refined explicit estimate of exponential decay rate of underdamped Langevin dynamics in $L^2$ distance, based on a framework developed in [1]. To achieve this, we first prove a Poincar\'{e}-type inequality with Gibbs measure in…
We introduce a novel framework for efficient sampling from complex, unnormalised target distributions by exploiting multiscale dynamics. Traditional score-based sampling methods either rely on learned approximations of the score function or…
In this paper we introduce a new sampling algorithm which has the potential to be adopted as a universal replacement to the Metropolis--Hastings algorithm. It is related to the slice sampler, and motivated by an algorithm which is…
The main objective of this paper is to present an efficient structure-preserving scheme, which is based on the idea of the scalar auxiliary variable approach, for solving the space fractional nonlinear Schr\"{o}dinger equation. First, we…