The Impact of Hydration Shell Inclusion and Chain Exclusion in the Efficacy of Reaction Coordinates for Homogeneous and Heterogeneous Ice Nucleation

Abstract

Ice nucleation plays a pivotal role in many natural and industrial processes, and molecular simulations play have proven vital in uncovering its kinetics and mechanisms. A fundamental component of such simulations is the choice of an order parameter (OP) that quantifies the progress of nucleation, with the efficacy of an OP typically measured by its ability to predict the committor probabilities. Here, we leverage a machine learning framework introduced in our earlier work (Domingues,~\emph{et al.}, \emph{J. Phys. Chem. Lett.}, 15, 1279, {\bf 2024}) to systematically investigate how key implementation details influence the efficacy of standard Steinhardt OPs in capturing the progress of both homogeneous and heterogeneous ice nucleation. Our analysis identify distance and $q_6$ cutoffs, as the primary determinants of OP performance, regardless of the mode of nucleation. We also examine the impact of two popular refinement strategies, namely chain exclusion and hydration shell inclusion, on OP efficacy. We find neither strategy to exhibit a universally consistent impact. Instead, their efficacy depends strongly on the chosen distance and $q_6$ cutoffs. Chain exclusion enhances OP efficacy when the underlying OP lacks sufficient selectivity, whereas hydration shell inclusion is beneficial for overly selective OPs. Consequently, we demonstrate that selecting optimal combinations of such cutoffs can eliminate the need for these refinement strategies altogether. These findings provide a systematic understanding of how to design and optimize OPs for accurately describing complex nucleation phenomena, offering valuable guidance for improving the predictive power of molecular simulations.

11 pages, 9 figures

References in corpus (12)

• Accurate determination of crystal structures based on averaged local bond order parameters
• Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations
• Rate of Homogeneous Crystal Nucleation in molten NaCl
• Ice Nucleation on Carbon Surface Supports the Classical Theory for Heterogeneous Nucleation
• Local order parameters for use in driving homogeneous ice nucleation with all-atom models of water
• Microscopic Mechanism and Kinetics of Ice Formation at Complex Interfaces: Zooming in on Kaolinite
• Studying Rare Events using Forward-Flux Sampling: Recent Breakthroughs and Future Outlook
• Suppression of Sub-surface Freezing in Free-Standing Thin Films of a Coarse-grained Model of Water
• How to Quantify and Avoid Finite Size Effects in Computational Studies of Crystal Nucleation: The Case of Heterogeneous Ice Nucleation
• Role of Nanoscale Interfacial Proximity in Contact Freezing in Water
• How to quantify and avoid finite size effects in computational studies of crystal nucleation: The case of homogeneous crystal nucleation
• LeaPP: Learning Pathways to Polymorphs through machine learning analysis of atomic trajectories

Cited by in corpus (2)

• Robustness of classical nucleation theory to chemical heterogeneity of crystal nucleating substrates
• Quantum Markov chain Monte Carlo method with programmable quantum simulators