arXiv preprint ·
Transient subtraction method: An alternative approach for variance reduction in dissipative particle dynamics
- Published
- Source
- arXiv preprint, posted before peer review
- Original
- Read it at arXiv
- Topics
Abstract
The efficient and accurate computation of transport coefficients remains a numerical challenge in dissipative particle dynamics (DPD), which can require substantial sampling to control statistical uncertainty. We extend recent results obtained in Langevin dynamics to consider the transient subtraction method for reducing the variance of transport coefficient estimations in DPD. The system is initialised by perturbing equilibrium samples, and its relaxation towards equilibrium is monitored under the reference dynamics. A reference trajectory, initialised at equilibrium and driven by the same Brownian noise, provides a control variate for the transient response. We apply the transient subtraction method to the computation of the mobility and the shear viscosity in a Lennard-Jones fluid in DPD. The shear viscosity is obtained from the time-integrated response of a Fourier observable associated with the longitudinal velocity. Numerical comparisons with the naive transient estimator show that the transient subtraction method reduces the variance at various integration times, and the reduction is more pronounced with a higher friction coefficient. These results indicate that the transient subtraction method is able to improve the statistical precision of transport coefficient estimations in DPD and reduce the number of trajectories needed to attain a prescribed standard error.
The authors’ abstract, as arXiv lists it. arXiv’s listings, abstracts included, are free for anyone to reuse.
Read it on ReadSearch
Each opens after you create an account, and brings you back to this paper.
More in Computational Physics
Machine Learning-Enabled Time-Skipping Method: Theory and Perturbation-Based Verification
ROSE: Deep Learning exchange-correlation functional from reduced density-matrix representations
Feature-preserving interface reconstruction using piecewise circles and cusps
Structure-Preserving Full Approximation Scheme for Real-Space Hartree-Fock Equations
Characterization of Phase Transitions with the Annealed Importance Sampling Algorithm
Get the few papers worth your time in Computational Physics every morning.
ReadSearch reads the new papers in your field each day and brings you the few worth reading, each with its key finding in a sentence.