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.

Relative Distribution of Water Clusters at Temperature (300-3000K) and Pressure (1-500MPa)

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

The Affection to Water Clusters on Temperatures and Electric Fields by Molecular Dynamics Simulation

Ab Initio Investigation of the Lower Energy Candidate Structures for (h2o)(5)+ Water Cluster

Ab Initio Investigation of Possible Candidate Structures and Properties of Water Cluster (H2O)7+ Via Particle Swarm Optimization Method

Water Clusters and Density Fluctuations in Liquid Water Based on Extended Hierarchical Clustering Methods

Dynamics of Water Clusters Confined in Ionic Liquid at an Elevated Pressure

High Pressure-Temperature Brillouin Study of Liquid Water: Evidence of the Structural Transition from Low-Density Water to High-Density Water

Pressure Dependence Of Hydrogen Bonding In A Novel H2o-H-2 Clathrate

High Temperatures and High Pressures Brillouin Scattering Studies of Liquid H2O+CO2 Mixtures

A Theoretical Study on the Formation, Binding Energy and Monomer Dipole Moment of Small Water Cluster Systems

Hydrogen bonds and dynamics of liquid water and alcohols

Water Cluster Distribution under Ambient Conditions and Dissociation of H+(H2O)(n) (n=4 similar to 16) Using Mass Spectrometry in Atmosphere

A computational study of H-bonded networks in cyclic water clusters, (H2O)n (n = 3–12)

Molecular Dynamics Simulation Studies of Cryoprotective Agent Solutions: the Relation Between Melting Temperature and the Ratio of Hydrogen Bonding Acceptor to Donor Number

Water Clusters: Contributions of Binding Energy and Entropy to Stability

Dynamics and structural changes of small water clusters on ionization

Infrared Spectroscopy of Neutral Water Clusters at Finite Temperature: Evidence for a Noncyclic Pentamer.

Structural and Energetic Stability of the Lowest Equilibrium Structures of Water Clusters

Accurate Relative Energies and Binding Energies of Large Ice-Liquid Water Clusters and Periodic Structures

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