Abstract:Due to their low density in atmosphere, theoretical simulations of the nucleation of gas-phase molecules are computationally very expensive. In this study, neural network potential (NNP) is combined with enhanced sampling techniques to effectively investigate the nucleation of water clusters in gas phase. The neural network potential is trained based on water-cluster energies and forces from density functional theory (DFT). The problem that the binding between water molecules is too weak in the previous empirical force field model has been solved in the NNP. This NNP potential is then applied to Monte Carlo simulations in grand canonical ensemble with enhanced sampling methods such as aggregation-volume-bias Monte Carlo (AVBMC) and transition-matrix Monte Carlo (TMMC) to realize a random walk among different cluster sizes. Probability distribution of water cluster sizes and the corresponding Gibbs free energies can then be obtained. Subsequently, the evaporation rates of water clusters can be calculated via umbrella sampling Monte Carlo simulations in canonical ensemble combined with variational transition state theory (VTST). We observe a big change of free energy and evaporation rate from tetramer to pentamer. A statistical analysis of the number of hydrogen bonds suggests that more hydrogen bonds are required to be broken in the evaporation reaction of tetramer compared to that of trimer and pentamer. Structure analysis indicates that, although the ground state of the pentamer has a two-dimensional ring structure, three-dimensional hydrogen bond network begins to form in pentamer at finite temperature. Therefore, it is a two-dimensional to three-dimensional transition from tetramer to pentamer. The fact that the most probable configuration of pentamer is different from the lowest energy configuration demonstrates the importance of molecular simulations. Simply finding the lowest energy configuration via global geometry optimization and then calculating the free energy within a harmonic approximation of vibrations are not a universal protocol for cluster systems. Methods used in this study are expected to be applicable for more complicated multicomponent systems, which opens an avenue for the research of particulate matter formation in atmosphere.

An Explicit Cluster Model for Binary Nuclei in Water–alcohol Systems

On the Cluster Compositions in the Classical Binary Nucleation Theory

Ab Initio Investigation Of Structure, Stability, Thermal Behavior, Bonding, And Infrared Spectra Of Ionized Water Cluster (H2o)(6)(+) (Vol 145, 154307, 2016)

Clustering and Micro-immiscibility in Alcohol-Water Mixtures: Evidence from Molecular Dynamics Simulations

Investigation of liquor microstructure (ethanol-water clusters): Molecular dynamics simulation and density functional theory

A simple model for density anomaly and excess thermodynamic properties of alcohol-like in water-like mixture

The Self-Consistency Correction to Homogeneous Nucleation: Extension to Binary Systems

Measurement of the Molecular Content of Binary Nuclei. II. Use of the Nucleation Rate Surface for Water–ethanol

Temperature of maximum density and excess properties of short-chain alcohol aqueous solutions: A simplified model simulation study

Site-site interaction model for alcohol models in two-dimensions

Mathematical models of water nucleation process for the Direct Simulation Monte Carlo method

Binary Homogeneous Nucleation in Water-Succinic Acid and Water-Glutaric Acid Systems

Strange Predictions by Binary Heterogeneous Nucleation Theory Compared with A Quantitative Experiment

Ternary nucleation: kinetics and application to water–ammonia–hydrochloric acid system

Assessment of Anharmonicities in Clusters: Developing and Validating a Minimum-Information Partition Function

Density functional theory (DFT) study of water autoionization in solvated clusters

Nucleation of Water Clusters in Gas Phase: A Computational Study Based on Neural Network Potential and Enhanced Sampling

Dynamics and structural changes of small water clusters on ionization

A New Molecular-Thermodynamic Model Based on Lattice Fluid Theory: Application to Pure Fluids and Their Mixtures

Core-softened water-alcohol mixtures: the solute-size effects

Two state model for the ML-BOP potential