Rotational excitation of C2(X1Σg +) by para- and ortho-H2.

A new four dimensional (4D) potential energy surface for the C2(X1Σg +)-H2 van der Waals system is generated. The potential was obtained from a multi-reference internally contracted configuration-interaction method including the Davidson correction (MRCI+Q). The four atoms were described using the augmented correlation-consistent quadruple zeta (aug-cc-pVQZ) basis sets. Both molecules were treated as rigid rotors. Close-coupling calculations of the inelastic integral cross sections of C2 in collisions with para-H2(j H2 = 0) and ortho-H2(j H2 = 1) were also carried out at low energies. After Boltzmann thermal averaging, rate coefficients were obtained for temperatures ranging from 5 to 100 K. The rate coefficients for collisions with ortho-H2 are significantly larger than the rate coefficients for collisions with para-H2. This journal is © The Royal Society of Chemistry.

Introduction

The C2 molecule is one of the best known constituents of the interstellar medium. It was first detected by Souza and Lutz.[1] Subsequently, many absorption lines have been detected in the near-infrared region for the A1Πu–X1Σg+ Phillips system,[2-8] and in the UV region for the D1Σu–X1Σg+ Mulliken system[9] and the F1Πu–X1Σg+ band.[10,11] The C2 molecule is used as a diagnostic probe of the physical condition in diffuse molecular clouds. The populations of its rotational levels are determined by competition between collisional excitation and radiative excitation through absorption into excited electronic states followed by fluorescence to various vibrational levels of X1Σg+.[2] The analysis of rotationally excited C2 provides useful information about the density of collision partners, the kinetic temperature and the strength of the radiation field permeating the gas.[12] In order to perform a more accurate analysis of C2 excitation, rotational rate coefficients of C2 excited by its principal interstellar collisional partner, H2, are required. In ref. 13, the full dimensional PES was calculated by the CCSD(T)-F12a/cc-pCVTZ-F12 method. However, we suppose that there is a lack of static correlation because of the multi configurational character of the PES for interacting C2 and H2. In this paper, we present a new complete four dimensional (4D) PES for the ground electronic state of the C2–H2 collisional system calculated at the MRCI+Q level. Then, this accurate PES was used to perform rotational excitation of C2 with both para-H2 and ortho-H2.

The activation energy of the C2 + H2 ⇌ C2H + H reaction being about 900 cm−1,[14] inelastic cross sections of C2 in collisions with H2(jH = 0) and H2(jH = 1) are reported for total energies of up to 1000 cm−1 for the first 9 even rotational levels, yielding converged rate coefficients up to 100 K.

Excitation rate coefficients of C2 excited by para-H2(jH = 0) were provided by Najar et al.[15] In this study, the PES was derived by averaging only over three specific H2 orientations. These calculations gave a first estimate of C2–para-H2 collisional rate coefficients. However, rotational excitation rate coefficients for C2 in collision with ortho-H2 are also needed for a complete description of C2 excitation.[16]

The paper is organized as follows: Section 2 describes the calculation of the potential energy surface, Section 3 contains a concise description of the scattering calculations, and in Section 4, we present and discuss our results.

Potential energy surface

As mentioned in previous work,[17] the ground electronic state of the C2–H2 system is described by two dominant configurations. Because of this multi configurational character, the potential energy surface (PES) of this van der Waals complex was accomplished by first performing full valence complete active space (CASSCF) calculations[18,19] followed by multi-reference internally contracted configuration-interaction (MRCI) calculations,[20,21] including all single and double excitations from the CASSCF reference function. The contributions of higher-order excitations were estimated using the normalized multi-reference Davidson correction.[22,23] In the CASSCF-MRCI calculations, the active space is defined by distributing 10 electrons in 10 molecular orbitals. The 1s2 orbitals of the carbon atoms were kept frozen. To achieve a good description of weak interactions, the calculations were performed using a rather large aVQZ basis set of Dunning.[24-26] All ab initio calculations were carried out using the MOLPRO program package.[27]

The body-fixed coordinate system presented in Fig. 1 was used. The center of the coordinates is placed in the middle of the bond of the C2 molecule. Vector R is directed along the z-axis and connects the C2 and H2 centers of mass. The rotation of the C2 molecule is defined by angle θC and the rotation of the H2 molecule is defined by angle θH, while ϕ is the dihedral angle.

Fig. 1

Geometry describing the C2–H2 complex.

The ab initio calculations were performed using the approximation of rigid interacting molecules with the H2 bond length corresponding to the ground vibrational state: rH–H = 1.449 Bohr.[28] For C2, we used the experimental equilibrium distance rC–C = 2.348 Bohr,[29] which differs by less than 1% from the calculated value obtained by averaging over the ν = 0 wavefunction. We note that employing state-averaged geometries is a reliable approximation for including zero-point vibrational effects within a rigid-rotor PES.[30-32]

In order to determine the interaction potential V(R,θC,θH,ϕ), the basis set superposition error (BSSE) was corrected for all geometries with the Boys and Bernardi counterpoise scheme:[33]where the energies of the C2 and H2 subsystems are computed using a full basis set of the complex. The size inconsistency was corrected for all geometries by adding the energy −33.75 cm−1.

At each point R the calculations were carried out for 122 angular orientations with θC varying from 0° to 180°, θH varying from 0° to 90°, and ϕ varying from 0° to 90°. R-distances were varied from 4.8 to 20 Bohr, giving 29 grid points.

For the solution of the close-coupling scattering equations, it is most convenient to expand the interaction potential V(R,θC,θH,ϕ) into the angular functions at each value of R. For the scattering of two linear rigid rotors, we used the expansion of Green:[34]

The basis functions A(θC,θH,ϕ) are products of the associated Legendre polynomials P:where 〈…∣…〉 is a Clebsch–Gordan coefficient. The P functions are related to coupled spherical harmonics through Y(θ,ϕ) = P(θ)e(im. Here l1 and l2 are associated respectively with the rotational motion of C2 and H2, and l = |l1 − l2|, …, |l1 + l2|. In eqn (2), the homonuclear symmetry of the interacting molecules forces the indexes l1 and l2 to be even.

To obtain the expansion described above, at each point of the radial grid R we fitted the ab initio points to expression (2) using a linear least squares method, with lmax1 = 6 and lmax2 = 4. The choice of such expansion parameters limits the results to a total of 28 radial expansion coefficients v(R). The root-mean-square error (rmse) of the fits in the long and (attractive) medium range was of the order of 10−3 to 10−2 cm−1. The medium repulsive wall (energies of the order of 102 to 103 cm−1) was reproduced with rmse of 1–10 cm−1.

Each radial coefficient was interpolated using cubic splines for 4.8 < R ≤ 16a0. For R > 16a0 the energy (not the coefficients) was extrapolated using the c5/R5 + c6/R6 functional form based on its values at R = 16a0 and R = 18a0. For R < 4.8a0 the energy was extrapolated using the functional form A exp(BR) based on its values at R = 4.8a0 and R = 5a0. Thus, our PES covers the whole range of R. However, we suggest to use it in the range from 3.0 to 100a0.

The global minimum of the 4D PES corresponds to a T-shaped structure with θC = 0°, θH = 90°, ϕ = 0° and R = 6.92a0 with V = −109.22 cm−1. In Fig. 2–4 we present the contour plots of the interaction energy surface for fixed equilibrium parameters. It is seen that there is a weak anisotropy of the PES for the Jacoby angle ϕ and a strong anisotropy for the angle θC. The local minimum of the PES occurs for another T-shaped structure with θC = 90°, θH = 0°, ϕ = 0° and R = 6.70a0. The interaction energy at this point is V = −101.63 cm−1.

Scattering study of excitation of C2 by para- and ortho-H2

In the present work, we consider collisions of 12C2 with para- and ortho-H2. Only even rotational levels j were considered here because, in homonuclear molecules with zero nuclear spins like C2, rotational levels with odd values are not allowed. On the other hand, the H2 molecule can exist in two forms denoted para- and ortho-H2. In the former, the nuclear spins of the hydrogen atoms are opposed resulting in total I = 0, whereas, in the latter, the nuclear spins are aligned and I = 1. As the nuclear wavefunction must be asymmetric, in para-H2, only even rotational states, jH = 0, 2,… are allowed, and in ortho-H2, only odd values, jH = 1, 3,… are allowed.

Full close-coupling calculations[35] were carried out for total energies up to 1000 cm−1, using C2 and H2 as rigid rotors with rotational and centrifugal distortion constants BC = 1.820 cm−1 and DC = 6.92 × 10−6 cm−1 (ref. 36), and BH = 60.853 cm−1 and DH = 3.3 × 10−6 cm−1,[29] respectively. Transitions among C2 rotational levels up to j = 16 were computed for collisions involving para-H2(jH = 0) and ortho-H2(jH = 1). For convergence of the cross sections, the rotational states jH = 0, 2 and jH = 1, 3 were included in the rotational basis sets for para- and ortho-H2, respectively. The inclusion of the jH = 4 channel in para-H2 affects the cross sections by less than 1 percent. For the C2 molecule, the rotational basis set was extended up to jmax = 24 at the largest total energy considered (1000 cm−1). The coupled scattering equations were solved by the MOLSCAT code[37] using the propagator of Manolopoulos.[38] For the C2–H2 system the reduced mass is μ = 1.8597 amu. The integration parameters were chosen to ensure convergence of the cross sections in this range. The minimum and maximum integration distances chosen were Rmin = 3 Bohr and Rmax = 50 Bohr, respectively. The maximum value of the total angular momentum J used in the calculations was chosen according to a convergence criterion of 0.01 Å2. The energy grid was adjusted to describe all the details of the resonances of the inelastic cross sections. The thermal rate coefficients at temperature T are obtained by averaging the rotational inelastic cross sections, σ, over a Boltzmann distribution of collision energy Ec:with , where kB is the Boltzmann constant.

Results

First we compare the σ0→2 cross sections for the excitation of C2 with para-H2(jH = 0, 2) calculated with the 4D PES and including jH = 0, 2 in the rotational basis set of H2 with the cross sections for the excitation of C2 with para-H2 calculated previously by Najar et al.[15] with a 2D effective PES averaged over three specific H2 orientations using equal weights of 1/3 for the three geometries considered. As we can see in Fig. 5, the two cross sections are close, having the same form of resonances with a shift and small difference in amplitudes. The difference is appreciable at high energies, approaching 35 percent. To understand where the difference comes from, we present in Table 1 cross sections at different energies for C2-paraH2 calculations with a PES averaged over 3 orientations, C2-paraH2 4D PES calculations with just the basis jH = 0 and C2-paraH2 4D PES calculations with the rotational basis involving jH = 0 and jH = 2 levels (jH = 0, 2). We note that at low energies the major difference (∼15%) is due to lack of anisotropy in the PES averaged over 3 orientations. At high energies the inaccuracies of the PES have only a minor influence on the cross sections, and the differences stem mainly from neglect of the jH = 2 rotational state. This result was also confirmed by Lique et al.[39] in the SiS–H2 system.

Fig. 5

Cross sections as a function of collision energy for 0 → 2 rotational excitation of C2 with para-H2 obtained by a 2D effective PES (red curve) and by the 4D PES including (jH = 2) in the rotational basis set of H2 (black curve).

Cross sections σ0→2 (Å2) for C2–paraH2 obtained using a 2D effective PES, a 4D PES with basis jH = 0 and a 4D PES with the basis jH = 0, 2

Energy (cm−1)	2D effective PES	4D PES, jH2 = 0	4D PES, jH2 = 0, 2
50	11.249	13.307	14.160
100	5.613	6.580	7.087
300	2.622	2.650	2.515
500	2.978	2.819	2.379
1000	3.540	3.298	2.567

Fig. 6 shows the collision energy dependence of the 0 → 2 inelastic integral cross sections of C2 in collision with para-H2(jH = 0) (including jH = 0, 2) and ortho-H2(jH = 1) (including jH = 1, 3). The figure illustrates many resonances in the cross sections for collision energies below 100 cm−1, which are a consequence of the temporary trapping of the H2 molecule in the attractive potential well, involving the formation of quasi-bound states of the van der Waals C2–H2 complex. We note a difference in amplitude and shape of the resonances between the two cross sections. The collisions with ortho-H2 are more strongly inelastic. This is expected because, unlike with ortho-H2, the long-range interactions with para-H2(jH = 0) do not involve the permanent quadrupole moment of the H2 molecule; only induction and dispersion forces participate.[40] This behavior was also found in H2–H2 (ref. 41) and SO2–H2 (ref. 42) collisions.

Fig. 6

Cross sections for rotational excitation of C2 by para-H2(jH = 0, 2) and ortho-H2(jH = 1, 3) for the 0 → 2 transition as a function of collision energy.

Rate coefficients for rotational excitation of C2 by para-H2(jH = 0) and ortho-H2(jH = 1) have been computed from cross sections generated from the grid of collision energies. The representative variation with temperature is illustrated in Fig. 7. The curves in Fig. 7 and the results in Table 2 show that the rates for ortho-H2(jH = 1) are significantly larger by a factor of 5 on average. This result, also found previously for some systems like CO–H2,[43] H2O–H2 (ref. 44) and SiS–H2,[45] confirms the expectation that ortho-H2(jH = 1) is a more efficient perturbant than para-H2(jH = 0).

Fig. 7

Temperature dependence of rotational de-excitation rate coefficients of C2 excited by H2 for Δj = 2 transitions. Dashed lines: para-H2(jH = 0); solid lines: ortho-H2(jH = 1).

Rate coefficients (cm3 s−1) for rotational de-excitation in C2 induced by para-H2(jH = 0) and ortho-H2(jH = 1). The upper entry for each transition is for T = 50 K, and the lower entry is for T = 100 K

Transition	Para-H2	Ortho-H2
2 → 0	0.234(-10)	1.122(-10)
	0.188(-10)	1.100(-10)
4 → 2	0.322(-10)	1.191(-10)
	0.273(-10)	1.216(-10)
6 → 4	0.192(-10)	0.870(-10)
	0.201(-10)	0.949(-10)
8 → 6	0.110(-10)	0.610(-10)
	0.141(-10)	0.704(-10)

In the paper of Lavendy et al.,[46] for construction of the PES the authors considered only 5 relative orientations of C2 and H2 for each R value. Later, Phillips[47] used this PES to compute the rate coefficients using the close-coupling approach. In our previous work,[15] we constructed a 2D PES based on nearly 900 ab initio points. For each (R,θ) we considered 3 angular orientations of the H2 molecule, the energies of which were then averaged. In that work,[15] we compared our collision rates for C2–para-H2(jH = 0) to those of Phillips[47] for the temperature range T < 100 K. We found that the rates calculated by Phillips are greater by a factor of 5. This difference is due to the fact that the interaction potential used by Phillips overestimates the anisotropy of the potential. In addition, Phillips found that para-H2(jH = 0) and ortho-H2(jH = 1) are nearly identical colliders, which doesn’t agree with our findings that there is a notable difference. This can also be explained by the fact that the potential of Phillips based on the ab initio energies of Lavendy et al. was constructed by neglecting the terms responsible for the quadrupole–quadrupole interactions.

Conclusions

We have applied full close-coupling calculations to investigate rotational excitation in the C2 molecule for collisions involving para-H2(jH = 0) and ortho-H2(jH = 1). The calculations are performed using a highly accurate 4D potential energy surface. The inclusion of jH = 2 in the rotational basis set of H2 makes an appreciable contribution at high energies to the para-H2 cross sections, approaching 35 percent. Rate coefficients for transitions among the levels of C2 up to j = 16 were determined for the range of temperatures 5 ≤ T ≤ 100 K. Rate coefficients with ortho-H2 are larger by a factor of 5 on average than rate coefficients with para-H2.

Conflicts of interest

There are no conflicts to declare.