Can molecular simulations reliably compare homogeneous and heterogeneous ice nucleation?
arXiv:2201.06338 · doi:10.1063/5.0085750
Abstract
In principle, the answer to the posed titular question is undoubtedly 'yes.' But in practice, requisite reference data for homogeneous systems have been obtained with a treatment of intermolecular interactions that is different from that typically employed for heterogeneous systems. In this article, we assess the impact of the choice of truncation scheme when comparing water in homogeneous and inhomogeneous environments. Specifically, we use explicit free energy calculations and a simple mean field analysis to demonstrate that using the 'cut-and-shift' version of the Lennard-Jones potential (common to most simple point charge models of water) results in a systematic increase in the melting temperature of ice I. In addition, by drawing an analogy between a change in cutoff and a change in pressure, we use existing literature data for homogeneous ice nucleation at negative pressures to suggest that enhancements due to heterogeneous nucleation may have been overestimated by several orders of magnitude.
Main article: 9 pages, 7 figures. Supplementary Material: 10 pages, 4 figures
References in corpus (12)
- Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations
- Determination of phase diagrams via computer simulation: Methodology and applications to water, electrolytes and proteins
- Application of Ewald summations to long-range dispersion forces
- Homogeneous ice nucleation evaluated for several water models
- Microscopic Mechanism and Kinetics of Ice Formation at Complex Interfaces: Zooming in on Kaolinite
- Ice Formation on Kaolinite: Insights from Molecular Dynamics Simulations
- The Microscopic Features of Heterogeneous Ice Nucleation May Affect the Macroscopic Morphology of Atmospheric Ice Crystals
- Quantum-mechanical exploration of the phase diagram of water
- Computing the free energy of molecular solids by the Einstein molecule approach: Ices XIII and XIV, hard-dumbbells and a patchy model of proteins
- Stabilization of AgI's polar surfaces by the aqueous environment, and its implications for ice formation
- Finite field formalism for bulk electrolyte solutions
- Communication: Truncated non-bonded potentials can yield unphysical behavior in molecular dynamics simulations of interfaces