Time-domain calculations of the polarized Raman spectra, the transient infrared absorption anisotropy, and the extent of delocalization of the OH stretching mode of liquid water

The polarized Raman spectrum and the time dependence of the transient infrared (TRIR) absorption anisotropy are calculated for the OH stretching mode of liquid water (neat liquid H2O) by using time-domain formulations, which include the effects of both the diagonal frequency modulations (of individual oscillators) induced by the interactions between the dipole derivatives and the intermolecular electric field, and the off-diagonal (intermolecular) vibrational coupling described by the transition dipole coupling (TDC) mechanism. The IR spectrum of neat liquid H2O and the TRIR anisotropy of a liquid mixture of H2O/HDO/D2O are also calculated. It is shown that the calculated features of these optical signals, including the temperature dependence of the polarized Raman and IR spectra, are in reasonable agreement with the experimental results, indicating that the frequency separation between the isotropic and anisotropic components of the polarized Raman spectrum and the rapid decay (approximately 0.1 ps) of the TRIR anisotropy of the OH stretching mode of neat liquid H2O are mainly controlled by the resonant intermolecular vibrational coupling described by the TDC mechanism. Comparing with the time evolution of vibrational excitations, it is suggested that the TRIR anisotropy decays in the time needed for the initially localized vibrational excitations to delocalize over a few oscillators. It is also shown that the enhancement of the dipole derivatives by the interactions with surrounding molecules is an important factor in generating the spectral profiles of the OH stretching Raman band. The time-domain behavior of the molecular motions that affect the spectroscopic features is discussed.     

Local hydrogen bonding dynamics and collective reorganization in water: ultrafast infrared spectroscopy of HOD/D(2)O

We present an investigation into hydrogen bonding dynamics and kinetics in water using femtosecond infrared spectroscopy of the OH stretching vibration of HOD in D(2)O. Infrared vibrational echo peak shift and polarization-selective pump-probe experiments were performed with mid-IR pulses short enough to capture all relevant dynamical processes. The experiments are self-consistently analyzed with a nonlinear response function expressed in terms of three dynamical parameters for the OH stretching vibration: the frequency correlation function, the lifetime, and the second Legendre polynomial dipole reorientation correlation function. It also accounts for vibrational-relaxation-induced excitation of intermolecular motion that appears as heating. The long time, picosecond behavior is consistent with previous work, but new dynamics are revealed on the sub-200 fs time scale. The frequency correlation function is characterized by a 50 fs decay and 180 fs beat associated with underdamped intermolecular vibrations of hydrogen bonding partners prior to 1.4 ps exponential relaxation. The reorientational correlation function observes a 50 fs librational decay prior to 3 ps diffusive reorientation. Both of these correlation functions compare favorably with the predictions from classical molecular dynamics simulations. The time-dependent behavior can be separated into short and long time scales by the 340 fs correlation time for OH frequency shifts. The fast time scales arise from dynamics that are mainly local: fluctuations in hydrogen bond distances and angles within relatively fixed intermolecular configurations. On time scales longer than the correlation time, dephasing and reorientations reflect collective reorganization of the liquid structure. Since the OH transition frequency and dipole are only weakly sensitive to these collective coordinates, this is a kinetic regime which gives an effective rate for exchange of intermolecular structures.     

The anharmonic vibrational potential and relaxation pathways of the amide I and II modes of N-methylacetamide

We investigate the influence of isotopic substitution and solvation of N-methylacetamide (NMA) on anharmonic vibrational coupling and vibrational relaxation of the amide I and amide II modes. Differences in the anharmonic potential of isotopic derivatives of NMA in D2O and DMSO-d6 are quantified by extraction of the anharmonic parameters and the transition dipole moment angles from cross-peaks in the two-dimensional infrared (2D-IR) spectra. To interpret the effects of isotopic substitution and solvent interaction on the anharmonic potential, density functional theory and potential energy distribution calculations are performed. It is shown that the origin of anharmonic variation arises from differing local mode contributions to the normal modes of the NMA isotopologues, particularly in amide II. The time domain manifestation of the coupling is the coherent exchange of excitation between amide modes seen as the quantum beats in femtosecond pump-probes. The biphasic behavior of population relaxation of the pump-probe and 2D-IR experiments can be understood by the rapid exchange of strongly coupled modes within the peptide backbone, followed by picosecond dissipation into weakly coupled modes of the bath.     

Isotopic dilution effect on the IR absorption bandshape of CH2Cl2

The shape of the 1265 cm⁻¹ absorption band in liquid CH₂Cl₂ changes on dilution with the deuterated compound. This is discussed in terms of intermolecular vibrational relaxation. The important influence of molecular reorientation on vibrational relaxation is pointed out.     

A theoretical study of the chiroptical properties of molecules with isotopically engendered chirality

There has been a considerable interest in the chiroptical properties of molecules whose chirality is exclusively due to an isotopic substitution and numerous examples for the electronic circular dichroism (CD) spectra of isotopically chiral systems have been reported in literature. Four different explanations have been proposed for the mechanism as to how the isotopic substitution induces a chiral perturbation of the otherwise achiral electronic wave function; however, up to now no conclusive answer has been given about the dominating effect responsible for the experimental observations. In this study we will present, for the first time, fully quantum-mechanical calculations of the CD spectra of three different molecular systems with isotopically engendered chirality. As examples, we consider the spectra of organic molecules with ketone and alpha-diketone carbonyl and diene chromophores. The effect of vibronic couplings for the reorientation of the electric and magnetic transition dipole moments is taken into account within the Herzberg-Teller approximation. The ground and excited state geometries and vibrational normal modes are obtained with (time-dependent) density functional theory [(TD)DFT], while the vibronic coupling effects are calculated at the TDDFT and density functional theory/multireference configuration interaction (DFT/MRCI) levels of theory. Generally, the band shapes of the experimental CD spectra are reproduced very well, and also the absolute CD intensities from the simulations are of the right order of magnitude. The sign and the intensity of the CD band are determined by a delicate balance of the contributions of a large number of individual vibronic transitions, and it is found that the vibrational normal modes with a large displacement are dominant. The separation of the calculated CD spectrum into the different contributions due to the overlap of the in-plane and out-of-plane components (regarding the symmetry plane of the unsubstituted molecule) of the electric and magnetic transition dipole moments yields information about the influence of the vibronic coupling effects for the reorientation of the corresponding transition dipole moments. In conclusion, the calculations clearly show that vibronic effects are responsible or at least dominant for the chiroptical properties of isotopically chiral organic molecules.     

Infrared and molecular-dynamics studies of the rotational dynamics of water highly diluted in supercritical CO2

Far-infrared (FIR) and mid-infrared (MIR) profiles of D2O infinitely dilute in supercritical CO2 have been studied using molecular-dynamics simulations. For this purpose, we have proposed an intermolecular potential model taking implicitly into account electron donor-acceptor (EDA) interactions between water and CO2 evaluated from ab initio calculations of the intermolecular potential-energy surface (IPS). Interaction-induced dipole mechanisms have been also taken into account in addition to the water permanent dipole to evaluate the simulated FIR profiles of water and CO2 polarizable molecules. They were found to play a minor role in the genesis of the FIR profiles of water/CO2 under supercritical conditions. The analysis of the reorientational dynamics of D2O shows that the rotational dynamics of water is weakly anisotropic due to the EDA interactions which affect more specifically the reorientational motions of the C2 symmetry axis of solute. These results have been used to assess the contribution of the vibrational relaxation in the experimental mid-infrared profiles associated with the nu1 symmetric and nu3 antisymmetric stretching and nu2 bending modes of D2O. It was found that the rotational dynamics mainly contribute to the broadening of the infrared (IR) profiles. Nevertheless, the vibrational processes play a role in the frequency shifts of the band centers and the relative intensity enhancements of the nu1 and nu3 modes of D2O. In particular, the EDA interactions between water and CO2 lead to the appearance of a well-defined IR band of the nu1 mode of D2O. Finally, a comparison with another model taking only into account dipole-quadrupole electrostatic interactions between water and CO2 molecules clearly reveals that EDA interactions have to be considered to reproduce both MIR and FIR measurements. From this point of view CO2 can be classified on a hydrophilic solvent scale based upon the solubility criterion as an intermediate solvent between "inert" xenon and carbon tetrachloride.     

Evolution of intermolecular structure and dynamics in supercritical carbon dioxide with pressure: an ab initio molecular dynamics study

The effect of pressure on supercritical carbon dioxide (scCO2) has been characterized by using Car-Parrinello molecular dynamics simulations. Structural and dynamical properties along an isotherm of 318.15 K and at pressures ranging from 190 to 5000 bar have been obtained. Intermolecular pair correlation functions and three-dimensional atomic probability density map calculations indicate that the local environment of a central CO2 molecule becomes more structured with increasing pressure. The closest neighbors are predominantly oriented in a distorted T-shaped geometry while neighbors separated by larger distances are likely oriented in a slipped parallel arrangement. The structure of scCO2 at high densities has been compared with that of crystalline CO2. The probability distributions of intramolecular distances narrow down with increasing pressure. A marginal but non-negligible effect of pressure on the instantaneous intramolecular OCO angle is observed, lending credence to the idea that intermolecular interactions between CO2 molecules in an inhomogeneous near neighbor environment could contribute to the observed instantaneous molecular dipole moment. The extent of deviation from a perfect linear geometry of the carbon dioxide molecule decreases with increasing pressure. Time constants derived from reorientational time correlation functions of the molecular backbone compare well with experimental data. Within the range of thermodynamic conditions explored here, no significant changes are observed in the frequencies of intramolecular vibrational modes. However, a blue shift is observed in the low-frequency cage rattling mode with increasing pressure.     

Infrared and Raman line shapes for ice Ih. II. H2O and D2O

We present a theoretical study of infrared and Raman line shapes of polycrystalline and single crystal ice Ih, for both water and heavy water, at 1, 125, and 245 K. Our calculations involve a mixed quantum/classical approach, a new water simulation model with explicit three-body interactions, transition frequency and dipole maps, and intramolecular and intermolecular vibrational coupling maps. Our theoretical spectra are in reasonable agreement with experimental spectra (available only near the two higher temperatures). We trace the origins of the different spectral peaks to weak and strong intermolecular couplings. We also discuss the delocalization of the vibrational eigenstates in terms of the competing effects of disorder and coupling.     

Direct Vibrational Energy Transfer in Monomeric Water Probed with Ultrafast Two Dimensional Infrared Spectroscopy

Vibrational relaxation dynamics of monomeric water molecule dissolved in d-chloroform solution were revisited using the two dimensional Infrared (2D IR) spectroscopy. The vibrational lifetime of OH bending in monomeric water shows a bi-exponential decay. The fast component (T₁=(1.2±0.1) ps) is caused by the rapid population equilibration between the vibrational modes of the monomeric water molecule. The slow component (T₂=(26.4±0.2) ps) is mainly caused by the vibrational population decay of OH bending mode. The reorientation of the OH bending in monomeric water is determined with a time constant of τ=(1.2±0.1) ps which is much faster than the rotational dynamics of water molecules in the bulk solution. Furthermore, we are able to reveal the direct vibrational energy transfer from OH stretching to OH bending in monomeric water dissolved in d-chloroform for the first time. The vibrational coupling and relative orientation of transition dipole moment between OH bending and stretching that effect their intra-molecular vibrational energy transfer rates are discussed in detail.     

Multidimensional anharmonic calculation of the vibrational frequencies and intensities for the trans and cis isomers of HONO with the use of normal coordinates

The equilibrium geometry and the potential energy and dipole moment surfaces have been determined for the cis and trans isomers of the HONO molecule by an ab initio Moller–Plesset (MP2) calculation with a wide set of atomic orbitals. The multidimensional anharmonic vibrational Schrodinger equations are solved using the variational method with the Hamiltonian and wave functions written in the normal coordinates of cis and trans isomers. All one- and two-dimensional and a number of three-dimensional vibrational problems are solved to obtain the energy levels and vibrational eigenfunctions. The frequencies and intensities for the fundamental, overtone and some combination bands are determined in good agreement with the available experimental results. The calculation shows the strength of coupling between different vibrational modes and reveals the presence of strong resonances between the (v₁, v₃, v₆) and (v₁, v₃−1, v₆+2) states of cis-HONO. This fact may be important for understanding the energy redistribution between the intermolecular degrees of freedom. The magnitude and direction of vibrationally averaged ground-state dipole moment of both isomers, as well as the direction of transition dipole moments, are in good agreement with the experimental findings. The changes in the values of dipole moment and some geometrical parameters of cis- and trans-HONO on vibrational excitation are also computed.     

Isotropic and anisotropic Raman scattering study in liquid p-dioxane

The isotropic and anisotropic profiles of the 835 and 2965 cm⁻¹ Raman lines of p-dioxane in the neat liquid and in solution have been studied as a function of temperature and concentration. From the correlation functions obtained by Fourier inversion of band contours, the possible interaction responsible for the vibrational dephasing of the oscillators and their reorientational relaxation are considered. It is shown that the p-dioxane molecule tumbling about the C₂ axis in the molecular plane perpendicular to the oxygen-oxygen direction proceeds by small-step Brownian diffusion associated with an Arrhenius activation energy of 9.0 kJ mol⁻¹. The vibrational relaxation mechanism of the two modes is interpreted in terms of pure dephasing due to weak collisions.     

Pressure and temperature dependence of the hydrogen bonding in supercritical ethanol: a computer simulation study

As a step toward deeper insight on the "hydrogen bonding" in supercritical ethanol (scEtOH), we carried out NVT molecular dynamics simulations of the fluid over a wide range of temperatures and pressures. The fluid was studied at SC conditions for which thermodynamic and spectroscopic (NMR, infrared, Raman, dielectric) data are available. The various site-site pair distribution functions (pdf's) were calculated, and their temperature and pressure dependence was obtained. It was found that over the thermodynamic conditions investigated here, scEtOH remains highly structured. Moreover, the characteristic behavior of the first peaks in H-H, O-O, and H-O pdf's reveals that hydrogen bonds still exist in scEtOH. The analysis focuses also on the reorientational dynamics of the bond unit vectors O-H, C-O, and of the permanent dipole moment of the molecules as well as the total dipole moment of the sample. The corresponding Legendre time correlation functions were discussed in connection to the "hydrogen bonding" in the fluid and in the context of experimental results. Specifically, the behavior of the O-H dynamics exhibits the well-known associative nature of the molecules in the system. A further analysis of the hydrogen bonds was carried out, and the degree of aggregation (average number of H-bonds per molecule) was obtained and compared with results from NMR chemical shift studies. Also the estimated monomer and free O-H groups in the fluid were compared with results from IR and Raman vibrational spectroscopy. The percentage analysis fi of the liquid and scEtOH molecules, with i = 0, 1, 2, 3, ... hydrogen bonds per molecule, has been obtained. The results show the existence of small, linear-chain oligomers formed mainly by two molecules, whereas the number of the three body oligomers, and specifically that of four body oligomers in the sample, is relatively small.     

Ab initio vibrational transition dipole moments and intensities of formaldehyde

Vibrational transition dipole moments and absorption band intensities for the ground state of formaldehyde, including the deuterated isotopic forms, are calculated. The analysis is based on ab initio SCF and CI potential energy and dipole moment surfaces. The formalism derives from second-order perturbation theory and involves the expansion of the dipole moment in terms of normal coordinates, as well as the incorporation of point group symmetry in the selection of the dipole moment components for the allowed transitions. Dipole moment expansion coefficients for the three molecule-fixed Cartesian coordinates of formaldehyde are calculated for internal and normal coordinate representations. Transition dipole moments and absorption band intensities of the fundamental, first overtone, combination, and second overtone transitions are reported. The calculated intensities and dipole moment derivatives are compared to experiment and discussed in the context of molecular orbital and bond polarization theory.     

Intermolecular coupling of the CO stretching vibrations in electrolyte solutions of carbonyl compounds

The polarized Raman scattering and infrared spectra of perchlorate solutions in acetone have been investigated in the CO stretching band region. The cation-dependent separation between the isotropic and anisotropic maxima of Raman band is interpreted in terms of intermolecular coupling of the CO vibrators in the solvation shell of cation. The linear correlation between the isotropic-anisotropic separation and integral intensities of the IR band indicates that the induced dipole mechanism of the coupling dominates.     

Vibrational coupling in collision-induced raman scattering

The observed splitting of the collision-induced Raman band v₃ in liquid SF₆ is interpreted as vibrational exciton line splitting The bandwidth of the collision-induced Raman bands v₃ and v₆, in liquid SF₆ and v₂ and v₃ in liquid CO₂ and CS₂ can be explained by contributions from reonentational motion and from transition dipole-transition dipole vibrational coupling.     

Computer simulation of the system N2 in fluid argon-correlation functions and relaxation times

The correlation functions of the dipole moment, P₂[u(0) · u(t)], angular and linear velocity, and bond forces have been calculated from computer simulated data for four different density-temperature states of N₂ in fluid argon. From these functions infrared and Raman line shapes, NMR relaxation times, and rotational and classical vibrational relaxation times have been computed.     

Comparison between infrared and raman band profiles in the absence of rotational relaxation in the liquid phase: ν(OH) of methanol in acetonitrile

In the case of weakly hydrogen bonded systems OH…B it has been shown that the ν(OH) IR and Raman band profiles are not identical. The origin of this difference is discussed in terms of changes in the dipole moment and polarizability derivatives with intermolecular forces     

Angular correlation functions of spherical-top carbonyl complexes in gaseous N2 influenced by coriolis coupling

The Fourier transform of measured line shapes of the triply degenerate vibrations of carbonyl complexes in gaseous N₂ shows strong deviations from the angular correlation functions calculated by Steele for the rigid spherical-top molecules. These deviations are explained theoretically with the aid of Coriolis coupling and vibrational relaxation. For molecules in the gas phase it is possible to separate the Coriolis coupling from the vibrational relaxation.     

Vibrational dephasing and hydrodynamic effects on vibrational relaxation rates in acetophenone: Raman bandshape analysis

The isotropic component of Raman band for C=O stretching mode of acetophenone in solution was analyzed by estimating the correlation coefficient with reference to Lorentzian lineshape. In the intermediate region of solute/solvent concentration there is a sharp decrease in the correlation coefficient and there appears to be a transition from non-Lorentzian to Lorentzian lineshape. The vibrational relaxation rates have been estimated from the isotropic component of Raman band in different solvents. The rate is shown to be dependent on several macroscopic as well as microscopic properties of the solute-solvent system and intermolecular interactions. The hydrodynamic and dispersion forces appear to play a major role in determining the vibrational relaxation rate and the broadening of the bands.     

Studies of molecular motions and vibrational relaxation VII. Rotational diffusion model to investigate the band shapes of the ν5 and ν8 ban

Reorientational correlation functions G_rot^(l)(t)) for the degenerate (E) bands of liquid acetonitrile (CH₃CN) have been computed us NMR spin-lattice relaxation data (for CD₃CN) and gas phase Raman band profiles, assuming that the rotational diffusion model is valid. The effects both anisotropic rotational motion and of Coriolis coupling are included. The predicted correlation functions along with those calculated using ther cl “free” rotor equations, have been compared with those obtained from the υ_s (Raman) and υ_s (IR and Raman) experimental band profiles. It shown that, despite the simplicity of the model and obvious (understandable) discrepancies at short times, sensible conclusions may be drawn. This work a starting point for the testing of more complicated models for reorientational motion in dense phases.