# Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae

Francesca Tonolo, Ernesto Quintas-Sánchez, Adrian Batista-Planas, Richard

Dawes, François Lique

► To cite this version:

Francesca Tonolo, Ernesto Quintas-Sánchez, Adrian Batista-Planas, Richard Dawes, François Lique. Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae. Journal of Physical Chemistry A, 2025, 129 (41), pp.9583-9590. 10.1021/acs.jpca.5c05448. http://dx.doi.org/10.1021/jp05448a

HAL Id: hal-05296021

https://hal.science/hal-05296021v1

Submitted on 18 Feb 2026

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
Collisional excitation of HCN by CO to refine the modeling
of cometary comae

Francesca Tonolo,*,† Ernesto Quintas-Sánchez,‡ Adrian Batista-Planas,‡ Richard Dawes,‡ and
François Lique*,†

†Univ. Rennes, CNRS, IPR (Institut de Physique de Rennes), UMR 6251, Rennes, F-35000,
France

‡Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri
65409, USA

E-mail: francesca.tonolo.1@univ-rennes.fr; francois.lique@univ-rennes.fr

Abstract

We present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5–50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called “thermalized” rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

Introduction

Cometary ices are among the most pristine materi-
als of the Solar System, as they preserve the chem-
ical conditions that existed in its earliest stages. ¹
Thanks to recent advances in observational tech-
niques, we can now observe with very high spa-
tial and spectral resolutions the molecular con-

tent of cometary comae, i.e., the ensemble of gases
which sublimate from cometary ices by solar radia-
tion, forming an atmosphere around the cometary
nucleus. 2–5 However, the accurate interpretation
of these spectra requires sophisticated radiative
transfer models. A key challenge in this regard
is describing how a molecule’s energy is distributed
among its rotational levels. At the typical low den-
sities of cometary comae (1 < *n* < 10⁻¹⁰ cm⁻³),
this distribution can strongly deviate from the lo-
cal thermodynamic equilibrium (LTE). In these
cases, accurate datasets of collisional coefficients
with the most abundant perturbing species are cru-
cial (as an example, the reader is referred to Loreau
et al.⁶), and usually rely on state-of-the-art theo-
retical calculations.

However, when dealing with heavy projectiles
such as H₂O, CO and CO₂ (the most abundant
gases in cometary comae), collisional calculations
notably complicate: the high densities of states of
the projectiles dramatically increase the number of
channels to be considered, which makes full quan-
tum scattering calculations unaffordable in terms
of computational time and memory requirements.
To address this, the statistical adiabatic channel
model (SACM) was recently tested as an alterna-
tive. Developed by Quack and Troe⁷ and subse-
quently refined by Loreau et al.⁸, this approach
can generate reliable sets of collisional coefficients
in favorable cases at significantly reduced compu-
tational cost. In terms of accuracy, the SACM is
particularly effective for systems that form a stable
intermediate complex, often showing discrepancies
of less than a factor 2 with respect to the exact
close coupling (CC) calculations.⁸⁻¹¹ For cometary applications, the SACM currently offers the most practical solution to overcome the prohibitive costs associated with heavy projectiles and to deal with the low temperatures (~5-50 K) that rule out the application of (quasi)classical methods *. For the CS-CO collisional system, it was recently shown that it can produce state-to-state outcomes differing by less than an order of magnitude from full quantum calculations.¹⁴

In order to refine the modeling of cometary comae, we address in this work the collisional excitation of HCN by CO. HCN is one of the most abundant nitrogen-bearing molecules in cometary comae.¹⁵ Detected in both millimeter and IR wavelengths, HCN primarily originates from direct sublimation of ices on the comet nucleus. Hence, it serves as a valuable tracer of the thermal and chemical history of pre-cometary material.¹⁶,¹⁷ Understanding the abundance of HCN in comets compared to that in the interstellar medium provides insights into the early conditions of the Solar System and the possible delivery of prebiotic molecules to Earth.¹⁸,¹⁹ Moreover, the abundance of HCN relative to other nitrogen-bearing species (e.g., NH₃, N₂) can indicate the temperature at which cometary ices formed and the extent and nature of the chemical process they underwent.²⁰⁻²² However, the only collisional investigations of cometary interest regarding HCN concern its excitation processes induced by H₂O.²³,²⁴ These data are not sufficient to model cometary comae at intermediate and large heliocentric distances, where the solar radiation is not intense enough to sublimate the water from the nucleus. Distant cometary comae are instead generally dominated by CO.²⁵⁻²⁷ These comets are particularly relevant for study of the composition of the early Solar System because they typically contain materials not previously altered by solar radiation. Moreover, the low densities that characterize cometary comae lead to strong deviations from LTE, making accurate collisional rate coefficients essential for reliable radiative transfer modeling.

This work aims to bridge this gap by providing the first dataset of collisional coefficients for the rotational (de)-excitation of HCN by CO. The methodology and results are presented as follows:

<figure><img src="image_2.png" /><figcaption>Figure 1: Coordinates used to describe the HCN–CO interaction. See the text for details.</figcaption></figure>

the first Section describes the employed computational strategy. There, in the context of the Born-Oppenheimer approximation, the characterization of the interaction potential and the scattering calculations used to solve the nuclear Schrödinger equation are presented in two different subsections. The subsequent Section presents and discusses the results. Finally, in the last Section the main conclusions are drawn.

## Computational Details

### Potential Energy Surface

The Jacobi coordinates used to represent the four-dimensional (4D) HCN–CO interactions are depicted in Figure 1. Note that for the representation of the PES, CO is monomer 1, and HCN is monomer 2. $\vec{R}$ is the vector between the centers of mass of the two fragments, and $\vec{r}_1$ and $\vec{r}_2$ are vectors aligned with each molecule. Coordinate $R$ is the length of vector $\vec{R}$, while coordinates $\theta_1$ and $\theta_2$ represent (respectively) the angles between $\vec{R}$ and the vectors $\vec{r}_1$ and $\vec{r}_2$. The fourth coordinate is the dihedral (out of plane) torsional angle, labeled $\phi$, which is the angle between the vectors $\vec{R} \times \vec{r}_1$ and $\vec{R} \times \vec{r}_2$.

In line with previous works on other van der Waals (vdW) dimers composed of two linear molecules,²⁸⁻³⁶ the analytical representation of the PES was constructed using an automated interpolating moving least squares (IMLS) methodology, freely available as a software package under the name AUTOSURF.³⁷ The fitting basis and other aspects of the IMLS procedure are similar to those employed for other previously treated

<footnote>*The mixed quantum/classical theory (MQCT) offers a promising compromise in this direction, but it currently remains insufficient to accurately characterize collisional processes at temperatures below 50 K.¹²,¹³</footnote>
<figure><img src="image_3.png" /><figcaption>Figure 2: R-optimized contour plot of the PES as a function of the extended angles $\hat{\theta}_1$ and $\hat{\theta}_2$ for planar configurations ($\varphi = 0^\circ$ and $\varphi = 180^\circ$). For each pair of angles, the energy (given in cm$^{-1}$) is optimized with respect to the center-of-mass distance R. See text for details.</figcaption></figure>

systems that have been described in detail elsewhere.$^{37-40}$ The shortest intermonomer center-of-mass distance considered is $R = 2.4$ Å. The short-range part of the PES was restricted by excluding regions with repulsive energies above a maximum of 6 kcal/mol ($\sim 2000$ cm$^{-1}$) relative to the separated monomers asymptote. The *ab initio* data coverage in the fitted PES extended to $R = 25$ Å, while the zero of energy was set at infinite center-of-mass separation between the monomers. To guide the placement of high-level *ab initio* data, a lower-level guide surface was constructed using 3000 points at the explicitly correlated CCSD(T)-F12a/VDZ-F12 level,$^{41}$ distributed using a Sobol sequence$^{42}$ biased to sample the short range region more densely. Since the two fragments are not identical, and each has only $C_{\infty v}$ symmetry, there is no additional symmetry to respect/exploit in the PES beyond that of the angle $\varphi$, for which energies were only computed in the reduced angular range: $0 < \varphi < 180^\circ$, as is the case in this treatment for all 4D systems.

The final high-level PES was computed using the CCSD(T)-F12b method, extrapolated to the complete basis set (CBS) limit. The CBS extrapolation was performed using the VTZ-F12 and VQZ-F12 bases$^{43}$ and the $l^{-3}$ formula.$^{44}$ All *ab initio* calculations were performed using the Molpro electronic structure code package.$^{45}$

<table><caption>Table 1: Geometric parameters and potential energy for stable structures in the PES. Energies are given relative to the asymptote. Units are Å, degrees, and cm$^{-1}$.</caption><thead><tr><th></th><th>GM <sup>1</sup></th><th>LM <sup>1</sup></th></tr></thead><tbody><tr><td><i>R</i></td><td>4.826</td><td>4.532</td></tr><tr><td><i>θ</i><sub>1</sub></td><td>180.0</td><td>0.0</td></tr><tr><td><i>θ</i><sub>2</sub></td><td>0.0</td><td>0.0</td></tr><tr><td><i>ϕ</i></td><td>—</td><td>—</td></tr><tr><td><i>V</i></td><td>-604.13</td><td>-363.32</td></tr></tbody></table>

¹ Defined in the text.

For construction of the PES, both monomers were held rigid. The bond distance for CO was fixed at $r_{CO} = 1.128206$ Å, the vibrationally-averaged bond distance for the ground rotational state of CO, consistent with its rotational constant of $B(CO) = 1.9225$ cm$^{-1}$.$^{46}$ For HCN, the equilibrium geometry is linear and the structural parameters employed in this study ($r_{HC} = 1.0655$ Å and $r_{CN} = 2.2187$ Å) are consistent with the experimental rotational constant, $B(HCN) = 1.4782$ cm$^{-1}$.$^{47}$ Masses of 15.9949146221, 12, 1.007825032, and 14.0030740052 u were used for $^{16}$O, $^{12}$C, $^{1}$H and $^{14}$N, respectively.

To represent the long range (out to arbitrary distances), the PES switches smoothly to an analytic expression representing electrostatic, induction, and dispersion interactions between the fragments (truncated at 8th order). For the high-level PES, 3841 points were used, resulting in an estimated global root-mean-squared fitting error of 0.8 cm$^{-1}$ (excluding the long range where the fitting error is extremely small). The estimated error is 0.1 cm$^{-1}$ for geometries corresponding to energies in the wells (below the asymptotic zero). The analytical representation of the PES is available from the authors upon request.

Figure 3 shows radial cuts through the PES for nine different angular orientations. The points represent *ab initio* data (not included in the fit), while the lines plot the fitted PES. This helps one appreciate the anisotropy of the PES as reflected in the different values of $R$ at which the onset of a steep repulsive wall begins, as well as the accuracy of the fit.

The PES is characterized by two minima, the global minimum (GM), a collinear arrangement
<figure><img src="image_4.png" /><figcaption>Figure 3: Various radial cuts defined by different orientational orientations of the monomers. In all cases, energies are in cm⁻¹, lines represent the fitted PES, and points represent *ab initio* calculations (not used in the fit).</figcaption></figure>

where the C-atom of CO approaches the H-atom of HCN, and a local minimum (LM), also collinear, where the O-atom of CO approaches the H-atom of HCN. These minima are seen across the middle of Figure 2, with GM appearing twice in the plot due to the use of extended angles. The energies and structural parameters are given in Table 1 where it is seen that GM (*E* = −604.1 cm⁻¹) is much more stable than LM (*E* = −363.3 cm⁻¹), although it is LM with a significantly shorter separation between monomers, due to the location of the CO monomer’s center-of-mass.

## Scattering Calculations

The constructed state-of-the-art interaction potential was subsequently implemented in the scattering equations to solve the nuclear part of the Schrödinger equation and characterize the dynamics of the collisional system.

The most accurate method to solve scattering equations and to obtain the so-called scattering matrix, which contains all the collisional observables of the system, is the full quantum CC approach. However, for systems characterized by high densities of rotational states, combined with (a) deep well(s) of the intermediate complex (such as HCN interacting with CO), the number of coupled equations required for convergence becomes

computationally prohibitive.

To address this, as introduced above, we decided to adopt the SACM approach. This approximation requires to compute only the adiabatic curves that are obtained by excluding the nuclear kinetic term and diagonalizing the interaction Hamiltonian in a basis of rotational functions for the two colliders. This noticeably simplifies the scattering calculations, as each adiabatic curve is independent of the energy and hence needs to be computed only once for each value of $J_{\text{tot}}$. The calculations were performed by using the HIBRIDON code,⁴⁸ which contains an implementation that permits one to retrieve only the adiabatic states. The scattering matrices were then rebuilt by including all the open adiabatic channels at each total energy, assuming that they occur with the same probability. This yields to state-to-state collisional cross sections, which can be subsequently averaged over a Maxwell-Boltzmann distribution of energies to obtain the corresponding rate coefficients.

As mentioned in the Introduction, the SACM method assumes that the intermediate complex formed during each collision lives long enough for energies to be statistically redistributed. This is a good approximation for the HCN and CO system, where the hydrogen bonding interaction with HCN (2.97D) stabilizes the complex significantly. In fact, the well depth of the HCN–CO complex is nearly three times greater than that of CS–CO, for which the SACM has already shown good accuracy.¹⁴ However, HCN has a lower density of rotational states compared to CS, and is therefore less efficiently stabilized in the potential well during collisions. On the other hand, this characteristic reduces the computational demand and enables the exploration of a broader energy range.

The cross sections were calculated for transitions between rotational levels up to $j_2 = 9$ for HCN and $j_1 = 8$ for CO, which are the most populated states up to 50 K. At this temperature, the adiabatic channels associated with energy levels up to ~500 cm⁻¹ would typically be included in the calculations. In this energy range, the accessible rotational levels of HCN and CO are $j_2^{\text{max}} = 18$ and $j_1^{\text{max}} = 15$, respectively. However, the inclusion of the full rotational basis led to a substantial memory demand, reaching over 15000 open channels even at only $J_{\text{tot}} = 2$. In order to ease the calculations from a computational point of view, the rotational basis was further truncated to $j_2^{\text{max}} = 12$ and $j_1^{\text{max}} = 10$. The accuracy of this truncation
was ensured by confirming that the convergence of the rate coefficients at $J_{tot} = 2$ is within 20% with respect to the full rotational basis over the 5–50 K temperature range. Since the main contributions from higher values of $J_{tot}$ are expected to increase the centrifugal barrier for each collisional process, the resulting deviations can reasonably be considered within the uncertainty range of the fully converged calculations.

Using the reduced rotational basis, a partial wave expansion from $J_{tot} = 0$ up to $J_{tot} = 142$ was performed, ensuring convergence of the cross sections better than 5%. Starting from the cross-sections for each $j_1j_2 \rightarrow j'_1j'_2$ transition ($\sigma_{j_1j_2 \rightarrow j'_1j'_2}$), the (de)-excitation rate coefficients ($k_{j_1j_2 \rightarrow j'_1j'_2}$) from 5 to 50 K were derived by thermal averaging over the collisional energy ($E_c$):

$$k_{j_1j_2 \rightarrow j'_1j'_2}(T) = \left(\frac{8}{\pi\mu k_B^3 T^3}\right)^{1/2} \times \int_0^\infty \sigma_{j_1j_2 \rightarrow j'_1j'_2}(E_c) E_c \exp(-E_c/k_B T) dE_c, \quad (1)$$

where $k_B$ is the Boltzmann constant.

In addition to the dataset of state-to-state collisional coefficients, thermalized rate coefficients for the (de)-excitation of HCN were also retrieved. These coefficients were obtained by summing the $k_{j_1j_2 \rightarrow j'_1j'_2}(T)$ coefficients over all possible final states of CO and averaging over the initial rotational states of CO:

$$\bar{k}_{j_2 \rightarrow j'_2}(T) = \sum_{j_1} n_{j_1}(T) \sum_{j'_1} k_{j_1j_2 \rightarrow j'_1j'_2}(T), \quad (2)$$

where $n_{j_1}(T)$ is the relative population of the projectile:

$$n_{j_1}(T) = \frac{(2j_1 + 1) \exp\left(-\frac{E_{j_1}}{k_B T}\right)}{\sum_{j'_1} (2j'_1 + 1) \exp\left(-\frac{E_{j'_1}}{k_B T}\right)}, \quad (3)$$

with $(2j_1 + 1)$ being the degeneracy of the rotational levels of CO.

Thermalized rate coefficients are commonly used in radiative transfer models⁴⁹ to address collisional excitation of the target molecule, under the assumption that the rotational temperature of the projectile is equal to the kinetic temperature of the surrounding gas. Therefore, this dataset can be applied to model the abundance of HCN in cometary comae thermalized with respect to CO.

To assess the accuracy of the SACM method, we compared a subset of rate coefficients with the results obtained with full quantum CC calculations. Due to the deep potential well of the system, a direct comparison of the inelastic cross sections could indeed be biased by the presence of resonances extending over a wide range of kinetic energies. To minimize this effect, instead of the inelastic cross sections, we directly compared the state-to-state rate coefficients, computed in an energy interval between 3 and 200 cm⁻¹. The 3–20 cm⁻¹ range was sampled with very fine steps (0.2 cm⁻¹), which gradually increased up to 2 cm⁻¹ above 100 cm⁻¹. To optimize the rotational basis for both HCN and CO, we evaluated how the inclusion of different rotational levels affected the convergence of the inelastic cross sections and we selected the basis sets that ensured convergence within 2%. Specifically, from 2 and 50 cm⁻¹, we included the first 22 rotational states of HCN and the first 19 states of CO. These were progressively increased up to 28 and 26, respectively, at 200 cm⁻¹. The angular part of the nuclear Schrödinger equation was restricted to a single partial wave $J_{tot} = 0$ to keep the computational cost affordable in the considered energy range. The radial part was subsequently included by numerical propagation over the $R = 2.9–42.0$ Å range. We computed CC rate coefficients in the 5–50 K temperature range and accounting for (de)-excitations among the rotational levels of HCN-CO below 80 cm⁻¹. These were compared with their counterparts computed within the SACM method, using the same rotational basis used for data production but accounting only for the adiabatic channels up to 200 cm⁻¹ and for $J_{tot} = 0$. The comparison for $T = 10, 30$ and 50 K is reported in Fig. 4.

At these temperatures, most of the rate coefficients are reproduced within a factor of 2, and deviations do not exceed one order of magnitude, confirming the method’s reliability. This was somewhat expected, as the agreement is consistent with the findings of Godard Palluet et al.¹⁴ for the CS and CO system.

The agreement between the CC and SACM approaches significantly improves when comparing thermalized rate coefficients, showing deviations of less than a factor of 2, as illustrated in Fig. 5. Thermal averaging over the rotational states of CO effectively smooths out the discrepancies between the two computational methods, making the SACM approach particularly accurate in capturing the overall collisional behavior of the system.
<figure><img src="image_5.png" /><figcaption>Figure 4: Comparison between the state-to-state rate coefficients computed with the CC and SACM approaches at $J_{tot} = 0$ for 10, 30 and 50 K.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>Figure 5: Comparison between the thermalized rate coefficients for the (de)-excitation among the rotational levels of the HCN and CO system computed with the CC and SACM approaches at $J_{\text{tot}} = 0$ for 10 K, 30 K and 50 K.</figcaption></figure>

## Results and Discussion

We derived the first dataset of (de)-excitation rate coefficients between the lowest ten rotational levels of HCN in collision with the first nine rotational levels of CO. While the full dataset will be made available through the EMAA†, LAMDA⁵⁰ and BASECOL⁵¹ databases, we summarize here the main outcomes and trends observed.

Fig. 6 shows the thermalized rate coefficients of de-excitation from the first ten rotational levels of HCN. For each initial $j_2$ level (indicated by a different color), the rate coefficients are plotted as a function of the $j'_2$ final state. The weak dependence of these trends with respect to the kinetic temperature stands out: although absolute values of the rate coefficients are slightly higher at low temperature, the relative tendency remains largely consistent across the different thermal contributions. The same behavior (although not shown here) is observed for the state-to-state rate coefficients resolved for the rotational levels of the CO projectile. In general, the thermalized rate coefficients tend to decrease with increasing $\Delta j_2$, reflecting the growing energy gap that they address. However, exceptions to this trend emerge when the $j'_2$ final state is equal to 1 or 4. In these cases, regardless of the initial state, a slight increase of the rate coefficient value is observed, suggesting a favor for these particular

transitions. At present, the authors have not identified a clear physical rationale for this behavior.

We now address the influence of the rotational state of the CO projectile, which is illustrated in Fig. 7. Here, the values of the state-to-state rate coefficients for de-excitation from the $j_2 = 9$ rotational state of HCN and for the transitions involving the lowest nine $j_1$ rotational states of CO are reported. The values are mapped with respect to $\Delta j_2$ (x-axis) and $-\Delta j_1$ (y-axis). Hence, the diagonal elements correspond to transitions where a de-excitation of HCN is matched by an excitation of CO of the same magnitude. Likewise, the elements of the first superdiagonal and subdiagonal represent transitions where the energetic jump between the levels of the two molecules differs by $\pm 1$, and so on. The largest values of the state-to-state rate coefficients are predominantly concentrated in this region, particularly for transition where $\Delta j_2 = -\Delta j_1$, $\Delta j_2 = -\Delta j_1 \pm 1$, and $\Delta j_2 = -\Delta j_1 \pm 2$. This behavior is a signature of near resonant energy transfer regimes,⁵²⁻⁵⁴ where the system tends to conserve its total internal energy during the collision. This phenomenon is observed when the colliding partners, such as HCN and CO, have similar rotational constants, allowing their energy levels to align closely in resonance. Noteworthy, the near resonant energy transfer is efficiently captured by the statistical SACM method, highlighting its robustness in describing even purely quantum effects.

We expect these results to be relevant to model HCN in cometary comae, where non-LTE radiative transfer is key for determining molecular abundances. As demonstrated in the previous section, the SACM approach provides sufficiently accurate rate coefficients, proving that the deep potential well and dense rotational structure of the system allow for a statistically valid redistribution of its energy after the formation of the intermediate complex.

## Conclusions

The present study aims to support the interpretation of current and future observations of HCN in cometary comae. In these environments, HCN is largely abundant and widely spread, making it an important tracer of the thermal and chemical history of precometary material. Furthermore, HCN is considered a fundamental building block in prebiotic chemistry, and its presence in comets offers

<footnote>†https://emaa.osug.fr//</footnote>
<figure><img src="image_7.png" /><figcaption>Figure 6: Variation of the thermalized de-excitation rate coefficients as a function of $j_2'$, starting from $j_2 = 0-9$ and for T = 10 K, 30 K and 50 K.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>Figure 7: Variation of some state-to-state rate coefficients as a function of Δj₂ (x axis) and −Δj₁ (y axis), starting from j₂ = 9 and j₁ = 0 − 8 and for T= 30 K.</figcaption></figure>
valuable insights into the processes that may have
contributed to the origin of life in the primordial
Earth. To contribute to this understanding, we
have computed the first dataset of inelastic rate
coefficients involving the ten lowest rotational lev-
els of HCN in collision with the lowest nine states
of CO, over a temperature range of 5–50 K.

We characterized the interaction potential using
coupled-cluster electronic structure theory and an
accurate fitting method. Subsequently, we per-
formed scattering calculations by means of the
SACM statistical approach. In order to facili-
tate direct application in non-LTE radiative trans-
fer models, the computed state-to-state rate coeffi-
cients were thermalized over the rotational states of
CO, assuming equilibrium between the rotational
population of CO and the kinetic temperature.

The accuracy of the SACM method was assessed
by comparison with a subset of full quantum CC
calculations at J<sub>tot</sub> = 0. Given the complexity
of the system and the lack of existing data, the
SACM showed an encouraging agreement, repro-
ducing state-to-state rate coefficients within a fac-
tor of 10 and thermalized rate coefficients within a
factor of 2.

With the contribution of this work, accurate col-
lisional datasets for HCN interacting with both
H₂O²³,²⁴ and CO are now available. Given the
ubiquity of HCN in cometary comae, these datasets
offer a valuable test case to investigate the sen-
sitivity of non-LTE effects under diverse physical
conditions, such as different heliocentric distances.
Worth to note, also CO₂ is quite abundant in many
cometary comae, especially at intermediate helio-
centric distances. However, its high density of ro-
tational states (resulting from a very low rotational
constant) makes its collisional characterization sig-
nificantly more demanding from a computational
point of view. Although it is not expected to be the
dominant collisional partner, computing rate coef-
ficients with CO₂ could nonetheless provide addi-
tional insights. All of these data will play a central
role in driving the non-LTE distribution of popu-
lation of HCN in cometary comae, bringing new
knowledge to probe the chemical and physical con-
ditions of the early Solar System.

**Acknowledgement** This work has been sup-
ported by Region Bretagne. The authors acknowl-
edge the "Programme National de Planétologie"
(PNP) under the responsibility of INSU, CNRS
(France). FL acknowledges the Institut Universi-

taire de France. RD, EQS, and ABP are supported
by the U.S. Department of Energy under Award
DE-SC0025420.

References

(1) Mumma, M. J.; Charnley, S. B. The chemical composition of comets—emerging taxonomies and natal heritage. *Annual Review of Astronomy and Astrophysics* **2011**, *49*, 471–524.

(2) Bockelée-Morvan, D. Cometary science with ALMA. *Astrophysics and Space Science* **2008**, *313*, 183–189.

(3) Roth, N. X.; Milam, S. N.; Cordiner, M. A.; Bockelée-Morvan, D.; Biver, N.; Boissier, J.; Lis, D. C.; Remijan, A. J.; Charnley, S. B. Leveraging the ALMA Atacama Compact Array for Cometary Science: An Interferometric Survey of Comet C/2015 ER61 (PanSTARRS) and Evidence for a Distributed Source of Carbon Monosulfide. *The Astrophysical Journal* **2021**, *921*, 14.

(4) Biver, N.; Boissier, J.; Bockelée-Morvan, D.; Crovisier, J.; Cottin, H.; Cordiner, M.; Roth, N.; Moreno, R. Observations of comet C/2020 F3 (NEOWISE) with IRAM telescopes. *Astronomy & Astrophysics* **2022**, *668*, A171.

(5) Cordiner, M.; Roth, N.; Milam, S.; Vil-
lanueva, G.; Bockelée-Morvan, D.; Remi-
jan, A.; Charnley, S.; Biver, N.; Lis, D.; Qi, C.
et al. Gas sources from the coma and nu-
cleus of comet 46P/Wirtanen observed us-
ing ALMA. *The Astrophysical Journal* **2023**,
*953*, 59.

(6) Loreau, J.; Faure, A.; Lique, F. The effect of water and electron collisions in the rotational excitation of HF in comets. *Monthly Notices of the Royal Astronomical Society* **2022**, *516*, 5964–5971.

(7) Quack, M.; Troe, J. Complex formation in reactive and inelastic scattering: Statistical adiabatic channel model of unimolecular processes III. *Berichte der Bunsengesellschaft für physikalische Chemie* **1975**, *79*, 170–183.
(8) Loreau, J.; Lique, F.; Faure, A. An efficient statistical method to compute molecular collisional rate coefficients. *The Astrophysical Journal Letters* **2018**, *853*, L5.

(9) Loreau, J.; Faure, A.; Lique, F. Scattering of CO with H₂O: Statistical and classical alternatives to close-coupling calculations. *The Journal of Chemical Physics* **2018**, *148*, 244308.

(10) Balança, C.; Scribano, Y.; Loreau, J.; Lique, F.; Feautrier, N. Inelastic rate coefficients for collisions of N₂H⁺ with H₂. *Monthly Notices of the Royal Astronomical Society* **2020**, *495*, 2524–2530.

(11) Pirlot Jankowiak, P.; Lique, F. Collisional Energy Transfer in the Highly Reactive OH⁺–H₂ System. *Physical Review Letters* **2025**, *134*, 253002.

(12) Joy, C.; Mandal, B.; Bostan, D.; Dubernet, M.-L.; Babikov, D. Mixed quantum/classical theory (MQCT) approach to the dynamics of molecule-molecule collisions in complex systems. *Faraday Discussions* **2024**, *251*, 225–248.

(13) Mandal, B.; Bostan, D.; Joy, C.; Babikov, D. MQCT 2024: A program for calculations of inelastic scattering of two molecules (new version announcement). *Computer Physics Communications* **2024**, *294*, 108938.

(14) Godard Palluet, A.; Dawes, R.; Quintas-Sánchez, E.; Batista-Planas, A.; Lique, F. A promising statistical approach for studying the collisional excitation induced by CO: Application to the CS–CO system. *The Journal of Chemical Physics* **2025**, *162*, 241101.

(15) Bockelée-Morvan, D. An overview of comet composition. *Proceedings of the International Astronomical Union* **2011**, *7*, 261–274.

(16) Drahus, M.; Küppers, M.; Jarchow, C.; Paganini, L.; Hartogh, P.; Villanueva, G. The HCN molecule as a tracer of the nucleus rotation of comet 73P-C/Schwassmann-Wachmann 3. *Astronomy & Astrophysics* **2010**, *510*, A55.

(17) Villanueva, G.; Magee-Sauer, K.; Mumma, M. Modeling of nitrogen compounds in cometary atmospheres: Fluorescence models of ammonia (NH₃), hydrogen cyanide (HCN), hydrogen isocyanide (HNC) and cyanoacetylene (HC₃N). *Journal of Quantitative Spectroscopy and Radiative Transfer* **2013**, *129*, 158–168.

(18) Matthews, C. N.; Minard, R. D. Hydrogen cyanide polymers connect cosmochemistry and biochemistry. *Proceedings of the International Astronomical Union* **2008**, *4*, 453–458.

(19) Jung, S. H.; Choe, J. C. Mechanisms of prebiotic adenine synthesis from HCN by oligomerization in the gas phase. *Astrobiology* **2013**, *13*, 465–475.

(20) Mumma, M. J.; Weissman, P. R.; Stern, S. Comets and the origin of the solar system-Reading the Rosetta Stone. Protostars and planets III. 1993; pp 1177–1252.

(21) Magee-Sauer, K.; Mumma, M. J.; DiSanti, M. A.; Dello Russo, N. Hydrogen cyanide in comet C/1996 B2 Hyakutake. *Journal of Geophysical Research: Planets* **2002**, *107*, 6–1.

(22) Biver, N.; Russo, N. D.; Opitom, C.; Rubin, M. Chemistry of comet atmospheres. Comets III. University of Arizona Press. 2022.

(23) Dubernet, M.-L.; Quintas-Sánchez, E. First quantum study of the rotational excitation of HCN by para-H₂O: Convergence of quantum results, influence of the potential energy surface, and approximate rate coefficients of interest for cometary atmospheres. *Molecular Astrophysics* **2019**, *16*, 100046.

(24) Zoltowski, M.; Lique, F.; Zuchowski, P.; Klos, J.; Loreau, J.; Kedziera, D. Collisional excitation of HCN by H₂O. *Monthly Notices of the Royal Astronomical Society* **2025**, *540*, 626–632.

(25) Ootsubo, T.; Kawakita, H.; Hamada, S.; Kobayashi, H.; Yamaguchi, M.; Usui, F.; Nakagawa, T.; Ueno, M.; Ishiguro, M.; Sekiguchi, T. et al. AKARI near-infrared spectroscopic survey for CO₂ in 18 comets. *The Astrophysical Journal* **2012**, *752*, 15.
(26) Yang, B.; Jewitt, D.; Zhao, Y.; Jiang, X.; Ye, Q.; Chen, Y.-T. Discovery of carbon monoxide in distant comet C/2017 K2 (PANSTARRS). *The Astrophysical Journal Letters* **2021**, *914*, L17.

(27) Ejeta, C.; Gibb, E.; DiSanti, M. A.; Kawakita, H.; Bonev, B. P.; Russo, N. D.; Roth, N. X.; Khan, Y.; McKay, A. J.; Combi, M. R. et al. Infrared Compositional Measurements in Comet C/2017 K2 (PAN-STARRS) at Heliocentric Distances Beyond 2.3 au. *The Astronomical Journal* **2025**, *169*, 102.

(28) Castro-Juárez, E.; Wang, X.-G.; Carrington Jr, T.; Quintas-Sánchez, E.; Dawes, R. Computational study of the ro-vibrational spectrum of CO–CO₂. *The Journal of chemical physics* **2019**, *151*, 084307.

(29) Quintas-Sánchez, E.; Dawes, R.; Wang, X.-G.; Carrington, T. Computational study of the rovibrational spectrum of CO₂–N₂. *Physical Chemistry Chemical Physics* **2020**, *22*, 22674–22683.

(30) Gancewski, M.; Jóźwiak, H.; Quintas-Sánchez, E.; Dawes, R.; Thibault, F.; Wcisło, P. Fully quantum calculations of O₂–N₂ scattering using a new potential energy surface: Collisional perturbations of the oxygen 118 GHz fine structure line. *The Journal of Chemical Physics* **2021**, *155*, 124307.

(31) Desrousseaux, B.; Quintas-Sánchez, E.; Dawes, R.; Marinakis, S.; Lique, F. Collisional excitation of interstellar PN by H₂: New interaction potential and scattering calculations. *The Journal of Chemical Physics* **2021**, *154*, 034304.

(32) Quintas-Sánchez, E.; Dawes, R.; Denis-Alpizar, O. Theoretical study of the HCS⁺–H₂ van der Waals complex: potential energy surface, rovibrational bound states, and rotationally inelastic collisional cross sections. *Molecular Physics* **2021**, *119*, 1980234.

(33) Zadrożny, A.; Jóźwiak, H.; Quintas-Sánchez, E.; Dawes, R.; Wcisło, P. Ab initio quantum scattering calculations for the CO–O₂ system and a new CO–O₂ potential energy surface: O₂ and air broadening of the

R (0) line in CO. *The Journal of Chemical Physics* **2022**, *157*.

(34) Ajili, Y.; Quintas-Sánchez, E.; Mehnen, B.; Żuchowski, P. S.; Brzek, F.; El-Kork, N.; Gacesa, M.; Dawes, R.; Hochlaf, M. Theoretical study of the CO₂–O₂ van der Waals complex: potential energy surface and applications. *Physical Chemistry Chemical Physics* **2022**, *24*, 28984–28993.

(35) Olejnik, A.; Jóźwiak, H.; Gancewski, M.; Quintas-Sánchez, E.; Dawes, R.; Wcisło, P. Ab initio quantum scattering calculations and a new potential energy surface for the HCl(X¹Σ⁺)–O₂(X³Σ⁻ₗ) system: Collision-induced line shape parameters for O₂-perturbed R(0) 0–0 line in H³⁵Cl. *The Journal of Chemical Physics* **2023**, *159*, 134301.

(36) Bostan, D.; Mandal, B.; Joy, C.; Żółtowski, M.; Lique, F.; Loreau, J.; Quintas-Sánchez, E.; Batista-Planas, A.; Dawes, R.; Babikov, D. Mixed quantum/classical calculations of rotationally inelastic scattering in the CO + CO system: a comparison with fully quantum results. *Physical Chemistry Chemical Physics* **2024**, *26*, 6627–6637.

(37) Quintas-Sánchez, E.; Dawes, R. AUTOSURF: A Freely Available Program To Construct Potential Energy Surfaces. *J. Chem. Inf. Model.* **2019**, *59*, 262–271.

(38) Majumder, M.; Ndengué, S. A.; Dawes, R. Automated construction of potential energy surfaces. *Mol. Phys.* **2016**, *114*, 1–18.

(39) Dawes, R.; Quintas-Sánchez, E. *Reviews in Computational Chemistry* vol. 31; John Wiley & Sons, 2018; Chapter 5, pp 199–264.

(40) Quintas-Sánchez, E.; Dawes, R. Spectroscopy and scattering studies using interpolated ab initio potentials. *Annual Review of Physical Chemistry* **2021**, *72*, 399–421.

(41) Werner, H.-J.; Knizia, G.; Manby, F. R. Explicitly correlated coupled cluster methods with pair-specific geminals. *Molecular Physics* **2011**, *109*, 407–417.

(42) Sobol, I. M. Uniformly distributed sequences with an additional uniform property. *USSR
*Computational Mathematics and Mathematical Physics* **1976**, *16*, 236–242.

(43) Peterson, K. A.; Adler, T. B.; Werner, H.-J. Systematically convergent basis sets for explicitly correlated wavefunctions: The atoms H, He, B–Ne, and Al–Ar. *J. Chem. Phys.* **2008**, *128*, 084102.

(44) Feller, D.; Peterson, K. A.; Crawford, T. D. Sources of error in electronic structure calculations on small chemical systems. *The Journal of chemical physics* **2006**, *124*, 054107.

(45) Werner, H.-J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M. Molpro: a general-purpose quantum chemistry program package. *Wiley Interdisciplinary Reviews: Computational Molecular Science* **2012**, *2*, 242–253.

(46) Winnewisser, G.; Belov, S. P.; Klaus, T.; Schieder, R. Sub-Doppler Measurements on the Rotational Transitions of Carbon Monoxide. *Journal of Molecular Spectroscopy* **1997**, *184*, 468–472.

(47) Ahrens, V.; Lewen, F.; Takano, S.; Winnewisser, G.; Urban, Š.; Negirev, A.; Koroliev, A. Sub-Doppler saturation spectroscopy of HCN up to 1 THz and detection of $J = 3 \rightarrow 2 (4 \rightarrow 3)$ emission from TMC1. *Zeitschrift für Naturforschung A* **2002**, *57*, 669–681.

(48) Alexander, M.; Dagdigian, P. J.; Werner, H.-J.; Kłos, J.; Desrousseaux, B.; Raffy, G.; Lique, F. Hibridon: A program suite for time-independent non-reactive quantum scattering calculations. *Computer Physics Communications* **2023**, *289*, 108761.

(49) Żółtowski, M.; Lique, F.; Loreau, J.; Faure, A.; Cordiner, M. The excitation of CO in CO-dominated cometary comae. *Monthly Notices of the Royal Astronomical Society* **2023**, *520*, 3887–3894.

(50) Schöier, F. L.; van der Tak, F. F.; van Dishoeck, E. F.; Black, J. H. An atomic and molecular database for analysis of submillimetre line observations. *Astronomy & Astrophysics* **2005**, *432*, 369–379.

(51) Dubernet, M.; Boursier, C.; Denis-Alpizar, O.; Ba, Y.; Moreau, N.; Zwölf, C.; Amor, M.; Babikov, D.; Balakrishnan, N.; Balança, C. et al. BASECOL2023 scientific content. *Astronomy & Astrophysics* **2024**, *683*, A40.

(52) Santos, S. F. d.; Balakrishnan, N.; Lepp, S.; Quéméner, G.; Forrey, R. C.; Hinde, R. J.; Stancil, P. C. Quantum dynamics of rovibrational transitions in H₂-H₂ collisions: Internal energy and rotational angular momentum conservation effects. *J. Chem. Phys.* **2011**, *134*, 214303–214303.

(53) Lique, F.; Alexander, M. H. Spin-orbit quenching of Cl($^2$P$_{1/2}$) by H$_2$. *J. Chem. Phys.* **2012**, *136*, 124312–124312.

(54) Lanza, M.; Kalugina, Y.; Wiesenfeld, L.; Lique, F. Near-resonant rotational energy transfer in HCl-H₂ inelastic collisions. *J. Chem. Phys.* **2014**, *140*, 064316.

TOC Graphic

<figure><img src="image_9.png" /></figure>