Does twist angle affect the properties of water confined inside twisted bilayer graphene?

Graphene nanoslit pores are used for nanofluidic devices, such as, in water desalination, ion-selective channels, ionic transistors, sensing, molecular sieving, blue energy harvesting, and protein sequencing. It is a strenuous task to prepare nanofluidic devices, because a small misalignment leads to a significant alteration in various properties of the devices. Here, we focus on the rotational misalignment between two parallel graphene sheets. Using molecular dynamics simulation, we probe the structure and dynamics of monolayer water confined inside graphene nanochannels for a range of commensurate twist angles. With SPC/E and TIP4P/2005 water models, our simulations reveal the independence of the equilibrium number density- n ∼ 13 nm^-2 for SPC/E and n ∼ 11.5 nm^-2 for TIP4P/2005- across twists. Based on the respective densities of the water models, the structure and dielectric constant are invariant of twist angles. The confined water structure at this density shows square ice ordering for SPC/E water only. TIP4P/2005 shows ordering at the vicinity of a critical density (n ∼ 12.5 nm^-2). The average perpendicular dielectric constant of the confined water remains anomalously low (∼2 for SPC/E and ∼6 for TIP4P/2005) for the studied twist angles. We find that the friction coefficient of confined water molecules varies for small twist angles, while becoming independent for twists greater than 5.1°. Our results indicate that a small, angular misalignment will not impair the dielectric properties of monolayer water within a graphene slit-pore, but can significantly influence its dynamics.

Dynamical Manifestations of Supercooling in Amyloid Hydration

The effect of extreme temperature on amyloidogenic species remains sparsely explored. In a recent study (J. Phys. Chem. Lett., 2019, 10, (10)), we employed exhaustive molecular dynamics simulations to explore the cold thermal response of a putative small amyloid oligomer and to elicit the role of solvent modulation. Herein, we investigate the dynamical response of the hydration waters of the oligomer within the supercooled states. Using NMR-based formalism, we delineate the entropic response in terms of the side-chain conformational entropy that corroborates the weakening of the hydrophobic core with lowering of temperature. The translational dynamics of the protein and hydration waters reveal the coupling of protein dynamical fluctuations with solvent dynamics under supercooled conditions. Probing the translational motion as a space-time correlation indicates glassy dynamics exhibited by hydration waters in the supercooled regime. Caging of the water molecules with lowering of temperature and the resultant hopping dynamics are reflected in the longer β-relaxation timescales of translational motion. Furthermore, we utilized mode-coupling theory (MCT) and derived the ideal glass transition temperature from translational and rotational dynamics, around ∼196 and 209 K, respectively. Interestingly, rotational motion in the supercooled regime deviates from the MCT law, exhibits Arrhenius motion, and marks a fragile-to-strong crossover at 227 K. The low-frequency vibrational modes also coincide with the dynamical transition. This exposition lends dynamical insights into the hydration coupling of an amyloid aggregate under cryogenic conditions.

Amphiphilicity of Intricate Layered Graphene/g-C3N4 Nanosheets

The hybrid heterostructure of the tri-s-triazine form of graphitic carbon nitride (g-C₃N₄), a stable two-dimensional material, results from intricate layer formation with graphene. In this material, g-C₃N₄, an amphiphilic material, stabilizes Pickering emulsions as an emulsifier and can effectively disperse graphene. Due to the various technological applications of the hybrid nanosheets in an aqueous environment, it is essential to study the interaction of water molecules with graphene and g-C₃N₄ (Gr/g-C₃N₄)-combined heterostructure. Although few studies have been performed signifying the water orientation in the interfacial layer, we find that there is a lack of detailed studies using various dynamical and structural properties of the interfacial water molecules. The interface of the Gr/g-C₃N₄ hybrid structure, one of the rarely found amphiphilic interfaces (on the g-C₃N₃ side), is appropriate for exploring the water affinity due to the availability of heterogeneous interfacial aqueous interactions. We adopted classical molecular dynamics simulations using two models for water molecules to study the structure and dynamics of an aqueous interface. We have correlated the structural properties to dynamics and spectral properties to understand the overall behavior of the amphiphilic interface. Our results branch into two significant hydrogen bond (HB) properties in HB count and HB strength among the water molecules in the different layers. The HB counts in the different layers of water are correlated using the average distance distribution (P_rO4), tetrahedral order parameters, HB donor/acceptor count, and total HBs per water molecule. A conspicuous difference is found in the HB count and related dynamics of the system. The HB lifetime and diffusion coefficient hint at the equivalent strength of HBs in the different layers. All the findings conclude that the amphiphilicity of the Gr/g-C₃N₄ interface can help in understanding various interfacial physical and chemical processes.

Aqueous Affinity and Interfacial Dynamics of Anisotropic Buckled Black Phosphorous

The structure of black phosphorous (BP) is similar to the honeycomb arrangement of graphene, but the layered BP is found to be buckled and highly anisotropic. The buckled surface structure affects interfacial molecule mobility and plays a vital role in various nanomaterial applications. The BP is also known for wettability, droplet formation, stability, and hydrophobicity in the aqueous environment. However, there is a gap concerning the structural and dynamical behavior of water molecules, which is available in abundance for other monoatomic and polyatomic two-dimensional (2D) materials. Motivated by the technological importance, we try to bridge the gap by explaining the surface anisotropy-facilitated behavior of water molecules on bilayer BP using classical and first principles molecular dynamics (MD) simulations. From our classical MD study, we find three distinct layers of water molecules. The water layer closest to the interface is L1, followed by L2 and L3/bulk perpendicular to the BP surface. Water molecules in the L1 layer experience some structural disintegration in hydrogen bond (HB) phenomena compared to the bulk. There is a loss of HB donor-acceptor count per water molecule. The average HB count decreases because of an elevated rate of HB formation and deformation; this would affect the dynamic properties in terms of HB lifetime. Therefore, we observe the reduced lifetime of HB in the layer in close contact with BP, which again complements our finding on the diffusion coefficient of water molecules in distinct layers. Water diffuses relatively faster with diffusion coefficient 3.25 × 10^-9 m² s^-1 in L1, followed by L2 and L3. The BP layer shows moderate hydrophobic nature. Our results also indicate the anisotropic behavior as the diffusion along the x-direction is faster than that along the y-direction. The gap in the slope of the x and y components of mean-squared displacement (MSD) complements the pinning effect in an aqueous environment. We observe blue-shifted and red-shifted libration and O-H stretching modes from the calculated power spectra for the L1 water molecules compared to the L2 and L3 molecules from first principles MD simulations. Our analysis may help understand the physical phenomena that occur during the surface wetting of the predroplet formation process observed experimentally.

Dimensionality dependence of the Kauzmann temperature: A case study using bulk and confined water

The Kauzmann temperature (T_K) of a supercooled liquid is defined as the temperature at which the liquid entropy becomes equal to that of the crystal. The excess entropy, the difference between liquid and crystal entropies, is routinely used as a measure of the configurational entropy, whose vanishing signals the thermodynamic glass transition. The existence of the thermodynamic glass transition is a widely studied subject, and of particular recent interest is the role of dimensionality in determining the presence of a glass transition at a finite temperature. The glass transition in water has been investigated intensely and is challenging as the experimental glass transition appears to occur at a temperature where the metastable liquid is strongly prone to crystallization and is not stable. To understand the dimensionality dependence of the Kauzmann temperature in water, we study computationally bulk water (three-dimensions), water confined in the slit pore of the graphene sheet (two-dimensions), and water confined in the pore of the carbon nanotube of chirality (11,11) having a diameter of 14.9 Å (one-dimension), which is the lowest diameter where amorphous water does not always crystallize into nanotube ice in the supercooled region. Using molecular dynamics simulations, we compute the entropy of water in bulk and under reduced dimensional nanoscale confinement to investigate the variation of the Kauzmann temperature with dimension. We obtain a value of T_K (133 K) for bulk water in good agreement with experiments [136 K (C. A. Angell, Science 319, 582-587 (2008) and K. Amann-Winkel et al., Proc. Natl. Acad. Sci. U. S. A. 110, 17720-17725 (2013)]. However, for confined water, in two-dimensions and one-dimension, we find that there is no finite temperature Kauzmann point (in other words, the Kauzmann temperature is 0 K). Analysis of the fluidicity factor, a measure of anharmonicity in the oscillation of normal modes, reveals that the Kauzmann temperature can also be computed from the difference in the fluidicity factor between amorphous and ice phases.