Researcher profile

Gregory A. Voth

Gregory A. Voth contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2023arXiv

Understanding Dynamics in Coarse-Grained Models: I. Universal Excess Entropy Scaling Relationship

Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics.

preprint2023arXiv

Understanding Dynamics in Coarse-Grained Models: II. Coarse-Grained Diffusion Modeled Using Hard Sphere Theory

The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids.

preprint2022arXiv

Centroid Molecular Dynamics Can Be Greatly Accelerated Through Neural Network Learned Centroid Forces Derived from Path Integral Molecular Dynamics

For nearly the past 30 years, Centroid Molecular Dynamics (CMD) has proven to be a viable classical-like phase space formulation for the calculation of quantum dynamical properties. However, calculation of the centroid effective force remains a significant computational cost and limits the ability of CMD to be an efficient approach to study condensed phase quantum dynamics. In this paper we introduce a neural network-based methodology for first learning the centroid effective force from path integral molecular dynamics data, which is subsequently used as an effective force field to evolve the centroids directly with the CMD algorithm. This method, called Machine-Learned Centroid Molecular Dynamics (ML-CMD) is faster and far less costly than both standard on the fly CMD and ring polymer molecular dynamics (RPMD). The training aspect of ML-CMD is also straightforwardly implemented utilizing the DeepMD software kit. ML-CMD is then applied to two model systems to illustrate the approach: liquid para-hydrogen and water. The results show comparable accuracy to both CMD and RPMD in the estimation of quantum dynamical properties, including the self-diffusion constant and velocity time correlation function, but for significantly reduced overall computational cost.

preprint2021arXiv

Understanding the Essential Nature of the Hydrated Excess Proton Through Simulation and Interpretation of Recent Spectroscopic Experiments

Two-dimensional infrared spectroscopy experiments have presented new results regarding the dynamics of the hydrated excess proton (aka <q>hydronium</q> cation solvated in water). It has been suggested by these experiments that the hydrated excess proton has an anisotropy reorientation timescale of 2.5 ps, which can be viewed as being somewhat long lived. Through the use of both the reactive molecular dynamics Multistate-Empirical Valence Bond method and Experiment Directed Simulation Ab Initio Molecular Dynamics we show that timescales of the same magnitude are obtained that correspond to proton transport, while also involving structural reorientations of the hydrated proton structure that correspond to the so-called <q>special pair dance</q>. The latter is a process predicted by prior computational studies in which the central hydrated hydronium in a distorted Eigen cation (H<sub>9</sub>O<sub>4</sub><sup>+</sup>) structure continually switches special pair partners with its strongly hydrogen-bonded neighboring water molecules. These dynamics are further characterized through the time-evolution of instantaneous normal modes. It is concluded that the hydrated excess proton has a spectral signature unique from the other protons in the hydrated proton complex. However, the results call into question the use of a static picture based on a simple effective one dimensional potential well to describe the hydrated excess proton in water. Instead, they more conclusively point to a distorted and dynamic Eigen cation as the most prevalent hydrated proton species in acid solutions of dilute to moderate concentrations.

preprint2020arXiv

A Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation-based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation, will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. It is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

preprint2020arXiv

Density Functional Theory-based Quantum Mechanics/Coarse-grained Molecular Mechanics: Theory and Implementation

Quantum mechanics/molecular mechanics (QM/MM) is a standard computational tool for describing chemical reactivity in systems with many degrees of freedom, including polymers, enzymes, and reacting molecules in complex solvents. However, QM/MM is less suitable for systems with complex MM dynamics due to associated long relaxation times, the high computational cost of QM energy evaluations, and expensive long-range electrostatics. Recently, a systematic coarse-graining of the MM part was proposed to overcome these QM/MM limitations in the form of the quantum mechanics/coarse-grained molecular mechanics (QM/CG-MM) approach. Herein, we recast QM/CG-MM in the density functional theory formalism and, by employing the force-matching variational principle, access the method performance for two model systems: QM CCl4 in the MM CCl4 liquid and the reaction of tert-butyl hypochlorite with the benzyl radical in the MM CCl4 solvent. We find that DFT-QM/CG-MM accurately reproduces DFT-QM/MM radial distribution functions and 3-body correlations between QM and CG-MM subsystems. The free energy profile of the reaction is also described well, with an error < 1-2 kcal/mol. DFT-QM/CG-MM is a general, systematic, and computationally efficient approach to include chemical reactivity in coarse-grained molecular models.

preprint2020arXiv

Water Assisted Proton Transport in Confined Nanochannels

Hydrated excess protons under hydrophobic confinement are a critical component of charge transport behavior and reactivity in nanoporous materials and biomolecular systems. Herein excess proton confinement effects are computationally investigated for sub-2 nm hydrophobic nanopores by varying the diameters (d = 0.81, 0.95, 1.09, 1.22, 1.36, 1.63, and 1.90 nm), lengths (l ~3 and ~5 nm), curvature, and chirality of cylindrical carbon nanotube (CNT) nanopores. CNTs with a combination of different diameter segments are also explored. The spatial distribution of water molecules under confinement are diameter-dependent; however, proton solvation and transport is consistently found to occur in the water layer adjacent to the pore wall, showing an &#34;amphiphilic&#34; character of the hydrated excess proton hydronium-like structure. The proton transport free energy barrier also decreases significantly as the nanopore diameter increases and proton transport becomes almost barrierless in the d > 1 nm nanopores. Among the nanopores studied, the Zundel cation (${H_{5}O_{2}}^{+}$) is populated only in the d = 0.95 nm CNT (7,7) nanopore. The presence of the hydrated excess proton and $K^{+}$ inside the CNT (7,7) nanopore induces a water density increase by 40 and 20%, respectively. The $K^{+}$ transport through CNT nanopores is also consistently higher in free energy barrier than proton transport. Interestingly, the evolution of excess protonic charge defect distribution reveals a &#34;frozen&#34; single water wire configuration in the d = 0.81 nm CNT (6,6) nanopore (or segment), through which hydrated excess protons can only shuttle via the Grotthuss mechanism. Vehicular diffusion becomes relevant to proton transport in the &#34;flat&#34; free energy regions and in the wider nanopores, where protons do not primarily shuttle in the axial direction.

preprint2019arXiv

Adversarial-Residual-Coarse-Graining: Applying machine learning theory to systematic molecular coarse-graining

We utilize connections between molecular coarse-graining approaches and implicit generative models in machine learning to describe a new framework for systematic molecular coarse-graining (CG). Focus is placed on the formalism encompassing generative adversarial networks. The resulting method enables a variety of model parameterization strategies, some of which show similarity to previous CG methods. We demonstrate that the resulting framework can rigorously parameterize CG models containing CG sites with no prescribed connection to the reference atomistic system (termed virtual sites); however, this advantage is offset by the lack of a closed-form expression for the CG Hamiltonian at the resolution obtained after integration over the virtual CG sites. Computational examples are provided for cases in which these methods ideally return identical parameters as Relative Entropy Minimization (REM) CG but where traditional REM CG optimization equations are not applicable.