Tuning hydrogen bond network connectivity in the electric double layer with cations.

Hydrogen bond (H-bond) network connectivity in electric double layers (EDLs) is of paramount importance for interfacial HER/HOR electrocatalytic processes. However, it remains unclear whether the cation-specific effect on H-bond network connectivity in EDLs exists. Herein, we report simulation evidence from ab initio molecular dynamics that cations at Pt(111)/water interfaces can tune the structure and the connectivity of H-bond networks in EDLs. As the surface charge density σ becomes more negative, we show that the connectivity of the H-bond networks in EDLs of the Na+ and Ca2+ systems decreases markedly; in stark contrast, the connectivity of the H-bond networks in EDLs of the Mg2+ system increases slightly. Further analysis revealed that the interplay between the hydration of cations and the interfacial water structure plays a key role in the connectivity of H-bond networks in EDLs. These findings highlight the key roles of cations in EDLs and electrocatalysis.

The performance of OPC and OPC3 water models in predictions of 2D structures under nanoconfinement.

Nanoconfined water plays an important role in broad fields of science and engineering. Classical molecular dynamics (MD) simulations have been widely used to investigate water phases under nanoconfinement. The key ingredient of MD is the force field. In this study, we systematically investigated the performance of a recently introduced family of globally optimal water models, OPC and OPC3, and TIP4P/2005 in describing nanoconfined two-dimensional (2D) water ice. Our studies show that the melting points of the monolayer square ice (MSI) of all three water models are higher than the melting points of the corresponding bulk ice Ih. Under the same conditions, the melting points of MSI of OPC and TIP4P/2005 are the same and are ∼90 K lower than that of the OPC3 water model. In addition, we show that OPC and TIP4P/2005 water models are able to form a bilayer AA-stacked structure and a trilayer AAA-stacked structure, which are not the cases for the OPC3 model. Considering the available experimental data and first-principles simulations, we consider the OPC water model as a potential water model for 2D water ice MD studies.

Puckered Zigzag Monolayer Ice: Does a Confined Flat Four-Coordinated Monolayer Ice Always Have a Corresponding Puckered Phase?

We note that a flat, four-coordinated monolayer ice under confinement always has a corresponding puckered phase. Recently, a monolayer ice consisting of an array of zigzag water chains (ZZMI) predicted by first-principles calculations of water under confinement is a flat four-coordinated monolayer ice. Herein, to investigate whether puckered ZZMI exists stably, we perform molecular dynamics simulations of two-dimensional (2D) ice formation for water constrained in graphene nanocapillaries. We find a novel monolayer ice structure that can be viewed as the ZZMI puckered along the direction perpendicular to the zigzag chain (pZZMI). Unlike ZZMI that does not satisfy the ice rule, each water molecule in pZZMI can form four hydrogen bonds (HBs) via forming two stable intersublayer HBs and two intrasublayer HBs. This work provides a fresh perspective on 2D confined ice, highlighting the intrinsic connections between 2D confined ices.