Librational Band Research Articles

Distinguishing ice β-XV from deep glassy ice VI: Raman spectroscopy.

The nature of the hydrogen sublattice of an HCl-doped ice VI sample after cooling at 1.8 GPa has been a topic of recent interest. The samples are interpreted either as the new H-ordered ice phase ice β-XV with a thermodynamic stability region in the phase diagram [T. M. Gasser et al., Chem. Sci., 2018, 9, 4224], or alternatively as H-disordered, deep glassy ice VI [A. Rosu-Finsen and C. G. Salzmann, Chem. Sci., 2019, 10, 515]. Here we provide a comprehensive Raman spectroscopic study on ice β-XV, ice XV and ice VI, with the following key findings: (i) the Raman spectra of ice β-XV differ fundamentally from those of ice VI and ice XV, where the degree of H-order is even higher than in ice XV. (ii) Upon cooling ice VI there is competition between formation of ice XV and ice β-XV domains, where ice XV forms at 0.0 GPa, but ice β-XV at 1.8 GPa. Domains of ice β-XV are present in literature "ice XV" at 1.0 GPa. This result clarifies the puzzling earlier observation that the degree of H-order in ice XV apparently improves upon heating and recooling at ambient pressure. In reality, this procedure leaves the H-order in ice XV unaffected, but removes ice β-XV domains by transforming them to ice XV. (iii) Upon heating, the samples experience the transition sequence ice β-XV → ice XV → ice VI, i.e., an order-order transition at 103 K followed by an order-disorder transition at 129 K. The former progresses via a disordered transient state. (iv) D2O ice β-XV forms upon cooling DCl-doped D2O-ice VI, albeit at a much lower pace than in the hydrogenated case so that untransformed D2O ice VI domains are present even after slow cooling. The librational band at 380 cm-1 is identified to be the characteristic spectroscopic feature of deuterated ice β-XV. Taken together these findings clarify open questions in previous work on H-ordering in the ice VI lattice, rule out a glassy nature of ice β-XV and pave the way for a future neutron diffraction study to refine the crystal structure of D2O ice β-XV.

Ab initio spectroscopy of water under electric fields.

Whereas a broad range of literature exists on the spectroscopy of water in disparate conditions, infrared (IR) and Raman spectra of water subjected to electric fields have never extensively been investigated so far. Based on ab initio molecular dynamics simulations, here we present IR and Raman spectra of bulk liquid water under the effect of static electric fields. A contraction of the entire frequency range is recorded upon increasing the field intensity both in the IR and in the Raman spectra. Whilst the OH stretching band is progressively shifted toward lower frequencies - indicating a field-induced strengthening of the H-bond network - all the other bands are up-shifted by the field. Furthermore, an evident modification of the librational mode band appears in all the spectra. Finally, the order-maker action of the field emerges also from the increase of the water orientational tetrahedral order. Upon field exposure, the water structure becomes more "ice like".

Terahertz VRT Spectroscopy of the Water Hexamer-h12 Cage: Dramatic Libration-Induced Enhancement of Hydrogen Bond Tunneling Dynamics.

We report the assignment and analysis of 176 transitions belonging to a librational band of the (H2O)6 cage isomer near 525 cm-1(15 THz). From a fit of the transitions to an asymmetric top model, we observe both dramatic changes in the rotational constants relative to the ground state, indicating significant nonrigidity, and striking enhancement in the tunneling motions that break and reform the hydrogen bonds in the cluster. This is the fifth water cluster system to display such an enhancement in the 15 THz librational region, the details of which may help to elucidate the hydrogen bond dynamics occurring in bulk liquid water.

Hyper‐Raman investigation of intermolecular vibrations of water and ice

AbstractWe have constructed a hyper‐Raman optimized spectrometer using a high repetition rate (200 kHz) picosecond‐pulsed Nd:YAG laser. This system has been applied to measure liquid water and ice (H2O and D2O) to obtain their hyper‐Raman spectra with high signal‐to‐noise ratios. Unique vibrational spectra of water and ice, distinct either from infrared absorption and Raman scattering, have been obtained. The intermolecular bands in the low‐wavenumber region are observed in the hyper‐Raman spectra with intensities higher than those of the intramolecular ones. They show H/D isotope shifts fully consistent with their assignments to the librational and translational motions of assembled water. The librational bands are found to shift to higher wavenumber on going from liquid water to ice, indicating a more confining environment in ice than in water. Hyper‐Raman spectroscopy will shed new light, complementary to the traditional infrared and Raman vibrational spectroscopies, on intermolecular interactions in liquid water and ice.

Terahertz-infrared spectroscopy of Shewanella oneidensis MR-1 extracellular matrix.

Employing optical spectroscopy we have performed a comparative study of the dielectric response of extracellular matrix and filaments of electrogenic bacteria Shewanella oneidensis MR-1, cytochrome c, and bovine serum albumin. Combining infrared transmission measurements on thin layers with data of the terahertz spectra, we obtain the dielectric permittivity and AC conductivity spectra of the materials in a broad frequency band from a few cm-1 up to 7000cm-1 in the temperature range from 5 to 300K. Strong absorption bands are observed in the three materials that cover the range from 10 to 300cm-1 and mainly determine the terahertz absorption. When cooled down to liquid helium temperatures, the bands in Shewanella oneidensis MR-1 and cytochrome c reveal a distinct fine structure. In all three materials, we identify the presence of liquid bound water in the form of librational and translational absorption bands at ≈200 and ≈600cm-1, respectively. The sharp excitations seen above 1000cm-1 are assigned to intramolecular vibrations.

Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water

ABSTRACTUnderstanding the Raman spectroscopy at the atomistic level is important for the elucidation of dynamical processes in liquid water. Because the polarisability (or its time derivative) is often a highly nonlinear function of coordinates or/and momenta, we employ the linearized semiclassical initial value representation for quantum dynamical simulations of liquid water (and heavy water) under ambient conditions based on an ab initio-based, flexible, polarisable model (the POLI2VS force field). It is shown that quantum dynamical effects play a critical role in reproducing the peaks in the intermediate region between the librational and bending bands, those between the bending and stretching bands, and the double-peak in the stretching band in the experimental isotropic Raman spectrum. In contrast, quantum dynamical effects are important but less decisive in the anisotropic Raman spectrum. By selectively freezing either the intramolecular O–H stretching or H–O–H bending mode, we demonstrate that the peak in the intermediate region (2000–2400 cm−1) of the isotropic Raman spectrum arises from the interplay of the stretching and bending motions while a substantial part of the peak in the same intermediate region of the anisotropic Raman spectrum may be attributed to the combined motion of the bending and librational modes.

Simulations of the OKE Response in Simple Liquids Using a Polarizable and a Nonpolarizable Force Field.

Recently polarizable force fields are becoming increasingly popular for molecular dynamics simulations. As the signal obtained in the optical Kerr effect (OKE) experiment is due to the polarizability dynamics of the investigated system, a study is conducted in order to compare the experimental results with those obtained with the polarizable AMOEBA force field. The comparison is made in the frequency domain; however, time domain data are also included. The selected molecular systems are the isotropic carbon tetrachloride molecule, the anisotropic chloroform, carbon disulfide and acetone molecules, and the hydrogen-bonded water and methanol molecules. Different dipole-induced-dipole (DID) method variants are used for calculation of the OKE response, showing the importance of use of the all-atom approach with preoptimized atomic polarizabilities. In order to obtain a good intermolecular to intramolecular components amplitude ratio, the isotropic polarizability in the Thole correction needs to be updated between iterations. The convergence of the spectra calculated with different DID variants is also considered, and the approach that appears to be the best gives a very good approximation after three iterations. The comparison of the experimental and simulated spectra shows a rather good agreement for the non-hydrogen-bonded molecules, although the contribution of the reorientation of anisotropic molecules is overestimated. In the case of the hydrogen-bonded molecules, the theoretical spectra are far from the experimental ones. The highly overestimated librational bands indicate excessive polarizability anisotropy introduced by the potential model. Finally, in order to verify the significance of different components of the AMOEBA model, it is gradually simplified and compared with a simple reference potential model. Removal of polarizability shows a tremendous change in the case of hydrogen-bonded liquids, whereas for the other molecules it is of minor importance. The non-hydrogen-bonded liquids are, however, more sensitive to the presence of atomic multipoles in the model.

High-resolution synchrotron terahertz investigation of the large-amplitude hydrogen bond librational band of (HCN)2.

The high-resolution terahertz absorption spectrum of the large-amplitude intermolecular donor librational band ν of the homodimer (HCN)2 has been recorded by means of long-path static gas-phase Fourier transform spectroscopy at 207 K employing a highly brilliant electron storage ring source. The rovibrational structure of the ν band has the typical appearance of a perpendicular type band of a Σ-Π transition for a linear polyatomic molecule. The generated terahertz spectrum is analyzed employing a standard semi-rigid linear molecule Hamiltonian, yielding a band origin ν0 of 119.11526(60) cm-1 together with values for the excited state rotational constant B', the excited state quartic centrifugal distortion constant DJ' and the l-type doubling constant q for the degenerate state associated with the ν mode. The until now missing donor librational band origin enables the determination of an accurate experimental value for the vibrational zero-point energy of 2.50 ± 0.05 kJ mol-1 arising from the entire class of large-amplitude intermolecular modes. The spectroscopic findings are complemented by CCSD(T)-F12b/aug-cc-pV5Z (electronic energies) and CCSD(T)-F12b/aug-cc-pVQZ (force fields) electronic structure calculations, providing a (semi)-experimental value of 17.20 ± 0.20 kJ mol-1 for the dissociation energy D0 of this strictly linear weak intermolecular CHN hydrogen bond.

Solvation effects on the N–O and O–H stretching modes in hydrated NO3−(H2O)n clusters

NO3-(H2O)n clusters are a molecular model used to understand the solvation interaction between water and nitrate, an important anion in nature, industrial processes and biology. We demonstrate by ab initio molecular dynamics simulations that among the many isomeric structures at each cluster size with n = 1-6, thermal stability is an important consideration. The vibrational profile at a particular size, probed previously by infrared multiple photon dissociation (IRMPD) spectroscopy, can be accounted for by the isomers, which are both energetically and dynamically stable. Conversion and broadening due to the fluctuation of hydrogen bonds are important not only for the O-H stretching modes but also for the N-O stretching modes. Distinct patterns for the O-H stretching modes are predicted for the various solvation motifs. We also predict a surface structure for NO3-(H2O)n as n increases beyond 6, which can be verified by an early onset of strong libration bands for H2O in the IRMPD spectra and a flattening of the vertical detachment energy in the photoelectron spectra.

Development of non-bonded interaction parameters between graphene and water using particle swarm optimization.

New Lennard-Jones parameters have been developed to describe the interactions between atomistic model of graphene, represented by REBO potential, and five commonly used all-atom water models, namely SPC, SPC/E, SPC/Fw, SPC/Fd, and TIP3P/Fs by employing particle swarm optimization (PSO) method. These new parameters were optimized to reproduce the macroscopic contact angle of water on a graphene sheet. The calculated line tension was in the order of 10-11 J/m for the droplets of all water models. Our molecular dynamics simulations indicate the preferential orientation of water molecules near graphene-water interface with one OH bond pointing toward the graphene surface. Detailed analysis of simulation trajectories reveals the presence of water molecules with ≤∼1, ∼2, and ∼4 hydrogen bonds at the surface of air-water interface, graphene-water interface, and bulk region of the water droplet, respectively. Presence of water molecules with ≤∼1 and ∼2 hydrogen bonds suggest the existence of water clusters of different sizes at these interfaces. The trends observed in the libration, bending, and stretching bands of the vibrational spectra are closely associated with these structural features of water. The inhomogeneity in hydrogen bond network of water at the air-water and graphene-water interface is manifested by broadening of the peaks in the libration band for water present at these interfaces. The stretching band for the molecules in water droplet shows a blue shift as compared to the pure bulk water, which conjecture the presence of weaker hydrogen bond network in a droplet. © 2017 Wiley Periodicals, Inc.

Electrochemical ATR-SEIRAS Using Low-Cost, Micromachined Si Wafers

Thin, micromachined Si wafers, designed as internal reflection elements (IREs) for attenuated total reflectance infrared spectroscopy, are adapted to serve as substrates for electrochemical ATR surface enhanced infrared absorption spectroscopy (ATR-SEIRAS). The 500 μm thick wafer IREs with groove angles of 35° are significantly more transparent at long mid-IR wavelengths as compared to conventional large Si hemisphere IREs. The appeal of greater transparency is mitigated by smaller optical throughput at larger grazing angles and steeper angles of incidence at the reflecting plane that reduce the enhancement factor. Through use of the potential dependent adsorption of 4-methoxypyridine (MOP) as a test system, the microgroove IRE is shown to provide relatively strong electrochemical ATR-SEIRAS responses when the angle of incident radiation is between 50 and 55°, corresponding to refracted angles through the crystal of ∼40°. The higher than expected enhancement is attributed to attenuation of the reflection loss of p-polarized light and multiple reflections within the wafer-based IRE. The micromachined IREs are shown to outperform a 25 mm radius hemisphere in terms of S/N at wavenumbers less than ca. 1400 cm-1 despite the weaker signal enhancement derived from the steeper angle incident on the IRE/sample interface. The high optical transparency of the new IREs allows the spectral observation of displaced water libration bands at ca. 730 cm-1 upon solvent replacement by adsorbed MOP. The results are highly encouraging for the further development of low-cost, Si wafer-based IREs for electrochemical ATR-SEIRAS applications.

Influence of Polarization on Carbohydrate Hydration: A Comparative Study Using Additive and Polarizable Force Fields

Carbohydrates are known to closely modulate their surrounding solvent structures and influence solvation dynamics. Spectroscopic investigations studying far-IR regions (below 1000 cm(-1)) have observed spectral shifts in the libration band (around 600 cm(-1)) of water in the presence of monosaccharides and polysaccharides. In this paper, we use molecular dynamics simulations to gain atomistic insight into carbohydrate-water interactions and to specifically highlight the differences between additive (nonpolarizable) and polarizable simulations. A total of six monosaccharide systems, α and β anomers of glucose, galactose, and mannose, were studied using additive and polarizable Chemistry at HARvard Macromolecular Mechanics (CHARMM) carbohydrate force fields. Solvents were modeled using three additive water models TIP3P, TIP4P, and TIP5P in additive simulations and polarizable water model SWM4 in polarizable simulations. The presence of carbohydrate has a significant effect on the microscopic water structure, with the effects being pronounced for proximal water molecules. Notably, disruption of the tetrahedral arrangement of proximal water molecules was observed due to the formation of strong carbohydrate-water hydrogen bonds in both additive and polarizable simulations. However, the inclusion of polarization resulted in significant water-bridge occupancies, improved ordered water structures (tetrahedral order parameter), and longer carbohydrate-water H-bond correlations as compared to those for additive simulations. Additionally, polarizable simulations also allowed the calculation of power spectra from the dipole-dipole autocorrelation function, which corresponds to the IR spectra. From the power spectra, we could identify spectral signatures differentiating the proximal and bulk water structures, which could not be captured from additive simulations.

Quantum Local Monomer IR Spectrum of Liquid D2O at 300 K from 0 to 4000 cm(-1) Is in Near-Quantitative Agreement with Experiment.

The local monomer model is applied, with ab initio potential energy and dipole moment surfaces, to a calculation of the IR spectrum of liquid D2O at 300 K, over the spectral range 0 to 4000 cm(-1). The spectrum is an incoherent superposition of spectra of many monomers over snapshots of a molecular dynamics trajectory, where both intramolecular and intermolecular coupling in each monomer is treated. The comparison to experiment shows an unprecedented level of agreement for the stretch, bend, and libration bands and also the bend+libration and stretch+bend combination bands. This indicates that the incoherent approach captures much of the dynamics underlying the spectrum, provided monomer couplings are considered. The calculated spectrum is compared to the recently calculated IR spectrum of H2O, using the same method and potential energy and dipole moment surfaces, and shifts relative to that spectrum are presented and discussed.

The hydrogen-bond network of water supports propagating optical phonon-like modes.

The local structure of liquid water as a function of temperature is a source of intense research. This structure is intimately linked to the dynamics of water molecules, which can be measured using Raman and infrared spectroscopies. The assignment of spectral peaks depends on whether they are collective modes or single-molecule motions. Vibrational modes in liquids are usually considered to be associated to the motions of single molecules or small clusters. Using molecular dynamics simulations, here we find dispersive optical phonon-like modes in the librational and OH-stretching bands. We argue that on subpicosecond time scales these modes propagate through water's hydrogen-bond network over distances of up to 2 nm. In the long wavelength limit these optical modes exhibit longitudinal–transverse splitting, indicating the presence of coherent long-range dipole–dipole interactions, as in ice. Our results indicate the dynamics of liquid water have more similarities to ice than previously thought.

Spectroscopic identification of ethanol-water conformers by large-amplitude hydrogen bond librational modes

The far-infrared absorption spectra have been recorded for hydrogen-bonded complexes of water with ethanol embedded in cryogenic neon matrices at 2.8 K. The partial isotopic H/D-substitution of the ethanol subunit enabled by a dual inlet deposition procedure enables the observation and unambiguous assignment of the intermolecular high-frequency out-of-plane and the low-frequency in-plane donor OH librational modes for two different conformations of the mixed binary ethanol/water complex. The resolved donor OH librational bands confirm directly previous experimental evidence that ethanol acts as the O⋯HO hydrogen bond acceptor in the two most stable conformations. In the most stable conformation, the water subunit forces the ethanol molecule into its less stable gauche configuration upon dimerization owing to a cooperative secondary weak O⋯HC hydrogen bond interaction evidenced by a significantly blue-shift of the low-frequency in-plane donor OH librational band origin. The strong correlation between the low-frequency in-plane donor OH librational motion and the secondary intermolecular O⋯HC hydrogen bond is demonstrated by electronic structure calculations. The experimental findings are further supported by CCSD(T)-F12/aug-cc-pVQZ calculations of the conformational energy differences together with second-order vibrational perturbation theory calculations of the large-amplitude donor OH librational band origins.

Pressure Effect on the Boson Peak in Deeply Cooled Confined Water: Evidence of a Liquid-Liquid Transition.

The boson peak in deeply cooled water confined in nanopores is studied to examine the liquid-liquid transition (LLT). Below ∼180 K, the boson peaks at pressures P higher than ∼3.5 kbar are evidently distinct from those at low pressures by higher mean frequencies and lower heights. Moreover, the higher-P boson peaks can be rescaled to a master curve while the lower-P boson peaks can be rescaled to a different one. These phenomena agree with the existence of two liquid phases with different densities and local structures and the associated LLT in the measured (P, T) region. In addition, the P dependence of the librational band also agrees with the above conclusion.

Elucidation of the Water–DMSO Mixing Process Based on an IR Study

The IR spectra of water–DMSO mixtures were recorded in the range of 0 < x < 1 DMSO mole fraction in the spectral region of 400–6000 cm−1. The combination and libration bands of water and the S = O band of DMSO were used to monitor changes in the system. The combination and libration bands of water were analyzed in terms of different degrees of connectivity of water molecules. Two types of complexes differing by their water–DMSO ratio and the nature of water participating in their formation were found in the mixture. Complex I involves large water clusters stabilized by DMSO molecules and gradually diminishing in their size with increasing DMSO content. Water molecules separated from the bulk water interact with DMSO to form complex II. The water clusters were shown to exist in the mixture up to 90 mol% DMSO. The inertness of the inner molecules of the water cluster is the reason of the begining of DMSO accumulation at the notable water excess.

Proton Ordering of Cubic Ice Ic: Spectroscopy and Computer Simulations.

Several proton-disordered crystalline ice structures are known to proton order at sufficiently low temperatures, provided that the right preparation procedure is used. For cubic ice, ice Ic, however, no proton ordering has been observed so far. Here, we subject ice Ic to an experimental protocol similar to that used to proton order hexagonal ice. In situ FT-IR spectroscopy carried out during this procedure reveals that the librational band of the spectrum narrows and acquires a structure that is observed neither in proton-disordered ice Ic nor in ice XI, the proton-ordered variant of hexagonal ice. On the basis of vibrational spectra computed for ice Ic and four of its proton-ordered variants using classical molecular dynamics and ab initio simulations, we conclude that the features of our experimental spectra are due to partial proton ordering, providing the first evidence of proton ordering in cubic ice. We further find that the proton-ordered structure with the lowest energy is ferroelectric, while the structure with the second lowest energy is weakly ferroelectric. Both structures fit the experimental spectral similarly well such that no unique assignment of proton order is possible based on our results.

Quantum behaviour of water molecule in gemstone: terahertz fingerprints

We have shown that a weak interaction of a lone H2O molecule with the ‘‘walls’’ of nano-sized crystalline cage of gemstone (beryl) results in emergence of a rich set of molecular vibrational states. By analogy with translational and librational bands in liquid water and ice corresponding absorption bands are explained as due to translational (T) and librational (L) movements of the H2O molecule which is hydrogen bonded to the cage walls. In beryl crystal lattice, however, the six-fold symmetry of the cage brings about additional effect of splitting of the T and L bands into fine structure due to tunnelling within the six-well potential relief. The presented results will be of use for analysis of more complicated systems with confined water molecules like H2O chains in carbon nano-tubes, molecular clusters in e.g. zeolites, clays, silica gels and other natural or synthetic frameworks, as well as for interfacial water in biological systems.

The influence of charge on the structure and dynamics of water encapsulated in reverse micelles.

Hydrogen-bonded structure and relaxation dynamics of water entrapped inside reverse micelles (RMs) composed of surfactants with different charged head groups: sodium bis(2-ethylhexyl) sulfosuccinate (AOT) (anionic), didodecyldimethylammonium bromide (DDAB) (cationic) and Igepal CO-520 (Igepal) (nonionic) in cyclohexane (Cy) have been studied as a function of hydration (defined by ). Sub-diffusive slow (sub-ns) relaxation dynamics of water has been measured by the time resolved fluorescence spectroscopy (TRFS) technique using two fluorophores, namely 8-anilino-1-naphthalenesulfonic acid (ANS) and coumarin-343 (C-343). The hydrogen bonded connectivity network of water confined in these RMs has been investigated by monitoring the hydrogen bond stretching and libration bands of water using far-infrared FTIR spectroscopy. In addition, the ultrafast collective relaxation dynamics of water inside these RMs has been determined by dielectric relaxation in the THz region (0.2-2.0 THz) using THz time domain spectroscopy (THz-TDS). While TRFS measurements establish the retardation of water dynamics for all the RM systems, FTIR and THz-TDS measurements provide with signature of charge specificity.