Vibrational Spectroscopy as a Probe of Structure and Dynamics in Liquid Water

Entropy of Liquid Water as Predicted by the Two-Phase Thermodynamic Model and Data-Driven Many-Body Potentials

We investigate the entropy of liquid water at ambient conditions using the two-phase thermodynamic (2PT) model, applied to both common pairwise-additive water models and the MB-pol and MB-pol(2023) data-driven many-body potentials. Our simulations demonstrate that the 2PT model yields entropy values in semiquantitative agreement with experimental data when using MB-pol and MB-pol(2023). Additionally, our analyses indicate that the entropy values predicted by pairwise-additive water models may benefit from error compensation between the inherent approximations of the 2PT model and the known limitations of these water models in describing many-body interactions. Despite its approximate nature, the simplicity of the 2PT model makes it a valuable tool for estimating relative entropy changes of liquid water across various environments, especially when combined with water models that provide a consistently robust representation of the structural, thermodynamic, and dynamical properties of liquid water.

Tuning MOF/polymer interfacial pore geometry in mixed matrix membrane for upgrading CO2 separation performance

The current paradigm considers the control of the MOF/polymer interface mostly for achieving a good compatibility between the two components to ensure the fabrication of continuous mixed-matrix metal-organic framework (MMMOF) membranes. Here, we unravel that the interfacial pore shape nanostructure plays a key role for an optimum molecular transport. The prototypical ultrasmall pore AlFFIVE-1-Ni MOF was assembled with the polymer PIM-1 to design a composite with gradually expanding pore from the MOF entrance to the MOF/polymer interfacial region. Concentration gradient-driven molecular dynamics simulations demonstrated that this pore nanostructuring enables an optimum guided path for the gas molecules at the MOF/polymer interface that decisively leads to an acceleration of the molecular transport all along the MMMOF membrane. This numerical prediction resulted in the successful fabrication of a [001]-oriented nanosheets AlFFIVE-1-Ni/PIM-1 MMMOF membrane exhibiting an excellent CO2 permeability, better than many MMMs, and ideally associated with a sufficiently high CO2/CH4 selectivity that makes this membrane very promising for natural gas/biogas purification.

Positive and Negative Impacts of Interfacial Hydrogen Bonds on Photocatalytic Hydrogen Evolution

Understanding the behavior of water molecules at solid-liquid interfaces is crucial for various applications such as photocatalytic water splitting, a key technology for sustainable fuel production and chemical transformations. Despite extensive studies conducted in the past, the impact of the microscopic structure of interfacial water molecules on photocatalytic reactivity has not been directly examined. In this study, using real-time mass spectrometry and Fourier-transform infrared spectroscopy, we demonstrated the crucial role of hydrogen bond (H-bond) networks on the photocatalytic hydrogen evolution in thickness-controlled water adsorption layers on various TiO2 photocatalysts. Under controlled water vapor environments with relative humidity (RH) below 70%, we observed a monotonic increase in the H2 formation rate with increasing RH, indicating that reactive water molecules were present not only in the first adsorbed layer but also in several overlying layers. In contrast, at RH > 70%, when more than three water layers covered the catalyst surface, the H2 formation rate turned to decrease dramatically because of the structural rearrangement and hardening of the interfacial H-bond network induced during further water adsorption. This unique many-body effect of interfacial water was consistently observed for various TiO2 particles with different crystalline structures, including brookite, anatase, and a mixture of anatase and rutile. Our results demonstrated that depositing several water layers in a water vapor environment with RH ∼ 70% is optimal for photocatalytic hydrogen evolution rather than liquid-phase reaction conditions in aqueous solutions. This study provides molecular-level insights into designing interfacial water conditions to enhance photocatalytic performance.

Probing the structure of water in individual living cells

Water regulates or even governs a wide range of biological processes. Despite its fundamental importance, surprisingly little is known about the structure of intracellular water. Herein we employ a Raman micro-spectroscopy technique to uncover the composition, abundance and vibrational spectra of intracellular water in individual living cells. In three different cell types, we show a small but consistent population (~3%) of non-bulk-like water. It exhibits a weakened hydrogen-bonded network and a more disordered tetrahedral structure. We attribute this population to biointerfacial water located in the vicinity of biomolecules. Moreover, our whole-cell modeling suggests that all soluble (globular) proteins inside cells are surrounded by, on average, one full molecular layer (about 2.6 Angstrom) of biointerfacial water. Furthermore, relative invariance of biointerfacial water is observed among different single cells. Overall, our study not only opens up experimental possibilities of interrogating water structure in vivo but also provides insights into water in life.

The Dependence of Hydrophobic Interactions on the Shape of Solute Surface

According to our recent studies on hydrophobicity, this work is aimed at understanding the dependence of hydrophobic interactions on the shape of a solute's surface. It has been observed that dissolved solutes primarily affect the structure of interfacial water, which refers to the top layer of water at the interface between the solute and water. As solutes aggregate in a solution, hydrophobic interactions become closely related to the transition of water molecules from the interfacial region to the bulk water. It is inferred that hydrophobic interactions may depend on the shape of the solute surface. To enhance the strength of hydrophobic interactions, the solutes tend to aggregate, thereby minimizing their surface area-to-volume ratio. This also suggests that hydrophobic interactions may exhibit directional characteristics. Moreover, this phenomenon can be supported by calculated potential mean forces (PMFs) using molecular dynamics (MD) simulations, where different surfaces, such as convex, flat, or concave, are associated with a sphere. Furthermore, this concept can be extended to comprehend the molecular packing parameter, commonly utilized in studying the self-assembly behavior of amphiphilic molecules in aqueous solutions.

Nuclear quantum effects on the vibrational dynamics of the water-air interface

We have applied path-integral molecular dynamics simulations to investigate the impact of nuclear quantum effects on the vibrational dynamics of water molecules at the water-air interface. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time series analysis, while the time scales of vibrational spectral diffusion are determined from frequency-time correlation functions and joint probability distributions. We find that the inclusion of nuclear quantum effects leads not only to a redshift in the vibrational frequency distribution by about 120 cm-1 for both the bulk and interfacial water molecules but also to an acceleration of the vibrational dynamics at the water-air interface by as much as 35%. In addition, a blueshift of about 45 cm-1 is seen in the vibrational frequency distribution of interfacial water molecules compared to that of the bulk. Furthermore, the dynamics of water molecules beyond the topmost molecular layer was found to be rather similar to that of bulk water.

Monitoring water harvesting in metal-organic frameworks, one water molecule at a time

Metal-organic frameworks (MOFs) have gained prominence as potential materials for atmospheric water harvesting, a vital solution for arid regions and areas experiencing severe water shortages. However, the molecular factors influencing the performance of MOFs in capturing water from the air remain elusive. Among all MOFs, Ni2X2BTDD (X = F, Cl, Br) stands out as a promising water harvester due to its ability to adsorb substantial amounts of water at low relative humidity (RH). Here, we use advanced molecular dynamics simulations carried out with the state-of-the-art MB-pol data-driven many-body potential to monitor water adsorption in the three Ni2X2BTDD variants as a function of RH. Our simulations reveal that the type of halide atom in the three Ni2X2BTDD frameworks significantly influences the corresponding molecular mechanisms of water adsorption: while water molecules form strong hydrogen bonds with the fluoride atoms in Ni2F2BTDD, they tend to form hydrogen bonds with the nitrogen atoms of the triazolate linkers in Ni2Cl2BTDD and Ni2Br2BTDD. Importantly, the large size of the bromide atoms reduces the void volume in the Ni2Br2BTDD pores, which enable water molecules to initiate an extended hydrogen-bond network at lower RH. These findings not only underscore the prospect for precisely tuning structural and chemical modifications of the frameworks to optimize their interaction with water, but also highlight the predictive power of simulations with the MB-pol data-driven many-body potential. By providing a realistic description of water under different thermodynamic conditions and environments, these simulations yield unique, molecular-level insights that can guide the design and optimization of energy-efficient water harvesting materials.

Structure and dynamics of liquid water from ab initio simulations: adding Minnesota density functionals to Jacob's ladder

The accurate representation of the structural and dynamical properties of water is essential for simulating the unique behavior of this ubiquitous solvent. Here we assess the current status of describing liquid water using ab initio molecular dynamics, with a special focus on the performance of all the later generation Minnesota functionals. Findings are contextualized within the current knowledge on DFT for describing bulk water under ambient conditions and compared to experimental data. We find that, contrary to the prevalent idea that local and semilocal functionals overstructure water and underestimate dynamical properties, M06-L, revM06-L, and M11-L understructure water, while MN12-L and MN15-L overdistance water molecules due to weak cohesive effects. This can be attributed to a weakening of the hydrogen bond network, which leads to dynamical fingerprints that are over fast. While most of the hybrid Minnesota functionals (M06, M08-HX, M08-SO, M11, MN12-SX, and MN15) also yield understructured water, their dynamical properties generally improve over their semilocal counterparts. It emerges that exact exchange is a crucial component for accurately describing hydrogen bonds, which ultimately leads to corrections in both the dynamical and structural properties. However, an excessive amount of exact exchange strengthens hydrogen bonds and causes overstructuring and slow dynamics (M06-HF). As a compromise, M06-2X is the best performing Minnesota functional for water, and its D3 corrected variant shows very good structural agreement. From previous studies considering nuclear quantum effects (NQEs), the hybrid revPBE0-D3, and the rung-5 RPA (RPA@PBE) have been identified as the only two approximations that closely agree with experiments. Our results suggest that the M06-2X(-D3) functionals have the potential to further improve the reproduction of experimental properties when incorporating NQEs through path integral approaches. This work provides further proof that accurate modeling of water interactions requires the inclusion of both exact exchange and balanced (non-local) correlation, highlighting the need for higher rungs on Jacob's ladder to achieve predictive simulations of complex biological systems in aqueous environments.

Unified picture of vibrational relaxation of OH stretch at the air/water interface

The elucidation of the energy dissipation process is crucial for understanding various phenomena occurring in nature. Yet, the vibrational relaxation and its timescale at the water interface, where the hydrogen-bonding network is truncated, are not well understood and are still under debate. In the present study, we focus on the OH stretch of interfacial water at the air/water interface and investigate its vibrational relaxation by femtosecond time-resolved, heterodyne-detected vibrational sum-frequency generation (TR-HD-VSFG) spectroscopy. The temporal change of the vibrationally excited hydrogen-bonded (HB) OH stretch band (ν=1→2 transition) is measured, enabling us to determine reliable vibrational relaxation (T1) time. The T1 times obtained with direct excitations of HB OH stretch are 0.2-0.4 ps, which are similar to the T1 time in bulk water and do not noticeably change with the excitation frequency. It suggests that vibrational relaxation of the interfacial HB OH proceeds predominantly with the intramolecular relaxation mechanism as in the case of bulk water. The delayed rise and following decay of the excited-state HB OH band are observed with excitation of free OH stretch, indicating conversion from excited free OH to excited HB OH (~0.9 ps) followed by relaxation to low-frequency vibrations (~0.3 ps). This study provides a complete set of the T1 time of the interfacial OH stretch and presents a unified picture of its vibrational relaxation at the air/water interface.

Raman and IR spectra of water under graphene nanoconfinement at ambient and extreme pressure-temperature conditions: a first-principles study

The nanoconfinement of water can result in dramatic differences in its physical and chemical properties compared to bulk water. However, a detailed molecular-level understanding of these properties is still lacking. Vibrational spectroscopy, such as Raman and infrared, is a popular experimental tool for studying the structure and dynamics of water, and is often complemented by atomistic simulations to interpret experimental spectra, but there have been few theoretical spectroscopy studies of nanoconfined water using first-principles methods at ambient conditions, let alone under extreme pressure-temperature conditions. Here, we compute the Raman and IR spectra of water nanoconfined by graphene at ambient and extreme pressure-temperature conditions using ab initio simulations. Our results revealed alterations in the Raman stretching and low-frequency bands due to the graphene confinement. We also found spectroscopic evidence indicating that nanoconfinement considerably changes the tetrahedral hydrogen bond network, which is typically found in bulk water. Furthermore, we observed an unusual bending band in the Raman spectrum at ∼10 GPa and 1000 K, which is attributed to the unique molecular structure of confined ionic water. Additionally, we found that at ∼20 GPa and 1000 K, confined water transformed into a superionic fluid, making it challenging to identify the IR stretching band. Finally, we computed the ionic conductivity of confined water in the ionic and superionic phases. Our results highlight the efficacy of Raman and IR spectroscopy in studying the structure and dynamics of nanoconfined water in a large pressure-temperature range. Our predicted Raman and IR spectra can serve as a valuable guide for future experiments.

Hydrophobic hydration of the hydrocarbon adamantane in amorphous ice

Hydrophobic molecules are by definition difficult to hydrate. Previous studies in the area of hydrophobic hydration have therefore often relied on using amphiphilic molecules where the hydrophilic part of a molecule enabled the solubility in liquid water. Here, we show that the hydrophobic adamantane (C10H16) molecule can be fully hydrated through vapour codeposition with water onto a cryogenic substrate at 80 K resulting in the matrix isolation of adamantane in amorphous ice. Using neutron diffraction in combination with the isotopic substitution method and the empirical potential structure refinement technique, we find that the first hydration shell of adamantane is well structured consisting of a hydrogen-bonded cage of 28 water molecules that is also found in cubic structure II clathrate hydrates. The four hexagonal faces of the 51264 cage are situated above the four methine (CH) groups of adamantane whereas the methylene (CH2) groups are positioned below the edges of two adjoining pentagonal faces. The oxygen atoms of the 28 water molecules can be categorised on the basis of symmetry equivalences as twelve A, twelve B and four C oxygens. The water molecules of the first hydration shell display orientations consistent with those expected for a clathrate-hydrate-type cage, but also unfavourable ones with respect to the hydrogen bonding between the water molecules. Annealing the samples at 140 K, which is just below the crystallisation temperature of the matrix, removes the unfavourable orientations and leads to a slight increase in the structural order of the first hydration shell. The very closest water molecules display a tendency for their dipole moments to point towards the adamantane which is attributed to steric effects. Other than this, no significant polarisation effects are observed which is consistent with weak interactions between adamantane and the amorphous ice matrix. FT-IR spectroscopy shows that the incorporation of adamantane into amorphous ice leads to a weakening of the hydrogen bonds. In summary, the matrix-isolation of the highly symmetric adamantane in amorphous ice provides an interesting test case for hydrophobic hydration. Studying the structure and spectroscopic properties of water at the interface with hydrophobic hydrocarbons is also relevant for astrophysical environments, such as comets or the interstellar medium, where amorphous ice and hydrocarbons have been shown to coexist in large quantities.

Raman Spectra of Electrified Si-Water Interfaces: First-Principles Simulations

We investigate the Raman spectra of liquid water in contact with a semiconductor surface using first-principles molecular dynamics simulations. We focus on a hydrogenated silicon-water interface and compute the Raman spectra from time correlation functions of the polarizability. We establish a relationship between Raman spectral signatures and structural properties of the liquid at the interface, and we identify the vibrational impacts of an applied electric field. We show that negative bias leads to a reduction of the number of hydrogen bonds (HBs) formed between the surface and the topmost water layer and an enhancement of the HB interactions between water molecules. Instead, positive bias leads to an enhancement of both the HB interactions between water and the surface and between water molecules, creating a semi-ordered interfacial layer. Our work provides molecular-level insights into electrified semiconductor/water interfaces and the identification of specific structural features through Raman spectroscopy.

Predicting the Raman Spectra of Liquid Water with a Monomer-Field Model

The Raman spectrum of liquid water is quite complex, reflecting its strong sensitivity to the local environment of the individual waters. The OH-stretch region of the spectrum, which captures the influence of hydrogen bonding, has only just begun to be unraveled. Here we develop a model for predicting the Raman spectra of the OH-stretch region by considering how local electric fields distort the energy surface of each water monomer. We find that our model is capable of reproducing the bimodal nature of the main peak, with the shoulder at 3250 cm-1 resulting almost entirely from Fermi resonance. Furthermore, we capture the temperature and polarization dependence of the shoulder, which has proven to be difficult to obtain with previous methods, and analyze the origin of this dependence. We expect our model to be generally useful for understanding and predicting how Raman spectra change under different conditions and with different probe reporters beyond water.

Theoretical Study of the Two-Dimensional Vibrational Sum Frequency Generation Spectroscopy of the Air-Water Interface at Varying Temperature and Its Connections to the Interfacial Structure and Dynamics

We performed a theoretical study of the temperature variation of two-dimensional vibrational sum frequency generation (2D-VSFG) spectra of the OH stretch modes at air-water interfaces in the mid-IR region. The calculations are performed at four different temperatures from 250 to 325 K by using a combination of techniques involving response function formalism of nonlinear spectroscopy, electronic structure calculations, and molecular dynamics simulations. Also, the calculations are performed for isotopically dilute solutions so that the intra- and intermolecular coupling between the vibrational modes of interest can be ignored. We have established the connections of temperature variation of various frequency- and time-dependent features of the calculated spectra to the changes in the underlying structure and dynamics of the interfaces. The results reveal that interfacial water is dynamically more heterogeneous than bulk water, with three dominant dynamical processes exhibiting their corresponding time-dependent features in the 2D-VSFG spectrum. These are the spectral diffusion of hydrogen-bonded OH groups at the interface, conversion of an initially hydrogen-bonded OH group to a dangling OH which is a stable state for surface water, unlike the bulk water, and the third one, which involves the conversion of an initially free or dangling OH group to its hydrogen-bonded state at the interface. The temporal appearance of the cross peaks corresponding to interconversion of the hydrogen-bonded state to the dangling state or vice versa of an interfacial OH group is found to take place at a slower rate than the dynamics of spectral diffusion of hydrogen-bonded molecules at the interface, which, in turn, is slower than the corresponding spectral diffusion of bulk water molecules. The temperature variation of these dynamic processes can be linked to the decay of appropriate hydrogen-bond and non-hydrogen-bond time correlation functions of interfacial water molecules for the different air-water systems studied in this work.

Anion Recognition in Solution: Insights from Thermodynamics and Ultrafast Structural Dynamics

Anion recognition through noncovalent interactions stands as an emerging field in supramolecular chemistry, exerting a profound influence on the regulation of biological functions. Herein, the thermodynamics of complexation between sodium cyanate (NaOCN) and calix[4]pyrrole was systematically investigated by linear and nonlinear IR spectroscopy, highlighting enthalpy changes as the dominant driving force. The overall orientational relaxation of bound anion can be described by an Arrhenius-type activated process, yielding an activation energy of 15.0 ± 1.0 kJ mol-1. The structural dynamics of contact ion pairs (CIPs) formed between Na+ and OCN- in solution showed a negligible temperature effect, suggesting entropy changes as the principal governing factor. Further analysis revealed that anion recognition in solution is mediated by conformational changes of the receptor and collective rearrangement of hydrogen bond dynamics. This study, framed within the paradigms of thermodynamics and ultrafast structural dynamics, substantially advances our comprehension of the microscopic mechanisms underlying anion recognition in the realm of supramolecular chemistry.

Simulating the isotropic Raman spectra of O-H stretching mode in liquid H2O based on a machine learning potential: the influence of vibrational couplings

In this study, we trained a deep potential (DP) for H2O, an accurate machine learning (ML) potential. We performed molecular dynamics (MD) simulations of liquid water using the DP model (or DeePMD simulations). Our results showed that the DP model exhibits DFT-level accuracy, and the DeePMD simulation is a promising approach for modeling the structural properties of liquid water. Based on the DeePMD simulation trajectories, we calculated the isotropic Raman spectra of the O-H stretching mode using the surface-specific velocity-velocity correlation function (ssVVCF), showing that the DeePMD/ssVVCF approach can correctly capture the bimodal characteristics of the experimental Raman spectra, with one peak located near 3400 cm-1 and the other near 3250 cm-1. The success of the DeePMD/ssVVCF approach should be credited to (1) the DFT-level accuracy of the DP model for H2O, (2) the ssVVCF formulation considering the coupling between vibrational modes, and (3) non-Condon effects. Furthermore, the DeePMD simulations revealed that the anharmonic interactions between the coupled water molecules in the first and second hydration shells should play an essential role in the strong mixing of the H-O-H bending mode and the O-H stretching mode, leading to the delocalization of the O-H stretching band. In particular, increasing the strength of hydrogen bonds would enhance the bend-stretch coupling, leading to the red-shifting of the O-H vibrational spectra and the increase in the intensity of the shoulder around 3250 cm-1.

Elucidating the Competitive Adsorption of H2O and CO2 in CALF-20: New Insights for Enhanced Carbon Capture Metal-Organic Frameworks

In light of the pressing need for efficient carbon capture solutions, our study investigates the simultaneous adsorption of water (H2O) and carbon dioxide (CO2) as a function of relative humidity in CALF-20, a highly scalable and stable metal-organic framework (MOF). Advanced computer simulations reveal that due to their similar interactions with the framework, H2O and CO2 molecules compete for the same binding sites, occupying similar void regions within the CALF-20 pores. This competition results in distinct thermodynamic and dynamical behaviors of H2O and CO2 molecules, depending on whether one or both guest species are present. Notably, the presence of CO2 molecules forces the H2O molecules to form more connected hydrogen-bond networks within smaller regions, slowing water reorientation dynamics and decreasing water entropy. Conversely, the presence of water speeds up the reorientation of CO2 molecules, decreases the CO2 entropy, and increases the propensity for CO2 to be adsorbed within the framework due to stronger water-mediated interactions. Due to the competition for the same void spaces, both H2O and CO2 molecules exhibit slower diffusion when molecules of the other guest species are present. These findings offer valuable strategies and insights into enhancing the differential affinity of H2O and CO2 for MOFs specifically designed for carbon capture applications.

Water bend-libration as a cellular Raman imaging probe of hydration

Water is a ubiquitous and vital component of living systems. Hydration, which is the interaction between water and intracellular biomolecules, plays an important role in cellular processes. However, it is technically challenging to study water structure within cells directly. Here, we demonstrate the utility and power of the water bend-libration combination band as a unique Raman spectral imaging probe of cellular hydration. Hydration maps reveal distinct water environments within subcellular compartments (e.g., nucleolus and lipid droplet) due to the spectral sensitivity of this coupled vibrational band. Spectroscopic studies using the water bend-libration are broadly applicable, offering the potential to capture the chemical complexity of hydration in numerous systems.

Counteracting Effects of Trimethylamine N-Oxide against Urea in Aqueous Solutions: Insights from Theoretical Two-Dimensional Infrared Spectroscopy

The study of small osmolytes in their aqueous solutions has gained significant attention because of their relevance to structure and thermodynamics of proteins in aqueous media. Special attention has been given to the binary and ternary aqueous solutions of urea and trimethylamine N-oxide (TMAO). Urea is a well-known protein denaturant, while TMAO protects proteins in their native states. Interestingly, TMAO counteracts urea's ability to denature proteins when present in solutions with approximately half of the concentration of urea. Vibrational spectroscopy can improve our understanding of the molecular origin of this counteracting effect because of its sensitivity to local structure and dynamics. We present results of theoretical linear vibrational and two-dimensional infrared (2DIR) spectroscopy of water in the binary and ternary aqueous solutions of TMAO and urea. The 2DIR spectra are calculated using the electronic structure/molecular dynamics approach. The non-Condon effects in spectral transitions are incorporated in the theoretical calculations of 2DIR spectra. It is found that TMAO disrupts the local structure of water, while urea leaves it essentially unaffected. The 2DIR results show that both TMAO and urea slow down the dynamics of spectral diffusion of water. The extent of slowing down is found to be particularly significant for both hydration and bulk water in the presence of TMAO which can be attributed to strong hydrogen bonds between the water and TMAO molecules. The water molecules present in the hydration layer of the solutes in the ternary solutions are found to relax at even slower rates compared to that in their binary solutions in water. The hydrogen bonds between TMAO and urea are found to be not stable. Thus, the counteracting effect of TMAO against urea is seen to take place mainly through water-mediated interactions. Such TMAO-induced effects giving rise to more structured and slower hydrogen-bonded network are successfully captured through 2DIR spectroscopic calculations.

Exploring the Origins of the Intensity of the OH Stretch-HOH Bend Combination Band in Water

While the intensity of the OH stretching fundamental transition is strongly correlated to hydrogen-bond strength, the intensity of the corresponding transition to the state with one quantum of excitation in both the OH stretching and HOH bending vibrations in the same water molecule shows a much weaker sensitivity to the hydrogen-bonding environment. The origins of this difference are explored through analyses of the contributions of terms in the expansion of the dipole moment to the calculated intensity. It is found that the leading contribution to the stretch-bend intensity involves the second derivative of the dipole moment with respect to the OH bond length and HOH angle. While this is not surprising, the insensitivity of this derivative to the hydrogen-bonding environment is unexpected. Possible contributions of mode mixing are also explored. While mode mixing leads to splittings of the energies of nearly degenerate excited states, it does not result in significant changes in the sum of the intensities of these transitions. Analysis of changes in the partial charges on the hydrogen atoms upon displacement of the HOH angles shows that these charges generally increase with increasing HOH angle. This effect is partially canceled by a decrease in the charge of the hydrogen atom when a hydrogen bond is broken. The extent of this cancellation increases with the hydrogen bond strength, which is reflected in the observed insensitivity of the intensity of the stretch-bend transition to hydrogen-bond strength.

Chain-like Structures Facilitate Li+ Transport in Concentrated Aqueous Electrolytes: Insights from Ultrafast Infrared Spectroscopy and Molecular Dynamics Simulations

Highly concentrated aqueous electrolytes have attracted attention due to their unique applications in lithium ion batteries (LIBs). However, the solvation structure and transport mechanism of Li+ cations at concentrated concentrations remain largely unexplored. To address this gap in knowledge, we employ ultrafast infrared spectroscopy and molecular dynamics (MD) simulations to reveal the dynamic and spatial structural heterogeneity in aqueous lithium chloride (LiCl) solutions. The coupling between the reorientation dynamics of the extrinsic probe and the macroscopic viscosity in aqueous LiCl solutions was analyzed using the Stokes-Einstein-Debye (SED) equations. MD simulations reveal that the Cl- and Li+ form chain-like structures through electrostatic interactions, supporting the vehicular migration of Li+ through the chain-like structure. The concentration dependent conductivity of the LiCl solution is well reproduced, where Li(H2O)2+ and Li(H2O)3+ are the dominant species that contribute to the conduction of Li+. This study is expected to establish correlations between ion pair structures and macroscopic properties.

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.

Anharmonic Assignment of the Water Octamer Spectrum in the OH Stretch Region

We interface the quasi-classical trajectory approach with an ab initio potential energy surface for water to assign the vibrational spectroscopical features of the OH stretch region of the water octamer cluster, which is considered to be a precursor of ice. An attempt by Li et al. to assign their recent reference experiment involved lower-level calculations based on an ad hoc scaled harmonic approach. Differently from the conclusions of this previous assignment, which invoked the contribution of 5 conformers and a solvated form of the water heptamer in the spectrum, we find out that the spectroscopic features can be related to the 4 conformers of the octamer lying lower in energy.

Manipulating hydrogen bond dissociation rates and mechanisms in water dimer through vibrational strong coupling

The vibrational strong coupling (VSC) between molecular vibrations and cavity photon modes has recently emerged as a promising tool for influencing chemical reactivities. Despite numerous experimental and theoretical efforts, the underlying mechanism of VSC effects remains elusive. In this study, we combine state-of-art quantum cavity vibrational self-consistent field/configuration interaction theory (cav-VSCF/VCI), quasi-classical trajectory method, along with the quantum-chemical CCSD(T)-level machine learning potential, to simulate the hydrogen bond dissociation dynamics of water dimer under VSC. We observe that manipulating the light-matter coupling strength and cavity frequencies can either inhibit or accelerate the dissociation rate. Furthermore, we discover that the cavity surprisingly modifies the vibrational dissociation channels, with a pathway involving both water fragments in their ground vibrational states becoming the major channel, which is a minor one when the water dimer is outside the cavity. We elucidate the mechanisms behind these effects by investigating the critical role of the optical cavity in modifying the intramolecular and intermolecular coupling patterns. While our work focuses on single water dimer system, it provides direct and statistically significant evidence of VSC effects on molecular reaction dynamics.

Transient Terahertz Stark Effect: A Dynamic Probe of Electric Interactions in Polar Liquids

Electric forces acting on molecules in liquids at ambient temperature fluctuate at terahertz (THz) frequencies with a direct impact on their electronic and optical properties. We introduce the transient THz Stark effect to modify the electronic absorption spectra of dye molecules and, thus, elucidate and determine the underlying molecular interactions and dynamics. Picosecond electric fields of megavolts/cm induce a nonequilibrium response of the prototypical Betaine-30 in polar solution that is probed via transient absorption changes. The field-induced broadening of the absorption band follows the THz intensity in time, with a minor impact of solvent dynamics. The ground and excited state dipole energies in the THz field govern this response, allowing for a quantification of electric forces in a structurally frozen molecular environment.

Combinational Vibration Modes in H2O/HDO/D2O Mixtures Detected Thanks to the Superior Sensitivity of Femtosecond Stimulated Raman Scattering

Overtones and combinational modes frequently play essential roles in ultrafast vibrational energy relaxation in liquid water. However, these modes are very weak and often overlap with fundamental modes, particularly in isotopologues mixtures. We measured VV and HV Raman spectra of H₂O and D₂O mixtures with femtosecond stimulated Raman scattering (FSRS) and compared the results with calculated spectra. Specifically, we observed the mode at around 1850 cm^-1 and assigned it to H-O-D bend + rocking libration. Second, we found that the H-O-D bend overtone band and the OD stretch + rocking libration combination band contribute to the band located between 2850 and 3050 cm^-1. Furthermore, we assigned the broad band located between 4000 and 4200 cm^-1 to be composed of combinational modes of high-frequency OH stretching modes with predominantly twisting and rocking librations. These results should help in a proper interpretation of Raman spectra of aqueous systems as well as in the identification of vibrational relaxation pathways in isotopically diluted water.

Theoretical Two Dimensional Infrared Spectroscopy of Aqueous Solutions of tert-Butyl Alcohol: Variation of the Dynamics of Spectral Diffusion along the Percolation Transition

Binary mixtures of water and tert-butyl alcohol (TBA) are known to exhibit the so-called percolation transition where small clusters of TBA molecules span into large aggregates beyond a threshold concentration of the alcohol. In the present study, we have investigated the linear and two-dimensional infrared spectral features of aqueous solutions of TBA for varying concentration of the alcohol along the percolation transition. The percolation transition is characterized through calculations of intermolecular radial distribution functions and average size of the largest cluster of TBA molecules. It is found that, with variation of alcohol concentration, the radial distribution functions of the central carbon atoms of TBA molecules show a nonmonotonic change in the height of the first peak and also the size of the largest cluster of TBA molecules show a jump in the increase of its size for TBA mole fraction between 0.04 and 0.06 corresponding to a transition from smaller clusters to larger spanning aggregates. However, it is found that the linear infrared spectrum of water does not exhibit any noticeable changes on variation of TBA concentration along the percolation transition. Subsequently, two-dimensional infrared (2DIR) spectra and vibrational frequency time correlation function of water are calculated for all the TBA-water solutions considered in this study. The spectral diffusion of water calculated from 2DIR is found to slow down with increase of the TBA concentration. The time scales of spectral diffusion of water, as characterized by the relaxation of frequency time correlation function, 2DIR metric of central line slope, and also the hydrogen bond time correlation functions, are found to exhibit a noticeable jump along the percolation transition. The hydrophilic group of TBA is found to retard the water dynamics more effectively than the hydrophobic groups. Also, the jump in the dynamical slowdown along the percolation transition is found to be more significant for water molecules at the hydrophilic sites.

Dynamic Heterogeneity and Kovacs' Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. We used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water's structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

Is Unified Understanding of Vibrational Coupling of Water Possible? Hyper-Raman Measurement and Machine Learning Spectra

The impact of the vibrational coupling of the OH stretch mode on the spectra differs significantly between IR and Raman spectra of water. Unified understanding of the vibrational couplings is not yet achieved. By using a different class of vibrational spectroscopy, hyper-Raman (HR) spectroscopy, together with machine-learning-assisted HR spectra calculation, we examine the impact of the vibrational couplings of water through the comparison of isotopically diluted H₂O and pure H₂O. We found that the isotopic dilution reduces the HR bandwidths, but the impact of the vibrational coupling is smaller than in the IR and parallel-polarized Raman. Machine learning HR spectra indicate that the intermolecular coupling plays a major role in broadening the bandwidth, while the intramolecular coupling is negligibly small, which is consistent with the IR and Raman spectra. Our result clearly demonstrates a limited impact of the intramolecular vibration, independent of the selection rules of vibrational spectroscopies.

Temperature Dependence of Non-Condon Effects in Two-Dimensional Vibrational Spectroscopy of Water

Non-Condon effects in vibrational spectroscopy refers to the dependence of a molecule's vibrational transition dipole and polarizability on the coordinates of the surrounding environment. Earlier studies have shown that such effects can be pronounced for hydrogen-bonded systems like liquid water. Here, we present a theoretical study of two-dimensional vibrational spectroscopy under the non-Condon and Condon approximations at varying temperatures. We have performed calculations of both two-dimensional infrared and two-dimensional vibrational Raman spectra to gain insights into the temperature dependence of non-Condon effects in nonlinear vibrational spectroscopy. The two-dimensional spectra are calculated for the OH vibration of interest in the isotopic dilution limit where the coupling between the oscillators is ignored. Generally, both the infrared and Raman line shapes undergo red shifts with decrease in temperature due to strengthening of hydrogen bonds and decrease in the fraction of OH modes with weaker or no hydrogen bonds. The infrared line shape is further red-shifted under the non-Condon effects at a given temperature, while the Raman line shape does not show any such red shift due to non-Condon effects. The spectral dynamics becomes slower on decrease of temperature due to slower hydrogen bond relaxation and, for a given temperature, the spectral diffusion occurs at a faster rate upon inclusion of non-Condon effects. The time scales of spectral diffusion extracted from different metrics agree well with each other and also with experiments. The changes in the spectrum due to non-Condon effects are found to be more significant at lower temperatures.

Water model for hydrophobic cavities: structure and energy from quantum-chemical calculations

This ab initio study aims to design a series of large water clusters having a hollow clathrate-like cage able to host hydrophobic solutes of various sizes. Starting from the (H₂O)_n (n = 18, 20, 24 and 28) hollow cages, water layers have been added in a stepwise manner in order to model the configuration of water molecules beyond the primary shell. The large (H₂O)₁₀₀, (H₂O)₁₂₀ and (H₂O)₁₄₀ clusters complete the hydrogen bonding network of the cage with optimal and regular tiling of the do-, tetra-decahedron and hexa-decahedron, respectively. This study is corroborated by an investigation of dense water clusters up to the (H₂O)₁₂₃ one, being highly consistent with experimental data on ice concerning the electronic and zero-point energies for aggregate formation at 0 K and enthalpy and entropy at 273 K. The cavity creation profoundly alters the orientation of water molecules compared with those found in dense clusters. Nevertheless, such a large reorganization is necessary to maximize the water-water attraction by making it similar to the one found in dense clusters. The cage formation is an endothermic process; however, the computed values are large compared with previous reports for hydrocarbon aqueous solutions. Larger clusters are required for a more fruitful comparison.

Direct calculation of the temperature dependence of 2D-IR spectra: Urea in water

A method for directly calculating the temperature derivative of two-dimensional infrared (2D-IR) spectra from simulations at a single temperature is presented. The approach is demonstrated by application to the OD stretching spectrum of isotopically dilute aqueous (HOD in H₂O) solutions of urea as a function of concentration. Urea is an important osmolyte because of its ability to denature proteins, which has motivated significant interest in its effect on the structure and dynamics of water. The present results show that the temperature dependence of both the linear IR and 2D-IR spectra, which report on the underlying energetic driving forces, is more sensitive to urea concentration than the spectra themselves. Additional physical insight is provided by calculation of the contributions to the temperature derivative from different interactions, e.g., water-water, water-urea, and urea-urea, present in the system. Finally, it is demonstrated how 2D-IR spectra at other temperatures can be obtained from only room temperature simulations.

Spectroscopy from Machine Learning by Accurately Representing the Atomic Polar Tensor

Vibrational spectroscopy is a key technique to elucidate microscopic structure and dynamics. Without the aid of theoretical approaches, it is, however, often difficult to understand such spectra at a microscopic level. Ab initio molecular dynamics has repeatedly proved to be suitable for this purpose; however, the computational cost can be daunting. Here, the E(3)-equivariant neural network e3nn is used to fit the atomic polar tensor of liquid water a posteriori on top of existing molecular dynamics simulations. Notably, the introduced methodology is general and thus transferable to any other system as well. The target property is most fundamental and gives access to the IR spectrum, and more importantly, it is a highly powerful tool to directly assign IR spectral features to nuclear motion─a connection which has been pursued in the past but only using severe approximations due to the prohibitive computational cost. The herein introduced methodology overcomes this bottleneck. To benchmark the machine learning model, the IR spectrum of liquid water is calculated, indeed showing excellent agreement with the explicit reference calculation. In conclusion, the presented methodology gives a new route to calculate accurate IR spectra from molecular dynamics simulations and will facilitate the understanding of such spectra on a microscopic level.

Low-temperature vibronic spectroscopy of condensed chromophore exhibiting inhomogeneous distribution of vibrational frequencies in a mixed quantum-classical environment

This work has been motivated by the recent paper by the author [M. Toutounji, Phys. Chem. Chem. Phys., 2021, 23, 21981] whereby a mixed quantum-classical Liouville equation was used to probe the spectroscopy and dynamics of a spin-boson system. A mixed quantum-classical Liouville equation treats the system of interest quantum mechanically, the bath classically, and the coupling term mixed quantum-classical mechanically. This paper offers a two-fold advantage: correcting the treatment of the electronic transition decay (width in frequency domain) and assessing the local heterogeneous vibrational structure. The homogeneous linear absorption spectrum of a chromophore embedded in a mixed quantum-classical environment at low temperature is composed of a sharp peak called a zero-phonon line (ZPL) and a broad phonon sideband (PSB), whereby the ZPL and the PSB are assimilated by a Lorentzian function and Voigt profiles, respectively. The PSB, in this case, is characterized by a local heterogeneous structure due to a dispersive medium of vibrations, modeled by vibrational Gaussian distributions to represent the arising inhomogeneous broadening and Lorentzians to model the homogeneous vibrations. This description seems to model proteins and amorphous solids exhibiting a local heterogeneous structure as both electronic and vibrational inhomogeneous broadening seems to be large in these media. This work provides a derivation of linear absorption lineshape and vibronic transition dipole moment time correlation functions, both of which account for pure electronic dephasing (ZPL width) the Voigt profile description of the phonon profiles (PSB) in dispersive media.

Non-negligible Axial Ligand Effect on Electrocatalytic CO2 Reduction with Iron Porphyrin Complexes

Iron(III) porphyrin complexes have been demonstrated as one of the efficient molecular catalysts for the electrochemical reduction of CO₂. However, the role of axial ligands coordinated with a metal center in the complex on the electrochemical CO₂ reduction activity has not been fully explored yet. Herein, iron(III) tetraphenylporphyrin thiocyanate (FeTPP-SCN) is synthesized from a commercially available catalyst of FeTPP-Cl by a counteranion exchanging reaction. Cyclic voltammetry measurements showed that the catalytic activity of FeTPP-SCN is noticeably suppressed in the DMF solutions. The structural dynamics of the axial ligand in FeTPP-SCN are further examined by the FTIR and ultrafast IR spectroscopies, where the SCN ligand is employed as the local vibrational probe. Vibrational relaxation measurements showed that the reorientational dynamics of SCN ligands was strongly restricted in DMF solution, suggesting that the subtle electrostatic interaction between the ligands and metal center in the complex can have a non-negligible effect on its catalytic activity.

Observation of Strong Synergy in the Interfacial Water Response of Binary Ionic and Nonionic Surfactant Mixtures

Interfacial vibrational footprints of the binary mixture of sodium dodecyl sulfate (SDS) and hexaethylene glycol monododecyl ether (C₁₂E₆) were probed using heterodyne detected vibrational sum frequency generation (HDVSFG). Our results show that in the presence of C₁₂E₆ at CMC (70 μM) the effect of SDS on the orientation of interfacial water molecules is enhanced 10 times compared to just pure surfactants. The experimental results contest the traditional Langmuir adsorption model predictions. This is also evidenced by our molecular dynamics simulations that show a remarkable restructuring and enhanced orientation of the interfacial water molecules upon DS^- adsorption to the C₁₂E₆ surface. The simulations show that the adsorption free energy of DS^- ions to a water surface covered with C₁₂E₆ is an enthalpy-driven process and more attractive by ∼10 k_BT compared to the adsorption energy of DS^- to the surface of pure water.

Dynamic Community Detection Decouples Multiple Time Scale Behavior of Complex Chemical Systems

Although community or cluster identification is becoming a standard tool within the simulation community, traditional algorithms are challenging to adapt to time-dependent data. Here, we introduce temporal community identification using the Δ-screening algorithm, which has the flexibility to account for varying community compositions, merging and splitting behaviors within dynamically evolving chemical networks. When applied to a complex chemical system whose varying chemical environments cause multiple time scale behavior, Δ-screening is able to resolve the multiple time scales of temporal communities. This computationally efficient algorithm is easily adapted to a wide range of dynamic chemical systems; flexibility in implementation allows the user to increase or decrease the resolution of temporal features by controlling parameters associated with community composition and fluctuations therein.

Hydrogen-Bond Configurations of Hydration Water around Glycerol Investigated by HOH Bending and OH Stretching Analysis

Toward a comprehensive understanding of the mechanism of glycerol as a moisturizer, studies on the hydrogen-bond (HB) structure of hydration water, which is known to be disordered by glycerol, are insufficient. To this aim, we evaluated the HB configurations based on the HOH bending and OH stretching spectra of the hydration water from those of glycerol/water mixtures by subtracting the contributions of bulk water and glycerol using dielectric relaxation spectroscopy. Analysis of the HOH bending band showed that hydration water-donating HBs lose the intermolecular bending coupling with increasing glycerol by replacing the water-water HBs with water-glycerol HBs. The OH stretching band provided more detailed insight into the HB configuration, indicating that the double-donor double-acceptor and double-donor single-acceptor configurations in bulk water change to a predominantly double-donor single-acceptor configuration in hydration water around glycerol. The formation of more donor HBs than acceptor HBs may be due to the steric constrains by glycerol and/or differences in the partial charge on the oxygen atom between water and glycerol.

Van der Waals interactions regulating the hydration of 2-methacryloyloxyethyl phosphorylcholine, the constructing monomer of biocompatible polymers

Van der Waals (VDW) interactions provide fantastic properties for biological systems that function at room temperature. The VDW interaction, which primarily contributes to weak hydrogen bonding, is expected to play a key role in regulating hydrophobic hydration to express the biologically inert biocompatible function of polymerized MPCs (2-methacryloyloxyethyl phosphorylcholine). This report explores at the molecular level the biologically inert function of polymerized MPCs through an array of vibrational spectroscopic and computational characterization of MPC monomers, as temperature-dependent change of intramolecular weak hydrogen bonding. Synchrotron Fourier transform infrared microspectroscopy and terahertz time-domain spectroscopy were used to investigate temperature-dependent spectral changes in the low frequency vibrations of the MPC over the temperature range from cryogenic to room temperature, and the results were analysed by highly reliable well-established density functional theory (DFT) calculations. Complicated spectral features in the low frequency energy region and the uncertain conformations of the MPC in the amorphous powder state are clearly resolved under a polarizable continuum model and dispersion correction to pure DFT calculations.

Coupled Cluster Molecular Dynamics of Condensed Phase Systems Enabled by Machine Learning Potentials: Liquid Water Benchmark

Coupled cluster theory is a general and systematic electronic structure method, but in particular the highly accurate "gold standard" coupled cluster singles, doubles and perturbative triples, CCSD(T), can only be applied to small systems. To overcome this limitation, we introduce a framework to transfer CCSD(T) accuracy of finite molecular clusters to extended condensed phase systems using a high-dimensional neural network potential. This approach, which is automated, allows one to perform high-quality coupled cluster molecular dynamics, CCMD, as we demonstrate for liquid water including nuclear quantum effects. The machine learning strategy is very efficient, generic, can be systematically improved, and is applicable to a variety of complex systems.

Isotope effects observed in diluted D2O/H2O mixtures identify HOD-induced low-density structures in D2O but not H2O

Normal and heavy water are solvents most commonly used to study the isotope effect. The isotope effect of a solvent significantly influences the behavior of a single molecule in a solution, especially when there are interactions between the solvent and the solute. The influence of the isotope effect becomes more significant in D₂O/H₂O since the hydrogen bond in H₂O is slightly weaker than its counterpart (deuterium bond) in D₂O. Herein, we characterize the isotope effect in a mixture of normal and heavy water on the solvation of a HOD molecule. We show that the HOD molecule affects the proximal solvent molecules, and these disturbances are much more significant in heavy water than in normal water. Moreover, in D₂O, we observe the formation of low-density structures indicative of an ordering of the solvent around the HOD molecule. The qualitative differences between HOD interaction with D₂O and H₂O were consistently confirmed with Raman spectroscopy and NMR diffusometry.

Water Structure in the Electrical Double Layer and the Contributions to the Total Interfacial Potential at Different Surface Charge Densities

The electric double layer governs the processes of all charged surfaces in aqueous solutions; however, elucidating the structure of the water molecules is challenging for even the most advanced spectroscopic techniques. Here, we present the individual Stern layer and diffuse layer OH stretching spectra at the silica/water interface in the presence of NaCl over a wide pH range using a combination of vibrational sum frequency generation spectroscopy, heterodyned second harmonic generation, and streaming potential measurements. We find that the Stern layer water molecules and diffuse layer water molecules respond differently to pH changes: unlike the diffuse layer, whose water molecules remain net-oriented in one direction, water molecules in the Stern layer flip their net orientation as the solution pH is reduced from basic to acidic. We obtain an experimental estimate of the non-Gouy-Chapman (Stern) potential contribution to the total potential drop across the insulator/electrolyte interface and discuss it in the context of dipolar, quadrupolar, and higher order potential contributions that vary with the observed changes in the net orientation of water in the Stern layer. Our findings show that a purely Gouy-Chapman (Stern) view is insufficient to accurately describe the electrical double layer of aqueous interfaces.

Probing electrolyte–silica interactions through simulations of the infrared spectroscopy of nanoscale pores

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D2O confined in hydroxylated amorphous silica slit pores of width 1–6 nm and pH [Formula: see text]. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump–probe anisotropy and two-dimensional IR measurements, are investigated. The aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Solvent effects on the NMR shieldings of stacked DNA base pairs

Stacking effects are among the most important effects in DNA. We have recently studied their influence in fragments of DNA through the analysis of NMR magnetic shieldings, firstly in vacuo. As a continuation of this line of research we show here the influence of solvent effects on the shieldings through the application of both explicit and implicit models. We found that the explicit solvent model is more appropriate for consideration due to the results matching better in general with experiments, as well as providing clear knowledge of the electronic origin of the value of the shieldings. Our study is grounded on a recently developed theoretical model of our own, by which we are able to learn about the magnetic effects of given fragments of DNA molecules on selected base pairs. We use the shieldings of the atoms of a central base pair (guanine-cytosine) of a selected fragment of DNA molecules as descriptors of physical effects, like π-stacking and solvent effects. They can be taken separately and altogether. The effect of π-stacking is introduced through the addition of some pairs above and below of the central base pair, and now, the solvent effect is considered including a network of water molecules that consist of two solvation layers, which were fixed in the calculations performed in all fragments. We show that the solvent effects enhance the stacking effects on the magnetic shieldings of atoms that belong to the external N-H bonds. The net effect is of deshielding on both atoms. There is also a deshielding effect on the carbon atoms that belong to CO bonds, for which the oxygen atom has an explicit hydrogen bond (HB) with a solvent water molecule. Solvent effects are found to be no higher than a few percent of the total value of the shieldings (between 1% and 5%) for most atoms, although there are few for which such an effect can be higher. There is one nitrogen atom, the acceptor of the HB between guanine and cytosine, that is more highly shielded (around 15 ppm or 10%) when the explicit solvent is considered. In a similar manner, the most external nitrogen atom of cytosine and the hydrogen atom that is bonded to it are highly deshielded (around 10 ppm for nitrogen and around 3 ppm for hydrogen).

From Networked to Isolated: Observing Water Hydrogen Bonds in Concentrated Electrolytes with Two-Dimensional Infrared Spectroscopy

Superconcentrated electrolytes have emerged as a promising class of materials for energy storage devices, with evidence that high voltage performance is possible even with water as the solvent. Here, we study the changes in the water hydrogen bonding network induced by the dissolution of lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) in concentrations ranging from the dilute to the superconcentrated regimes. Using time-resolved two-dimensional infrared spectroscopy, we observe the progressive disruption of the water-water hydrogen bond network and the appearance of isolated water molecules interacting only with ions, which can be identified and spectroscopically isolated through the intermolecular cross-peaks between the water and the TFSI^- ions. Analyzing the vibrational relaxation of excitations of the H₂O stretching mode, we observe a transition in the dominant relaxation path as the bulk-like water vanishes and is replaced by ion-solvation water with the rapid single-step relaxation of delocalized stretching vibrations into the low frequency modes being replaced by multistep relaxation through the intramolecular H₂O bend and into the TFSI^- high frequency modes prior to relaxing to the low frequency structural degrees of freedom. These results definitively demonstrate the absence of vibrationally bulk-like water in the presence of high concentrations of LiTFSI and especially in the superconcentrated regime, while additionally revealing aspects of the water hydrogen bond network that have been difficult to discern from the vibrational spectroscopy of the neat liquid.