Hydrogen-Bonding in Liquid Water at Multikilobar Pressures

Generating interstitial water within the persisting tetrahedral H-bond network explains density increase upon compressing liquid water

Despite its ubiquitous nature, the atomic structure of water in its liquid state is still controversially debated. We use a combination of X-ray Raman scattering spectroscopy in conjunction with ab initio and path integral molecular dynamics simulations to study the local atomic and electronic structure of water under high pressure conditions. Systematically increasing fingerprints of non-hydrogen-bonded H[Formula: see text]O molecules in the first hydration shell are identified in the experimental and computational oxygen K-edge excitation spectra. This provides evidence for a compaction mechanism in terms of a continuous collapse of the second hydration shell with increasing pressure via generation of interstitial water within locally tetrahedral hydrogen-bonding environments.

Solvation Properties of Neutral Gold Species in Supercritical Water Studied By THz Spectroscopy

Supercritical water provides distinctly different solvation properties compared to what is known from liquid water. Despite its prevalence deep in the Earth's crust and its role in chemosynthetic ecosystems in the vicinity of hydrothermal vents, molecular insights into its solvation mechanisms are still very scarce compared to what is known for liquid water. Recently, neutral metal particles have been detected in hydrothermal fluids and proposed to explain the transport of gold species to ore deposits on Earth. Using ab initio molecular dynamics, we elucidate the solvation properties of small gold species at supercritical conditions. The neutral metal clusters themselves contribute enormous THz intensity not because of their intramolecular vibrations, but due to their pronounced electronic polarization coupling to the dynamical supercritical solvent, leading to a continuum absorption up to about 1000 cm-1. On top, long-lived interactions between the gold clusters and solvation water leads at these supercritical conditions to a sharp THz resonance that happens to be close to the one due to H-bonding in liquid water at ambient conditions. The resulting distinct resonances can be used to analyse the solvation properties of neutral metal particles in supercritical aqueous solutions.

Terahertz spectroscopy of thick and diluted water solutions

While bright terahertz sources are used to perform nonlinear experiments, they can be advantageous for high-precision linear measurements of opaque samples. By placing the sample away from the focus, nonlinearities can be suppressed, and sizeable amounts of transmitted radiation detected. Here, this approach is demonstrated for a 0.5 mm thick layer of liquid water in a static sample holder. Variations of the index of refraction as small as (7 ± 2) · 10-4 were detected at 0.58 THz for an aqueous salt solution containing ten millimoles of sodium chloride. To my knowledge, this precision is unprecedented in time-domain spectroscopy studies of diluted aqueous systems or other optically thick and opaque materials.

Characterization of the Hydrogen-Bond Network in High-Pressure Water by Deep Potential Molecular Dynamics

The hydrogen-bond (H-bond) network of high-pressure water is investigated by neural-network-based molecular dynamics (MD) simulations with first-principles accuracy. The static structure factors (SSFs) of water at three densities, i.e., 1, 1.115, and 1.24 g/cm3, are directly evaluated from 512 water MD trajectories, which are in quantitative agreement with the experiments. We propose a new method to decompose the computed SSF and identify the changes in the SSF with respect to the changes in H-bond structures. We find that a larger water density results in a higher probability for one or two non-H-bonded water molecules to be inserted into the inner shell, explaining the changes in the tetrahedrality of water under pressure. We predict that the structure of the accepting end of water molecules is more easily influenced by the pressure than by the donating end. Our work sheds new light on explaining the SSF and H-bond properties in related fields.

Water structure modification by d-(+)-glucose at different concentrations and temperatures-effect of mutarotation

Water structure modification by carbohydrates is essential both in chemistry and life processes and in particular, molecular level interaction of glucose with water is very important. With a view to developing a fundamental knowledge base, thermodynamic parameters derived from measurements of density, viscosity, and refractive index have been analyzed to investigate how d-(+)-glucose alters the structure of water at various concentrations and temperatures. The nature and extent of the interactions have been investigated using apparent molar volume, Jones-Dole constants, changes in free energy (ΔG), changes in entropy (ΔS), and changes in enthalpy (ΔH) for viscous flow. Using measurements from dynamic light scattering (DLS), the sizes of the aggregates were studied. The kinetics of mutarotation have been investigated using polarimetry and the structural effect on water during mutarotation between α-d-glucose and β-d-glucose with time has been explored by near-infrared (NIR) spectroscopy. The spectroscopic results were examined using difference spectroscopy and two-dimensional correlation spectroscopy (2DCOS). The absorption bands of water shift to a higher wavenumber irrespective of the concentration of the solution with time due to the enhancement of the cleavage of hydrogen bonding in water. At high temperatures, three bands in the region 7100-7350 cm^-1 are attributed to the first overtones of the hydrogen-bonded -O-H stretching vibration. Refractive index values indicate an increase in the density of the anomer solutions with time, suggesting an increase in free water concentration. These results provide evidence for more than one water molecule being involved in the mechanism of mutarotation and propose a concerted mechanism for proton transfer.

Life in Multi-Extreme Environments: Brines, Osmotic and Hydrostatic Pressure─A Physicochemical View

Elucidating the details of the formation, stability, interactions, and reactivity of biomolecular systems under extreme environmental conditions, including high salt concentrations in brines and high osmotic and high hydrostatic pressures, is of fundamental biological, astrobiological, and biotechnological importance. Bacteria and archaea are able to survive in the deep ocean or subsurface of Earth, where pressures of up to 1 kbar are reached. The deep subsurface of Mars may host high concentrations of ions in brines, such as perchlorates, but we know little about how these conditions and the resulting osmotic stress conditions would affect the habitability of such environments for cellular life. We discuss the combined effects of osmotic (salts, organic cosolvents) and hydrostatic pressures on the structure, stability, and reactivity of biomolecular systems, including membranes, proteins, and nucleic acids. To this end, a variety of biophysical techniques have been applied, including calorimetry, UV/vis, FTIR and fluorescence spectroscopy, and neutron and X-ray scattering, in conjunction with high pressure techniques. Knowledge of these effects is essential to our understanding of life exposed to such harsh conditions, and of the physical limits of life in general. Finally, we discuss strategies that not only help us understand the adaptive mechanisms of organisms that thrive in such harsh geological settings but could also have important ramifications in biotechnological and pharmaceutical applications.

Aqueous TMAO solution under high hydrostatic pressure

Trimethylamine N-oxide (TMAO) is a well known osmolyte in nature, which is used by deep sea fish to stabilize proteins against High Hydrostatic Pressure (HHP). We present a combined ab initio molecular dynamics, force field molecular dynamics, and THz absorption study of TMAO in water up to 12 kbar to decipher its solvation properties upon extreme compression. On the hydrophilic oxygen side of TMAO, AIMD simulations at 1 bar and 10 kbar predict a change of the coordination number from a dominating TMAO·(H₂O)₃ complex at ambient conditions towards an increased population of a TMAO·(H₂O)₄ complex at HHP conditions. This increase of the TMAO-oxygen coordination number goes in line with a weakening of the local hydrogen bond network, spectroscopic shifts and intensity changes of the corresponding intermolecular THz bands. Using a pressure-dependent HHP force field, FFMD simulations predict a significant increase of hydrophobic hydration from 1 bar up to 4-5 kbar, which levels off at higher pressures up to 10 kbar. THz spectroscopic data reveal two important pressure regimes with spectroscopic inflection points of the dominant intermolecular modes: The first regime (1.5-2 kbar) is barely recognizable in the simulation data. However, it relates well with the observation that the apparent molar volume of solvated TMAO is nearly constant in the biologically relevant pressure range up to 1 kbar as found in the deepest habitats on Earth in the ocean. The second inflection point around 4-5 kbar is related to the amount of hydrophobic hydration as predicted by the FFMD simulations. In particular, the blueshift of the intramolecular CNC bending mode of TMAO at about 390 cm^-1 is the spectroscopic signature of increasingly pronounced pressure-induced changes in the solvation shell of TMAO. Thus, the CNC bend can serve as local pressure sensor in the multi-kbar pressure regime.

Pressure-Dependent Structure of Methanol-Water Mixtures up to 1.2 GPa: Neutron Diffraction Experiments and Molecular Dynamics Simulations

Total scattering structure factors of per-deuterated methanol and heavy water, CD₃OD and D₂O, have been determined across the entire composition range as a function of pressure up to 1.2 GPa, by neutron diffraction. The largest variations due to increasing pressure were observed below a scattering variable value of 5 Å⁻¹, mostly as shifts in terms of the positions of the first and second maxima. Molecular dynamics computer simulations, using combinations of all-atom potentials for methanol and various water force fields, were conducted at the experimental pressures with the aim of interpreting neutron diffraction results. The peak-position shifts mentioned above could be qualitatively reproduced by simulations, although in terms of peak intensities, the accord between neutron diffraction and molecular dynamics was much less satisfactory. However, bearing in mind that increasing pressure must have a profound effect on repulsive forces between neighboring molecules, the agreement between experiment and computer simulation can certainly be termed as satisfactory. In order to reveal the influence of changing pressure on local intermolecular structure in these “simplest of complex” hydrogen-bonded liquid mixtures, simulated structures were analyzed in terms of hydrogen bond-related partial radial distribution functions and size distributions of hydrogen-bonded cyclic entities. Distinct differences between pressure-dependent structures of water-rich and methanol-rich composition regions were revealed.

Dynamical Collectivity and Nuclear Quantum Effects on the Intermolecular Stretching Mode of Liquid Water

This study investigated the broadband terahertz and low-frequency Raman spectroscopy of liquid water (H₂O, D₂O, and H₂¹⁸O) over 2 decades of frequency to address long-standing challenges regarding the interpretation of the intermolecular stretching mode at around 5 THz. We experimentally demonstrated that the intermolecular stretching mode of liquid water obtained via terahertz spectroscopy is significantly redshifted and broadened compared with that via Raman. This result was rationalized by the enhanced dynamical collectivity probed by terahertz spectroscopy, although both have a common origin in the kinetic motion. Their temperature and isotope dependences emphasize the significance of oscillation mass in determining the intermolecular stretching lineshape, while quantum effects cannot be overlooked in both terahertz and low-frequency Raman spectra.

An isolated water droplet in the aqueous solution of a supramolecular tetrahedral cage

Water under nanoconfinement at ambient conditions has exhibited low-dimensional ice formation and liquid-solid phase transitions, but with structural and dynamical signatures that map onto known regions of water's phase diagram. Using terahertz (THz) absorption spectroscopy and ab initio molecular dynamics, we have investigated the ambient water confined in a supramolecular tetrahedral assembly, and determined that a dynamically distinct network of 9 ± 1 water molecules is present within the nanocavity of the host. The low-frequency absorption spectrum and theoretical analysis of the water in the Ga₄L₆ ^12- host demonstrate that the structure and dynamics of the encapsulated droplet is distinct from any known phase of water. A further inference is that the release of the highly unusual encapsulated water droplet creates a strong thermodynamic driver for the high-affinity binding of guests in aqueous solution for the Ga₄L₆ ^12- supramolecular construct.

Supercritical Water is not Hydrogen Bonded

Thinking about water is inextricably linked to hydrogen bonds, which are highly directional in character and determine the unique structure of water, in particular its tetrahedral H-bond network. Here, we assess if this common connotation also holds for supercritical water. We employ extensive ab initio molecular dynamics simulations to systematically monitor the evolution of the H-bond network mode of water from room temperature, where it is the hallmark of its fluctuating three-dimensional network structure, to supercritical conditions. Our simulations reveal that the oscillation period required for H-bond vibrations to occur exceeds the lifetime of H-bonds in supercritical water by far. Instead, the corresponding low-frequency intermolecular vibrations of water pairs as seen in supercritical water are found to be well represented by isotropic van-der-Waals interactions only. Based on these findings, we conclude that water in its supercritical phase is not a H-bonded fluid.

How Can Protons Migrate in Extremely Compressed Liquid Water?

Compression of liquid water up to multi-kbar pressures is known to perturb dramatically its local structure required for charge defects to migrate as topological defects in the hydrogen-bonded network. Our ab initio simulations show that the migration of excess protons is not much affected at 10 kbar, whereas that of proton holes is significantly reduced. Non-Markovian analyses show that this is not due to modifying the free energy barriers of both charge transfer and migration. It is rather pressure-induced modifications of the population of activated states, depending on interstitial water, which rules charge migration at extreme compression.

Evidence of a Low-High Density Turning Point in Liquid Water at Ordinary Temperature under Pressure: A Molecular Dynamics Study

Water has a fundamental role in important processes spanning a wide range of pressure and temperature conditions. Knowledge of structural, dynamic and thermodynamic properties of water at nonstandard conditions is a primary concern since interest in astronomical, geological, and technological processes is continuously growing. Molecular dynamics simulations allow us to study thermodynamic conditions that require sophisticated techniques and instruments, while at the same time offering the interpretation of properties at the atomic level. It is established that the behavior of water is strongly affected by the temperature and pressure conditions, determining the existence of low and high density regimes. For the first time, a thermodynamic property, isothermal compressibility, has been adopted to detect the low-high density turning point at ambient temperature in liquid water due to pressure. Molecular dynamics simulations have been performed with five three-site models, allowing us to characterize the complexity of water nature at these conditions at the atomic level.