886
Views
0
CrossRef citations to date
0
Altmetric
Jürgen Troe Special Issue

Generalised Landau–Teller fitting of the vibrational relaxation rates in atom–diatom collisions: O2–Ar system as a case study

, ORCID Icon &
Article: e1938266 | Received 02 Mar 2021, Accepted 25 May 2021, Published online: 18 Jun 2021

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

1. Introduction

The experimental values of vibration relaxation rate coefficients k10 (T) for v=1v=0 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 k10LT(T) reads (1) k10LT(T)exp324π2ω2μαLT2kBT1/3(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 exp(αLTR), 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, k10CT(T), can be conveniently expressed through two characteristic temperatures – the LT temperature and the vibrational (v) temperature: (2) TLT=4π2ω2μ/αLT2kB,Tv=ω/kB(2)

In this way (3) k10CT(T)C10SSH(T)exp32TLTT1/3+Tv2T(3) For this model, the temperature dependence of C10SSH(T) is also known, C10SSH(T) is proportional to T1/3. It comes from the interplay of the temperature dependence of the collision number, which is proportional to T1/2, and the temperature dependence of the effective energy width in the thermal averaging, T1/6. Expression (3) is widely used for the interpretation of experimental data [Citation7–10].

Within this approach, the question of relation of TLT 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 F10CTEM through the dividing surface S* created by revolution of the generatrix R=R(T,θ) (4) k10CTEM(T)=SF10CTEM(T,θ)dS,dS(θ)=2πR(R)2+(dR/dθ)2R=R(T,θ)dθ(4) The mean configuration-dependent flux F10CTEM(T,θ) 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 F10CTEM(T,θ) 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 RC(E,θ),θ on the equipotential line Uc(R,θ)=E or equivalently by E,θ. In the CTEM approximation, the flux is proportional to the exponential factor (5) F10CTEM(E,θ)exp(2ωτCTEM(E,θ))(5) where the exponent contains the characteristic CTEM collision time τCTEM(E,θ) (6) τCTEM(E,θ)=12qt(E,θ)qs2μEM(E,θ)dqUCT(q,θ)E(6)

The quantity τCTEM(E,θ) refers to the driving mode coordinate q that starts at the point RC(E,θ),θ (turning point qt) and runs in the classically forbidden region of the common potential UCT(q;θ) along gradU(R,θ) 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 RC(E,θ),θ. In this way, the two-dimensional potential U(R+ΔR,θ+Δθ) is then replaced by a one-dimensional potential UCT(q;θ) along gradU(R,θ). As in our earlier work [Citation27], we approximate UCT(q;θ) by an exponential function of q with the characteristic slope parameter αCT defined by the gradient of U(R,θ) as (7) αCT(E,θ)=|gradU(R,θ)|RC(E,θ),θ/E(7) Then, the expression in Equation (6) becomes (8) τCTEM(E,θ)=πμEM(E,θ)/αCT(E,θ)2E(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) μEM(E,θ)μ=1+μMre2RCθ21(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 |RC/RCθ|, is not too large, i.e. |RC/RCθ|<<1. 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 μRC2/Mre2 is large enough.

The thermal averaging of the flux that transforms F10CTEM(E,θ) into F10CTEM(T,θ) is accomplished by the steepest descent (SD) approximation, yielding (10) P¯10CTEM(T,θ)exp(2ωτCTEM(E,θ)E/kBT)(10) where the optimal energy E=E(T,θ) is found from the SD condition (11) 2ωdτCTEM(E,θ)/dE|E=E+1/kBT=0(11)

The solution of Equation (11) with account taken for Equations (7)–(9) yields the function E=E(T,θ) and thus the exponent in Equation (10) for a fixed θ. A minimum of E(T,θ) with respect to θ determines the optimal configuration R*,θ*. We also note that the equation E(T,θ)=U(R,θ) defines the curve R=R(T,θ) which serves as a generatrix of the ‘dividing surface’ in Equation (4) (vide supra).

For an isotropic exponential interaction, with U(R,θ)=U(R)exp(α¯R) and E being the radial collision energy, the time τCTEM coincides with the characteristic time for a collinear collision discussed by Landau and Teller [Citation1]: (12) τLT=πμ/α¯2E(12) where the energy dependence of τLT appears explicitly. Then equation in Equation (11) can be solved explicitly to give (13) ELT=kBTLT=kBT2TLTT1/3(13) with TLT defined by Equation (2).

In general case, it is convenient to represent the exponential in Equation (10) as (14) P¯10CTEM(T,θ)exp32TLTTΘ(T,θ)1/3(14) Here TLT is introduced as a reference quantity from Equation (2), and Θ(T,θ) is a correction function. With a reasonably defined steepness parameter α¯, the function Θ(T,θ) is expected to be not too different from unity across the whole range of θ (0θπ) and the temperature range of interest.

If the preexponential factor in the expression in Equation (12) is known, i.e. the probability P¯10CTEM(T,θ) is represented as (15) P¯10CTEM(T,θ)=A10CTEM(T,θ)exp32TLTTΘ(T,θ)1/3+Tv2T(15) it can be used in Equation (4) for calculation of k10CTEM(T). With the exponential accuracy, k10CTEM(T) from Equation (15) can be estimated as (16) k10CTEM(T)exp32TLTT1/3+Tv2T,(16) where (17) TLT=TLTΘ(T)(17) with Θ(T,θ) being the minimal value of Θ(T,θ) with respect to θ which is attained at θ = θ*: (18) Θ(T)=minθ{Θ(T,θ)}=Θ(T,θ)|θ=θ(18)

If, in addition, Θ(T), in the temperature range of interest, weakly depends on T, and one replaces Θ(T) by a certain constant Θ¯=Θ(T¯) at the reference temperature T¯, expression in Equation (16) becomes (19) k10CTEM(T)k¯10LT(T)exp32T¯LTT1/3+Tv2T(19) with T¯LT=TLTΘ¯.

In a more accurate treatment, one can assume the standard θ dependence of A10CTEM(T,θ) which follows from the SSH model, A10CTEM(T,θ)cos2θ. Then, Equation (3) transformed into the form that contains standard LT exponential, assumes the form (20) k10CTEM(T)=B(T,T¯)C10LT(T)exp32T¯LTT1/3+Tv2T(20) where the preexponential B(T,T¯) includes contributions from θ and T dependence of the CTEM exponents. In this way, the factor B(T,T¯) is written as (21) B(T,T¯)=b1(T)b2(T,T¯),b1(T)=0π(3/2)cos2θ×exp32TLT(Θ(T,θ)Θ(T))T1/3sinθdθ,b2(T,T¯)=exp32TLT[Θ(T)Θ¯]T1/3(21)

The factor b1(T) roughly corresponds to the range of θ, Δθ that contributes to the vibrational transition, Δθ/π<1. It is expected to have a weak positive temperature dependence due to the broadening of the effective integration range over θ. The factor b2(T,T¯) compensates the inconsistency in replacing the temperature-dependent quantity TLT in Equation (16) by the temperature-independent quantity T¯LT in Equation (19). If Θ(T) has a positive temperature dependence (as is the case for O2–Ar relaxation, see Section 3), the factor b2(T,T¯) has a negative temperature dependence and passes through unity at T=T¯. The expression in Equation (20) formally looks like a corrected (correction factor B) SSH expression. However, now the LT temperature T¯LT 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 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=1v=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.

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.

Shown here are several equipotential contours (red full lines) with green dashed curve for E=E(T,θ)|T=3000K. superimposed on them. The full portion of the driving mode coordinate q qualitatively indicates the range 1/αCT(E,θ)|T=3000K.

The main features of the vibrational relaxation rate coefficient k10CTEM(T) for O2 at collision with Ar are summarised below in items (i)–(v) related to results of Section 2.

(i). The calculated quantities μEM(E,θ)/μ (left part) and αCT(E,θ) (right part) from Equations (7) and (9) are shown in Figure  as plots vs. θ for different values of E. Since plots are symmetrical with respect to θ=π/2, only halves of them are shown.

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  and ). All curves μEM(E,θ) 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.

A noticeable increase of αCT(E,θ) with decreasing E reflects the effect of weak attraction at large interpartner distances. A slight variation of αCT(E,θ) 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 μEM/μ comes from the interplay of the large value of the mean frequencies ratio for the rotation of the complex to that of the molecule, μR/Mre, and the small interaction anisotropy characterised by the parameter |dR/Rdθ|. A small value of this parameter, |dR/Rdθ|<<1, serves as a condition of applicability of the EM expression in the form given by Equation (9). The dependence of μEM(E,θ) on E is so weak that it remains indiscernible.

(ii). The optimal energy E=E(T,θ) which determines the CTEM exponent in Equation (10) was found from numerical solutions of Equation (11) in conjunction with definitions of αCT(E,θ) and μEM(E,θ) by Equations (7) and (9). The surface plot of E(T,θ) 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.

Figure 4. Plots of the optimal energies E* vs. T for different angles θ.

Figure 4. Plots of the optimal energies E* vs. T for different angles θ.

Note that the relation E>>kBT 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 E(T,θ)/kB with increasing T indicates at virtual breakdown of the SD approximation at very high temperatures.

(iii). The CTEM correction function Θ(T,θ) was calculated from the CTEM exponent in Equation (10) by factoring out a reference LT temperature with α¯=1.9a01 as taken from the highly repulsive portion of PES reported in [Citation22]. Plots of Θ(T,θ) vs. θ for different temperatures are shown in Figure  within the range 0θπ/2.

Figure 5. Plots of the correction CTEM function Θ(T,θ) vs. θ for different temperatures.

Figure 5. Plots of the correction CTEM function Θ(T,θ) vs. θ for different temperatures.

For a fixed T, the global θ dependence of Θ(T,θ) reflects that of μEM(E,θ)/μ with negligible dependence of μEM(E,θ) on E. For a fixed θ, the dependence of Θ(T,θ) on T reflects increase of αCT(E,θ) with decreasing E*.

(iv). The optimal configuration of the collision complex O2–Ar for the v=1v=0 transition at T=3000K corresponds to the coordinates R=R(T,θ)=4.27a0,θ=θ=0.35π, 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 B(T,T¯) 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 T¯LT.

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 = µ).

For this case, theoretical value of T¯LT, calculated for T¯=2000K, is 3.3106K. This would correspond to expression in Equation (2) where ratio μEM/μ=0.7 (i.e. µ is replaced by μ~=0.7μ=12.5a.m.u.) and α by α~=2.10a01. Also shown is LT plot for the same value of α and for μEM/μ=1.0, what is equivalent to T¯LT=4.7106K (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 μ/M 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 R,θ 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 k10(T) 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