Computing three-dimensional densities from force densities improves statistical efficiency

Abstract

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 collecting histograms, with each discrete voxel in three-dimensional space needing to be visited significantly. This paper presents an extension of a previous methodology for the extraction of 3D solvent number densities with a reduced variance principle [Borgis et al., Mol. Phys. 111, 3486-3492 (2013)] to other 3D densities such as charge and polarization densities. The approach is also generalized to cover molecular solvents with structures described using rigid geometrical constraints, which include in particular popular water models such as SPC/E and TIPnP class of models. The noise reduction is illustrated for the microscopic hydration structure of a small molecule, in various simulation conditions, and for a protein. The method has large applicability to simulations of solvation in many fields, for example around biomolecules, nanoparticles, or within porous materials.

arXiv commentMain paper 13 pages 8 figures. Supporting information 2 pages 4 figures. Published in the Journal of Chemical Physics

References in corpus (7)

• Canonical sampling through velocity-rescaling
• The electric double layer has a life of its own
• Molecular Density Functional Theory of Water
• Molecular dynamics simulation of behaviour of water in nano-confined ionic liquid--water mixtures
• Efficient molecular density functional theory using generalized spherical harmonics expansions
• Better than counting: Density profiles from force sampling
• A Molecular Density Functional Theory Approach to Electron Transfer Reactions

Cited by in corpus (11)

• Use the force! Reduced variance estimators for densities, radial distribution functions and local mobilities in molecular simulations
• Sampling mobility profiles of confined fluids with equilibrium molecular dynamics simulations
• Gauge Invariance of Equilibrium Statistical Mechanics
• Hyperforce balance via thermal Noether invariance of any observable
• Reduced variance analysis of molecular dynamics simulations by linear combination of estimators
• Why hyperdensity functionals describe any equilibrium observable
• Custom Flow in Molecular Dynamics
• Noether invariance theory for the equilibrium force structure of soft matter
• Reduced-variance orientational distribution functions from torque sampling
• Why gauge invariance applies to statistical mechanics
• Robust Kirkwood-Buff inversion in complex mixtures via reciprocal-space methods