.png){kind=icon}
.png){kind=icon}
Return
From Clusters to Bulk: Searching for the Stability Crossover between Ih and Ic via AVBMC Simulations
.png){kind=icon}
OA
AI
B
DOI：10.1021/acs.jpcc.6c03068.png)
.png){kind=icon}
Original request for help.png){kind=icon}
Share.png){kind=icon}
Save
Abstract
En 中文
Among the more than 20 known crystalline phases of ice, hexagonal (Ih) and cubic (Ic) ice are the dominant polymorphs at ambient pressure. Their relative thermodynamic stability and competition have been the subject of extensive investigation. This work presents a systematic study of the free energy landscape spanning from small clusters to the bulk for both polymorphs using a lattice-based aggregation-volume-bias Monte Carlo (LB-AVBMC) approach. Simulations were performed with the TIP4P/2005 water model at 300 K. Although this temperature exceeds the model’s melting point, the lattice constraint stabilizes crystalline clusters even at small sizes. It was found that the relative stability of Ih and Ic exhibits nontrivial size dependence, with their free energies interweaving in the small-cluster regime. Nucleation free energies computed for large clusters were analyzed within the Tolman-corrected classical nucleation theory framework to extrapolate bulk chemical potentials. While Ih is more stable than Ic, both crystalline phases remain thermodynamically metastable relative to the liquid phase at 300 K. By combining bulk chemical potentials with enthalpy data obtained from additional NpT simulations, the melting temperature is estimated to be around 250 K, in good agreement with previous studies. The enthalpy difference between Ih and Ic is approximately 20 J/mol at 300 K and decreases with decreasing temperature.
Keywords:
Cluster chemistry
Liquids
Nucleation
Surface tension
Thermodynamic properties
Journal
T
IF:
3.2
Papers:
1.2K
Citations:
4
Organization
Cited Papers