Abstract:Photothermal methods permit measurement of molecular volume changes of solvated molecules over nanosecond timescales. Such experiments are an important tool in investigating complex biophysical phenomena including identifying transient species in solution. Developing a microscopic understanding of the origin of volume changes in the condensed phase is needed to complement the experimental measurements. A molecular dynamics (MD) method exploiting available simulation methodology is demonstrated here that both mimics experimental measurements and provides microscopic resolution to the thermodynamic measurements. To calculate thermodynamic volume changes over time, isothermal-isobaric (NPT) MD is performed on a solution for a chosen length of time and the volume of the system is thus established. A further simulation is then performed by "plucking" out a solute molecule of interest to determine the volume of the system in its absence. The difference between these volumes is the thermodynamic volume of the solute molecule. NPT MD allows the volume of the system to fluctuate over time and this results in a statistical uncertainty in volumes that are calculated. It is found in the systems investigated here that simulations lasting a few nanoseconds can discern volume changes of approximately 1.0 ml/mole. This precision is comparable to that achieved empirically, making the experimental and theoretical techniques synergistic. The technique is demonstrated here on model systems including neat water, both charged and neutral aqueous methane, and an aqueous beta-sheet peptide.

A molecular dynamics method for simulations in the canonical ensemble

Rigid Multibody Molecular Dynamics Algorithm in Canonical Ensemble

Born-Oppenheimer Molecular Dynamics with a Linear Combination of Atomic Orbitals and Hybrid Functionals for Condensed Matter Simulations Made Possible. Theory and Performance for the Microcanonical and Canonical Ensembles

Lattice Models for Magnetic Fluids: Correlations Between Order Parameters

Theory and Monte Carlo simulation of the ideal gas with shell particles in the canonical, isothermal-isobaric, grand canonical, and Gibbs ensembles

Origin of the Canonical Ensemble: Thermalization with Decoherence

Grand canonical steady-state simulation of nucleation

Unified Efficient Thermostat Scheme for the Canonical Ensemble with Holonomic or Isokinetic Constraints via Molecular Dynamics

Constant-temperature molecular-dynamics algorithms for mixed hard-core/continuous potentials

Canonical mean-field molecular dynamics derived from quantum mechanics

A Unified Efficient Thermostat Schemefor the Canonical Ensemblewith Holonomic or Isokinetic Constraints via Molecular Dynamics.

A Conjecture about Molecular Dynamics

Gaussian concentration bound and Ensemble equivalence in generic quantum many-body systems including long-range interaction

A Unified Thermostat Scheme for Efficient Configurational Sampling for Classical/quantum Canonical Ensembles Via Molecular Dynamics.

Molecular Dynamics Simulation on Cavitation Bubble Formation in Canonical Ensemble

Asymptotic methods for confined fluids

Grand Canonical Adaptive Resolution Simulation for Molecules with Electrons: A Theoretical Framework based on Physical Consistency

Virtual Molecular Dynamics

Generalized Gibbs canonical ensemble: A possible physical scenario

A molecular dynamics method for calculating molecular volume changes appropriate for biomolecular simulation

On a direct method of calculating pressure in the canonical ensemble