![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
The common trajectory effective mass (CTEM) approach is applied to the calculation of the temperature dependence of the vibrational relaxation rate coefficient k10 in collisions of a diatomic molecule with chemically inert atom on an anisotropic PES. As a case study, we consider the O2–Ar system for which experimental data across the wide temperature range and ab initio results for the PES are available. It is shown that temperature dependence of the vibrational relaxation rate coefficient k10 is described by the CTEM approach with non-empirical value of the effective Landau–Teller characteristic temperature that indicates the contribution of O2 rotation in inducing vibrational transition.
GRAPHICAL ABSTRACT
![](/cms/asset/7759f705-8e96-4e6f-b1d9-6afb508cbf2b/tmph_a_1938266_uf0001_oc.jpg)
1. Introduction
The experimental values of vibration relaxation rate coefficients k10 (T) for transitions in diatomic molecules in collision with atom are often represented in the form of the Landau–Teller (LT) temperature dependence or its later generalisations. The original LT paper [Citation1] gives only the exponential factor in the expression of k10 (T). A respective rate coefficient
reads
(1)
(1)
This expression refers to the one-dimensional (1D) atom–diatom collision partners of the reduced mass µ moving under the action of a repulsive exponential potential proportional to , where R denotes collision coordinate. In LT model, diatomic molecule is represented by harmonic oscillator with frequency ω. The 3D generalisation of the 1D expression (1) was suggested by Schwartz, Slawsky and Herzfeld (SSH) [Citation2,Citation3] for isotropic exponential interaction with a model coupling between two vibrational states [Citation4,Citation5]. A further refinement of the LT model comes from the first WKB correction to the classical LT exponent [Citation6] which is important in restoring the detailed balance relation between rate coefficients for vibrational deactivation and activation. The respective expression for the rate coefficient known now as a common trajectory (CT) rate coefficient,
, can be conveniently expressed through two characteristic temperatures – the LT temperature and the vibrational (v) temperature:
(2)
(2)
In this way
(3)
(3) For this model, the temperature dependence of
is also known,
is proportional to
. It comes from the interplay of the temperature dependence of the collision number, which is proportional to
, and the temperature dependence of the effective energy width in the thermal averaging,
. Expression (3) is widely used for the interpretation of experimental data [Citation7–10].
Within this approach, the question of relation of to an anisotropic PES remains open. Of course, with a modern computational technique for scattering calculations and available PES that depend on three coordinates (the atom–diatom distance R, the angle θ between the molecular axis and the collision axis and the interatomic distance r in a diatom) the relaxation rate coefficient can be calculated from the state-to-state vibrational–rotational inelastic cross sections. However, the subsequent thermal averaging over initial translational and rotational degrees of freedom is accompanied by tremendous loss of fine dynamical details, such that one is tempted to develop simplified methods that can shed light on some particular features of the relaxation rate coefficients. One of these features is the temperature dependence of the relaxation rate coefficients and its relation to the properties of PES. Here we are encouraged by the remark by Landau and Teller [Citation1] that the main temperature dependence of the relaxation rate coefficients can be deduced from the properties of the atom–diatom interaction for the equilibrium value of r. In other words, the vibrationally inelastic event in this way is related to the vibrationally elastic motion of colliding partners across PES for a fixed value of the vibrational coordinate of the diatom. As a case study, we consider here the relaxation of O2 in Ar for which there exist a wealth of experimental data from classical works [Citation11–17] discussed at lengths in reviews [Citation8, Citation9] and supplemented also by recent study [Citation18–20]. Here, we benefit also from information on ab initio PES for this system [Citation21, Citation22]. We acknowledge also the results of full dynamical calculations of vibrational relaxation of O2 in Ar [Citation23–26]. In this respect, our study can provide simple (though quite restrictive) insight into the results of complicated numerical calculations of the relaxation rate coefficient.
2. Configuration-dependent transition probability.
In the following, we adopt the common trajectory effective mass (CTEM) approach which represents an approximation to the semiclassical effective mass (SCEM) method used earlier in the interpretation of N2 relaxation in He [Citation27]. The main difference of the CTEM approach from the standard 3D generalisation of the classical LT model (such as the SSH model) is that here the reduced mass µ is replaced by the configuration-dependent effective mass µEM and that the global isotropic repulsive potential is replaced by the local configuration-dependent repulsive potential.
In line with our previous work [Citation27], the semiclassical CTEM rate coefficient is represented as the mean probability-weighted flux through the dividing surface S* created by revolution of the generatrix
(4)
(4) The mean configuration-dependent flux
in Equation (4) corresponds to the one-dimensional relative motion of the partners along the driving mode coordinate q through the dividing surface S*. The flux
is proportional, under near-adiabatic conditions of vibrational transition, to the exponential factor that originates from the product of two exponential factors, one coming from the transition probability at a fixed energy of the driving mode and the other from the Boltzmann distribution.
In what follows, we concentrate on the calculation of the mean exponential factor that determines main features of the temperature dependence of the rate coefficient. According to [Citation27], we introduce the configuration-dependent transition probability which is parameterised by the point on the equipotential line
or equivalently by E,θ. In the CTEM approximation, the flux is proportional to the exponential factor
(5)
(5) where the exponent contains the characteristic CTEM collision time
(6)
(6)
The quantity refers to the driving mode coordinate q that starts at the point
(turning point qt) and runs in the classically forbidden region of the common potential
along
until the integral over q can be considered as converged (the upper limit qs). The driving mode is a specific combination of the coordinate displacements ΔR and Δθ that describes the motion in a small region centred at
. In this way, the two-dimensional potential
is then replaced by a one-dimensional potential
along
. As in our earlier work [Citation27], we approximate
by an exponential function of q with the characteristic slope parameter αCT defined by the gradient of
as
(7)
(7) Then, the expression in Equation (6) becomes
(8)
(8)
The effective mass µEM that enters into Equations (6) and (8) takes into accounts that motion along q near qt is a superposition of the atom–diatom radial motion and the rotational motion of the diatom with respect to the collision axis. Explicitly, the expression for µEM (E,θ) reads [Citation27]
(9)
(9) where M and re are the reduced mass of the diatom and its equilibrium interatomic distance, respectively. This expression is valid when the anisotropy of PES, characterised by the ratio
, is not too large, i.e.
. We also note that foregoing condition permits one to regard the dividing surface S* in Equation (4) as a sphere. Nonetheless, the second term in the parentheses in r.h.s. of Equation (9) can be noticeable provided the ratio
is large enough.
The thermal averaging of the flux that transforms into
is accomplished by the steepest descent (SD) approximation, yielding
(10)
(10) where the optimal energy
is found from the SD condition
(11)
(11)
The solution of Equation (11) with account taken for Equations (7)–(9) yields the function and thus the exponent in Equation (10) for a fixed θ. A minimum of
with respect to θ determines the optimal configuration R*,θ*. We also note that the equation
defines the curve
which serves as a generatrix of the ‘dividing surface’ in Equation (4) (vide supra).
For an isotropic exponential interaction, with and E being the radial collision energy, the time
coincides with the characteristic time for a collinear collision discussed by Landau and Teller [Citation1]:
(12)
(12) where the energy dependence of
appears explicitly. Then equation in Equation (11) can be solved explicitly to give
(13)
(13) with TLT defined by Equation (2).
In general case, it is convenient to represent the exponential in Equation (10) as
(14)
(14) Here TLT is introduced as a reference quantity from Equation (2), and
is a correction function. With a reasonably defined steepness parameter
, the function
is expected to be not too different from unity across the whole range of θ
and the temperature range of interest.
If the preexponential factor in the expression in Equation (12) is known, i.e. the probability is represented as
(15)
(15) it can be used in Equation (4) for calculation of
. With the exponential accuracy,
from Equation (15) can be estimated as
(16)
(16) where
(17)
(17) with
being the minimal value of
with respect to θ which is attained at θ = θ*:
(18)
(18)
If, in addition, , in the temperature range of interest, weakly depends on T, and one replaces
by a certain constant
at the reference temperature
, expression in Equation (16) becomes
(19)
(19) with
.
In a more accurate treatment, one can assume the standard θ dependence of which follows from the SSH model,
. Then, Equation (3) transformed into the form that contains standard LT exponential, assumes the form
(20)
(20) where the preexponential
includes contributions from θ and T dependence of the CTEM exponents. In this way, the factor
is written as
(21)
(21)
The factor b1(T) roughly corresponds to the range of θ, Δθ that contributes to the vibrational transition, . It is expected to have a weak positive temperature dependence due to the broadening of the effective integration range over θ. The factor
compensates the inconsistency in replacing the temperature-dependent quantity
in Equation (16) by the temperature-independent quantity
in Equation (19). If
has a positive temperature dependence (as is the case for O2–Ar relaxation, see Section 3), the factor
has a negative temperature dependence and passes through unity at
. The expression in Equation (20) formally looks like a corrected (correction factor B) SSH expression. However, now the LT temperature
is not a fitting parameter, but rather a quantity that directly related to the parameters of PES, reduced mass of the partners and the reduced mass of the diatom.
3. Application to O2–Ar system
In this section, we illustrate the general features described in Section 2 by the vibrational relaxation of O2 in Ar. For this purpose, we use ab initio PES from [Citation22] and experimental data from recent work [Citation18–20]. The contour map of the potential for O2–Ar system with equilibrium distance between O–O atoms is presented in Figure .
Figure 1. The contour map of the potential for O2–Ar system (red full lines). Equipotential contours are shown for values of
equal to 500 K (a), 1200 K (b), 2500 K (c), 5000 K (d), 10,000 K (e), 20,000 K (f) and 40,000 K (g). Dashed green line corresponds to the optimal energy
for
transition at
. Optimal configuration at
is marked by a small circle. The driving mode corresponds to the coordinate q, and its full part indicates the range
.
![Figure 1. The contour map of the potential for O2–Ar system (red full lines). Equipotential contours RC(E,θ) are shown for values of TE=E/kB equal to 500 K (a), 1200 K (b), 2500 K (c), 5000 K (d), 10,000 K (e), 20,000 K (f) and 40,000 K (g). Dashed green line corresponds to the optimal energy E∗(T,θ) for v=1→v=0 transition at T=3000K. Optimal configuration at θ∗=0.35π is marked by a small circle. The driving mode corresponds to the coordinate q, and its full part indicates the range 1/αCT(E∗,θ∗)|T=3000K.](/cms/asset/77844582-5cd1-4c0c-aeed-3f01f68ece7b/tmph_a_1938266_f0001_oc.jpg)
Shown here are several equipotential contours (red full lines) with green dashed curve for . superimposed on them. The full portion of the driving mode coordinate q qualitatively indicates the range
.
The main features of the vibrational relaxation rate coefficient for O2 at collision with Ar are summarised below in items (i)–(v) related to results of Section 2.
(i). The calculated quantities (left part) and
(right part) from Equations (7) and (9) are shown in Figure as plots vs. θ for different values of
. Since plots are symmetrical with respect to
, only halves of them are shown.
Figure 2. Plots of and
vs. θ for different values of E. The range of E shown roughly corresponds to that of the optimal energy E* (see Figures and ). All curves
collapse into a single one.
![Figure 2. Plots of αCT(E,θ) and μEM(E,θ) vs. θ for different values of E. The range of E shown roughly corresponds to that of the optimal energy E* (see Figures 3 and 4). All curves μEM(E,θ) collapse into a single one.](/cms/asset/67a8d693-1408-45f2-9182-8b4769e7ecd9/tmph_a_1938266_f0002_oc.jpg)
A noticeable increase of with decreasing E reflects the effect of weak attraction at large interpartner distances. A slight variation of
with θ (for a fixed E) is presumably due to two maxima of the electron density of O2 that come from the outer MO orbital of O2. The angular dependence of the ratio
comes from the interplay of the large value of the mean frequencies ratio for the rotation of the complex to that of the molecule,
, and the small interaction anisotropy characterised by the parameter
. A small value of this parameter,
, serves as a condition of applicability of the EM expression in the form given by Equation (9). The dependence of
on E is so weak that it remains indiscernible.
(ii). The optimal energy which determines the CTEM exponent in Equation (10) was found from numerical solutions of Equation (11) in conjunction with definitions of
and
by Equations (7) and (9). The surface plot of
vs. T and θ is shown in Figure in which the optimal configuration is indicated (see below, item iv). Sections of this surface plot for several values of θ are shown in Figure .
Figure 3. 2D plot of the scaled optimal energy for the relaxing O2–Ar system vs. T and θ. Red full line corresponds to the optimal configurations and a small circle to that in Figure .
![Figure 3. 2D plot of the scaled optimal energy for the relaxing O2–Ar system vs. T and θ. Red full line corresponds to the optimal configurations and a small circle to that in Figure 1.](/cms/asset/437edf33-a39e-450e-aacf-ff89495da87d/tmph_a_1938266_f0003_oc.jpg)
Note that the relation which is fulfilled for all graphs is the condition of applicability of the SD approximation for the calculation of the CTEM exponent. The drop of the ratio
with increasing T indicates at virtual breakdown of the SD approximation at very high temperatures.
(iii). The CTEM correction function was calculated from the CTEM exponent in Equation (10) by factoring out a reference LT temperature with
as taken from the highly repulsive portion of PES reported in [Citation22]. Plots of
vs. θ for different temperatures are shown in Figure within the range
.
For a fixed T, the global θ dependence of reflects that of
with negligible dependence of
on E. For a fixed θ, the dependence of
on T reflects increase of
with decreasing E*.
(iv). The optimal configuration of the collision complex O2–Ar for the transition at
corresponds to the coordinates
, it is marked by a circle in Figures and .
(v). Within the interval where the high precision experimental data are available, the temperature dependence of is quite weak such that the overall temperature dependence of the CTEM rate coefficient is close to that of SSH rate coefficient provided that the LT temperature, TLT, of the latter is identified with
.
The comparison of CTEM rate coefficient from Equations (20)–(21) (full line) with experimental data is presented in Figure . Previous experimental results taken from Refs. [Citation13,Citation16,Citation17] (open symbols) are shown only for temperature range where recently measured high precision data ([Citation18–20], red filled symbols) are not available.
Figure 6. Comparison of CTEM rate coefficient from Equations (20), (21) (full line) with experiment. The CTEM rate coefficient was fitted (by a scaling prefactor) to experimental data at T = 4000 K. Experimental data measured recently with high precision are shown by filled red symbols: squares [Citation18], triangles [Citation19] and circles [Citation20]. Earlier experimental results are presented by open symbols: green triangles [Citation13], blue squares [Citation16], and magenta rhombs [Citation17]. The dotted line corresponds to the case when the rotation of O2 is ignored (µEM = µ).
![Figure 6. Comparison of CTEM rate coefficient from Equations (20), (21) (full line) with experiment. The CTEM rate coefficient was fitted (by a scaling prefactor) to experimental data at T = 4000 K. Experimental data measured recently with high precision are shown by filled red symbols: squares [Citation18], triangles [Citation19] and circles [Citation20]. Earlier experimental results are presented by open symbols: green triangles [Citation13], blue squares [Citation16], and magenta rhombs [Citation17]. The dotted line corresponds to the case when the rotation of O2 is ignored (µEM = µ).](/cms/asset/f89aa84c-c54f-4057-b2b2-369babeea416/tmph_a_1938266_f0006_oc.jpg)
For this case, theoretical value of , calculated for
, is
. This would correspond to expression in Equation (2) where ratio
(i.e. µ is replaced by
) and α by
. Also shown is LT plot for the same value of α and for
, what is equivalent to
(dotted line). In this case common trajectory model is used, but rotation of molecule in vibrational relaxation event is neglected.
4. Conclusion
The common trajectory effective mass (CTEM) approach describes the collision-induced 1–0 vibrational transition in a diatom as induced by the classical translational–rotational motion of an atom–diatom system across the PES for the atom–diatom interaction at the equilibrium value of the vibrational coordinate of the diatom. The CTEM rate coefficient is calculated up to a temperature-independent scaling coefficient which incorporates the coupling of diatom vibration to the translation and rotation of the collision system. The application of the theoretical results to the vibrational relaxation of O2 in Ar allows one to rationalise the experimental temperature dependence of the rate coefficient in terms of the effective Landau–Teller temperature calculated from ab initio PES. The latter indicates at the noticeable participation of the rotation of O2 in the inducing vibrational transition which shows in the substantial decrease of the effective mass in the optimal configuration of the collision complex. This is at variance with our earlier findings for N2–He system where the rotational contribution was found to be small [Citation27]. The explanation goes back to different values of the ratio which is larger than unity for O2–Ar and smaller than unity for N2–He.
In summary, we regard the common trajectory, CT (motion along a classical trajectory running across a common potential that simulates two vibrationally adiabatic potentials originating from two vibrational states of a diatom) effective mass, EM (simulation of a two-dimensional trajectory in the space by an one-dimensional trajectory in the motion along the driving mode coordinate) approach (CTEM) as a reasonable approximation in the problem of VRT energy transfer under conditions of high temperatures and not very strong anisotropy (see comments in Section 3) for interpretation of experimental temperature dependence of the relaxation rate coefficient
in terms of the properties of the PES for an atom and a non–vibrating diatom.
Acknowledgements
This paper acknowledges many stimulating discussions with Professor J. Troe and his continuing interest in our work.
Disclosure statement
No potential conflict of interest was reported by the author(s).
References
- L. Landau and E. Teller, Phys. Z. Sowjetunion. 10, 34 (1936).
- R.N. Schwartz, Z. Slawsky and K.F. Herzfeld, J. Chem. Phys. 20, 1591 (1952). doi:https://doi.org/10.1063/1.1700221
- R.N. Schwartz and K.F. Herzfeld, J. Chem. Phys. 22, 767 (1954). doi:https://doi.org/10.1063/1.1740190
- T.L. Cottrell and J.C. McCoubrey, Molecular Energy Transfer in Gases (Butterworths, London, 1961).
- E.E. Nikitin, Theory of Elementary Atomic and Molecular Processes in Gases (Clarendon Press, Oxford, 1974).
- E.E. Nikitin and J. Troe, Phys. Chem. Chem. Phys. 8, 2012 (2006). doi:https://doi.org/10.1039/B517843F
- B.F. Gordiets, A.I. Osipov and L.A. Shelepin, Kinetic Processes in Gases and Molecular Lasers (Gordon and Breach, New York London Paris, 1988).
- E.E. Nikitin, A.I. Osipov and S.Y. Umansky, in Reviews of Plasma Chemistry, edited by B.M. Smirnov (Consultants Bureau, New York London, 1994).
- M. Capitelli, C.M. Ferreira, B.F. Gordiets and A.I. Osipov, Plasma Kinetics in Atmospheric Gases (Springer, Berlin Heidelberg, 2000).
- R. Brun, High Temperature Phenomena in Shock Waves (Springer, Berlin Heidelberg, 2012).
- M. Camac, J. Chem. Phys. 34, 448 (1961). doi:https://doi.org/10.1063/1.4757208
- K.I. Wray, J.Chem. Phys. 37, 1254 (1962). doi:https://doi.org/10.1063/1.1733273
- D.R. White and R.C. Millikan, J. Chem. Phys. 39, 1807 (1963). doi:https://doi.org/10.1063/1.1734533
- R.C. Millikan and D.R. White, J. Chem. Phys. 39, 3209 (1963). doi:https://doi.org/10.1063/1.1734182
- N.A. Generalov and S.A. Losev, J. Quan. Spectrosc. Radiat. Transfer. 6, 101–125 (1966). doi:https://doi.org/10.1016/0022-4073(66)90066-5
- C. Chiang and G.B. Skinner, J. Quan. Spectrosc. Radiat. Transfer. 24, 525–528 (1980). doi:https://doi.org/10.1016/0022-4073(80)90021-7
- V.S. Rao and G.B. Skinner, J. Chem. Phys. 81, 775–779 (1984). doi:https://doi.org/10.1063/1.447710
- J.W. Streicher, A. Krish, R.K. Hanson, K.M. Hanquist, R.S. Chaudhry and I.D. Boyd, Phys. Fluids. 32, 076103 (2020). doi:https://doi.org/10.1063/5.0012426
- J. W. Streicher, A. Krish, S. Wang, D. F. Davidson, and R. K. Hanson, AIAA Report 2019–0795, AIAA Science and Technology Forum and Exposition – Non–Equilibrium Flows I, San Diego, CA (2019).
- K.G. Owen, D.F. Davidson and R.K. Hanson, J. Thermophys. Heat Transfer. 30, 791–798 (2016). doi:https://doi.org/10.2514/1.T4505
- F.A. Gianturco and A. Storozhev, J. Chem. Phys. 101, 9624 (1994). doi:https://doi.org/10.1063/1.467927
- S.M. Cybulski, R.A. Kendal, G. Chalasinski, M.W. Severson and M. Szczesniak, J. Chem. Phys. 106, 7731 (1997). doi:https://doi.org/10.1063/1.473798
- I. S. Ulusoy, D. A. Andrienko, I. D. Boyd, and R. Hernandez, AIAA Report 2017–0661, AIAA Science and Technology Forum and Exposition – Aerothermodynamics II, 55th AIAA Aerospace Sciences Meeting, Grapevine, Texas (2017). https://doi.org/https://doi.org/10.2514/6.2017-0661
- I.S. Ulusoy, D.A. Andrienko, I.D. Boyd and R. Hernandez, J. Chem. Phys. 144, 234311 (2016). J. Chem. Phys. 145, 239902 (2016) [Erratum]. doi:https://doi.org/10.1063/1.4954041
- J.G. Kim and I.D. Boyd, Chem. Phys. 446, 76–85 (2015). doi:https://doi.org/10.1016/j.chemphys.2014.11.009
- K. M. Hanquist, R. S. Chaudhry, I. D. Boyd, J. W. Streicher, A. Krish, and R. K. Hanson, AIAA Report 2020–3275, AIAA Aviation Forum – Hypersonics, Virtual Event (2020). https://doi.org/https://doi.org/10.2514/6.2020-3275
- E.I. Dashevskaya, I. Litvin, E.E. Nikitin and J. Troe, J. Chem. Phys. 125, 154315 (2006). doi:https://doi.org/10.1063/1.2357951