-
The need to implement FAIR principles in biomolecular simulations
Authors:
Rommie Amaro,
Johan Åqvist,
Ivet Bahar,
Federica Battistini,
Adam Bellaiche,
Daniel Beltran,
Philip C. Biggin,
Massimiliano Bonomi,
Gregory R. Bowman,
Richard Bryce,
Giovanni Bussi,
Paolo Carloni,
David Case,
Andrea Cavalli,
Chie-En A. Chang,
Thomas E. Cheatham III,
Margaret S. Cheung,
Cris Chipot,
Lillian T. Chong,
Preeti Choudhary,
Gerardo Andres Cisneros,
Cecilia Clementi,
Rosana Collepardo-Guevara,
Peter Coveney,
Roberto Covino
, et al. (101 additional authors not shown)
Abstract:
This letter illustrates the opinion of the molecular dynamics (MD) community on the need to adopt a new FAIR paradigm for the use of molecular simulations. It highlights the necessity of a collaborative effort to create, establish, and sustain a database that allows findability, accessibility, interoperability, and reusability of molecular dynamics simulation data. Such a development would democra…
▽ More
This letter illustrates the opinion of the molecular dynamics (MD) community on the need to adopt a new FAIR paradigm for the use of molecular simulations. It highlights the necessity of a collaborative effort to create, establish, and sustain a database that allows findability, accessibility, interoperability, and reusability of molecular dynamics simulation data. Such a development would democratize the field and significantly improve the impact of MD simulations on life science research. This will transform our working paradigm, pushing the field to a new frontier. We invite you to support our initiative at the MDDB community (https://mddbr.eu/community/)
△ Less
Submitted 30 August, 2024; v1 submitted 23 July, 2024;
originally announced July 2024.
-
Understanding Dynamics in Coarse-Grained Models: IV. Connection of Fine-Grained and Coarse-Grained Dynamics with the Stokes-Einstein and Stokes-Einstein-Debye Relations
Authors:
Jaehyeok Jin,
Gregory A. Voth
Abstract:
Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations…
▽ More
Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes-Einstein and Stokes-Einstein-Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work. Additional theoretical links to our previous work are also established. First, we demonstrate that the effective hard sphere radius determined by the classical perturbation theory can approximate the complex hydrodynamic radius value reasonably well. Also, we present a simple derivation of an excess entropy scaling relationship for viscosity by estimating the elliptical integral of molecules. In turn, since the translational and rotational motions at the FG level are correlated to each other, we conclude that the "entropy-free" CG diffusion only depends on the shape of the reference molecule. Our results and analyses impart an alternative way of recovering the FG diffusion from the CG description by coupling the translational and rotational motions at the hydrodynamic level.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
Coarse-Graining with Equivariant Neural Networks: A Path Towards Accurate and Data-Efficient Models
Authors:
Timothy D. Loose,
Patrick G. Sahrmann,
Thomas S. Qu,
Gregory A. Voth
Abstract:
Machine learning has recently entered into the mainstream of coarse-grained (CG) molecular modeling and simulation. While a variety of methods for incorporating deep learning into these models exist, many of them involve training neural networks to act directly as the CG force field. This has several benefits, the most significant of which is accuracy. Neural networks can inherently incorporate mu…
▽ More
Machine learning has recently entered into the mainstream of coarse-grained (CG) molecular modeling and simulation. While a variety of methods for incorporating deep learning into these models exist, many of them involve training neural networks to act directly as the CG force field. This has several benefits, the most significant of which is accuracy. Neural networks can inherently incorporate multi-body effects during the calculation of CG forces, and a well-trained neural network force field outperforms pairwise basis sets generated from essentially any methodology. However, this comes at a significant cost. First, these models are typically slower than pairwise force fields even when accounting for specialized hardware which accelerates the training and integration of such networks. The second, and the focus of this paper, is the need for the considerable amount of data needed to train such force fields. It is common to use 10s of microseconds of molecular dynamics data to train a single CG model, which approaches the point of eliminating the CG models usefulness in the first place. As we investigate in this work, this data-hunger trap from neural networks for predicting molecular energies and forces can be remediated in part by incorporating equivariant convolutional operations. We demonstrate that for CG water, networks which incorporate equivariant convolutional operations can produce functional models using datasets as small as a single frame of reference data, while networks without these operations cannot.
△ Less
Submitted 1 November, 2023; v1 submitted 17 September, 2023;
originally announced September 2023.
-
Understanding Dynamics in Coarse-Grained Models: III. Roles of Rotational Motion and Translation-Rotation Coupling in Coarse-Grained Dynamics
Authors:
Jaehyeok Jin,
Eok Kyun Lee,
Gregory A. Voth
Abstract:
This paper series aims to establish a complete correspondence between fine-grained (FG) and coarse-grained (CG) dynamics by way of excess entropy scaling (introduced in Paper I). While Paper II successfully captured translational motions in CG systems using a hard sphere mapping, the absence of rotational motions in single-site CG models introduces differences between FG and CG dynamics. In this t…
▽ More
This paper series aims to establish a complete correspondence between fine-grained (FG) and coarse-grained (CG) dynamics by way of excess entropy scaling (introduced in Paper I). While Paper II successfully captured translational motions in CG systems using a hard sphere mapping, the absence of rotational motions in single-site CG models introduces differences between FG and CG dynamics. In this third paper, our objective is to faithfully recover atomistic diffusion coefficients from CG dynamics by incorporating rotational dynamics. By extracting FG rotational diffusion, we unravel, for the first time reported to our knowledge, a universality in excess entropy scaling between the rotational and translational diffusion. Once the missing rotational dynamics are integrated into the CG translational dynamics, an effective translation-rotation coupling becomes essential. We propose two different approaches for estimating this coupling parameter: the rough hard sphere theory with acentric factor (temperature-independent) or the rough Lennard-Jones model with CG attractions (temperature-dependent). Altogether, we demonstrate that FG diffusion coefficients can be recovered from CG diffusion coefficients by (1) incorporating "entropy-free" rotational diffusion with translation-rotation coupling and (2) recapturing the missing entropy. Our findings shed light on the fundamental relationship between FG and CG dynamics in molecular fluids.
△ Less
Submitted 31 July, 2023;
originally announced July 2023.
-
Orientation of Finite Reynolds Number Anisotropic Particles Settling in Turbulence
Authors:
Anubhab Roy,
Stefan Kramel,
Udayshankar Menon,
Greg A. Voth,
Donald L. Koch
Abstract:
We present experimental and computational results for the orientation distributions of slender fibers and ramified particles settling in an isotropic turbulent flow. The rotational dynamics of the particles is modeled using a slender-body theory that includes the inertial torque due to sedimentation that tends to rotate the particles toward a broadside orientation. The particles are assumed to rot…
▽ More
We present experimental and computational results for the orientation distributions of slender fibers and ramified particles settling in an isotropic turbulent flow. The rotational dynamics of the particles is modeled using a slender-body theory that includes the inertial torque due to sedimentation that tends to rotate the particles toward a broadside orientation. The particles are assumed to rotate due to viscous forces associated with the turbulent velocity gradients occurring on the particle length scale. In the simulations, the turbulence is obtained from a stochastic model of the velocity gradient in a Lagrangian reference frame. In the experiments, the turbulence is generated by active jets in a vertical water tunnel. It is well known that axisymmetric particles rotate according to Jeffery's solution for the rotation of a spheroidal particle if one adopts an appropriate effective aspect ratio. We show that the same result applies to a ramified particle consisting of three coplanar fibers connected with equal angles at a central point which rotates like a thin oblate spheroid. The orientation statistics can be quantified with a single non-dimensional parameter, the settling factor $S_F$, defined as the ratio of rotations due to sedimentation and turbulent shear. For low values of $S_F$, we observe nearly isotropically oriented particles, whereas particles become strongly aligned near the horizontal plane for high values of $S_F$. The variance of the angle away from horizontal scales as $S_F^{-2}$ for $S_F \gg 1$, but the orientation distribution is non-Gaussian due to turbulent intermittency in this limit.
△ Less
Submitted 30 March, 2023;
originally announced March 2023.
-
Utilizing Machine Learning to Greatly Expand the Range and Accuracy of Bottom-Up Coarse-Grained Models Through Virtual Particles
Authors:
Patrick G. Sahrmann,
Timothy D. Loose,
Aleksander E. P. Durumeric,
Gregory A. Voth
Abstract:
Coarse-grained (CG) models parameterized using atomistic reference data, i.e., 'bottom up' CG models, have proven useful in the study of biomolecules and other soft matter. However, the construction of highly accurate, low resolution CG models of biomolecules remains challenging. We demonstrate in this work how virtual particles, CG sites with no atomistic correspondence, can be incorporated into…
▽ More
Coarse-grained (CG) models parameterized using atomistic reference data, i.e., 'bottom up' CG models, have proven useful in the study of biomolecules and other soft matter. However, the construction of highly accurate, low resolution CG models of biomolecules remains challenging. We demonstrate in this work how virtual particles, CG sites with no atomistic correspondence, can be incorporated into CG models within the context of relative entropy minimization (REM) as latent variables. The methodology presented, variational derivative relative entropy minimization (VD-REM), enables optimization of virtual particle interactions through a gradient descent algorithm aided by machine learning. We apply this methodology to the challenging case of a solvent-free CG model of a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid bilayer and demonstrate that introduction of virtual particles captures solvent-mediated behavior and higher-order correlations which REM alone cannot capture in a more standard CG model based only on the mapping of collections of atoms to the CG sites.
△ Less
Submitted 8 December, 2022;
originally announced December 2022.
-
Centroid Molecular Dynamics Can Be Greatly Accelerated Through Neural Network Learned Centroid Forces Derived from Path Integral Molecular Dynamics
Authors:
Timothy D. Loose,
Patrick G. Sahrmann,
Gregory A. Voth
Abstract:
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 introd…
▽ More
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.
△ Less
Submitted 13 September, 2022; v1 submitted 16 August, 2022;
originally announced August 2022.
-
Understanding Dynamics in Coarse-Grained Models: II. Coarse-Grained Diffusion Modeled Using Hard Sphere Theory
Authors:
Jaehyeok Jin,
Kenneth S. Schweizer,
Gregory A. Voth
Abstract:
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 bo…
▽ More
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.
△ Less
Submitted 17 January, 2023; v1 submitted 2 August, 2022;
originally announced August 2022.
-
Understanding Dynamics in Coarse-Grained Models: I. Universal Excess Entropy Scaling Relationship
Authors:
Jaehyeok Jin,
Kenneth S. Schweizer,
Gregory A. Voth
Abstract:
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…
▽ More
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.
△ Less
Submitted 17 January, 2023; v1 submitted 29 July, 2022;
originally announced August 2022.
-
Explaining classifiers to understand coarse-grained models
Authors:
Aleksander Evren Paetzold Durumeric,
Gregory A. Voth
Abstract:
Bottom-up coarse-grained molecular dynamics models are parameterized using complex effective Hamiltonians. These models are typically optimized to approximate high dimensional data from atomistic simulations. In contrast, human validation of these models is often limited to low dimensional statistics that do not necessarily differentiate between the CG model and said atomistic simulations. We prop…
▽ More
Bottom-up coarse-grained molecular dynamics models are parameterized using complex effective Hamiltonians. These models are typically optimized to approximate high dimensional data from atomistic simulations. In contrast, human validation of these models is often limited to low dimensional statistics that do not necessarily differentiate between the CG model and said atomistic simulations. We propose that explainable machine learning can directly convey high-dimensional error to scientists and use Shapley additive explanations do so in two coarse-grained protein models.
△ Less
Submitted 15 September, 2021;
originally announced September 2021.
-
Understanding the Essential Nature of the Hydrated Excess Proton Through Simulation and Interpretation of Recent Spectroscopic Experiments
Authors:
Paul B. Calio,
Chenghan Li,
Gregory A. Voth
Abstract:
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 molecu…
▽ More
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.
△ Less
Submitted 5 January, 2021;
originally announced January 2021.
-
The Hopping Mechanism of the Hydrated Excess Proton and Its Contribution to Proton Diffusion in Water
Authors:
Christopher Arnsten,
Chen Chen,
Paul B. Calio,
Chenghan Li,
Gregory A. Voth
Abstract:
In this work a series of analyses are performed on ab initio molecular dynamics (AIMD) simulations of a hydrated excess proton in water to quantify the relative occurrence of concerted hopping events and <span>rattling</span> events, and thus to further elucidate the hopping mechanism of proton transport in water. Contrary to results reported in certain earlier papers, the new analysis finds that…
▽ More
In this work a series of analyses are performed on ab initio molecular dynamics (AIMD) simulations of a hydrated excess proton in water to quantify the relative occurrence of concerted hopping events and <span>rattling</span> events, and thus to further elucidate the hopping mechanism of proton transport in water. Contrary to results reported in certain earlier papers, the new analysis finds that concerted hopping events do occur in all simulations, but that the majority of events are the product of proton rattling, where the excess proton will rattle between two or more waters. The results are consistent with the proposed <span>special-pair dance</span> model of the hydrated excess proton, wherein the acceptor water molecule for the proton transfer will quickly change (resonate between three equivalent special pairs) until a decisive proton hop occurs. To remove the misleading effect of simple rattling, a filter was applied to the trajectory such that hopping events that were followed by back hops to the original water are not counted. A steep reduction in the number of multiple hopping events is found when the filter is applied, suggesting that many multiple hopping events that occur in the unfiltered trajectory are largely the product of rattling, contrary to prior suggestions. Comparing the continuous correlation function of the filtered and unfiltered trajectories, we find agreement with experimental values for the proton hopping time and Eigen-Zundel interconversion time, respectively.
△ Less
Submitted 21 December, 2020;
originally announced December 2020.
-
Density Functional Theory-based Quantum Mechanics/Coarse-grained Molecular Mechanics: Theory and Implementation
Authors:
Alexander V. Mironenko,
Gregory A. Voth
Abstract:
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 expe…
▽ More
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.
△ Less
Submitted 20 July, 2020;
originally announced July 2020.
-
Lord Kelvin's isotropic helicoid
Authors:
Darci Collins,
Rami J. Hamati,
Fabien Candelier,
Kristian Gustavsson,
Bernhard Mehlig,
Greg A. Voth
Abstract:
Nearly 150 years ago, Lord Kelvin proposed the isotropic helicoid, a particle with isotropic yet chiral interactions with a fluid, so that translation couples to rotation. An implementation of his design fabricated with a three-dimensional printer is found experimentally to have no detectable translation-rotation coupling, although the particle point-group symmetry allows this coupling. We explain…
▽ More
Nearly 150 years ago, Lord Kelvin proposed the isotropic helicoid, a particle with isotropic yet chiral interactions with a fluid, so that translation couples to rotation. An implementation of his design fabricated with a three-dimensional printer is found experimentally to have no detectable translation-rotation coupling, although the particle point-group symmetry allows this coupling. We explain these results by demonstrating that in Stokes flow, the chiral coupling of such isotropic helicoids made out of non-chiral vanes is due only to hydrodynamic interactions between these vanes. Therefore it is small. In summary, Kelvin's predicted isotropic helicoid exists, but only as a weak breaking of a symmetry of non-interacting vanes in Stokes flow.
△ Less
Submitted 21 July, 2021; v1 submitted 15 June, 2020;
originally announced June 2020.
-
A Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water
Authors:
Paul B. Calio,
Glen M. Hocky,
Gregory A. Voth
Abstract:
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 reproduc…
▽ More
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.
△ Less
Submitted 2 June, 2020;
originally announced June 2020.
-
Water Assisted Proton Transport in Confined Nanochannels
Authors:
Xinyou Ma,
Chenghan Li,
Alex B. F. Martinson,
Gregory A. Voth
Abstract:
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, a…
▽ More
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 "amphiphilic" 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 "frozen" 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 "flat" free energy regions and in the wider nanopores, where protons do not primarily shuttle in the axial direction.
△ Less
Submitted 2 June, 2020;
originally announced June 2020.
-
Using deformable particles for single particle measurements of velocity gradient tensors
Authors:
Bardia Hejazi,
Michael Krellenstein,
Greg A. Voth
Abstract:
We measure the deformation of particles made of several slender arms in a two-dimensional (2D) linear shear and a three-dimensional (3D) turbulent flow. We show how these measurements of arm deformations along with the rotation rate of the particle allow us to extract the velocity gradient tensor of the flow. The particles used in the experiments have three symmetric arms in a plane (triads) and a…
▽ More
We measure the deformation of particles made of several slender arms in a two-dimensional (2D) linear shear and a three-dimensional (3D) turbulent flow. We show how these measurements of arm deformations along with the rotation rate of the particle allow us to extract the velocity gradient tensor of the flow. The particles used in the experiments have three symmetric arms in a plane (triads) and are fabricated using 3D printing of a flexible polymeric material. Deformation measurements of a particle free to rotate about a fixed axis in a 2D simple shear flow are used to validate our model relating particle deformations to the fluid strain. We then examine deformable particles in a 3D turbulent flow created by a jet array in a vertical water tunnel. Particle orientations and deformations are measured with high precision using four high speed cameras and have an uncertainty on the order of $10^{-4}$ radians. Measured deformations in 3D turbulence are small and only slightly larger than our orientation measurement uncertainty. Simulation results for triads in turbulence show deformations similar to the experimental observations. Deformable particles offer a promising method for measuring the full local velocity gradient tensor from measurements of a single particle where traditionally a high concentration of tracer particles would be required.
△ Less
Submitted 4 June, 2019;
originally announced June 2019.
-
Adversarial-Residual-Coarse-Graining: Applying machine learning theory to systematic molecular coarse-graining
Authors:
Aleksander E. P. Durumeric,
Gregory A. Voth
Abstract:
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 metho…
▽ More
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.
△ Less
Submitted 27 June, 2019; v1 submitted 1 April, 2019;
originally announced April 2019.
-
Quantum Mechanics / Coarse-Grained Molecular Mechanics (QM/CG-MM)
Authors:
Anton V. Sinitskiy,
Gregory A. Voth
Abstract:
Numerous molecular systems, including solutions, proteins, and composite materials, can be modeled using mixed-resolution representations, of which the quantum mechanics/molecular mechanics (QM/MM) approach has become the most widely used. However, the QM/MM approach often faces a number of challenges, including the slow sampling of the large configuration space for the MM part, the high cost of r…
▽ More
Numerous molecular systems, including solutions, proteins, and composite materials, can be modeled using mixed-resolution representations, of which the quantum mechanics/molecular mechanics (QM/MM) approach has become the most widely used. However, the QM/MM approach often faces a number of challenges, including the slow sampling of the large configuration space for the MM part, the high cost of repetitive QM computations for changing coordinates of atoms in the MM surroundings, and a difficulty in providing a simple, qualitative interpretation of numerical results in terms of the influence of the molecular environment upon the active QM region. In this paper, we address these issues by combining QM/MM modeling with the methodology of "bottom-up" coarse-graining (CG) to provide the theoretical basis for a systematic quantum-mechanical/coarse-grained molecular mechanics (QM/CG-MM) mixed resolution approach. A derivation of the method is presented based on a combination of statistical mechanics and quantum mechanics, leading to an equation for the effective Hamiltonian of the QM part, a central concept in the QM/CG-MM theory. A detailed analysis of different contributions to the effective Hamiltonian from electrostatic, induction, dispersion and exchange interactions between the QM part and the surroundings is provided, serving as a foundation for a potential hierarchy of QM/CG-MM methods varying in their accuracy and computational cost. A relationship of the QM/CG-MM methodology to other mixed resolution approaches is also discussed.
△ Less
Submitted 27 September, 2017;
originally announced September 2017.
-
Emergent scar lines in chaotic advection of passive directors
Authors:
Bardia Hejazi,
Bernhard Mehlig,
Greg A. Voth
Abstract:
We examine the spatial field of orientations of slender fibers that are advected by a two-dimensional fluid flow. The orientation field of these passive directors are important in a wide range of industrial and geophysical flows. We introduce emergent scar lines as the dominant coherent structures in the orientation field of passive directors in chaotic flows. Previous work has identified the exis…
▽ More
We examine the spatial field of orientations of slender fibers that are advected by a two-dimensional fluid flow. The orientation field of these passive directors are important in a wide range of industrial and geophysical flows. We introduce emergent scar lines as the dominant coherent structures in the orientation field of passive directors in chaotic flows. Previous work has identified the existence of scar lines where the orientation rotates by π over short distances, but the lines that were identified disappeared as time progressed. As a result, earlier work focused on topological singularities in the orientation field which we find to play a negligible role at long times. We use the standard map as a simple time-periodic two-dimensional (2D) flow that produces Lagrangian chaos. This class of flows produce persistent patterns in passive scalar advection, and we find that a different kind of persistent pattern develops in the passive director orientation field. We identify the mechanism by which emergent scar lines grow to dominate these patterns at long times in complex flows. Emergent scar lines form where the recent stretching of the fluid element is perpendicular to earlier stretching. Thus these scar lines can be labeled by their age, defined as the time since their stretching reached a maximum.
△ Less
Submitted 4 December, 2017; v1 submitted 22 June, 2017;
originally announced June 2017.
-
Non-uniqueness of quantum transition state theory and general dividing surfaces in the path integral space
Authors:
Seogjoo Jang,
Gregory A. Voth
Abstract:
Despite the fact that quantum mechanical principles do not allow the establishment of an exact quantum analogue of the classical transition state theory (TST), the development of a quantum TST (QTST) with a proper dynamical justification, while recovering the TST in the classical limit, has been a long standing theoretical challenge in chemical physics. One of the most recent efforts of this kind…
▽ More
Despite the fact that quantum mechanical principles do not allow the establishment of an exact quantum analogue of the classical transition state theory (TST), the development of a quantum TST (QTST) with a proper dynamical justification, while recovering the TST in the classical limit, has been a long standing theoretical challenge in chemical physics. One of the most recent efforts of this kind was put forth by Hele and Althorpe (HA) [ J. Chem. Phys. 138 , 084108 (2013)], which can be specified for any cyclically invariant dividing surface defined in the space of the imaginary time path integral. The present work revisits the issue of the non-uniqueness of QTST and provides a detailed theoretical analysis of HA-QTST for a general class of such path integral dividing surfaces. While we confirm that HA-QTST reproduces the result based on the ring polymer molecular dynamics (RPMD) rate theory for dividing surfaces containing only a quadratic form of low frequency Fourier modes, we find that it produces different results for those containing higher frequency imaginary time paths which accommodate greater quantum fluctuations. This result confirms the assessment made in our previous work [J. Chem. Phys. 144, 084110 (2016)] that HA-QTST does not provide a derivation of RPMD-TST in general, and points to a new ambiguity of HA-QTST with respect to its justification for general cyclically invariant dividing surfaces defined in the space of imaginary time path integrals. Our analysis also offers new insights into similar path integral based QTST approaches.
△ Less
Submitted 11 April, 2017;
originally announced April 2017.
-
Improved Ab Initio Molecular Dynamics by Minimal Biasing with Experimental Data
Authors:
Andrew D. White,
Chris Knight,
Glen M. Hocky,
Gregory A. Voth
Abstract:
Accounting for electrons and nuclei simultaneously is a powerful capability of ab initio molecular dynamics (AIMD). However, AIMD is often unable to accurately reproduce properties of systems such as water due to inaccuracies in the underlying electronic density functionals. This shortcoming is often addressed by added empirical corrections and/or increasing the simulation temperature. We present…
▽ More
Accounting for electrons and nuclei simultaneously is a powerful capability of ab initio molecular dynamics (AIMD). However, AIMD is often unable to accurately reproduce properties of systems such as water due to inaccuracies in the underlying electronic density functionals. This shortcoming is often addressed by added empirical corrections and/or increasing the simulation temperature. We present here a maximum-entropy approach to directly incorporate limited experimental data via a minimal bias. Biased AIMD simulations of water and an excess proton in water are shown to give significantly improved properties both for observables which were biased to match experimental data and for unbiased observables. This approach also yields new physical insight into inaccuracies in the underlying density functional theory as utilized in the unbiased AIMD.
△ Less
Submitted 24 January, 2017; v1 submitted 1 July, 2016;
originally announced July 2016.
-
Preferential rotation of chiral dipoles in isotropic turbulence
Authors:
Stefan Kramel,
Saskia Tympel,
Federico Toschi,
Greg A. Voth
Abstract:
Particles in the shape of chiral dipoles show a preferential rotation in three dimensional homogeneous isotropic turbulence. A chiral dipole consists of a rod with two helices of opposite handedness, one at each end. We can use 3d printing to fabricate these particles with length in the inertial range and track their rotations in a turbulent flow between oscillating grids. High aspect ratio chiral…
▽ More
Particles in the shape of chiral dipoles show a preferential rotation in three dimensional homogeneous isotropic turbulence. A chiral dipole consists of a rod with two helices of opposite handedness, one at each end. We can use 3d printing to fabricate these particles with length in the inertial range and track their rotations in a turbulent flow between oscillating grids. High aspect ratio chiral dipoles will align with the extensional eigenvectors of the strain rate tensor and the helical ends will respond to the strain field by spinning around its long axis. The mean of the measured spinning rate is non-zero and reflects the average stretching the particles experience. We use Stokesian dynamics simulations of chiral dipoles in pure strain flow to quantify the dependence of spinning on particle shape. Based on the known response to pure strain, we build a model that gives the spinning rate of small chiral dipoles using Lagrangian velocity gradients from high resolution direct numerical simulations. The statistics of chiral dipole spinning determined with this model show surprisingly good agreement with the measured spinning of much larger chiral dipoles in the experiments.
△ Less
Submitted 24 February, 2016;
originally announced February 2016.
-
Can quantum transition state theory be defined as an exact t=0+ limit?
Authors:
Seogjoo Jang,
Gregory A. Voth
Abstract:
The definition of the classical transition state theory (TST) as a t= 0+ limit of the flux-side time correlation function relies on the assumption that simultaneous measurement of population and flux is a well defined physical process. However, the noncommutativity of the two measurements in quantum mechanics makes the extension of such a concept to the quantum regime impossible. For this reason,…
▽ More
The definition of the classical transition state theory (TST) as a t= 0+ limit of the flux-side time correlation function relies on the assumption that simultaneous measurement of population and flux is a well defined physical process. However, the noncommutativity of the two measurements in quantum mechanics makes the extension of such a concept to the quantum regime impossible. For this reason, quantum TST (QTST) has been generally accepted as any kind of quantum rate theory reproducing the TST in the classical limit, and there has been a broad consensus that no unique QTST retaining all the properties of TST can be defined. Hele and Althorpe (HA) [J. Chem. Phys. 138, 084108 (2013)] recently suggested that a true QTST can be defined as the exact t=0+ limit of a certain kind of quantum flux-side time correlation function and that it is equivalent to the ring polymer molecular dynamics (RPMD) TST. This work seeks to question and clarify certain assumptions underlying these suggestions and their implications. First, the time correlation function used by HA as a starting expression is not related to the kinetic rate constant by virtue of linear response theory, which is the first important step in relating a t=0+ limit to a physically measurable rate. Second, a theoretical analysis calls into question a key step in HA's proof which appears not to rely on an exact quantum mechanical identity. The correction of this makes the true t=0+ limit of HA's QTST different from the RPMD-TST rate expression, but rather equal to the well-known path integral quantum transition state theory rate expression for the case of centroid dividing surface. An alternative path integral quantum rate expression is then formulated based on the linear response theory. It is shown that the t=0+ limit of the new rate expression vanishes in the exact quantum limit.
△ Less
Submitted 10 February, 2016; v1 submitted 31 October, 2015;
originally announced November 2015.
-
Measurements of the coupling between the tumbling of rods and the velocity gradient tensor in turbulence
Authors:
Rui Ni,
Stefan Kramel,
Nicholas T. Ouellette,
Greg A. Voth
Abstract:
We present simultaneous experimental measurements of the dynamics of anisotropic particles transported by a turbulent flow and the velocity gradient tensor of the flow surrounding them. We track both rod-shaped particles and small spherical flow tracers using stereoscopic particle tracking. By using scanned illumination, we are able to obtain a high enough seeding density of tracers to measure the…
▽ More
We present simultaneous experimental measurements of the dynamics of anisotropic particles transported by a turbulent flow and the velocity gradient tensor of the flow surrounding them. We track both rod-shaped particles and small spherical flow tracers using stereoscopic particle tracking. By using scanned illumination, we are able to obtain a high enough seeding density of tracers to measure the full velocity gradient tensor near the rod. The alignment of rods with the vorticity and the eigenvectors of the strain rate show agreement with numerical simulations. A full description of the tumbling of rods in turbulence requires specifying a seven-dimensional joint probability density function (PDF) of five scalars characterizing the velocity gradient tensor and two scalars describing the relative orientation of the rod. If these seven parameters are known, then Jeffery's equation specifies the rod tumbling rate and any statistic of rod rotations can be obtained as a weighted average over the joint PDF. To look for a lower-dimensional projection to simplify the problem, we explore conditional averages of the mean-squared tumbling rate. The conditional dependence of the mean-squared tumbling rate on the magnitude of both the vorticity and the strain rate is strong, as expected, and similar. There is also a strong dependence on the orientation between the rod and the vorticity, since a rod aligned with the vorticity vector tumbles due to strain but not vorticity. When conditioned on the alignment of the rod with the eigenvectors of the strain rate, the largest tumbling rate is obtained when the rod is oriented at a certain angle to the eigenvector that corresponds to the smallest eigenvalue, because this particular orientation maximizes the contribution from both the vorticity and strain.
△ Less
Submitted 21 November, 2014;
originally announced November 2014.
-
Measurements of the Solid-body Rotation of Anisotropic Particles in 3D Turbulence
Authors:
Guy G. Marcus,
Shima Parsa,
Stefan Kramel,
Rui Ni,
Greg A. Voth
Abstract:
We introduce a new method to measure Lagrangian vorticity and the rotational dynamics of anisotropic particles in a turbulent fluid flow. We use 3D printing technology to fabricate crosses (two perpendicular rods) and jacks (three mutually perpendicular rods). Time-resolved measurements of their orientation and solid-body rotation rate are obtained from stereoscopic video images of their motion in…
▽ More
We introduce a new method to measure Lagrangian vorticity and the rotational dynamics of anisotropic particles in a turbulent fluid flow. We use 3D printing technology to fabricate crosses (two perpendicular rods) and jacks (three mutually perpendicular rods). Time-resolved measurements of their orientation and solid-body rotation rate are obtained from stereoscopic video images of their motion in a turbulent flow between oscillating grids with $R_λ$=$91$. The advected particles have a largest dimension of 6 times the Kolmogorov length, making them a good approximation to anisotropic tracer particles. Crosses rotate like disks and jacks rotate like spheres, so these measurements, combined with previous measurements of tracer rods, allow experimental study of ellipsoids across the full range of aspect ratios. The measured mean square tumbling rate, $\langle \dot{p}_i \dot{p}_i \rangle$, confirms previous direct numerical simulations that indicate that disks tumble much more rapidly than rods. Measurements of the alignment of crosses with the direction of the solid-body rotation rate vector provide the first direct observation of the alignment of anisotropic particles by the velocity gradients of the flow.
△ Less
Submitted 27 July, 2014;
originally announced July 2014.
-
Alignment of vorticity and rods with Lagrangian fluid stretching in turbulence
Authors:
Rui Ni,
Nicholas T. Ouellette,
Greg A. Voth
Abstract:
Stretching in continuum mechanics is naturally described using the Cauchy-Green strain tensors. These tensors quantify the Lagrangian stretching experienced by a material element, and provide a powerful way to study processes in turbulent fluid flows that involve stretching such as vortex stretching and alignment of anisotropic particles. Analyzing data from a simulation of isotropic turbulence, w…
▽ More
Stretching in continuum mechanics is naturally described using the Cauchy-Green strain tensors. These tensors quantify the Lagrangian stretching experienced by a material element, and provide a powerful way to study processes in turbulent fluid flows that involve stretching such as vortex stretching and alignment of anisotropic particles. Analyzing data from a simulation of isotropic turbulence, we observe preferential alignment between anisotropic particles and vorticity. We show that this alignment arises because both of these quantities independently tend to align with the strongest Lagrangian stretching direction, as defined by the maximum eigenvector of the left Cauchy-Green strain tensor. In particular, anisotropic particles approach almost perfect alignment with the strongest stretching direction. The alignment of vorticity with stretching is weaker, but still much stronger than previously observed alignment of vorticity with the eigenvectors of the Eulerian strain rate tensor. The alignment of strong vorticity is almost the same as that of rods that have experienced the same stretching.
△ Less
Submitted 4 November, 2013;
originally announced November 2013.
-
Effects of Fluctuating Energy Input on the Small Scales in Turbulence
Authors:
Chen-Chi Chien,
Daniel B. Blum,
Greg A. Voth
Abstract:
In the standard cascade picture of 3D turbulent fluid flows, energy is input at a constant rate at large scales. Energy is then transferred to smaller scales by an intermittent process that has been the focus of a vast literature. However, the energy input at large scales is not constant in most real turbulent flows. We explore the signatures of these fluctuations of large scale energy input on sm…
▽ More
In the standard cascade picture of 3D turbulent fluid flows, energy is input at a constant rate at large scales. Energy is then transferred to smaller scales by an intermittent process that has been the focus of a vast literature. However, the energy input at large scales is not constant in most real turbulent flows. We explore the signatures of these fluctuations of large scale energy input on small scale turbulence statistics. Measurements were made in a flow between oscillating grids, with Re up to 271, in which temporal variations in the large scale energy input can be introduced by modulating the oscillating grid frequency. We find that the Kolmogorov constant for second order longitudinal structure functions depends on the magnitude of the fluctuations in the large scale energy input. We can quantitatively predict the measured change with a model based on Kolmogorov's refined similarity theory. The effects of fluctuations of the energy input can also be observed using structure functions conditioned on the instantaneous large scale velocity. A linear parameterization using the curvature of the iconditional structure functions provides a fairly good match with the measured changes in the Kolmogorov constant. Conditional structure functions are found to provide a more sensitive measure of the presence of fluctuations in the large scale energy input than inertial range scaling coefficients.
△ Less
Submitted 29 June, 2013;
originally announced July 2013.
-
Inertial Range Scaling in Rotations of Long Rods in Turbulence
Authors:
Shima Parsa,
Greg A. Voth
Abstract:
We measure the rotational statistics of neutrally buoyant rods with lengths $2.8 < l/η<72.9 $ in turbulence. For particles with length in the inertial range, we derive a scaling relationship for the mean square rotation rate, $< \dot{p}_i \dot{p}_i > \propto l^{-4/3}$ and show that measurements approach this scaling. Deviations from the proposed scaling are explained as the effect of dissipation r…
▽ More
We measure the rotational statistics of neutrally buoyant rods with lengths $2.8 < l/η<72.9 $ in turbulence. For particles with length in the inertial range, we derive a scaling relationship for the mean square rotation rate, $< \dot{p}_i \dot{p}_i > \propto l^{-4/3}$ and show that measurements approach this scaling. Deviations from the proposed scaling are explained as the effect of dissipation range scales. The correlation time of the Lagrangian autocorrelation of rod rotation rate scales as the turn over time of eddies of the size of the rod. Measuring rotational dynamics of single long rods provides a new way to access the spatial structure of the flow at different length scales.
△ Less
Submitted 19 April, 2013;
originally announced April 2013.
-
Rotation Rate of Rods in Turbulent Fluid Flow
Authors:
Shima Parsa,
Enrico Calzavarini,
Federico Toschi,
Greg A. Voth
Abstract:
The rotational dynamics of anisotropic particles advected in a turbulent fluid flow are important in many industrial and natural setting. Particle rotations are controlled by small scale properties of turbulence that are nearly universal, and so provide a rich system where experiments can be directly compared with theory and simulations. Here we report the first three-dimensional experimental meas…
▽ More
The rotational dynamics of anisotropic particles advected in a turbulent fluid flow are important in many industrial and natural setting. Particle rotations are controlled by small scale properties of turbulence that are nearly universal, and so provide a rich system where experiments can be directly compared with theory and simulations. Here we report the first three-dimensional experimental measurements of the orientation dynamics of rod-like particles as they are advected in a turbulent fluid flow. We also present numerical simulations that show good agreement with the experiments and allow extension to a wide range of particle shapes. Anisotropic tracer particles preferentially sample the flow since their orientations become correlated with the velocity gradient tensor. The rotation rate is heavily influenced by this preferential alignment, and the alignment depends strongly on particle shape.
△ Less
Submitted 1 May, 2012;
originally announced May 2012.
-
Effects of non-universal large scales on conditional structure functions in turbulence
Authors:
Daniel B. Blum,
Surendra Kunwar,
James Johnson,
Greg A. Voth
Abstract:
We report measurements of conditional Eulerian and Lagrangian structure functions in order to assess the effects of non-universal properties of the large scales on the small scales in turbulence. We study a 1m $\times$ 1m $\times$ 1.5m flow between oscillating grids which produces $R_λ=285$ while containing regions of nearly homogeneous and highly inhomogeneous turbulence. Large data sets of thr…
▽ More
We report measurements of conditional Eulerian and Lagrangian structure functions in order to assess the effects of non-universal properties of the large scales on the small scales in turbulence. We study a 1m $\times$ 1m $\times$ 1.5m flow between oscillating grids which produces $R_λ=285$ while containing regions of nearly homogeneous and highly inhomogeneous turbulence. Large data sets of three-dimensional tracer particle velocities have been collected using stereoscopic high speed cameras with real-time image compression technology. Eulerian and Lagrangian structure functions are measured in both homogeneous and inhomogeneous regions of the flow. We condition the structure functions on the instantaneous large scale velocity or on the grid phase. At all scales, the structure functions depend strongly on the large scale velocity, but are independent of the grid phase. We see clear signatures of inhomogeneity near the oscillating grids, but even in the homogeneous region in the center we see a surprisingly strong dependence on the large scale velocity that remains at all scales. Previous work has shown that similar correlations extend to very high Reynolds numbers. Comprehensive measurements of these effects in a laboratory flow provide a powerful tool for assessing the effects of shear, inhomogeneity and intermittency of the large scales on the small scales in turbulence.
△ Less
Submitted 1 December, 2009; v1 submitted 5 August, 2009;
originally announced August 2009.
-
Measurements of Accelerations of Large Neutrally-buoyant Particles in Intense Turbulence
Authors:
Rachel D. Brown,
Zellman Warhaft,
Greg A. Voth
Abstract:
We measure acceleration statistics of neutrally buoyant spherical particles with diameter 0.4 < d/eta <27 in intense turbulence (400< R_lambda <815). High speed cameras image polystyrene tracer particles in a flow between counter-rotating disks. The measurements of acceleration variance, <a^2>, clearly resolve the transition from the tracer like behavior of small particles to the much smaller ac…
▽ More
We measure acceleration statistics of neutrally buoyant spherical particles with diameter 0.4 < d/eta <27 in intense turbulence (400< R_lambda <815). High speed cameras image polystyrene tracer particles in a flow between counter-rotating disks. The measurements of acceleration variance, <a^2>, clearly resolve the transition from the tracer like behavior of small particles to the much smaller accelerations of large particles. For d>5 eta, <a^2> decreases with diameter as d^{-2/3} in agreement with inertial range scaling arguments. A model relating <a^2> to the pressure structure functions matches the transition from small to large particle behavior if the particles respond to pressure differences over (1.7 +- 0.3) d. A model relating <a^2> to the fluid acceleration averaged over the particle diameter predicts the transition with no free parameters, but does not show clean inertial range scaling in the size range studied. Consistent with earlier work, we find that the scaled acceleration probability density function shows very little dependence on particle size.
△ Less
Submitted 19 May, 2009;
originally announced May 2009.
-
Measurements of the steady streaming flow around oscillating spheres using 3D particle tracking velocimetry
Authors:
Florian Otto,
Emmalee K. Riegler,
Greg A. Voth
Abstract:
Granular particles vibrated in a fluid have been found to exhibit self-organization with attractive and repulsive interactions between the particles. These interactions have been attributed to the steady streaming flow around oscillating particles. Here we examine the steady streaming flow surrounding a vertically oscillating sphere using three dimensional particle tracking velocimetry. We prese…
▽ More
Granular particles vibrated in a fluid have been found to exhibit self-organization with attractive and repulsive interactions between the particles. These interactions have been attributed to the steady streaming flow around oscillating particles. Here we examine the steady streaming flow surrounding a vertically oscillating sphere using three dimensional particle tracking velocimetry. We present measurements of the flow with the sphere far from boundaries, close to the bottom wall of the tank, and near another oscillating sphere. The steady velocity field is found to disagree with available analytic calculations. When the sphere is oscillated near the bottom wall the entire topology of the flow changes, resulting in a larger repulsive region than expected. Previous experiments saw attraction between particles in the region where the flow around a single particle is repulsive. We conclude that advection in the streaming flow due to a single particle cannot explain the observed attractive and repulsive interactions, rather non-linear interactions between the flows around two or more spheres must be responsible.
△ Less
Submitted 26 September, 2007;
originally announced September 2007.
-
Ordered clusters and dynamical states of particles in a vibrated fluid
Authors:
Greg A. Voth,
B. Bigger,
M. R. Buckley,
W. Losert,
M. P. Brenner,
H. A. Stone,
J. P. Gollub
Abstract:
Fluid-mediated interactions between particles in a vibrating fluid lead to both long range attraction and short range repulsion. The resulting patterns include hexagonally ordered micro-crystallites, time-periodic structures, and chaotic fluctuating patterns with complex dynamics. A model based on streaming flow gives a good quantitative account of the attractive part of the interaction.
Fluid-mediated interactions between particles in a vibrating fluid lead to both long range attraction and short range repulsion. The resulting patterns include hexagonally ordered micro-crystallites, time-periodic structures, and chaotic fluctuating patterns with complex dynamics. A model based on streaming flow gives a good quantitative account of the attractive part of the interaction.
△ Less
Submitted 7 February, 2002;
originally announced February 2002.
-
Measurement of Particle Accelerations in Fully Developed Turbulence
Authors:
Greg A. Voth,
A. La Porta,
Alice M. Crawford,
Eberhard Bodenschatz,
Jim Alexander
Abstract:
We use silicon strip detectors (originally developed for the CLEO III high energy particle physics experiment) to measure fluid particle trajectories in turbulence with temporal resolution of up to 70,000 frames per second. This high frame rate allows the Kolmogorov time scale of a turbulent water flow to be fully resolved for 140 <= R_lambda <= 970. Particle trajectories exhibiting acceleration…
▽ More
We use silicon strip detectors (originally developed for the CLEO III high energy particle physics experiment) to measure fluid particle trajectories in turbulence with temporal resolution of up to 70,000 frames per second. This high frame rate allows the Kolmogorov time scale of a turbulent water flow to be fully resolved for 140 <= R_lambda <= 970. Particle trajectories exhibiting accelerations up to 16,000 m\s^2 (40 times the rms value) are routinely observed. The probability density function of the acceleration is found to have Reynolds number dependent stretched exponential tails. The moments of the acceleration distribution are calculated. The scaling of the acceleration component variance with the energy dissipation is found to be consistent with the results for low Reynolds number direct numerical simulations, and with the K41 based Heisenberg-Yaglom prediction for R_lambda >= 500. The acceleration flatness is found to increase with Reynolds number, and to exceed 60 at R_lambda = 970. The coupling of the acceleration to the large scale anisotropy is found to be large at low Reynolds number and to decrease as the Reynolds number increases, but to persist at all Reynolds numbers measured. The dependence of the acceleration variance on the size and density of the tracer particles is measured. The autocorrelation function of an acceleration component is measured, and is found to scale with the Kolmogorov time tau_eta.
△ Less
Submitted 3 December, 2001; v1 submitted 8 October, 2001;
originally announced October 2001.
-
Fluid Particle Accelerations in Fully Developed Turbulence
Authors:
A. La Porta,
Greg A. Voth,
Alice M. Crawford,
Jim Alexander,
Eberhard Bodenschatz
Abstract:
The motion of fluid particles as they are pushed along erratic trajectories by fluctuating pressure gradients is fundamental to transport and mixing in turbulence. It is essential in cloud formation and atmospheric transport, processes in stirred chemical reactors and combustion systems, and in the industrial production of nanoparticles. The perspective of particle trajectories has been used suc…
▽ More
The motion of fluid particles as they are pushed along erratic trajectories by fluctuating pressure gradients is fundamental to transport and mixing in turbulence. It is essential in cloud formation and atmospheric transport, processes in stirred chemical reactors and combustion systems, and in the industrial production of nanoparticles. The perspective of particle trajectories has been used successfully to describe mixing and transport in turbulence, but issues of fundamental importance remain unresolved. One such issue is the Heisenberg-Yaglom prediction of fluid particle accelerations, based on the 1941 scaling theory of Kolmogorov (K41). Here we report acceleration measurements using a detector adapted from high-energy physics to track particles in a laboratory water flow at Reynolds numbers up to 63,000. We find that universal K41 scaling of the acceleration variance is attained at high Reynolds numbers. Our data show strong intermittency---particles are observed with accelerations of up to 1,500 times the acceleration of gravity (40 times the root mean square value). Finally, we find that accelerations manifest the anisotropy of the large scale flow at all Reynolds numbers studied.
△ Less
Submitted 8 November, 2000;
originally announced November 2000.