Combined electronic structure/molecular dynamics approach for ultrafast infrared spectroscopy of dilute HOD in liquid H2O and D2O

We present a new approach that combines electronic structure methods and molecular dynamics simulations to investigate the infrared spectroscopy of condensed phase systems. This approach is applied to the OH stretch band of dilute HOD in liquid D2O and the OD stretch band of dilute HOD in liquid H2O for two commonly employed models of water, TIP4P and SPC/E. Ab initio OH and OD anharmonic transition frequencies are calculated for 100 HOD x (D2O)n and HOD x(H2O)n (n = 4-9) clusters randomly selected from liquid water simulations. A linear empirical relationship between the ab initio frequencies and the component of the electric field from the solvent along the bond of interest is developed. This relationship is used in a molecular dynamics simulation to compute frequency fluctuation time-correlation functions and infrared absorption line shapes. The normalized frequency fluctuation time-correlation functions are in good agreement with the results of previous theoretical approaches. Their long-time decay times are 0.5 ps for the TIP4P model and 0.9 ps for the SPC/E model, both of which appear to be somewhat too fast compared to recent experiments. The calculated line shapes are in good agreement with experiment, and improve upon the results of previous theoretical approaches. The methods presented are simple, and transferable to more complicated systems.     

Theoretical study of the infrared and Raman line shapes of liquid methanol

Using a combined electronic structure and molecular dynamics simulation method, we calculated the infrared and Raman spectra for the OH vibrations in liquid CH(3)OH. The vibrational frequencies, transition dipole moments, and transition polarizabilities are obtained from density functional theory calculations and then mapped into an empirical relation to the electric field on the H atom along the OH bond. Vibrational couplings between OH chromophores on different molecules are treated using transition dipole interactions. The simulated infrared and Raman line shapes are in good agreement with experimental observations. We have also shown that the vibrations of non-hydrogen-bonded OH groups contribute significantly to the difference between the IR and Raman line shapes.     

Collective solvent coordinates for the infrared spectrum of HOD in D2O based on an ab initio electrostatic map

An ab initio MP2 vibrational Hamiltonian of HOD in an external electrostatic potential parametrized by the electric field and its gradient-tensor is constructed. By combining it with the fluctuating electric field induced by the D(2)O solvent obtained from molecular dynamics simulations, we calculate the infrared absorption of the O-H stretch. The resulting solvent shift and infrared line shape for three force fields (TIP4P, SPC/E, and SW) are in good agreement with the experiment. A collective coordinate response for the solvent effect is constructed by identifying the main electrostatic field and gradient components contributing to the line shape. This allows a realistic stochastic Liouville equation simulation of the line shapes which is not restricted to Gaussian frequency fluctuations.     

Multidimensional infrared spectroscopy of water. I. Vibrational dynamics in two-dimensional IR line shapes

In this and the following paper, we describe the ultrafast structural fluctuations and rearrangements of the hydrogen bonding network of water using two-dimensional (2D) infrared spectroscopy. 2D IR spectra covering all the relevant time scales of molecular dynamics of the hydrogen bonding network of water were studied for the OH stretching absorption of HOD in D2O. Time-dependent evolution of the 2D IR line shape serves as a spectroscopic observable that tracks how different hydrogen bonding environments interconvert while changes in spectral intensity result from vibrational relaxation and molecular reorientation of the OH dipole. For waiting times up to the vibrational lifetime of 700 fs, changes in the 2D line shape reflect the spectral evolution of OH oscillators induced by hydrogen bond dynamics. These dynamics, characterized through a set of 2D line shape analysis metrics, show a rapid 60 fs decay, an underdamped oscillation on a 130 fs time scale induced by hydrogen bond stretching, and a long time decay constant of 1.4 ps. 2D surfaces for waiting times larger than 700 fs are dominated by the effects of vibrational relaxation and the thermalization of this excess energy by the solvent bath. Our modeling based on fluctuations with Gaussian statistics is able to reproduce the changes in dispersed pump-probe and 2D IR spectra induced by these relaxation processes, but misses the asymmetry resulting from frequency-dependent spectral diffusion. The dynamical origin of this asymmetry is discussed in the companion paper.     

HOD分子在(~C)1B1态的光解动力学:产物OH/OD的分支比以及OD键解离能

Photodissociation of jet-cooled HOD via the ? state around 124 nm has been studied using the H(D)-atom Rydberg tagging time-of-flight technique. Rotational state resolved action spectrum and the product translational energy distribution spectra have been recorded for both D+OH and H+OD dissociation channels. Product channel OH/OD branching ratios for the individual ?- X rotational transition have been determined. A comparison is also given with the B- X and ?- X transitions. In addition, the dissociation energy of the OD bond in HOD has been determined accurately to be 41751.3±5 cm^-1.     

Rotational state specific dissociation dynamics of HOD → H + OD via two-photon excitation to the C̃ electronic state

The dissociation dynamics of HOD via two-photon excitation to the C? state have been investigated using the H-atom Rydberg tagging time-of-flight (TOF) technique. The H-atom action spectrum for the C? ← X? transition shows resolved rotational structure. Product translational energy distributions and angular distributions have also been recorded for the H + OD channel for three excited levels each with k(a)′ = 2. From these distributions, quantum state distributions and angular anisotropy parameters (β2 and β4) for the OD product were determined. These results are consistent with the nonadiabatic predissociation picture illustrated in the one-photon dissociation process for H2O. The heterogeneous dissociation pathway via Coriolis coupling is the dominant dissociation process in the present study. A high proportion of the total available energy is deposited into the rotational energy of the OD product. The anisotropic recoil distributions reveal the distinctive contributions from the alignment of the excited states and the dissociation process. Comparisons are also made between the results for HOD and H2O via the equivalent rotational transitions. The OH bond energy, D(o)(H?OD), of the HOD molecule is also determined to be 41283.0 ± 5 cm(-1).     

Vibrational spectra of Na, K, Mn2+, Ni2+ and Zn2+ salts of 1,2,4,5-benzenetetracarboxylic (pyromellitic) acid--a short hydrogen bond evidence

Five salts of 1,2,4,5-benzenetetracarboxylic acid (pyromellitic acid), [C6H2(COO)4H4], have been synthesized and investigated by infrared and Raman spectroscopy and by single crystal X-ray diffraction methods: sodium salt [Na2(H2O)2][C6H2(COO)4H2], potassium salt [K(H2O)3][C6H2(COO)4H3] and transition metal salts [M(H2O)6][C6H2(COO)4H2], which M = Mn, Ni and Zn. Crystal structures of all five compounds show short intramolecular asymmetric hydrogen bonds (SHB) between adjacent carboxyl groups with O...O distance average of 2.40 A. The Raman and infrared spectra reported indicate the presence of short hydrogen bonds in all salts, in agreement with the X-ray data. The O-H stretching mode [nu(OH)] had been observed at about 2500 cm(-1). Deuterated analogues were synthesized and their Raman spectra show that nu(OH)/nu(OD) ratio average is about unit. The symmetric [nu(sym)(O..H..O)] and asymmetric [nu(asym)(O..H..O)] stretching modes have been attributed about 300 and 870 cm(-1), respectively, in all salts, and for deuterated analogues, the ratio nu(OH)/nu(OD) to nu(sym)(O..H..O, O..D..O) is close to unit like it occurs in nu(OH). The vibrational modes, mainly SHB modes, are tentatively assigned by molecular orbital ab initio calculations of pyromellitic acid and anions [C6H2(COO)4H3]- and [C6H2(COO)4H2]2-. Geometry optimizations showed a good agreement with experimental data. Frequency calculation confirms the assignment of specific vibrational modes. Ab initio calculations show that nu(C=O) and nu(sym)(COO) are strongly coupled with in plane OH bending [delta(OH)]. In Raman spectra of deuterated analogues is observed a frequency shift of these bands.     

The vibrational Stokes shift of water (HOD in D2O)

The vibrational Stokes shift of the OH stretching transition nu(OH) of water is the shift between the ground-state absorption and the excited-state (v=1) emission. A recent measurement on HOD in D(2)O solvent [S. Woutersen and H. J. Bakker, Phys. Rev. Lett. 83, 2077 (1999)] of a 70 cm(-1) redshift, and a subsequent calculation of a 57 cm(-1) redshift using equilibrium molecular dynamics simulations [C. P. Lawrence and J. L. Skinner, J. Chem. Phys. 117, 8847 (2002)] were in good agreement. We now report extensive measurements of the vibrational Stokes shift in HOD/D(2)O using an ultrafast IR pump, Raman probe method. The vibrational Stokes shift is seen to depend on the pump pulse frequency and on time delay; by varying these parameters it can be made to range from 112 to -32 cm(-1) (negative values indicate a blueshift in the excited state). The equilibrium vibrational Stokes shift is actually a negative rather than a positive quantity. Possible reasons for the disagreement between experiment and theory are briefly discussed.     

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.     

Raman,infrared and force field studies of K2(12)C2O4 x H2O and K2(13)C2O4 x H2O in the solid state and in aqueous solution,and of (NH4)2(12)C2O4 x H2O and (NH4)2(13)C2O4 x H2O in the solid state

The Raman and infrared spectra of solid K2(12)C2O4 x H2O are reported together with, for the first time, the corresponding Raman and infrared spectra of solid K2(13)C2O4 x H2O. Raman spectra of aqueous solutions of both isotopomers are also reported. In the solid state the oxalate anion is planar with D2h symmetry in this salt, whereas in aqueous solution the Raman spectra of the anion are best interpreted on the basis of D2d symmetry. The Raman spectra of solid (NH4)2(12)C2O4 x H2O and (NH4)2(13)C2O4 x H2O, in which the oxalate anion is twisted from planarity by 28 degrees about the CC bond, have also been recorded. Several reassignments have been made. The harmonic force field for the oxalate anion in the D2h, D2 and D2d geometries has been determined in part, and approximate values of key valence force constants determined. All the observed band wavenumbers and 12C/13C isotopic shifts are well reproduced by the force fields. The potential energy distribution of the totally symmetric normal modes of planar oxalate indicates that each mode consists of extensively mixed symmetry corrdinates and that the labels previously used for the bands seen here at 475 and 879 cm(-1) would better be described as v(CC) and deltaS(CO2), respectively, putting them in the same wavenumber order as v(NN) and deltaS(NO2) for the isoelectronic and isostructural molecule N2O4. The stretching force constants of N2O4 and planar C2O4(2-) are established to be in the order f(NN) < f(CC) and f(NO) > f(CO), consistent with the known relative bond lengths.     

Two-dimensional infrared spectroscopy of isotope-diluted ice Ih

We present experimental 2D IR spectra of isotope diluted ice Ih (i.e., the OH stretch mode of HOD in D(2)O and the OD stretch mode of HOD in H(2)O) at T = 80 K. The main spectral features are the extremely broad 1-2 excited state transition, much broader than the corresponding 0-1 groundstate transition, as well as the presence of quantum beats. We do not observe any inhomogeneous broadening that might be expected due to proton disorder in ice Ih. Complementary, we perform simulations in the framework of the Lippincott-Schroeder model, which qualitatively reproduce the experimental observations. We conclude that the origin of the observed line shape features is the coupling of the OH-vibrational coordinate with crystal phonons and explain the beatings as a coherent oscillation of the O···O hydrogen bond degree of freedom.     

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.     

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.     

Structure and dynamics of water confined in dimethyl sulfoxide.

We study the structure and dynamics of hydrogen-bonded complexes of H2O/D2O and dimethyl sulfoxide (DMSO) by infrared spectroscopy, NMR spectroscopy and ab initio calculations. We find that single water molecules occur in two configurations. For one half of the water monomers both OH/OD groups form strong hydrogen bonds to DMSO molecules, whereas for the other half only one of the two OH/OD groups is hydrogen-bonded to a solvent molecule. The H-bond strength between water and DMSO is in the order of that in bulk water. NMR deuteron relaxation rates and calculated deuteron quadrupole coupling constants yield rotational correlation times of water. The molecular reorientation of water monomers in DMSO is two-and-a-half times slower than in bulk water. This result can be explained by local structure behavior.     

Simultaneous Sensing of H2O,D2O and HOD through Peroxo Vibrations

Detection of HOD simultaneously in the presence of a mixture of H₂O and D₂O is still an experimental challenge. Till date, there is no literature report of simultaneous detection of H₂O, D₂O and HOD based on vibrational spectra. Herein we report simultaneous quantitative detection of H₂O, D₂O and HOD in the same reaction mixture with the help of bridged polynuclear peroxo complex in absence and presence of Au nanoparticles on the basis of a peroxide vibrational mode in resonance Raman and surface enhanced resonance Raman spectrum. We synthesize bridged polynuclear peroxo complex in different solvent mixture of H₂O and D₂O. Due to the formation of different nature of hydrogen bonding between peroxide and solvent molecules (H₂O, D₂O and HOD), vibrational frequency of peroxo bond is significantly affected. Mixtures of different H₂O and D₂O concentrations produce different HOD concentrations and that lead to different intensities of peaks positioned at 897, 823 and 867 cm⁻¹ indicating H₂O, D₂O and HOD, respectively. The lowest detection limits (LODs) were 0.028 mole fraction of D₂O in H₂O and 0.046 mole faction of H₂O in D₂O. In addition, for the first time the results revealed that the cis-peroxide forms two hydrogen bonds with solvent molecules.     

Full-dimensional quantum dynamics study of the H + H2O and H + HOD exchange reactions

The initial state-selected time-dependent wave packet approach is employed to study the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H, HOH' + D exchange reactions with both OH bonds in the H(2)O reactant and OH(D) bond in the HOD reactant treated as reactive bonds. The total reaction probabilities for different partial waves, as well as the integral cross sections, which are the exact CC (coupled-channel) results, are first obtained in this study for the H(2)O(HOD) reactant initially in the ground rovibrational state. Because of the shallow C(3v) minimum along the reaction path, the reaction probabilities for the three reactions present several resonance peaks, with one dominant resonance peak just above the threshold. The cross sections for the H' + HOD → HOH' + D reaction are substantially smaller than those for the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H reactions, indicating that the H'/H exchange reactions are much more favored. In the CC calculations, the resonance peaks in the reaction probabilities diminish quickly with the increase in total angular momenta J, resulting in the existence of a clear step-like feature just above the threshold in the cross sections for the title reactions, which manifests the signature of shape resonances in these reactions. In the CS calculations, the resonance peaks on reaction probabilities persist in many partial waves, and thus the resonance structures can no longer survive the partial-wave summation and are washed out completely in the CS cross sections for the title reactions.     

Raman spectroscopic study of the uranyl minerals vanmeersscheite U(OH)4[(UO2)3(PO4)2(OH)2].4H2O and arsenouranylite Ca(UO2)[(UO2)3(AsO4)2(OH)2].(OH)2.6H2O

Raman and infrared spectra of secondary uranyl phosphate vanmeersscheite and Raman spectrum of secondary uranyl arsenate arsenuranylite were recorded and interpreted, and the spectra related to the structure of the minerals. Observed bands were attributed to the stretching and bending vibrations of uranyl, phosphate and/or arsenate units and OH (H(2)O and OH(-)) units. Phosphuranylite sheet topology is characteristic for both minerals. U-O bond lengths in uranyl were calculated from the spectra and compared with those inferred for vanmeersscheite from the X-ray single crystal structure analysis. O-H...O hydrogen bonds in both minerals were also inferred using the Libowitzky empirical relation.     

Mass spectrometric separation and quantitation of overlapping isotopologues. H2O/HOD/D2O and H2Se/HDSe/D2Se mixtures

Three conceptually different mathematical methods are presented for accurate mass spectrometric determination of H2O/HOD/D2O and H2Se/HDSe/D2Se concentrations from mixtures. These are alternating least-squares, weighted two-band target entropy minimization, and a statistical mass balance model. The otherwise nonmeasurable mass spectra of partially deuterated isotopologues (HOD and HDSe) are mathematically constructed. Any recorded isotopologue mixture mass spectra are then deconvoluted by least-squares into their components. This approach is used to study the H2O/D2O exchange reaction, and is externally validated gravimetrically. The H2O/D2O exchange equilibrium constant is also measured from the deconvoluted 70 eV electron impact GC/MS data (K = 3.85 +/- 0.03).     

高分辨偏正拉曼光谱对丙酸酐C=O振动模的拉曼光谱非一致效应研究

根据C=O振动的各向同性和各向异性拉曼光谱和红外光谱特点讨论研究了丙酸酐分子的局部有序排列以及振动耦合机理. 利用三级联共聚焦拉曼光谱仪测定了不同浓度丙酸酐的各向同性与各向异性拉曼光谱图, 分别采集了丙酸酐在四氯化碳和甲醇中的光谱以及不同极性溶剂中的光谱, 具体分析了丙酸酐C=O振动模的浓度效应、 溶剂效应以及拉曼光谱非一致效应(NCE). 结果表明, 丙酸酐C=O振动模的NCE效应随着浓度的降低而减小; 随着溶剂极性的减小而增加. 利用密度泛函理论的B3LYP-D3/31-311G(d,p)基组计算了丙酸酐单体和二聚体的几何稳定构型, 用聚集态理论模型解释了丙酸酐分子的NCE效应、 浓度效应与溶剂效应. 理论计算结果与实验结果相吻合.