Vibrational relaxation of OH and CH fundamentals of polar and nonpolar molecules in the condensed phase

Studies of vibrational energy flow in various polar and nonpolar molecules that follows the ultrafast excitation of the CH and OH stretch fundamentals, modeled using semiclassical methods, are reviewed. Relaxation rates are calculated using Landau-Teller theory and a time-dependent method, both of which consider a quantum mechanical solute molecule coupled to a classical bath of solvent molecules. A wide range of decay rates are observed, ranging from 1 ps for neat methanol to 50 ps for neat bromoform. In order to understand the flow rates, it is argued that an understanding of the subtle mixing between the solute eigenstates is needed and that solute anharmonicities are critical to facilitating condensed phase vibrational relaxation. The solvent-assisted shifts of the solute vibrational energy levels are seen to play a critical role of enhancing or decreasing lifetimes.     

70 years of Landau-Teller theory for collisional energy transfer. Semiclassical three-dimensional generalizations of the classical collinear model

This article, in historical retrospective, describes the development of the celebrated Landau-Teller (LT) model of 1936 for vibrational-translational energy exchange in collisions of an atom with a diatomic molecule. We discuss semiclassical generalizations of the classical LT model and generalizations of the collinear LT model to account for the effects of rotation of the diatom on the vibrational relaxation rate. The former is based on the recovery of the Landau semiclassical exponent from the classical LT encounter time, and the latter on the definition of a 1-D driving mode within the manifold of the translational and rotational degrees of freedom of the colliding partners. The utility of generalized LT models is illustrated by three case studies that exemplify weak and strong effects of the rotation as well as the efficiencies of different driving modes in the vibrational relaxation of highly asymmetric diatoms.     

Vibrational energy relaxation of the OH(D) stretch fundamental of methanol in carbon tetrachloride

The lifetimes of the hydroxyl stretch fundamentals of two methanol isotopomers, MeOH and MeOD, in carbon tetrachloride solvent are calculated through the use of the perturbative Landau-Teller and fluctuating Landau-Teller methods. Examination of these systems allows for insight into the nature of the vibrational couplings that lead to intramolecular vibrational energy transfer. While both systems display energy transfer to nearly degenerate modes, MeOD also displays strong coupling to an off-resonant vibration. The relaxation of MeOH and MeOD occurs through transitions involving a total change in the vibrational quanta of 4 and 3, respectively. We calculate vibrational energy relaxation lifetimes of 4-5 ps for MeOH and 2-3 ps for MeOD that agree well with the experimentally determined values.     

Comparison between the Landau-Teller and flux-flux methods for computing vibrational energy relaxation rate constants in the condensed phase

The calculation of vibrational energy relaxation (VER) rate constants in the condensed phase is usually based on the Landau-Teller formula, which puts them in terms of the Fourier transform, at the vibrational frequency, of the autocorrelation function of the force exerted on the relaxing mode by the bath modes. An alternative expression for the VER rate constant puts it in terms of the autocorrelation function of the vibrational energy flux. In this paper, we compare the predictions obtained via those two methods in the case of iodine in liquid xenon. We find that the computational cost underlying both methods is comparable and that they predict similar VER rates. However, while the calculation of the VER rate via the Landau-Teller formula is somewhat more direct, the predictions obtained via the flux-flux formula are in somewhat better agreement with the VER rates obtained from nonequilibrium molecular dynamics simulations.     

Time scales and pathways of vibrational energy relaxation in liquid CHBr3 and CDBr3

Molecular dynamics simulations are used in conjunction with Landau-Teller, fluctuating Landau-Teller, and time-dependent perturbation theories to investigate energy flow out of various vibrational states of liquid CHBr3 and CDBr3. The CH stretch overtone is found to relax with a time scale of about 1 ps compared to the 50 ps rate for the fundamental. The relaxation pathways and rates for the CD stretch decay in CDBr3 are computed in order to understand the changes arising from deuteration. While the computed relaxation rate agrees well with experiments, the pathway is found to be more complex than anticipated. In addition to the above channels for CH(D) stretch relaxation that involve only the hindered translations and rotations of the solvent, routes involving off-resonant and resonant excitations of solvent vibrational modes are also examined. Finally, the decay of energy from low frequency states to near-lying solute states and solvent vibrations are studied.     

Vibrational energy transfer in H2 + He

A semiclassical approach is developed to study vibrational energy transfer in H₂ + He by use of the a priori interaction potential including all nonzero impact parameter collisions. The calculated values of the rate coefficient are found to be in excellent agreement with experimental data which are available in the temperature ranges 60–450 K and 1350-3000 K. The temperature dependence is shown to seriously deviate from the Landau-Teller prediction below 1000 K. The calculation was carried out over the temperature range of 30 to 10000 K.     

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 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 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.     

Semiclassical extension of the Landau-Teller theory of collisional energy transfer

A semiclassical version of the quantum coupled-states approximation for the vibrational relaxation of diatomic molecules in collisions with monatomic bath gases is presented. It is based on the effective mass approximation and a recovery of the semiclassical Landau exponent from the classical Landau-Teller collision time. For an interaction with small anisotropy, the Landau exponent includes first order corrections with respect to the orientational dependence of the collision time and the effective mass. The relaxation N(2)(v=1)-->N(2)(v=0) in He is discussed as an example. Employing the available vibrationally elastic potential, the semiclassical approach describes the temperature dependence of the rate constant k(10)(T) over seven orders of magnitude across the temperature range of 70-3000 K in agreement with experimental data and quantum coupled-states calculations. For this system, the hierarchy of corrections to the Landau-Teller conventional treatment in the order of importance is the following: quantum effects in the energy release, dynamical contributions of the rotation of N(2) to the vibrational transition, and deviations of the interaction potential from a purely repulsive form. The described treatment provides significant simplifications over complete coupled-states calculations such that applications to more complex situations appear promising.     

Vibrational relaxation in the condensed phase

A modified Landau-Teller equation for vibrational relaxation in the condensed phase is proposed. This equation differs from previous approaches by accounting for the fluctuations of the energies of the vibrational levels that result from the interactions with the surroundings (bath). In the conventional approach the effects of the bath are only included in the coupling between the relaxing and accepting vibrational modes. It is shown that the additional inclusion of the fluctuations of the energy levels can lead to a dramatic change of the vibrational relaxation rate.     

Vibrational energy relaxation dynamics of Si---H stretching modes on stepped H/Si(111) 1×1 surfaces

The vibrational energy relaxation rates of excited Si---H stretching modes on the monohydride steps of miscut H/Si(111) 1×1 surfaces are calculated using Bloch-Redfield theory combined with classical molecular dynamics (MD) simulation. The structure and vibrational frequencies of the surface are first investigated using the Car-Parrinello ab initio MD method. The calculated Si---Si---H bending frequencies and relaxed structures are then used to refine the empirical potential for the classical MD simulations. The lifetime of the excited Si---H stretching mode at the step is found to be shorter than the modes on the terrace. Both the magnitude and the trend of the calculated results agree well with the experimental measurement on the 9° monohydride stepped surface. The vibrational relaxation rate of the Si---H stretching modes on the 15° monohydride stepped surface are also calculated and predicted to have a slightly shorter lifetime than for the 9° surface.     

Quantum mechanical calculation of energy dependence of OCl/OH product branching ratio and product quantum state distributions for the O(1D) + HCl reaction on all three contributing electronic state potential energy surfaces

OCl/OH product branching ratios are calculated as a function of total energy for the O( (1) D) + HCl reaction using quantum wavepacket methods. The calculations take account of reaction on all the three electronic state potential energy surfaces which correlate with both reactants and products. Our results show that reaction on the excited electronic state surfaces has a large effect on the branching ratio at higher energies and that these surfaces must therefore be fully taken into account. The calculations use the potential energy surfaces of Nanbu and co-workers. Product vibrational and rotational quantum state distributions are also calculated as a function of energy for both product channels. Inclusion of the excited electronic state potential energy surfaces improves the agreement of the predicted product vibrational quantum state distributions with experiment for the OH product channel. For OCl agreement between theory and experiment is retained for the vibrational quantum state distributions when the excited electronic state potential energy surfaces are included in the analysis. For the rotational state distributions good agreement between theory and experiment is maintained for energies at which experimental results are available. At higher energies, above 0.7 eV of total energy, the OCl rotational state distributions predicted using all three electronic state potential energy surfaces shift to markedly smaller rotational quantum numbers.     

Polymorphism of crystalline alpha-quaterthiophene and alpha-sexithiophene: ab initio analysis and comparison with inelastic neutron scattering response

Phonons in the alpha-quaterthiophene (4T) and alpha-sexithiophene (6T) polymorph phases are investigated using the direct method combined with density functional theory (DFT)-based total energy calculations. The simulation of inelastic neutron scattering spectra (INS) on the LT and HT polymorph phases of 4T and 6T enable the corresponding spectral signatures of these materials to be identified. In particular, there are two fingerprints: (i) the low-frequency vibrational modes (frequencies lower than 200 cm(-1)) and (ii) the vibrational modes in the 600-900 cm(-1) frequency range. The good agreement with the INS experimental data allows us to assign unambiguously the origin of all features (first-order and high-order processes) of these spectra and to predict that the LT phase is the phase measured experimentally both on the 4T and 6T materials. Moreover, the broad background in the 600-1400 cm(-1) frequency range and the well-defined features which appear around 940 cm(-1) in the calculated INS spectra of 4T/HT and 6T/HT are assigned to multiphonon contributions. This multiphonon contribution at 940 cm(-1), which is absent in the 4T/LT and 6T/LT INS spectra, also constitutes a fingerprint of the HT phases. Finally, the calculated dispersion curves of the two polymorph phases of 4T and 6T are given.     

Potential energy surface, van der Waals motions, and vibronic transitions in phenol-argon complex

The structure and intermolecular vibrational energy levels of the phenol-Ar complex are calculated from its potential energy surface. This surface is constructed from a large set of the interaction energy values computed using second-order Moller-Plesset perturbation theory with the augmented correlation consistent polarized valence double-zeta basis set. The global minimum in the potential energy surface corresponds to a cluster structure with Ar located over the geometric center of the phenol ring at a distance of 3.510 A and shifted by 0.1355 A towards oxygen. The calculated dissociation energy of 371 cm(-1) is in accordance with the experiment. Additional local minima higher in energy are with Ar placed in the phenol plane. However, they are too shallow to form the bound states corresponding to planar isomers. The deformation of the potential energy surface shape, created by the interaction of Ar with the phenolic oxygen, is responsible for a pronounced intermode mixing. As a result, a set of hybrid stretching-bending states appears which cannot be described in terms of the standard models. The intermode coupling is reflected in the vibronic structure of the S1-S0 electronic transition. The intensities of the vibronic bands are calculated from the electronic transition dipole moment surfaces determined using the ab initio single-excitation configuration interaction method. They allow us to correct and complete the assignment of the spectra observed in phenol-Ar, as well as in the analogous complexes of phenol with Kr and Xe.     

The molecular potential energy surface and vibrational energy levels of methyl fluoride. Part II

New analytical bending and stretching, ground electronic state, potential energy surfaces for CH(3)F are reported. The surfaces are expressed in bond-length, bond-angle internal coordinates. The four-dimensional stretching surface is an accurate, least squares fit to over 2000 symmetrically unique ab initio points calculated at the CCSD(T) level. Similarly, the five-dimensional bending surface is a fit to over 1200 symmetrically unique ab initio points. This is an important first stage towards a full nine-dimensional potential energy surface for the prototype CH(3)F molecule. Using these surfaces, highly excited stretching and (separately) bending vibrational energy levels of CH(3)F are calculated variationally using a finite basis representation method. The method uses the exact vibrational kinetic energy operator derived for XY(3)Z systems by Manson and Law (preceding paper, Part I, Phys. Chem. Chem. Phys., 2006, 8, DOI: 10.1039/b603106d). We use the full C(3v) symmetry and the computer codes are designed to use an arbitrary potential energy function. Ultimately, these results will be used to design a compact basis for fully coupled stretch-bend calculations of the vibrational energy levels of the CH(3)F system.     

Quantum state-resolved energy transfer dynamics at gas-liquid interfaces: IR laser studies of CO2 scattering from perfluorinated liquids

An apparatus for detailed study of quantum state-resolved inelastic energy transfer dynamics at the gas-liquid interface is described. The approach relies on supersonic jet-cooled molecular beams impinging on a continuously renewable liquid surface in a vacuum and exploits sub-Doppler high-resolution laser absorption methods to probe rotational, vibrational, and translational distributions in the scattered flux. First results are presented for skimmed beams of jet-cooled CO(2) (T(beam) approximately 15 K) colliding at normal incidence with a liquid perfluoropolyether (PFPE) surface at E(inc) = 10.6(8) kcal/mol. The experiment uses a tunable Pb-salt diode laser for direct absorption on the CO(2) nu(3) asymmetric stretch. Measured rotational distributions in both 00(0)0 and 01(1)0 vibrational manifolds indicate CO(2) inelastically scatters from the liquid surface into a clearly non-Boltzmann distribution, revealing nonequilibrium dynamics with average rotational energies in excess of the liquid (T(s) = 300 K). Furthermore, high-resolution analysis of the absorption profiles reveals that Doppler widths correspond to temperatures significantly warmer than T(s) and increase systematically with the J rotational state. These rotational and translational distributions are consistent with two distinct gas-liquid collision pathways: (i) a T approximately 300 K component due to trapping-desorption (TD) and (ii) a much hotter distribution (T approximately 750 K) due to "prompt" impulsive scattering (IS) from the gas-liquid interface. By way of contrast, vibrational populations in the CO(2) bending mode are inefficiently excited by scattering from the liquid, presumably reflecting much slower T-V collisional energy transfer rates.