# Reduced variance analysis of molecular dynamics simulations by linear combination of estimators.

@article{Coles2021ReducedVA, title={Reduced variance analysis of molecular dynamics simulations by linear combination of estimators.}, author={Samuel W. Coles and Etienne Mangaud and Daan Frenkel and Benjamin Rotenberg}, journal={The Journal of chemical physics}, year={2021}, volume={154 19}, pages={ 191101 } }

Building upon recent developments of force-based estimators with a reduced variance for the computation of densities, radial distribution functions, or local transport properties from molecular simulations, we show that the variance can be further reduced by considering optimal linear combinations of such estimators. This control variates approach, well known in statistics and already used in other branches of computational physics, has been comparatively much less exploited in molecular…

#### One Citation

Adaptive Brownian Dynamics.

- Medicine, PhysicsThe Journal of chemical physics
- 2021

A framework for performant Brownian Dynamics many-body simulations with adaptive timestepping is presented and it is shown that adaptive BD solves performance and stability issues of conventional BD, already outperforming the latter even in standard situations.

#### References

SHOWING 1-10 OF 32 REFERENCES

Partial differential equations and stochastic methods in molecular dynamics*

- Mathematics, Computer ScienceActa Numerica
- 2016

This review describes how techniques from the analysis of partial differential equations can be used to devise good algorithms and to quantify their efficiency and accuracy.

Accurate thermal conductivities from optimally short molecular dynamics simulations

- Physics, MedicineScientific Reports
- 2017

A new method to compute transport coefficients from optimally short molecular dynamics simulations, based on the Green-Kubo theory of linear response and the cepstral analysis of time series, finds that simulation times of one to a few hundred picoseconds are sufficient in these systems to achieve an accuracy of 10% on the estimated thermal conductivities.

Reformulation of Ensemble Averages via Coordinate Mapping.

- Computer Science, MedicineJournal of chemical theory and computation
- 2016

A general framework is established for reformulation of the ensemble averages commonly encountered in statistical mechanics. This "mapped-averaging" scheme allows approximate theoretical results that…

Computing three-dimensional densities from force densities improves statistical efficiency

- Mathematics, PhysicsThe Journal of Chemical Physics
- 2019

The extraction of inhomogeneous 3-dimensional densities around tagged solutes from molecular simulations is known to have a very high computational cost because this is traditionally performed by…

Computation of pair distribution functions and three-dimensional densities with a reduced variance principle

- Mathematics
- 2013

No fancy statistical objects here, we go back to the computation of one of the most basic and fundamental quantities in the statistical mechanics of fluids, namely the pair distribution functions.…

On efficient simulations of multiscale kinetic transport

- Computer Science, MedicinePhilosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
- 2013

It is shown that it is possible to derive evolution equations for the two parts from the governing kinetic equation, leading to a decomposition that is dynamically and automatically adaptive, and a multiscale method that seamlessly bridges the two descriptions without introducing any approximation.

Alternative ensemble averages in molecular dynamics simulation of hard spheres

- PhysicsMolecular Physics
- 2019

ABSTRACT We consider application to the hard sphere (HS) model of the mapped-averaging framework for generating alternative ensemble averages for thermodynamic properties. Specifically, we develop…

A Perturbative Approach to Control Variates in Molecular Dynamics

- Computer Science, MathematicsMultiscale Model. Simul.
- 2019

It is shown by a perturbative argument that a control variate computed for a simplified version of the model can provide an efficient control Variate for the actual problem at hand.

Control functionals for Monte Carlo integration

- Mathematics
- 2014

Summary
A non-parametric extension of control variates is presented. These leverage gradient information on the sampling density to achieve substantial variance reduction. It is not required that…

Efficient simulation of multidimensional phonon transport using energy-based variance-reduced Monte Carlo formulations

- Physics
- 2011

We present a Monte Carlo method for obtaining solutions of the Boltzmann equation to describe phonon transport in micro- and nanoscale devices. The proposed method can resolve arbitrarily small…