Display options
Share it on

Phys Chem Chem Phys. 2016 Feb 07;18(5):4102-11. doi: 10.1039/c5cp06535f.

Computational investigation of structure, dynamics and nucleation kinetics of a family of modified Stillinger-Weber model fluids in bulk and free-standing thin films.

Physical chemistry chemical physics : PCCP

Melisa M Gianetti, Amir Haji-Akbari, M Paula Longinotti, Pablo G Debenedetti

Affiliations

  1. DQIAQF/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Buenos Aires, Argentina.
  2. Department of Chemical and Biological Engineering, Princeton University, Princeton NJ 08544, USA. [email protected].

PMID: 26778494 DOI: 10.1039/c5cp06535f

Abstract

In recent years, computer simulations have found increasingly widespread use as powerful tools for studying phase transitions in wide variety of systems. In the particular and very important case of aqueous systems, the commonly used force-fields tend to offer quite different predictions with respect to a wide range of thermodynamic and kinetic properties, including the ease of ice nucleation, the propensity to freeze at a vapor-liquid interface, and the existence of a liquid-liquid phase transition. It is thus of fundamental and practical interest to understand how different features of a given water model affect its thermodynamic and kinetic properties. In this work, we use the forward-flux sampling technique to study the crystallization kinetics of a family of modified Stillinger-Weber (SW) potentials with energy (ε) and length (σ) scales taken from the monoatomic water (mW) model, but with different tetrahedrality parameters (λ). By increasing λ from 21 to 24, we observe the nucleation rate increases by 48 orders of magnitude at a supercooling of ζ = T/Tm = 0.845. Using classical nucleation theory, we are able to demonstrate that this change can largely be accounted for by the increase in |Δμ|, the thermodynamic driving force. We also perform rate calculations in freestanding thin films of the supercooled liquid, and observe a crossover from surface-enhanced crystallization at λ = 21 to bulk-dominated crystallization for λ ≥ 22.

Publication Types