Distinguishing Weak and Strong Hydrogen Bonds in Liquid Water-A Potential of Mean Force-Based Approach

Reorientation of Hydrogen Bonds Renders Unusual Enhancement in Thermal Transport of Water in Nanoconfined Environments

Liquid confined in a nanochannel or nanotube has exhibited a superfast transport phenomenon, providing an ideal heat and mass transfer platform to meet the increasingly stringent challenge of thermal management in developing high-power-density nanoelectronics and nanochips. However, understanding the thermal transport of confined liquid is currently lacking and is speculated to be fundamentally different from that of bulk counterparts due to the unprecedented thermodynamics of liquid in nanoconfined environments. Here, we report that the thermal conductivity of water confined in a silica nanotube is nearly 2-fold as that of bulk status. Further molecular dynamics simulations reveal that this unusual enhancement originates from the densification and reorientation of local hydrogen bonds close to the nanotubes. Thermal-confinement scaling law is established and quantitatively supported by comprehensive simulations with remarkable agreement. Our findings lay a theoretical foundation for designing nanofluidics-enabled cooling strategies and devices.

Second inflection point of supercooled water surface tension induced by hydrogen bonds: A molecular-dynamics study

Surface tension of supercooled water is a fundamental property in various scientific processes. In this study, we perform molecular dynamics simulations with the TIP4P-2005 model to investigate the surface tension of supercooled water down to 220 K. Our results show a second inflection point (SIP) in the surface tension at temperature TSIP ≈ 267.5 ± 2.3 K. Using an extended IAPWS-E functional fit for the water surface tension, we calculate the surface excess internal-energy and entropy terms of the excess Helmholtz free energy. Similar to prior studies [Wang et al., Phys. Chem. Chem. Phys. 21, 3360 (2019); Gorfer et al., J. Chem. Phys. 158, 054503 (2023)], our results show that the surface tension is governed by two driving forces: a surface excess entropy change above the SIP and a surface excess internal-energy change below it. We study hydrogen-bonding near the SIP because it is the main cause of water's anomalous properties. With decreasing temperature, our results show that the entropy contribution to the surface tension reaches a maximum slightly below the SIP and then decreases. This is because the number of hydrogen bonds increases more slowly below the SIP. Moreover, the strengths and lifetimes of the hydrogen bonds also rise dramatically below the SIP, causing the internal-energy term to dominate the excess surface free energy. Thus, the SIP in the surface tension of supercooled TIP4P-2005 water is associated with an increase in the strengths and lifetimes of hydrogen bonds, along with a decrease in the formation rate (#/K) of new hydrogen bonds.

Prediction of water anomalous properties by introducing the two-state theory in SAFT

Water is one of the most abundant substances on earth, but it is still not entirely understood. It shows unusual behavior, and its properties present characteristic extrema unlike any other fluid. This unusual behavior has been linked to the two-state theory of water, which proposes that water forms different clusters, one with a high density and one with a low density, which may even form two distinct phases at low temperatures. Models incorporating the two-state theory manage to capture the unusual extrema of water, unlike traditional equations of state, which fail. In this work, we have derived the framework to incorporate the two-state theory of water into the Statistical-Associating-Fluid-Theory (SAFT). More specifically, we have assumed that water is an ideal solution of high density water molecules and low density water molecules that are in chemical equilibrium. Using this assumption, we have generalized the association term SAFT to allow for the simultaneous existence of the two water types, which have the same physical parameters but different association properties. We have incorporated the newly derived association term in the context of the Perturbed Chain-SAFT (PC-SAFT). The new model is referred to as PC-SAFT-Two-State (PC-SAFT-TS). Using PC-SAFT-TS, we have succeeded in predicting the characteristic extrema of water, such as its density and speed of sound maximum, etc., without loss of accuracy compared to the original PC-SAFT. This new framework is readily extended to mixtures, and PC-SAFT-TS manages to capture the solubility minimum of hydrocarbons in water in a straightforward manner.

Monovalent ion-graphene oxide interactions are controlled by carboxylic acid groups: Sum frequency generation spectroscopy studies

Graphene oxide (GO) is a two-dimensional, mechanically strong, and chemically tunable material for separations. Elucidating GO-ion-water interactions at the molecular scale is highly important for predictive understanding of separation systems. However, direct observations of the nanometer region by GO surfaces under operando conditions are not trivial. Therefore, thin films of GO at the air/water interface can be used as model systems. With this approach, we study the effects of alkali metal ions on water organization near graphene oxide films at the air/water interface using vibrational sum frequency generation (SFG) spectroscopy. We also use an arachidic acid Langmuir monolayer as a benchmark for a pure carboxylic acid surface. Theoretical modeling of the concentration-dependent sum frequency signal from graphene oxide and arachidic acid surfaces reveals that the adsorption of monovalent ions is mainly controlled by the carboxylic acid groups on graphene oxide. An in-depth analysis of sum frequency spectra reveals at least three distinct water populations with different hydrogen bonding strengths. The origin of each population can be identified from concentration dependent variations of their SFG signal. Interestingly, an interfacial water structure seemed mostly insensitive to the character of the alkali cation, in contrast to similar studies conducted at the silica/water interface. However, we observed an ion-specific effect with lithium, whose strong hydration prevented direct interactions with the graphene oxide film.

Roles of hydrogen bonding interactions and hydrophobic effects on enhanced water structure in aqueous solutions of amphiphilic organic molecules

Insights into the solute-induced water structural transformations are essential to understand the role of water in biological and chemical reaction processes. Herein, the structural changes in water induced by amphiphilic organic molecules were investigated using concentration-dependent derivative Raman spectroscopy (DRS) combined with two-dimensional Raman correlation spectroscopy (2D Raman-COS). We shall restrict our attention in this work to binary mixtures of water with dimethyl sulfoxide (DMSO), acetone, and isopropanol (IPA), all of which have similar chemical structures. The spectral changes in O:H and OH stretching modes illustrate that the solute molecules induce an enhancement of the water structure in dilute solutions, where the enhanced degree of water structure is closely related to the size of the dipole moment of organic molecules. In addition, the transformations of solute-induced water-specific structures were evaluated by 2D Raman-COS, which shows that the strong hydrogen bond (H-bond) structure of water is more sensitive to organic molecules and induces a transition to the weak H-bond structure of water.

Exploring the Molecular Dynamics of a Lipid-A Vesicle at the Atom Level: Morphology and Permeation Mechanism

Lipid-A was previously shown to spontaneously aggregate into a vesicle via the hybrid particle field approach. We assess the validity of the proposed vesiculation mechanism by simulating the resulting lipid-A vesicle at the atom level. The spatial confinement imposed by the vesicle geometry on the conformation and packing of lipid-A induces significant heterogeneity of physical properties in the inner and outer leaflets. It also induces tighter molecular packing and lower acyl chain order compared to the lamellar arrangement. Around 5% of water molecules passively permeates the vesicle membrane inward and outward. The permeation is facilitated by interactions with water molecules that are transported across the membrane by a network of electrostatic interactions with the hydrogen bond donors/acceptors in the N-acetylglucosamine ring and upper region of the acyl chains of lipid-A. The permeation process takes place at low rates but still at higher frequencies than observed for the lamellar arrangement of lipid-A. These findings not only substantiate the proposed lipid-A vesiculation mechanism but also reveal the complex structural dynamics of an important nonlamellar arrangement of lipid-A.

Light, Water, and Melatonin: The Synergistic Regulation of Phase Separation in Dementia

The swift rise in acceptance of molecular principles defining phase separation by a broad array of scientific disciplines is shadowed by increasing discoveries linking phase separation to pathological aggregations associated with numerous neurodegenerative disorders, including Alzheimer’s disease, that contribute to dementia. Phase separation is powered by multivalent macromolecular interactions. Importantly, the release of water molecules from protein hydration shells into bulk creates entropic gains that promote phase separation and the subsequent generation of insoluble cytotoxic aggregates that drive healthy brain cells into diseased states. Higher viscosity in interfacial waters and limited hydration in interiors of biomolecular condensates facilitate phase separation. Light, water, and melatonin constitute an ancient synergy that ensures adequate protein hydration to prevent aberrant phase separation. The 670 nm visible red wavelength found in sunlight and employed in photobiomodulation reduces interfacial and mitochondrial matrix viscosity to enhance ATP production via increasing ATP synthase motor efficiency. Melatonin is a potent antioxidant that lowers viscosity to increase ATP by scavenging excess reactive oxygen species and free radicals. Reduced viscosity by light and melatonin elevates the availability of free water molecules that allow melatonin to adopt favorable conformations that enhance intrinsic features, including binding interactions with adenosine that reinforces the adenosine moiety effect of ATP responsible for preventing water removal that causes hydrophobic collapse and aggregation in phase separation. Precise recalibration of interspecies melatonin dosages that account for differences in metabolic rates and bioavailability will ensure the efficacious reinstatement of the once-powerful ancient synergy between light, water, and melatonin in a modern world.

Hydrogen-Bond Dynamics and Water Structure in Aqueous Ethylene Glycol Solution via Two-Dimensional Raman Correlation Spectroscopy

The hydrogen-bond (H-bond) dynamics and water structural transitions in aqueous ethylene glycol (EG) solution were investigated on the basis of concentration- and temperature-dependent two-dimensional Raman correlation spectroscopy (2D Raman-COS). At room temperature, EG-induced enhancement of the water structure when the EG/water molar ratio is less than 1:28 resulted from the hydrophobic effect around the methylene groups of EG. The decrease in the temperature caused an enhancement of the Raman peak at about 3200 cm^-1, representing an increase in the orderliness of water molecules. Further analysis of the water-specific structures by 2D Raman-COS reveals that the strong H-bond structure preferentially responds to external perturbations and induces a weak H-bond structural transition in water. Finally, EG-induced water structural transitions were calculated by the density functional theory (DFT). Hopefully, 2D Raman-COS combined with DFT calculations would advance the study of solute-induced water structural transitions in water-organic chemistry.

Nature of Cations Critically Affects Water at the Negatively Charged Silica Interface

Understanding the collective behavior of ions at charged surfaces is of paramount importance for geological and electrochemical processes. Ions screen the surface charge, and interfacial fields break the centro-symmetry near the surface, which can be probed using second-order nonlinear spectroscopies. The effect of electrolyte concentration on the nonlinear optical response has been semi-quantitatively explained by mean-field models based on the Poisson-Boltzmann equation. Yet, to explain previously reported ion-specific effects on the spectroscopic response, drastic ion-specific changes in the interfacial properties, including surface acidities and dielectric permittivities, or strong ion adsorption/desorption had to be invoked. Here, we use sum-frequency generation (SFG) spectroscopy to probe the symmetry-breaking of water molecules at a charged silica surface in contact with alkaline metal chloride solutions (LiCl, NaCl, KCl, and CsCl) at various concentrations. We find that the water response varies with the cation: the SFG response is markedly enhanced for LiCl compared to CsCl. We show that within mean-field models, neither specific ion-surface interactions nor a reduced dielectric constant of water near the interface can account for the variation of spectral intensities with cation nature. Molecular dynamics simulations confirm that the decay of the electrochemical potential only weakly depends on the salt type. Instead, the effect of different salts on the optical response is indirect, through the reorganization of the interfacial water: the salt-type-dependent alignment of water directly at the interface can explain the observations.

Multiple insights call for revision of modern thermodynamic models to account for structural fluctuations in water

Modern thermodynamic models incorporate the concept of association (hydrogen bonding) and they can describe very satisfactorily many properties of water containing mixtures. They have not been successful in representing water's anomalous properties and this work provides a possible explanation. We have analyzed and interpreted recent experimental data, molecular simulation results, and two-state theory approaches and compared against the predictions from thermodynamic models. We show that the dominance of the tetrahedral structure implemented in modern thermodynamic models may be the reason for their failure for describing water systems. While this study does not prove the two-state theories for water, it indicates that a high level of tetrahedral structure of water is not in agreement with water's anomalous properties when used in thermodynamic models.

Structural characteristics of low-density environments in liquid water

The existence of two structural forms in liquid water has been a point of discussion for a long time. A phase transition between these two forms of liquid water has been proposed based on evidence from molecular simulations, and experiments have also been very recently able to track the proposed transition of the low-density liquid form to the high-density liquid form. We propose to use the average angle an oxygen atom makes with its neighbors to describe the structural environment of a water molecule. The distribution of this order parameter is observed to have two peaks with one peak at ∼109.5^{∘}, corresponding to the internal angle of a regular tetrahedron, indicating tetrahedral arrangement. The other peak corresponds to an environment with a tighter arrangement of neighboring molecules. The distribution of O-O-O angles is decomposed into two skewed distributions to estimate the fractions of the two liquid forms in water. A good similarity is observed between the temperature and pressure trends of fractions of locally favored tetrahedral structure (LFTS) form estimated using the new order parameter and the reports in the literature, over a range of temperatures and pressures. We also compare the structural environments indicated by different order parameters and find that the order parameter proposed in this paper captures the structure of first solvation shell of the LFTS accurately.