A molecular theory of liquid phase vibrational energy relaxation (VER) [S. A. Adelman et al., Adv. Chem. Phys. 84, 73 (1993)] is applied to study the temperature T and density ρ dependencies of the VER rate constant k(T,ρ)=T1−1, where T1 is the energy relaxation time, of model Lennard-Jones systems that roughly simulate solutions of high-mass, low-frequency dihalogen solutes in rare gas solvents; specifically the I2/Xe, I2/Ar, and ICI/Xe solutions. For selected states of these systems, the theory’s assumptions are tested against molecular dynamics (MD) results. The theory is based on the expression T1=β−1(ωl), where ωl and β(ω) are, respectively, the solute’s liquid phase vibrational frequency and vibrational coordinate friction kernel. The friction kernel is evaluated as a cosine transform of the fluctuating force autocorrelation function of the solute vibrational coordinate, conditional that this coordinate is fixed at equilibrium. Additionally, the early-time decay of the force autocorrelation function is approximated by a Gaussian function which is exact to order t2. This Gaussian approximation permits evaluation of T1 in terms of integrals over equilibrium solute–solvent pair correlation functions. The pair correlation function formulas yield T1’s in semiquantitative agreement with those found by MD evaluations of the Gaussian approximation, but with three orders of magnitude less computational effort. For the isothermal ρ dependencies of k(T,ρ), the theory predicts for all systems that the Gaussian decay time τ is nearly independent of ρ. This in turn implies that k(T,ρ) factorizes into a liquid phase structural contribution and a gas phase dynamical contribution, yielding a first-principles form for k(T,ρ) similar to that postulated by the isolated binary collision model. Also, the theory predicts both “classical” superlinear rate isotherms, and “nonclassical” sublinear isotherms similar to those recently observed by Troe and co-workers for azulene relaxation in supercritical fluids. The isochoric T dependencies of k(T,ρ) are studied in the range 300 to 1000 K. For none of the solutions are the rate isochores found to accurately conform to either Arrhenius or Landau–Teller kinetics.
Read full abstract