First-Principles Molecular Dynamics Simulations of the Spontaneous Freezing Transition of 2D Water in a Nanoslit
Jian Jiang,Yurui Gao,Weiduo Zhu,Yuan Liu,Chongqin Zhu,Joseph S. Francisco,Xiao Cheng Zeng
DOI: https://doi.org/10.1021/jacs.1c03243
IF: 15
2021-05-19
Journal of the American Chemical Society
Abstract:As with bulk ices, two-dimensional (2D) ices exhibit diverse crystalline structures, and the majority of these 2D structures have been predicted based on classical molecular dynamics (MD) simulations. Here, the spontaneous freezing transition of 2D liquid water within hydrophobic nanoslits is demonstrated for the first time using first-principles MD simulations. Various 2D ices are observed under different lateral pressure and temperature conditions. Notably, the liquid water confined to a 6.0 Å-wide nanoslit can spontaneously freeze into a monolayer ice consisting of an array of zigzag water chains at 2.5 GPa and 250 K. Moreover, within an 8.0 Å-wide nanoslit and at 4.0 GPa and 300 K, a previously unreported bilayer ice forms spontaneously that has a structure resembling that of the double surface layers of bulk ice-VII. Both 2D crystalline ices do not obey the ice rule, suggesting first-principles simulation can access a certain phase space that is not easily approached using classical simulations.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/jacs.1c03243?goto=supporting-info">https://pubs.acs.org/doi/10.1021/jacs.1c03243</a>.Area density of water phases from BOMD simulations at T= 250 K or T = 300 K; time-dependent potential energy of the ZZMI formed at nanoslit width w = 6 Å and lateral pressure <i>P</i> = 2.5 GPa; snapshot of zigzag ice chain monolayer ice (ZZMI) obtained from BOMD simulation at nanoslit width of 6 Å, temperature of 250 K, and lateral pressure of 2.5 GPa; computed mean-squared displacements at different lateral pressures and at temperature of 250 or 300 K and at different nanoslit width w (6–8 Å; top and side views of BL-iVII ice structure (Figure S4); oxygen–oxygen radial distribution function, hydrogen bond O···H distance distributions, and O–H···O angle distributions of BL-iVII in BOMD simulations; oxygen–oxygen radial distribution function, hydrogen bond O···H distance distributions, and O–H···O angle distributions of ZZMI in BOMD simulations with BLYP-D3 and vdW-DF2 functionals; top view and side view of fRMI and ZZMI 2D ice phases; fitted Lennard-Jones 9-3 potential of water-graphene interaction (Figure S8) (<a class="ext-link" href="/doi/suppl/10.1021/jacs.1c03243/suppl_file/ja1c03243_si_001.pdf">PDF</a>)Top view of monolayer ice ZZMI formation through BOMD at 250 K, 2.5 GPa and 6 Å (<a class="ext-link" href="/doi/suppl/10.1021/jacs.1c03243/suppl_file/ja1c03243_si_002.mp4">MP4</a>)Top view of bilayer ice BL-iVII formation through BOMD at 300 K, 4.0 GPa and 8 Å (<a class="ext-link" href="/doi/suppl/10.1021/jacs.1c03243/suppl_file/ja1c03243_si_003.mp4">MP4</a>)This article has not yet been cited by other publications.
chemistry, multidisciplinary