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.

Molecular Dynamics Study on Self-Assembly of C36 Clusters

C72 (or C120) cluster growth from C36 (or C60) clusters

Isomers Of C-36 And Free Energy Criteria For Cluster Growth

Growth of C-30 and C-31 Clusters: Structures, Energetics and Dynamics

Theoretical investigation on isomer formation probability and free energy of small C clusters

Conditions for the self-assembling of cluster materials

Simulated Dynamics and Endohedral Relesase of Ne from Ne@C60 Clusters

Thermal Stability of Ag13- Clusters Studied by Ab Initio Molecular Dynamics Simulations

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

THEORY–EXPERIMENT RELATIONSHIP OF THE "SHRINKING HOT GIANT" ROAD OF DYNAMIC FULLERENE SELF-ASSEMBLY IN HOT CARBON VAPOR

Studies of the Formation of Carbon Clusters

Isomer Abundance of Small Carbon Clusters Formed in Buffer He Gas

Molecular dynamics simulations of lead clusters

Coalescence of nanoscale metal clusters: Molecular-dynamics study

Growth Patterns of Carbon Clusters Cn (n = 2-60) Identified Via ABCluster Searching and DFT Benchmarking.

Structure and dynamics of carbon dioxide clusters: A molecular dynamics study

Molecular Dynamics Simulation Of Phase Transition, Nucleation, Recrystalization Of Kf Clusters

Growth Patterns of Carbon Clusters C n (n = 2-60) Identified via ABCluster Searching and DFT Benchmarking

Isomer Spectrum of C-30 Cluster and the Growth Dynamics

Molecular Dynamics Simulation of Nucleation and Recrystalization of KCl Clusters

Investigating limiting factors for self-assembly of carbon nanotubes: A molecular dynamics simulation study