Test of classical nucleation theory on deeply supercooled high-pressure simulated silica

Abstract

We test classical nucleation theory (CNT) in the case of simulations of deeply supercooled, high density liquid silica, as modelled by the BKS potential. We find that at density $\rho =4.38$~g/cm${}^3$, spontaneous nucleation of crystalline stishovite occurs in conventional molecular dynamics simulations at temperature T=3000 K, and we evaluate the nucleation rate J directly at this T via "brute force" sampling of nucleation events. We then use parallel, constrained Monte Carlo simulations to evaluate $\Delta G(n)$, the free energy to form a crystalline embryo containing n silicon atoms, at T=3000, 3100, 3200 and 3300 K. We find that the prediction of CNT for the n-dependence of $\Delta G(n)$ fits reasonably well to the data at all T studied, and at 3300 K yields a chemical potential difference between liquid and stishovite that matches independent calculation. We find that $n^{\ast }$, the size of the critical nucleus, is approximately 10 silicon atoms at T=3300 K. At 3000 K, $n^{\ast }$ decreases to approximately 3, and at such small sizes methodological challenges arise in the evaluation of $\Delta G(n)$ when using standard techniques; indeed even the thermodynamic stability of the supercooled liquid comes into question under these conditions. We therefore present a modified approach that permits an estimation of $\Delta G(n)$ at 3000 K. Finally, we directly evaluate at T=3000 K the kinetic prefactors in the CNT expression for J, and find physically reasonable values; e.g. the diffusion length that Si atoms must travel in order to move from the liquid to the crystal embryo is approximately 0.2 nm. We are thereby able to compare the results for J at 3000 K obtained both directly and based on CNT, and find that they agree within an order of magnitude.

arXiv commentcorrected calculation, new figure, accepted in JCP

References in corpus (8)

• Static and Dynamic Properties of a Viscous Silica Melt Molecular Dynamics Computer Simulations
• Fragile-to-strong transition and polyamorphism in the energy landscape of liquid silica
• Interplay Between Time-Temperature-Transformation and the Liquid-Liquid Phase Transition in Water
• Computer simulations of liquid silica: equation of state and liquid-liquid phase transition
• Rate of Homogeneous Crystal Nucleation in molten NaCl
• Free energy and configurational entropy of liquid silica: fragile-to-strong crossover and polyamorphism
• Phase diagram of silica from computer simulation
• The Origin of the Fragile-to-Strong Crossover in Liquid Silica as Expressed by its Potential Energy Landscape

Cited by in corpus (12)

• The role of local structure in dynamical arrest
• Phase diagram of a two-dimensional system with anomalous liquid properties
• Crystal Nucleation in a Supercooled Liquid with Glassy Dynamics
• Nucleation barriers in tetrahedral liquids spanning glassy and crystallizing regimes
• A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water
• The barrier to ice nucleation in monatomic water
• Heterogeneous Nucleation in the Low Barrier Regime
• Crystallization of Lennard-Jones nanodroplets: from near melting to deeply supercooled
• Thermodynamics and kinetics of crystallisation in deeply supercooled Stillinger-Weber silicon
• Are nucleation bubbles in a liquid all independent?
• Competitive Heterogeneous Nucleation onto a Microscopic Impurity in a Potts Model
• A general relation between the largest nucleus and all nuclei distributions for free energy calculations