Two-dimensional infrared spectroscopy of neat ice Ih

Transferability of vibrational spectroscopic map from TIP4P to TIP4P-like water models

We computed the IR, Raman, and sum frequency generation spectra of water in the OH-stretch region by employing the quantum/classical mixed approach that consists of a vibrational spectroscopic map and molecular dynamics (MD) simulation. We carried out the MD simulation with the TIP4P, TIP4P/2005, and TIP4P/Ice models and applied the map designed for TIP4P by Skinner et al. to each MD trajectory. Although the map is not tuned for TIP4P-like models, TIP4P/2005 and TIP4P/Ice provide the best reproduction of the experimental vibrational spectra of liquid water and crystalline ice, respectively. This result demonstrates the transferability of the map from TIP4P to TIP4P/2005 and TIP4P/Ice, meaning that one can choose an appropriate TIP4P-like model to calculate the vibrational spectra of an aqueous system without rebuilding the map.

Quantum-chemical calculation of two-dimensional infrared spectra using localized-mode VSCF/VCI

Computational protocols for the simulation of two-dimensional infrared (2D IR) spectroscopy usually rely on vibrational exciton models which require an empirical parameterization. Here, we present an efficient quantum-chemical protocol for predicting static 2D IR spectra that does not require any empirical parameters. For the calculation of anharmonic vibrational energy levels and transition dipole moments, we employ the localized-mode vibrational self-consistent field (L-VSCF)/vibrational configuration interaction (L-VCI) approach previously established for (linear) anharmonic theoretical vibrational spectroscopy [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365-3377 (2014)]. We demonstrate that with an efficient expansion of the potential energy surface using anharmonic one-mode potentials and harmonic two-mode potentials, 2D IR spectra of metal carbonyl complexes and dipeptides can be predicted reliably. We further show how the close connection between L-VCI and vibrational exciton models can be exploited to extract the parameters of such models from those calculations. This provides a novel route to the fully quantum-chemical parameterization of vibrational exciton models for predicting 2D IR spectra.

Multiscale analysis of the hydrate based carbon capture from gas mixtures containing carbon dioxide

To reveal the kinetic performance of gas molecules in hydrate growth, hydrate formation from pure CO₂, flue gas, and biogas was measured using in-situ Raman and macroscopic methods at 271.6 K. In the in-situ Raman measurements, Raman peaks of gases in the hydrate phase were characterised and normalised by taking the water bands from 2800 to 3800 cm^-1 as a reference, whose line shapes were not found to have a noticeable change in the conversion from Ih ice to sI hydrate. The hydrate growth was suggested to start with the formation of unsaturated hydrate nuclei followed by gas adsorption. In hydrate formed from all tested gases, CO₂ concentrations in hydrate nuclei were found to be 23-33% of the saturation state. In the flue gas system, the N₂ concentration reached a saturation state once hydrate nuclei formed. In the biogas system, competitive adsorption of CH₄ and CO₂ molecules was observed, while N₂ molecules hardly evolved in hydrate formation. Combined with micro- and macroscopic analysis, small molecules such as N₂ and CO₂ were suggested to be more active in the formation of hydrate nuclei, and the preferential adsorption of CO₂ molecules took place in the subsequent gas adsorption process.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Icelike Vibrational Properties of Strong Hydrogen Bonds in Hydrated Lithium Nitrate

The hydrogen bond network accounts for many of the extraordinary physical properties of liquid water and ice. Its vibrational dynamics are quite complex in their entirety but can be accessed in detail by investigating small groups of only a few water molecules. Here, aqueous salt hydrates turned out to be an exceptional model system for water molecules arranged in well-defined geometrical structures that can be accessed by means of femtosecond spectroscopy of the OH stretching vibration. In this study, we find striking resemblance between the vibrational properties of three water molecules connected via strong hydrogen bonds in the trihydrate of LiNO₃ and those of ordinary ice I_h. As in ice, the vibrations of the hydrate water molecules show ultrafast excited state dynamics that are strongly accelerated when proceeding from deuterated to neat H₂O samples. The latter is analyzed by means of an additional relaxation channel that is due to Fermi resonance between the OH stretching vibration and the bend overtone accompanied by delocalization of the vibration over neighboring water molecules in the H₂O species. Moreover, in the hydrate and ice samples severe spectral broadening is examined when comparing fundamental and excited state absorption bands. Here, proton delocalization along the strong hydrogen bonds is given as a possible underlying mechanism.

Multiscale coupled Maxwell's equations and polarizable molecular dynamics simulation based on charge response kernel model

A computational scheme of coupled Maxwell's equations and polarizable molecular dynamics simulation has been developed based on a multi-scale model to describe the coupled dynamics of light electromagnetic waves and molecules in crystalline solids, where the charge response kernel model is employed to incorporate electronic polarization of the molecules. The method is applicable to electronically non-resonant light-matter interaction systems that involve atomic motions in spectroscopy and photonics. Since the scheme simultaneously traces the light propagation in a medium on a macroscopic scale and the microscopic molecular motion under the light electric field, this enables us to treat the experimental setup and mimic its measurement process. As the first applications, we demonstrate three numerical examples of basic spectroscopies of an ice crystalline solid: simulations of reflection and transmission of visible light, infrared absorption measurement, and stimulated Raman scattering measurement. These examples show the detailed behaviors of the interacting light fields and molecules in the spectroscopic processes.

Mapping the amide-I vibrations of model dipeptides with secondary structure sensitivity and amino acid residue specificity, and its application to amyloid β-peptide in aqueous solution

Vibrational spectroscopy has been known as particularly well-suited for deciphering the polypeptide's structure. To decode structural information encoded in IR spectra, we developed amide-I frequency maps on the basis of model dipeptides to correlate the amide-I frequency of interest to the combination of the calculated secondary structure dependent amide-I frequency by using DFT method and the electrostatic potentials that projected onto the amide unit from the micro-environment within molecular mechanics force field. The constructed maps were applied to model dipeptides and amyloid β-peptide fragment (Aβ_25-35). The dipeptide specified map (DS map) and the hybrid map (HYB map) predicted amide-I bands of Aβ_25-35 in solution satisfactorily reproduce experimental observation, and indicate the preference of forming β-sheet and random coil structure for Aβ_25-35 in D₂O just as the results of cluster analysis suggested. These maps with secondary structural sensitivity and amino acid residue specificity open up a way for the interpretation of amide-I vibrations and show their potentials in the understanding of molecular structure of polypeptides in solution.

Delocalization and stretch-bend mixing of the HOH bend in liquid water

Liquid water's rich sub-picosecond vibrational dynamics arise from the interplay of different high- and low-frequency modes evolving in a strong yet fluctuating hydrogen bond network. Recent studies of the OH stretching excitations of H₂O indicate that they are delocalized over several molecules, raising questions about whether the bending vibrations are similarly delocalized. In this paper, we take advantage of an improved 50 fs time-resolution and broadband infrared (IR) spectroscopy to interrogate the 2D IR lineshape and spectral dynamics of the HOH bending vibration of liquid H₂O. Indications of strong bend-stretch coupling are observed in early time 2D IR spectra through a broad excited state absorption that extends from 1500 cm^-1 to beyond 1900 cm^-1, which corresponds to transitions from the bend to the bend overtone and OH stretching band between 3150 and 3550 cm^-1. Pump-probe measurements reveal a fast 180 fs vibrational relaxation time, which results in a hot-ground state spectrum that is the same as observed for water IR excitation at any other frequency. The fastest dynamical time scale is 80 fs for the polarization anisotropy decay, providing evidence for the delocalized or excitonic character of the bend. Normal mode analysis conducted on water clusters extracted from molecular dynamics simulations corroborate significant stretch-bend mixing and indicate delocalization of δ_HOH on 2-7 water molecules.

Simple Quantum Dynamics with Thermalization

In this paper, we introduce two simple quantum dynamics methods. One is based on the popular surface-hopping method, and the other is based on rescaling of the propagation on the bath ground-state potential surface. The first method is special, as it avoids specific feedback from the simulated quantum system to the bath and can be applied for precalculated classical trajectories. It is based on the equipartition theorem to determine if hops between different potential energy surfaces are allowed. By comparing with the formally exact Hierarchical Equations Of Motion approach for four model systems we find that the method generally approximates the quantum dynamics toward thermal equilibrium very well. The second method is based on rescaling of the nonadiabatic coupling and also neglect the effect of the state of the quantum system on the bath. By the nature of the approximations, they cannot reproduce the effect of bath relaxation following excitation. However, the methods are both computationally more tractable than the conventional fewest switches surface hopping, and we foresee that the methods will be powerful for simulations of quantum dynamics in systems with complex bath dynamics, where the system-bath coupling is not too strong compared to the thermal energy.

Understanding water structure from Raman spectra of isotopic substitution H2O/D2O up to 573 K

The OH/OD stretch band features on Raman spectra of isotopic substitution H₂O/D₂O at temperatures up to 573 K are correlated with a multi-structure model that water has five dominant hydrogen bonding configurations: tetrahedral, deformed tetrahedral, single donor, single hydrogen bonded water and free water.

Differences in the Vibrational Dynamics of H(2)O and D(2)O: Observation of Symmetric and Antisymmetric Stretching Vibrations in Heavy Water

Water's ability to donate and accept hydrogen bonds leads to unique and complex collective dynamical phenomena associated with its hydrogen-bond network. It is appreciated that the vibrations governing liquid water's molecular dynamics are delocalized, with nuclear motion evolving coherently over the span of several molecules. Using two-dimensional infrared spectroscopy, we have found that the nuclear motions of heavy water, D2O, are qualitatively different than those of H2O. The nonlinear spectrum of liquid D2O reveals distinct O-D stretching resonances, in contrast to H2O. Furthermore, our data indicates that condensed-phase O-D vibrations have a different character than those in the gas phase, which we understand in terms of weakly delocalized symmetric and antisymmetric stretching vibrations. This difference in molecular dynamics reflects the shift in the balance between intra- and intermolecular couplings upon deuteration, an effect which can be understood in terms of the anharmonicity of the nuclear potential energy surface.