Reaction coordinates and rate constants for liquid droplet nucleation: quantifying the interplay between driving force and memory

Abstract

In this work we revisit the classic problem of homogeneous nucleation of a liquid droplet in a supersaturated vapor phase. We consider this at different extents of the driving force, which here is the extent of supersaturation, and calculate a reaction coordinate (RC) for nucleation as the driving force is varied. The RC is constructed as a linear combination of three order parameters, where one accounts for the number of liquid-like atoms, and the other two for local density fluctuations. The RC is calculated from all-atom biased and unbiased molecular dynamics (MD) simulations using the spectral gap optimization approach "SGOOP" [P. Tiwary and B. J. Berne, Proc. Natl. Acad. Sci. U. S. A. 113, 2839 (2016)]. Our key finding is that as the supersaturation decreases, the RC ceases to simply be the number of liquid-like atoms, and instead it becomes important to explicitly consider local density fluctuations that correlate with shape and density variations in the nucleus. All three order parameters are found to have similar barriers in their respective potentials of mean force, however, as the supersaturation decreases the density fluctuations decorrelate slower and thus carry longer memory. Thus at lower supersaturations density fluctuations are non-Markovian and can not be simply ignored from the RC by virtue of being noise. Finally, we use this optimized RC to calculate nucleation rates in the infrequent metadynamics framework, and show it leads to more accurate estimate of the nucleation rate with four orders of magnitude acceleration relative to unbiased MD.

arXiv comment10 pages, 5 figures

References in corpus (18)

• Canonical sampling through velocity-rescaling
• Escaping free-energy minima
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations
• A Novel Path Sampling Method for the Calculation of Rate Constants
• From Metadynamics to Dynamics
• Forward Flux Sampling for rare event simulations
• Simulating rare events in equilibrium or non-equilibrium stochastic systems
• Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
• A self-learning algorithm for biased molecular dynamics
• Investigating Rare Events by Transition Interface Sampling
• Homogeneous and heterogeneous nucleation of Lennard-Jones liquids
• Maximum Caliber: a general variational principle for dynamical systems
• Overcoming timescale and finite-size limitations to compute nucleation rates from small scale Well Tempered Metadynamics simulations
• Predicting reaction coordinates in energy landscapes with diffusion anisotropy
• How wet should be the reaction coordinate for ligand unbinding?
• Contour forward flux sampling: Sampling rare events along multiple collective variables
• Non-Markovian dynamics of clusters during nucleation

Cited by in corpus (8)

• Enhanced Sampling in the Age of Machine Learning: Algorithms and Applications
• Driving and characterizing nucleation of urea and glycine polymorphs in water
• Recent advances in describing and driving crystal nucleation using machine learning and artificial intelligence
• Manifold Learning in Atomistic Simulations: A Conceptual Review
• Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?
• Nucleation rates from small scale atomistic simulations and transition state theory
• Atomic scale insights into NaCl nucleation in nanoconfined environments
• Extending and validating bubble nucleation rate predictions in a Lennard-Jones fluid with enhanced sampling methods and transition state theory