Related papers: Stabilising complex Langevin simulations
This paper proposes a numerical method, termed PosDiv-CDG, that provably preserves both positivity and the globally divergence-free (DF) condition at arbitrarily high order in multiple dimensions. It resolves the fundamental structural…
We are interested in understanding stability (almost sure boundedness) of stochastic approximation algorithms (SAs) driven by a `controlled Markov' process. Analyzing this class of algorithms is important, since many reinforcement learning…
We demonstrate that modulation of the local strength of the cubic self-focusing (SF) nonlinearity in the two-dimensional (2D) geometry, in the form of a circle with contrast $\Delta g$ of the SF coefficient relative to the ambient medium…
An off-lattice Monte Carlo algorithm for solutions of equilibrium polymers (EP) is proposed. At low and moderate densities this is shown to reproduce faithfully the (static) properties found recently for flexible linear EP using a lattice…
This paper proposes a local dynamic model for large-eddy simulation (LES) without averaging in homogeneous directions. It is demonstrated that the widely-used dynamic Smagorinsky model (DSM) has a singular dynamic model constant if it is…
As an important Markov Chain Monte Carlo (MCMC) method, stochastic gradient Langevin dynamics (SGLD) algorithm has achieved great success in Bayesian learning and posterior sampling. However, SGLD typically suffers from slow convergence…
Approximate Thompson sampling with Langevin Monte Carlo broadens its reach from Gaussian posterior sampling to encompass more general smooth posteriors. However, it still encounters scalability issues in high-dimensional problems when…
In this paper, we revisit the recently established theoretical guarantees for the convergence of the Langevin Monte Carlo algorithm of sampling from a smooth and (strongly) log-concave density. We improve the existing results when the…
Subsurface datasets inherently possess big data characteristics such as vast volume, diverse features, and high sampling speeds, further compounded by the curse of dimensionality from various physical, engineering, and geological inputs.…
We revisit two basic Direct Simulation Monte Carlo Methods to model aggregation kinetics and extend them for aggregation processes with collisional fragmentation (shattering). We test the performance and accuracy of the extended methods and…
I report recent results from full QCD simulations by CP-PACS and JLQCD collaborations.
We show that the spectrum of the Dirac operator in complex Langevin simulations of QCD at non-zero chemical potential must behave in a way which is radically different from the one in simulations with ordinary non-complexified gauge fields:…
We propose a novel framework for learning stabilizable nonlinear dynamical systems for continuous control tasks in robotics. The key idea is to develop a new control-theoretic regularizer for dynamics fitting rooted in the notion of…
We re-examine the dynamical stability of the nakedly singular, static, spherical ly symmetric solutions of the Einstein-Klein Gordon system. We correct an earlier proof of the instability of these solutions, and demonstrate that there are…
We investigate the stabilization of a multidimensional system of coupled wave equations with only one Kelvin Voigt damping. Using a unique continuation result based on a Carleman estimate and a general criteria of Arendt Batty, we prove the…
We investigate the stability of a Sequential Monte Carlo (SMC) method applied to the problem of sampling from a target distribution on $\mathbb{R}^d$ for large $d$. It is well known that using a single importance sampling step one produces…
Generalized Langevin dynamics (GLD) arise in the modeling of a number of systems, ranging from structured fluids that exhibit a viscoelastic mechanical response, to biological systems, and other media that exhibit anomalous diffusive…
We theoretically analyze a scheme for fast stabilization of arbitrary qubit states with high fidelities, extending a protocol recently demonstrated experimentally [Lu et al., Phys. Rev. Lett. 119, 150502 (2017)]. That experiment utilized…
We discuss the design of state-of-the-art numerical methods for molecular dynamics, focusing on the demands of soft matter simulation, where the purposes include sampling and dynamics calculations both in and out of equilibrium. We discuss…
We discuss the interplay of low-x physics and QCD scaling violations by extending the unified approach describing inclusive structure functions and diffractive production in $\gamma* p$ interactions proposed in previous papers, to large…