Related papers: Efficient stochastic thermostatting of path integr…
We have discussed the dynamics of Langevin model subjected to colored noise, by using the functional-integral method (FIM) combined with equations of motion for mean and variance of the state variable. Two sets of colored noise have been…
In light of the recently developed complete GJ set of single random variable stochastic, discrete-time St{\o}rmer-Verlet algorithms for statistically accurate simulations of Langevin equations, we investigate two outstanding questions: 1)…
Path integral Monte Carlo (PIMC) simulations are used to calculate the momentum distribution of the homogeneous electron gas at finite temperature. This is done by calculating the off-diagonal elements of the real-space density matrix,…
We present a new algorithm for isothermal-isobaric molecular-dynamics simulation. The method uses an extended Hamiltonian with an Andersen piston combined with the Nos'e-Poincar'e thermostat, recently developed by Bond, Leimkuhler and Laird…
Measurement of the velocity field in thermal-hydraulic experiments is of great importance for phenomena interpretation and code validation. Direct measurement employing Particle Image Velocimetry (PIV) is challenging in some multiphase…
An exact and efficient new method to simulate dynamics in dissipative quantum systems is presented. A stochastic Liouville equation, deduced from Feynman and Vernon's path-integral expression of the reduced density matrix, is used to…
We present Path Integral Sampler~(PIS), a novel algorithm to draw samples from unnormalized probability density functions. The PIS is built on the Schr\"odinger bridge problem which aims to recover the most likely evolution of a diffusion…
Accounting for nuclear quantum effects (NQEs) can significantly alter material properties at finite temperatures. Atomic modeling using the path-integral molecular dynamics (PIMD) method can fully account for such effects, but requires…
Many coarse-grained models have been developed for equilibrium studies of lipid bilayer membranes. To achieve in simulations access to length-scales and time-scales difficult to attain in fully atomistic molecular dynamics, these…
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…
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 introduce a projected complex Langevin (CL) numerical sampling method -- a fictitious Langevin dynamics scheme that uses numerical projection to sample a constrained stationary distribution with highly oscillatory character. Despite the…
Nuclear quantum effects (NQEs) play an essential role in many atomistic systems, yet their explicit inclusion in molecular simulations remains challenging. Path-integral molecular dynamics (PIMD) provides a rigorous framework for…
Stick-slip is a ubiquitous phenomenon in many scientific fields, such as earthquake and glacier dynamics, acoustics, cell biology, interface science and tribology. As a fundamental mechanism of energy dissipation in nanofriction, it can be…
The evaluation of the path-integral representation for stochastic processes in the weak-noise limit shows that these systems are governed by a set of equations which are those of a classical dynamics. We show that, even when the noise is…
The particle-in-cell (PIC) method has been widely used for plasma simulation, because of its noise-reduction capability and moderate computational cost. The immersed finite element (IFE) method is efficient for solving interface problems on…
The Langevin equation accounts for unresolved bath degrees of freedom driving the system toward the bath temperature. Because of this, numerical solutions of the Langevin equation have a long history. Here, we recapitulate, combine, and…
We present a new path integral method to analyze stochastically perturbed ordinary differential equations with multiple time scales. The objective of this method is to derive from the original system a new stochastic differential equation…
Calibrating mathematical models of biological processes is essential for achieving predictive accuracy and gaining mechanistic insight. However, this task remains challenging due to limited and noisy data, significant biological…
Methods to efficiently determine the relative stability of polymorphs of organic crystals are highly desired in crystal structure predictions (CSPs). Current methodologies include use of static lattice phonons, quasi-harmonic approximation…