Use the force! Reduced variance estimators for densities, radial distribution functions and local mobilities in molecular simulations

Abstract

Even though the computation of local properties, such as densities or radial distribution functions, remains one of the most standard goals of molecular simulation, it still largely relies on straighforward histogram-based strategies. Here we highlight recent developments of alternative approaches leading, from different perspectives, to estimators with a reduced variance compared to conventional binning. They all make use of the force acting on the particles, in addition to their position, and allow to focus on the non-trivial part of the problem in order to alleviate (or even remove in some cases) the catastrophic behaviour of histograms as the bin size decreases. The corresponding computational cost is negligible for molecular dynamics simulations, since the forces are already computed to generate the configurations, and the benefit of reduced-variance estimators is even larger when the cost of generating the latter is high, in particular with ab initio simulations. The force sampling approach may result in spurious residual non-zero values of the density in regions where no particles are present, but strategies are available to mitigate this artefact. We illustrate this approach on number, charge and polarization densities, radial distribution functions and local transport coefficients, discuss the connections between the various perspectives and suggest future challenges for this promising approach.

References in corpus (9)

• The Electrostatic Screening Length in Concentrated Electrolytes Increases with Concentration
• First-order dynamical phase transition in models of glasses: an approach based on ensembles of histories
• The electric double layer has a life of its own
• Molecular Density Functional Theory of Water
• Revealing the three-dimensional structure of liquids using four-point correlation functions
• Efficient molecular density functional theory using generalized spherical harmonics expansions
• Zero-variance zero-bias quantum Monte Carlo estimators of the spherically and system-averaged pair density
• Nanoscale heterogeneity at the aqueous electrolyte-electrode interface
• Sampling mobility profiles of confined fluids with equilibrium molecular dynamics simulations

Cited by in corpus (33)

• Generalized Langevin Equation with a Non-Linear Potential of Mean Force and Non-Linear Memory Friction From a Hybrid Projection Scheme
• Power functional theory for many-body dynamics
• Perspective: How to overcome dynamical density functional theory
• Adaptive Brownian Dynamics
• Force density functional theory in- and out-of-equilibrium
• Noether-Constrained Correlations in Equilibrium Liquids
• Why neural functionals suit statistical mechanics
• Force balance in thermal quantum many-body systems from Noether's theorem
• Comparative study of force-based classical density functional theory
• Gauge Invariance of Equilibrium Statistical Mechanics
• Neural density functionals: Local learning and pair-correlation matching
• Hyperforce balance via thermal Noether invariance of any observable
• Reduced variance analysis of molecular dynamics simulations by linear combination of estimators
• Neural force functional for non-equilibrium many-body colloidal systems
• Why hyperdensity functionals describe any equilibrium observable
• Nuclear quantum effects in thermal conductivity from centroid molecular dynamics
• Fast quasi-centroid molecular dynamics for water and ice
• Local measures of fluctuations in inhomogeneous liquids: Statistical mechanics and illustrative applications
• Noether invariance theory for the equilibrium force structure of soft matter
• Normalizing flows as an enhanced sampling method for atomistic supercooled liquids
• Custom Flow in Molecular Dynamics
• Inhomogeneous steady shear dynamics of a three-body colloidal gel former
• Reduced-variance orientational distribution functions from torque sampling
• Why gauge invariance applies to statistical mechanics
• Neural Density Functional Theory in Higher Dimensions with Convolutional Layers
• Gauge invariance and hyperforce correlation theory for equilibrium fluid mixtures
• Multi-head committees enable direct uncertainty prediction for atomistic foundation models
• DL_POLY Quantum 2.1 software: A suite of real-time path integral methods for the simulation of dynamical properties and vibrational spectra
• h-CMD: An efficient hybrid fast centroid and quasi-centroid molecular dynamics method for the simulation of vibrational spectra
• Determining the chemical potential via universal density functional learning
• Robust Kirkwood-Buff inversion in complex mixtures via reciprocal-space methods
• Direct Boltzmann inversion method from particle configurations at arbitrary state points
• Connection between water's dynamical and structural properties: insights from ab initio simulations