Vibrational energy relaxation rates of H2 and D2 in liquid argon via the linearized semiclassical method

The vibrational energy relaxation (VER) rates for H2 and D2 in liquid argon (T=152 K, rho=1.45x1022 cm-3) are calculated using the linearized semiclassical (LSC) method (J. Phys. Chem. 2003, 107, 9059, 9070). The calculation is based on Fermi's golden rule. The VER rate constant is expressed in terms of the quantum-mechanical force-force correlation function, which is then estimated using the LSC method. A local harmonic approximation (LHA) is employed in order to compute the multidimensional Wigner integrals underlying the LSC approximation. The H2-Ar and D2-Ar interactions are described by the three-body potential of Bissonette et al. (J. Phys. Chem. A 1996, 105, 2639). The LHA-LSC-based VER rate constants for both D2 and H2 are found to be about 2-3 orders of magnitude slower than those obtained experimentally. However, their ratio agrees quantitatively with the corresponding experimental result. In contrast, the classical VER rate constants are found to be 8-9 orders of magnitude slower than those obtained experimentally, and their ratio is found to be qualitatively different from the corresponding experimental result.     

Vibrational energy relaxation of polyatomic molecules in liquid solution via the linearized semiclassical method

Vibrational energy relaxation (VER) of polyatomic, as opposed to diatomic, molecules can occur via different, often solvent assisted, intramolecular and/or intermolecular pathways. In this paper, we apply the linearized semiclassical (LSC) method for calculating VER rates in the prototypical case of a rigid, symmetrical and linear triatomic molecule (A-B-A) in a monatomic liquid. Starting at the first excited state of either the symmetric or asymmetric stretches, VER can occur either directly to the ground state or indirectly via intramolecular vibrational relaxation (IVR). The VER rate constants for the various pathways are calculated within the framework of the Landau-Teller formalism, where they are expressed in terms of two-time quantum-mechanical correlation functions. The latter are calculated by the LHA-LSC method, which puts them in a "Wignerized" form, and employs a local harmonic approximation (LHA) in order to compute the necessary multidimensional Wigner integrals. Results are reported for the LHL/Ar model of Deng and Stratt [J. Chem. Phys. 2002, 117, 1735], as well as for CO(2) in liquid argon and in liquid neon. The LHA-LSC method is shown to give rise to significantly faster VER and IVR rates in comparison to the classical treatment, particularly at lower temperatures. We also find that the type and extent of the quantum rate enhancement is strongly dependent on the particular VER pathway. Finally, we find that the classical and semiclassical treatments can give rise to opposite trends when it comes to the dependence of the VER rates on the solvent.     

Vibrational energy relaxation of a diatomic molecule in a room-temperature ionic liquid

Vibrational energy relaxation (VER) dynamics of a diatomic solute in ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate (EMI(+)PF(6) (-)) are studied via equilibrium and nonequilibrium molecular dynamics simulations. The time scale for VER is found to decrease markedly with the increasing solute dipole moment, consonant with many previous studies in polar solvents. A detailed analysis of nonequilibrium results shows that for a dipolar solute, dissipation of an excess solute vibrational energy occurs almost exclusively via the Lennard-Jones interactions between the solute and solvent, while an oscillatory energy exchange between the two is mainly controlled by their electrostatic interactions. Regardless of the anharmonicity of the solute vibrational potential, VER becomes accelerated as the initial vibrational energy increases. This is attributed primarily to the enhancement in variations of the solvent force on the solute bond, induced by large-amplitude solute vibrations. One interesting finding is that if a time variable scaled with the initial excitation energy is employed, dissipation dynamics of the excess vibrational energy of the dipolar solute tend to show a universal behavior irrespective of its initial vibrational state. Comparison with water and acetonitrile shows that overall characteristics of VER in EMI(+)PF(6) (-) are similar to those in acetonitrile, while relaxation in water is much faster than the two. It is also found that the Landau-Teller theory predictions for VER time scale obtained via equilibrium simulations of the solvent force autocorrelation function are in reasonable agreement with the nonequilibrium results.     

Vibrational energy relaxation in liquid N2CO mixtures

The vibrational energy relaxation rates of the liquid nitrogenCO system have been measured by optically pumping the collision-induced fundamental vibrational absorption band of liquid N₂ with the output of an HBr TEA laser. A radiatively dominated value of 56 ± 10 s is found for the intrinsic nitrogen relaxation time. The CO contribution to the decay rate is explained on the basis of a simple kinetic model and found also to be radiatively dominated at low CO concentrations. The importance of radiative trapping and energy transport in evaluating the lifetimes is demonstrated.     

Vibrational energy relaxation of small molecules and ions in liquids

A theoretical/computational framework for determining vibrational energy relaxation rates, pathways, and mechanisms, for small molecules and ions in liquids, is presented. The framework is based on the system—bath coupling approach, Fermi’s golden rule, classical time-correlation functions, and quantum correction factors. We provide results for three specific problems: relaxation of the oxygen stretch in neat liquid oxygen at 77 K, relaxation of the water bend in chloroform at room temperature, and relaxation of the azide ion anti-symmetric stretch in water at room temperature. In each case, our calculated lifetimes are in reasonable agreement with experiment. In the latter two cases, theory for the observed solvent isotope effects illuminates the relaxation pathways and mechanisms. Our results suggest several propensity rules for both pathways and mechanisms.     

Vibrational excitation and dissociation of a diatomic molecule in the diffusion approximation

A numeric solution is given for the excitation of vibration and dissociation in the diffusion approximation for an O₂-Ar (oxygen-argon) mixture. The time of vibrational relaxation is calculated, and also the rate constant for dissociation. The diffusion approximation describes the processes satisfactorily for a gas whose particles differ little in mass.I am indebted to B. V. Kuksenko, A. I. Osipov, and M. N. Safaryan for discussions and valuable advice.     

Vibrational energy relaxation of isotopically labeled amide I modes in cytochrome c: theoretical investigation of vibrational energy relaxation rates and pathways

With use of a time-dependent perturbation theory, vibrational energy relaxation (VER) of isotopically labeled amide I modes in cytochrome c solvated with water is investigated. Contributions to the VER are decomposed into two contributions from the protein and water. The VER pathways are visualized by using radial and angular excitation functions for resonant normal modes. Key differences of VER among different amide I modes are demonstrated, leading to a detailed picture of the spatial anisotropy of the VER. The results support the experimental observation that amide I modes in proteins relax with subpicosecond time scales, while the relaxation mechanism turns out to be sensitive to the environment of the amide I mode.     

Vibrational energy relaxation of the bend fundamental of dilute water in liquid chloroform and d-chloroform

The population lifetimes of the bend fundamental of dilute water in liquid chloroform (8.5 ps) and d-chloroform (28.5 ps) display an interesting solvent isotope effect. As the lowest excited vibrational state of the molecule, the water bend fundamental relaxes directly to the ground state with about 1600 cm-1 of energy released to the other degrees of freedom. The strong solvent isotope effect along with the large energy gap indicates the participation of solvent vibrational modes in this vibrational energy relaxation process. We calculate the vibrational energy relaxation rates of the water bend in chloroform and d-chloroform using the Landau-Teller formula with a new potential model developed and parametrized self-consistently to describe the chloroform-water interaction. The computed values are in reasonable agreement with the experimental results, and the trend for the isotope effect is correct. It is found that energy transfer to the solvent vibrations does indeed play an important role. Nevertheless, no single dominant solvent accepting mode can be identified; the relaxation appears to involve both the bend and the C-Cl stretches, and frequency changes of all of these modes upon deuteration contribute to the observed solvent isotope effect.     

Vibrational relaxation and line widths in liquids: a comparison of theory and experiment

Published experimental data on the Raman line widths of some a₁ bands of methyl iodide and benzene are compared with theoretical predictions of relaxation theories which assume either rapid or slow time scales for intermolecular fluctuations. The observed variations with density and with mole fraction are found to fit the rapid motion theory with an additional temperature independent core contribution.     

Linearized semiclassical initial value time correlation functions using the thermal Gaussian approximation: applications to condensed phase systems

The linearized approximation to the semiclassical initial value representation (LSC-IVR) has been used together with the thermal Gaussian approximation (TGA) (TGA/LSC-IVR) [J. Liu and W. H. Miller, J. Chem. Phys. 125, 224104 (2006)] to simulate quantum dynamical effects in realistic models of two condensed phase systems. This represents the first study of dynamical properties of the Ne(13) Lennard-Jones cluster in its liquid-solid phase transition region (temperature from 4 to 14 K). Calculation of the force autocorrelation function shows considerable differences from that given by classical mechanics, namely that the cluster is much more mobile (liquidlike) than in the classical case. Liquid para-hydrogen at two thermodynamic state points (25 and 14 K under nearly zero external pressure) has also been studied. The momentum autocorrelation function obtained from the TGA/LSC-IVR approach shows very good agreement with recent accurate path integral Monte Carlo results at 25 K [A. Nakayama and N. Makri, J. Chem. Phys. 125, 024503 (2006)]. The self-diffusion constants calculated by the TGA/LSC-IVR are in reasonable agreement with those from experiment and from other theoretical calculations. These applications demonstrate the TGA/LSC-IVR to be a practical and versatile method for quantum dynamics simulations of condensed phase systems.     

Vibrational relaxation and the noncoincidence effect in liquid N,N-dimethylacetamide

The Raman linewidth and the noncoincidence effect in N,N-dimethylacetamide have been studied and explained in terms of molecular attraction parameters. The transition dipole-transition dipole interaction is probably the predominant mechanism responsible for the noncoincidence effect. The line broadening of the isotropic component is explained on the basis of dispersion type interactions.     

Vibrational relaxation in liquid trimethylchlorogermane: an attempt to determine reliable vibrational correlation functions

Vibrational correlation functions and corresponding correlation times were determined by a variety of techniques for the v_s(GeCl) and sym-v_s(Ge(CH₃)₃) vibrations of trimethylchlorogermane from Raman bandshape analysis. The results clearly show that it is not possible to obtain reliable correlation functions and other related data for these vibrations due to the joint effects of isotopic splitting and Fermi resonances.     

Fluctuating energy level Landau-Teller theory: application to the vibrational energy relaxation of liquid methanol

State-to-state vibrational energy relaxation (VER) rates of the OH-stretch fundamental to select vibrational modes of liquid methanol are presented. The rates are calculated via a modified, fluctuating Landau-Teller (FLT) theory approach, which allow for dynamical vibrational energy level shifts. These rates are then compared to previously published results from Gulmen and Sibert [J. Phys. Chem. A 2004, 108, 2389] for the traditional Landau-Teller (LT) method as well as results calculated through time-dependent perturbation theory (TD), which naturally allow for the fluctuation. For the first time, this method is applied to a polyatomic molecular system, and the FLT theory greatly reduces the discrepancy between the LT and TD results or, at a minimum, is comparable to the LT approach with very little additional computational cost.     

Vibrational energy relaxation of the ND-stretching vibration of NH(2)D in liquid NH(3)

The vibrational energy relaxation from the first excited ND-stretching mode of NH(2)D dissolved in liquid NH(3) is studied using molecular dynamics simulations. The rate constants for inter- and intramolecular energy transfer are calculated in the framework of the quantum-classical Landau-Teller theory. At 273 K and an ammonia density of 0.642 g cm(-3) the calculated ND-stretch lifetime of τ = 9.1 ps is in good agreement with the experimental value of 8.6 ps. The main relaxation channel accounting for 52% of the energy transfer involves an intramolecular transition to the first excited state of the umbrella mode. The energy difference between both states is taken up by the near-resonant bending vibrations of the solvent. Less important for the ND-stretch lifetime are both the direct transition to the ground state and intramolecular relaxation via the NH(2)D bending modes contributing 23% each. Our calculations imply that the experimentally observed weak density dependence of τ is caused by detuning the resonance between the ND-stretch-umbrella energy gap and the solvent accepting modes which counteracts the expected linear increase of the relaxation rate with density.     

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.     

Test of the consistency of various linearized semiclassical initial value time correlation functions in application to inelastic neutron scattering from liquid para-hydrogen

The linearized approximation to the semiclassical initial value representation (LSC-IVR) is used to calculate time correlation functions relevant to the incoherent dynamic structure factor for inelastic neutron scattering from liquid para-hydrogen at 14 K. Various time correlations functions were used which, if evaluated exactly, would give identical results, but they do not because the LSC-IVR is approximate. Some of the correlation functions involve only linear operators, and others involve nonlinear operators. The consistency of the results obtained with the various time correlation functions thus provides a useful test of the accuracy of the LSC-IVR approximation and its ability to treat correlation functions involving both linear and nonlinear operators in realistic anharmonic systems. The good agreement of the results obtained from different correlation functions, their excellent behavior in the spectral moment tests based on the exact moment constraints, and their semiquantitative agreement with the inelastic neutron scattering experimental data all suggest that the LSC-IVR is indeed a good short-time approximation for quantum mechanical correlation functions.     

Delocalization corrections to the strictly localized molecular orbitals: A linearized SCF approximation

An explicit formula is derived for calculating the delocalization corrections (tails) to be added to the strictly localized bond orbitals. It was obtained by solving analytically the SCF problem for the interbond interactions in a linearized approximation. The model calculations at the CNDO/2 level show that this simple approach is sufficient to account for the molecular conformations.     

Applications of the RISM equation to diatomic fluids: the liquids nitrogen,oxygen and bromine

The reference interaction site model (RISM) integral equation is used to study the equilibrium pair correlation for one-component liquids composed of homonuclear diatomic molecules. The integral equation is first tested by comparing the results obtained from it with those of computer simulation calculations. An internal consistency test is developed which seems to provide an a priori measure of the accuracy of the RISM equation. Then the theory is used to interpret the neutron scattering structure factor data taken on the liquids nitrogen and oxygen. For both of these liquids, a satisfactory explanation of the data is obtained by assuming that the short ranged repulsive forces between molecules are mimicked by a two-site hard core model. However, this simple model does not provide a satisfactory explanation of the data taken from liquid bromine. But, it is shown that a slightly more sophisticated model for the short-ranged repulsion between Br₂ molecules does provide an adequate explanation. With the molecular models determined by fitting the neutron scattering data, the RISM equation provides a method for determining the atom, atom to center-of-mass, and center-of-mass to center-of-mass intermolecular distribution functions in the diatomic liquids. From these functions, the local structures in the three liquids are analyzed. While orientational pair correlations are nearly negligible in both the liquids nitrogen and oxygen, these correlations are fairly substantial in liquid bromine. Furthermore, even when the orientation of a molecule does not greatly influence the orientation of its neighbors, it is found that the orientation of a molecule does have an important effect on the location of its neighbors. Thus, the coupling of translational and orientational coordinates is significant, even in the liquids nitrogen and oxygen.     

Simultaneous flow measurement of specific heats and thermal expansion coefficients of liquids: Aqueoust-BuOH mixtures and neat alkanols and alkanediols at 25°C

A flow method is described for the simultaneous measurement of volumetric specific heat capacities c _p ^v and thermal expansion coefficients of liquids, using a Picker heat capacity microcalorimeter. the method involves a simple sequence of operations in which the calorimeter unit is used alternately to measure a difference in volumetric specific heats at constant flow rate and temperature, or a flow rate variation resulting from a temperature scan on a flow cell. The precision of the method in measurements is close to 1% which is comparable to that of current dilatometric methods; the overall time for the combined c _p ^v - measurement is less than 30 minutes. Typical uses and reliability of the method are illustrated through results for aqueous NaCl solutions, t-BuOH-water mixtures, and ,-diols (C₂–C₅) at 25°C. The and c _p ^v values in homologous series of normal alkanes, alcohols and diols appear well represented by equations comprising two contributions, one related to the hydrocarbon chain length, the other dependent on the number density of –OH groups. The general trends in these data and comparison with similar results for H₂O₂ and H₂O illustrate the magnitude of the c _p ^v and anomalies in liquid water.