Enhanced Configurational Entropy in High-Density Nanoconfined Bilayer Ice

Close-Packed Ices in Nanopores

Water molecules in any of the ice polymorphs organize themselves into a perfect four-coordinated hydrogen-bond network at the expense of dense packing. Even at high pressures, there seems to be no way to reconcile the ice rules with the close packing. Here, we report several close-packed ice phases in carbon nanotubes obtained from molecular dynamics simulations of two different water models. Typically they are in plastic states at high temperatures and are transformed into the hydrogen-ordered ice, keeping their close-packed structures at lower temperatures. The close-packed structures of water molecules in carbon nanotubes are identified with those of spheres in a cylinder. We present design principles of hydrogen-ordered, close-packed structures of ice in nanotubes, which suggest many possible dense ice forms with or without nonzero polarization. In fact, some of the simulated ices are found to exhibit ferroelectric ordering upon cooling.

Temperature-pressure phase diagram of confined monolayer water/ice at first-principles accuracy with a machine-learning force field

Understanding the phase behaviour of nanoconfined water films is of fundamental importance in broad fields of science and engineering. However, the phase behaviour of the thinnest water film - monolayer water - is still incompletely known. Here, we developed a machine-learning force field (MLFF) at first-principles accuracy to determine the phase diagram of monolayer water/ice in nanoconfinement with hydrophobic walls. We observed the spontaneous formation of two previously unreported high-density ices, namely, zigzag quasi-bilayer ice (ZZ-qBI) and branched-zigzag quasi-bilayer ice (bZZ-qBI). Unlike conventional bilayer ices, few inter-layer hydrogen bonds were observed in both quasi-bilayer ices. Notably, the bZZ-qBI entails a unique hydrogen-bonding network that consists of two distinctive types of hydrogen bonds. Moreover, we identified, for the first time, the stable region for the lowest-density [Formula: see text] monolayer ice (LD-48MI) at negative pressures (<-0.3 GPa). Overall, the MLFF enables large-scale first-principle-level molecular dynamics (MD) simulations of the spontaneous transition from the liquid water to a plethora of monolayer ices, including hexagonal, pentagonal, square, zigzag (ZZMI), and hexatic monolayer ices. These findings will enrich our understanding of the phase behaviour of the nanoconfined water/ices and provide a guide for future experimental realization of the 2D ices.

Two-dimensional partitioned square ice confined in graphene/graphite nanocapillaries

As one of the most fascinating confined water/ice phenomena, two-dimensional square ice has been extensively studied and experimentally confirmed in recent years. Apart from the unidirectional homogeneous square icing patterns considered in previous studies, the multidirectional partitioned square icing patterns are discovered in this study and characterized by molecular dynamics (MD) simulations. Square icing parameters are proposed to quantitatively distinguish the partitioned patterns from the homogeneous patterns and the liquid water. The number of graphene monolayers n is varied in this study, and the results show that it is more energetically favorable to form partitioned square icing patterns when the water molecules are confined between graphite sheets (n ≥ 2) compared to graphene (n = 1). This phenomenon is insensitive to n as long as n ≥ 2 because of the short-range nature of the interaction between water molecules and the carbon substrate. Moreover, it is energetically unfavorable to form partitioned square icing patterns for a single layer of water molecules even for n ≥ 2, verifying that the interaction between layers of water molecules is another dominant factor in the formation of partitioned structures. The conversion from partitioned structure to homogeneous square patterns is investigated by changing the pressure and the temperature. Based on the comprehensive MD simulations, this study unveils the formation mechanism of the partitioned square icing patterns.

Ab Initio Simulations of Water/Metal Interfaces

Structures and processes at water/metal interfaces play an important technological role in electrochemical energy conversion and storage, photoconversion, sensors, and corrosion, just to name a few. However, they are also of fundamental significance as a model system for the study of solid-liquid interfaces, which requires combining concepts from the chemistry and physics of crystalline materials and liquids. Particularly interesting is the fact that the water-water and water-metal interactions are of similar strength so that the structures at water/metal interfaces result from a competition between these comparable interactions. Because water is a polar molecule and water and metal surfaces are both polarizable, explicit consideration of the electronic degrees of freedom at water/metal interfaces is mandatory. In principle, ab initio molecular dynamics simulations are thus the method of choice to model water/metal interfaces, but they are computationally still rather demanding. Here, ab initio simulations of water/metal interfaces will be reviewed, starting from static systems such as the adsorption of single water molecules, water clusters, and icelike layers, followed by the properties of liquid water layers at metal surfaces. Technical issues such as the appropriate first-principles description of the water-water and water-metal interactions will be discussed, and electrochemical aspects will be addressed. Finally, more approximate but numerically less demanding approaches to treat water at metal surfaces from first-principles will be briefly discussed.

Nanoconfined Fluids: Uniqueness of Water Compared to Other Liquids

Nanoconfinement can drastically change the behavior of liquids, puzzling us with counterintuitive properties. It is relevant in applications, including decontamination and crystallization control. However, it still lacks a systematic analysis for fluids with different bulk properties. Here we address this gap. We compare, by molecular dynamics simulations, three different liquids in a graphene slit pore: (1) A simple fluid, such as argon, described by a Lennard-Jones potential; (2) an anomalous fluid, such as a liquid metal, modeled with an isotropic core-softened potential; and (3) water, the prototypical anomalous liquid, with directional HBs. We study how the slit-pore width affects the structure, thermodynamics, and dynamics of the fluids. All the fluids show similar oscillating properties by changing the pore size. However, their free-energy minima are quite different in nature: (i) are energy-driven for the simple liquid; (ii) are entropy-driven for the isotropic core-softened potential; and (iii) have a changing nature for water. Indeed, for water, the monolayer minimum is entropy driven, at variance with the simple liquid, while the bilayer minimum is energy driven, at variance with the other anomalous liquid. Also, water has a large increase in diffusion for subnm slit pores, becoming faster than bulk. Instead, the other two fluids have diffusion oscillations much smaller than water, slowing down for decreasing slit-pore width. Our results, clarifying that water confined at the subnm scale behaves differently from other (simple or anomalous) fluids under similar confinement, are possibly relevant in nanopores applications, for example, in water purification from contaminants.

Static and dynamic water structures at interfaces: A case study with focus on Pt(111)

An accurate atomistic treatment of aqueous solid-liquid interfaces necessitates the explicit description of interfacial water ideally via ab initio molecular dynamics simulations. Many applications, however, still rely on static interfacial water models, e.g., for the computation of (electro)chemical reaction barriers and focus on a single, prototypical structure. In this work, we systematically study the relation between density functional theory-derived static and dynamic interfacial water models with specific focus on the water-Pt(111) interface. We first introduce a general construction protocol for static 2D water layers on any substrate, which we apply to the low index surfaces of Pt. Subsequently, we compare these with structures from a broad selection of reference works based on the Smooth Overlap of Atomic Positions descriptor. The analysis reveals some structural overlap between static and dynamic water ensembles; however, static structures tend to overemphasize the in-plane hydrogen bonding network. This feature is especially pronounced for the widely used low-temperature hexagonal ice-like structure. In addition, a complex relation between structure, work function, and adsorption energy is observed, which suggests that the concentration on single, static water models might introduce systematic biases that are likely reduced by averaging over consistently created structural ensembles, as introduced here.

First-Principles Molecular Dynamics Simulations of the Spontaneous Freezing Transition of 2D Water in a Nanoslit

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.

Replica exchange MD simulations of two-dimensional water in graphene nanocapillaries: rhombic versus square structures, proton ordering, and phase transitions

The hydrogen bond patterns, proton ordering, and phase transitions of monolayer ice in two-dimensional hydrophobic confinement are fundamentally different from those found for bulk ice. To investigate the behavior of quasi-2D ice, we perform molecular dynamics simulations of water confined between fixed graphene plates at a distance of 0.65 nm. While experimental results are still limited and theoretical investigations are often based on a single, often empirically based force field model, this work presents a systematic study modeling the water-graphene interaction by effective Lennard-Jones potentials previously derived from high-level ab initio CCSD(T) calculations of water adsorbed on graphene [Phys. Chem. Chem. Phys., 2013, 15, 4995]. For the water-water interaction different water force fields, i.e. SPCE, TIP3P, TIP4P, TIP4P/ICE, and TIP5P, are used. The water occupancy of the graphene capillary at a pressure of 1000 MPa is determined to be between 13.5 and 13.9 water molecules per square nanometer, depending on the choice of the water force field. Based on these densities, we explore the structure and dynamics of quasi-2D water for temperatures ranging from 200 K to about 600 K for each of the five force fields. To ensure complete sampling of the configurational space and to overcome the barriers separating metastable structures, these simulations are based on the replica exchange molecular dynamics technique. We report different tetragonal hydrogen bond patterns, which are classified as nearly square or as rhombic. While many of these arrangements are essentially flat, in some cases puckered arrangements are found, too. Also the proton ordering of the quasi-2D water structures is considered, allowing us to identify them as ferroelectric, ferrielectric or antiferroelectric. For temperatures between 200 K and 400 K we find several second-order phase transitions from one ice structure to another, changing in many cases both the arrangements of the oxygen atoms and the proton ordering. For temperatures between 400 K and 600 K there are melting-like transitions from a monolayer of ice to a monolayer of liquid water. These first-order phase transitions have a latent heat between 3.4 and 4.0 kJ mol^-1. Both the values of the transition temperatures and of the latent heats display considerable model dependence for the five different water models investigated here.

Ionized water confined in graphene nanochannels

When confined between graphene layers, water in the presence of additional hydronium and hydroxide ions exhibits distinct properties such as ion layering structure determined by the channel size, disruption of the ion solvation shell, and slower ion recombination rate as compared to bulk water.

Phase behaviors of deeply supercooled bilayer water unseen in bulk water

Akin to bulk water, water confined to an isolated nanoslit can show a wealth of new 2D phases of ice and amorphous ice, as well as unusual phase behavior. Indeed, 2D water phases, such as bilayer hexagonal ice and monolayer square ice, have been detected in the laboratory, confirming earlier computational predictions. Herein, we report theoretical evidence of a hitherto unreported state, namely, bilayer very low density amorphous ice (BL-VLDA), as well as evidence of a strong first-order transition between BL-VLDA and the BL amorphous ice (BL-A), and a weak first-order transition between BL-VLDA and the BL very low density liquid (BL-VLDL) water. The diffusivity of BL-VLDA is typically in the range of 10-9 cm2/s to 10-10 cm2/s. Similar to bulk (3D) water, 2D water can exhibit two forms of liquid in the deeply supercooled state. However, unlike supercooled bulk water, for which the two forms of liquid can coexist and merge into one at a critical point, the 2D BL-VLDL and BL high-density liquid (BL-HDL) phases are separated by the highly stable solid phase of BL-A whose melting line exhibits the isochore end point (IEP) near 220 K in the temperature-pressure diagram. Above the IEP temperature, BL-VLDL and BL-HDL are indistinguishable. At negative pressures, the metastable BL-VLDL exhibits a spatially and temporally heterogeneous structure induced by dynamic changes in the nanodomains, a feature much less pronounced in the BL-HDL.

Simulations of water nano-confined between corrugated planes

Water confined to nanoscale widths in two dimensions between ideal planar walls has been the subject of ample study, aiming at understanding the intrinsic response of water to confinement, avoiding the consideration of the chemistry of actual confining materials. In this work, we study the response of such nanoconfined water to the imposition of a periodicity in the confinement by means of computer simulations, both using empirical potentials and from first-principles. For that we propose a periodic confining potential emulating the atomistic oscillation of the confining walls, which allows varying the lattice parameter and amplitude of the oscillation. We do it for a triangular lattice, with several values of the lattice parameter: one which is ideal for commensuration with layers of Ih ice and other values that would correspond to more realistic substrates. For the former, the phase diagram shows an overall rise of the melting temperature. The liquid maintains a bi-layer triangular structure, however, despite the fact that it is not favoured by the external periodicity. The first-principles liquid is significantly affected by the modulation in its layering and stacking even at relatively small amplitudes of the confinement modulation. Beyond some critical modulation amplitude, the hexatic phase present in flat confinement is replaced by a trilayer crystalline phase unlike any of the phases encountered for flat confinement. For more realistic lattice parameters, the liquid does not display higher tendency to freeze, but it clearly shows inhomogeneous behaviour as the strength of the rugosity increases. In spite of this expected inhomogeneity, the structural and dynamical response of the liquid is surprisingly insensitive to the external modulation. Although the first-principles calculations give a more triangular liquid than the one observed with empirical potentials (TIP4P/2005), both agree remarkably well for the main conclusions of the study.

Structures and thermodynamics of water encapsulated by graphene

Understanding phase behaviors of nanoconfined water has driven notable research interests recently. In this work, we examine water encapsulated under a graphene cover that offers an ideal testbed to explore its molecular structures and thermodynamics. We find layered water structures for up to ~1000 trapped water molecules, which is stabilized by the spatial confinement and pressure induced by interfacial adhesion. For monolayer encapsulations, we identify representative two-dimensional crystalline lattices as well as defects therein. Free energy analysis shows that the structural orders with low entropy are compensated by high formation energies due to the pressurized confinement. There exists an order-to-disorder transition for this condensed phase at ~480-490 K, with a sharp reduction in the number of hydrogen bonds and increase in the entropy. Fast diffusion of the encapsulated water demonstrates anomalous temperature dependence, indicating the solid-to-fluid nature of this structural transition. These findings offer fundamental understandings of the encapsulated water that can be used as a pressurized cell with trapped molecular species, and provide guidance for practical applications with its presence, for example, in the design of nanodevices and nanoconfined reactive cells.

CO Separation from H2 via Hydrate Formation in Single-Walled Carbon Nanotubes

Hydrogen is an alternative fuel without generating greenhouse gas or other harmful emissions. Industrial hydrogen production, however, always contains a small fraction of carbon monoxide (CO) (∼0.5-2%) that must be removed for use in fuel cells. Here, we present molecular dynamics simulation evidence on facile separation of CO from H₂ at ambient pressure via the formation of quasi-one-dimensional (Q1D) clathrate hydrates within single-walled carbon nanotubes (SW-CNTs). At ambient pressure, Q1D CO (or H₂) clathrates in SW-CNTs are formed spontaneously when the SW-CNTs are immersed in CO (or H₂) aqueous solution. More interestingly, for the CO/H₂ aqueous solution, highly preferential adsorption of CO over H₂ occurs within the octagonal or nonagonal ice nanotubes inside of SW-CNTs. These results suggest that the formation of Q1D hydrates within SW-CNTs can be a viable and safe method for the separation of CO from H₂, which can be exploited for hydrogen purification in fuel cells.

Buckling failure of square ice-nanotube arrays constrained in graphene nanocapillaries

Graphene confinement provides a new physical and mechanical environment with ultrahigh van der Waals pressure, resulting in new quasi-two-dimensional phases of few-layer ice. Polymorphic transition can occur in bilayer constrained water/ice system. Here, we perform a comprehensive study of the phase transition of AA-stacked bilayer water constrained within a graphene nanocapillary. The compression-limit and superheating-limit (phase) diagrams are obtained, based on the extensive molecular-dynamics simulations at numerous thermodynamic states. Liquid-to-solid, solid-to-solid, and solid-to-liquid-to-solid phase transitions are observed in the compression and superheating of bilayer water. Interestingly, there is a temperature threshold (∼275 K) in the compression-limit diagram, which indicates that the first-order and continuous-like phase transitions of bilayer water depend on the temperature. Two obviously different physical processes, compression and superheating, display similar structural evolution; that is, square ice-nanotube arrays (BL-VHDI) will bend first and then transform into bilayer triangular AA stacking ice (BL-AAI). The superheating limit of BL-VHDI exhibits local maxima, while that of BL-AAI increases monotonically. More importantly, from a mechanics point of view, we propose a novel mechanism of the transformation from BL-VHDI to BL-AAI, both for the compression and superheating limits. This structural transformation can be regarded as the "buckling failure" of the square-ice-nanotube columns, which is dominated by the lateral pressure.

Continuous melting through a hexatic phase in confined bilayer water

Liquid water is not only of obvious importance but also extremely intriguing, displaying many anomalies that still challenge our understanding of such an a priori simple system. The same is true when looking at nanoconfined water: The liquid between constituents in a cell is confined to such dimensions, and there is already evidence that such water can behave very differently from its bulk counterpart. A striking finding has been reported from computer simulations for two-dimensionally confined water: The liquid displays continuous or discontinuous melting depending on its density. In order to understand this behavior, we have analyzed the melting exhibited by a bilayer of nanoconfined water by means of molecular dynamics simulations. At high density we observe the continuous melting to be related to the phase change of the oxygens only, with the hydrogens remaining liquidlike throughout. Moreover, we find an intermediate hexatic phase for the oxygens between the liquid and a triangular solid ice phase, following the Kosterlitz-Thouless-Halperin-Nelson-Young theory for two-dimensional melting. The liquid itself tends to maintain the local structure of the triangular ice, with its two layers being strongly correlated yet with very slow exchange of matter. The decoupling in the behavior of the oxygens and hydrogens gives rise to a regime in which the complexity of water seems to disappear, resulting in what resembles a simple monoatomic liquid. This intrinsic tendency of our simulated water may be useful for understanding novel behaviors in other confined and interfacial water systems.

Structural and configurational properties of nanoconfined monolayer ice from first principles

Understanding the structural tendencies of nanoconfined water is of great interest for nanoscience and biology, where nano/micro-sized objects may be separated by very few layers of water. Here we investigate the properties of ice confined to a quasi-2D monolayer by a featureless, chemically neutral potential, in order to characterize its intrinsic behaviour. We use density-functional theory simulations with a non-local van der Waals density functional. An ab initio random structure search reveals all the energetically competitive monolayer configurations to belong to only two of the previously-identified families, characterized by a square or honeycomb hydrogen-bonding network, respectively. We discuss the modified ice rules needed for each network, and propose a simple point dipole 2D lattice model that successfully explains the energetics of the square configurations. All identified stable phases for both networks are found to be non-polar (but with a topologically non-trivial texture for the square) and, hence, non-ferroelectric, in contrast to previous predictions from a five-site empirical force-field model. Our results are in good agreement with very recently reported experimental observations.

AB-stacked square-like bilayer ice in graphene nanocapillaries

Water, when constrained between two graphene sheets and under ultrahigh pressure, can manifest dramatic differences from its bulk counterparts such as the van der Waals pressure induced water-to-ice transformation, known as the metastability limit of two-dimensional (2D) liquid.