Collisional excitation of NH(3{\Sigma}-) by Ar: A new ab initio 3D potential energy surface and scattering calculations

TL;DR

This study develops a new 3D potential energy surface for NH-Ar interactions, enabling detailed scattering calculations that improve understanding of molecular collisions relevant to atmospheric and astrophysical environments.

Contribution

We created a new ab initio 3D potential energy surface for NH-Ar, including NH bond vibration, and performed scattering calculations to determine collisional excitation rates.

Findings

Calculated cross sections for NH-Ar collisions up to 3000 cm-1.

Derived rate coefficients for temperatures up to 350 K.

Propensity rules align with similar systems, with some differences.

Abstract

Collisional excitation of light hydrides is important to fully understand the complex chemical and physical processes of atmospheric and astrophysical environments. Here, we focus on the NH(X3{\Sigma}-)-Ar van der Waals system. First, we have calculated a new three-dimensional Potential Energy Surface (PES), which explicitly includes the NH bond vibration. We have carried out the ab initio calculations of the PES employing the open-shell single- and double-excitation couple cluster method with noniterative perturbational treatment of the triple excitations. To achieve a better accuracy, we have first obtained the energies using the augmented correlation-consistent aug-cc-pVXZ (X = T, Q, 5) basis sets and then we have extrapolated the final values to the complete basis set limit. We have also studied the collisional excitation of NH(X3{\Sigma}-)-Ar at the close-coupling level, employing our new PES. We calculated collisional excitation cross sections of the fine-structure levels of NH by Ar for energies up to 3000 cm-1 . After thermal average of the cross sections, we have then obtained the rate coefficients for temperatures up to 350 K. The propensity rules between the fine-structure levels are in good agreement with those of similar collisional systems, even though they are not as strong and pronounced as for lighter systems, such as NH-He. The final theoretical values are also compared with the few available experimental data.

Tables

Table 1. Table 1: NH–Ar bound energy levels (in cm -1 ) obtained excluding the NH fine structure. Energies are relative to the ground-state energy of NH. All the levels correspond to the approximate quantum numbers N = 0 𝑁 0 N=0 . J 𝐽 J and l 𝑙 l correspond to the total and orbital angular momentum of the complex, respectively.

$J$	$l$	Energy (cm-1)
0	0	-73.1507
1	1	-72.9324
2	2	-72.4988
3	3	-71.8479

Table 2. Table 2: NH–Ar bound energy levels (in cm -1 ) obtained with the inclusion of the NH fine structure. Energies are relative to the ground-state energy of NH. All the levels correspond to the approximate quantum numbers N = 0 , F 1 𝑁 0 subscript 𝐹 1 N=0,F_{1} . J 𝐽 J and l 𝑙 l correspond to the total and orbital angular momentum of the complex, respectively.

$J$	$l$	Energy (cm-1)
1	0	-73.1519
0	1	-72.8964
1	1	-72.9507
2	1	-72.9305
1	2	-72.4804
2	2	-72.5169
3	2	-72.4947
2	3	-71.8333
3	3	-71.8658
4	3	-71.8426

Table 3. Table 3: Comparison between experimental and our theoretical cross sections at a collisional energy of 410 cm -1 and for transitions out of the N = 0 , F 1 𝑁 0 subscript 𝐹 1 N=0,F_{1} rotational level. All the values are normalized with respect to the cross section for the N = 0 , F 1 → N ′ = 1 , F 1 formulae-sequence formulae-sequence 𝑁 0 → subscript 𝐹 1 superscript 𝑁 ′ 1 subscript 𝐹 1 N=0,F_{1}\rightarrow N^{\prime}=1,F_{1} transition. Experimental error in parenthesis are in units of the last quoted digit.

	$F_1$		$F_2$		$F_3$
$N^{\prime }$	Exp111Ref. 14	Theory222Our work.	Exp111Ref. 14	Theory222Our work.	Exp111Ref. 14	Theory222Our work.
1	1.0	1.0	0.662(40)	0.560	0.255(54)	0.217
2	0.407(54)	1.906	0.284(46)	1.275	0.154(20)	0.833
3	0.068(15)	0.427	0.059(10)	0.321	0.047(08)	0.239
4	0.023(05)	0.061	-	-	-	-

Table 4. Table 4: Comparison between experimental and convolved theoretical cross sections at a collisional energy of 410 cm -1 and for transitions out of the N = 0 , F 1 𝑁 0 subscript 𝐹 1 N=0,F_{1} rotational level. All the values are normalized with respect to the cross section of the N = 0 , F 1 → N ′ = 1 , F 1 formulae-sequence formulae-sequence 𝑁 0 → subscript 𝐹 1 superscript 𝑁 ′ 1 subscript 𝐹 1 N=0,F_{1}\rightarrow N^{\prime}=1,F_{1} transition. Experimental error in parenthesis are in units of the last quoted digit. The convolution of the theoretical values is described in Sec. IV.

	$F_1$			$F_2$			$F_3$
$N^{\prime }$	Exp111Ref. 14	Theory222Our work. See text for details.	Exp111Ref. 14	Theory222Our work. See text for details.	Exp111Ref. 14	Theory222Our work. See text for details.
1	1.0	1.0	0.662(40)	0.842	0.255(54)	0.365
2	0.407(54)	0.954	0.284(46)	0.654	0.154(20)	0.411
3	0.068(15)	0.224	0.059(10)	0.172	0.047(08)	0.123
4	0.023(05)	0.037	-	-	-	-

Equations

$$V(R,\theta ,r)$$

$$E_{\mathrm{X}}=E_{\mathrm{CBS}}+A{e}^{-(X-1)}+B{e}^{-(X-1{)}^2},$$

$$V(R,\theta ,r)=\sum _{n=1}^N\sum _{l=1}^L{B}_{l,n}(R)(r-r_e{)}^{n-1}{d}_{m0}^{l+m-1}(\cos (\theta )),$$

$$\displaystyle B_{l,n}(R)=e^{-a_{l,n}(R-R_{l,n}^{(0)})}\Big{(}\sum_{i=0}^{2}b_{l,n}^{(i)}R^{i}\Big{)}$$

$$\displaystyle-\frac{1}{2}\Big{(}1+\tanh\frac{R-R_{l,n}^{(1)}}{R_{l,n}^{ref}}\Big{)}\sum_{j=6,8,10}\frac{c_{l,n}^{(j)}}{R^{j}}.$$

$$V_v(R,\theta )=⟨v(r)|V(R,\theta ,r)|v(r)⟩$$

$$|F_{1}jm⟩$$

$$|F_{2}jm⟩$$

$$|F_{3}jm⟩$$

$$\mathbf{j}=\mathbf{N}+\mathbf{S}$$

$$k_{F_{i}j\to F_i^{\prime }{j}^{\prime }}(T)$$

Full text

I Introduction

The study of inelastic collisions plays a relevant role in the understanding of important processes in different fields, such as atmospheric and astrophysical chemistry and physics. In particular open-shell molecules are crucial, being highly reactive compounds and intermediate in a large number of chemical reactions. A relevant chemical species is the NH radical. This compound serves as a prototype for other collisional studies involving open-shell molecules. Being diatomic it is also preferred for both experimental and theoretical scattering studies, owing to its large rotational energy level spacings. In addition, the magnetic moment of its ${}^{3}\Sigma^{-}$ electronic ground state, makes NH suitable for studies of ultracold molecules Friedrich and Doyle (2009); Egorov et al. (2004), because it can be easily thermalized at low temperatures through collision with cold buffer gas atoms. In the past, NH has been subject of many theoretical and experimental collisional studies in different electronic states and with a variety of perturbers, such as the rare gases He Alexander et al. (1991); Rinnenthal and Gericke (2002); Krems et al. (2003); Cybulski et al. (2005); Stoecklin (2009); Toboła et al. (2011); Dumouchel et al. (2012); Ramachandran et al. (2018) and Ne Rinnenthal and Gericke (2000); Kerenskaya et al. (2005); Bouhafs and Lique (2015).

In our work we focus on the calculation of a new ab initio 3D-averaged Potential Energy Surface (PES) and collisional excitation for the NH(${}^{3}\Sigma^{-}$)-Ar system. To our knowledge, there are no theoretical scattering studies for the fine-structure excitation of NH(${}^{3}\Sigma^{-}$) by Ar, while there is only one experimental work performed by Dagdigian 1989, employing a crossed beam apparatus. However, this experiment provides only relative collisional cross sections up to the rotational level $N$=4 and no rate coefficients are available.

The most recent PES is given by Kendall et al. 1998. They employed a combination of supermolecular and intermolecular unrestricted Møller-Plesset perturbation theory (UMPPT) Cybulski et al. (1995); Cybulski (1989) and a selection of monomer-centered basis sets augmented with bond functions. However, the NH bond length was kept frozen at 1.96 bohr. Recent studies Kalugina et al. (2014); Bouhafs and Lique (2015); Lique (2015); Ramachandran et al. (2018) have proven that the use of a 3D PES which takes into account molecular vibration leads to more accurate results when employed in collisional excitation studies of light hydrides by rare gases. Moreover, inclusion of the bond vibrational motion makes it possible to comprise excited vibrational states. Hence, we have computed a new ab initio PES for the NH(${}^{3}\Sigma^{-}$)–Ar van der Waals complex including the NH bond vibration.

Then, we present the first fully quantum close-coupling (CC) calculations of rotational inelastic cross sections for the NH(${}^{3}\Sigma^{-}$)–Ar collisional system. In addition, we have taken into account the spin-coupling splitting of the rotational levels and we have included the temperature dependence of the fine-structure resolved rate coefficients in the final results.

The paper is organized as follows: Sec. II covers the calculation of the new NH–Ar PES and informations about the bound states of the NH–Ar complex; in Sec. III, we present the scattering calculations, including the inelastic cross sections and rate coefficients. In Sec. IV, we compare the resulting cross sections with the available experimental data in Ref. 14. Conclusions are given in Sec. V.

II Potential energy surface

The two interacting species are considered in their ground electronic states NH(${}^{3}\Sigma^{-}$) and Ar(${}^{1}S$). The NH(${}^{3}\Sigma^{-}$)–Ar van der Waals system has 3A*′′* ground electronic state. In this work, we used the Jacobi coordinate system (see Fig. 1). The center of coordinates is placed in the NH center of mass (c.m.), and the vector R connects the NH c.m. with the Ar atom. The rotation of NH molecule is defined by the $\theta$ angle and the $r$ coordinate describes the NH bond length.

We performed the calculations for five NH bond lengths $r$ = [1.6, 1.8, 1.95, 2.15, 2.5] bohr which allows us to take into account vibrational motion of NH molecule up to $v=2$. We have carried out ab initio calculations of the PES of the NH–Ar van der Waals complex at the partially spin-restricted coupled cluster with single, double and perturbative triple excitations [RCCSD(T)] Hampel et al. (1992); Watts et al. (1993) level of theory, using MOLPRO 2015 package (mol, ). In order to determine the interaction potential, $V(R,\theta,r)$, the basis set superposition error (BSSE) was corrected at all geometries using the Boys and Bernardi counterpoise scheme Boys and Bernardi (1970):

[TABLE]

where the energies of the NH and Ar monomers are computed using the full basis set of the complex.

To achieve a good description of the charge-overlap effects we have performed the calculations in a rather large augmented correlation-consistent basis sets aug-cc-pVXZ (X = T, Q, 5) Dunning Jr (1989). Then, we have extrapolated the energies to the Complete Basis Set (CBS) limit using the following scheme Peterson et al. (1994):

[TABLE]

where $X$ is the cardinal number of the aug-cc-pVXZ basis set, $E_{\rm X}$ is the energy corresponding to aug-cc-pVXZ basis set, $E_{\rm CBS}$ is the energy extrapolated to CBS limit, $A$ and $B$ are the parameters to adjust. We have carried out the calculations for $\theta$ angle values from $0^{\circ}$ to $180^{\circ}$ in steps of $10^{\circ}$. $R$-distances were varied from 3.0 to 40.0  bohr, yielding 52 points for each angular orientation. Overall $\sim$5000 single point energies were calculated for the NH–Ar complex.

II.1 Analytical representation of the potential energy surface

The analytical expression employed for the interaction potential $V(R,\theta,r)$ has the following form Werner et al. (1988):

[TABLE]

where

[TABLE]

The basis functions $d_{m0}^{l+m-1}(\cos(\theta))$ are Wigner rotation functions, $N$ is the total number of $r$-distances, and $L$ is the total number of angles. The analytic potential was found to reproduce the calculated energies quite well: the mean difference between the analytic fit and the ab initio computed interaction energies is less than 2% over the entire grid.

Previous studies Kalugina et al. (2014) have shown that averaging of the PES over corresponding vibrational level $v$ leads to a better agreement with experimental results than using a purely two-dimension PES. The newly constructed PES, which takes into account the stretching of the NH molecule, can be averaged over any vibrational state, up to $v$ = 2. The averaging is done using the following formula:

[TABLE]

The NH vibrational wave functions $|v(r)\rangle$ were those computed in Bouhafs et al. Bouhafs and Lique (2015) that were evaluated using a discrete variable representation (DVR) method Colbert and Miller (1992) from ab initio calculations of the NH potential function using the internally contracted multireference configuration interaction (MRCI)Werner and Knowles (1988) level and a large aug-cc-pV5Z atomic basis set.

Figure 2 depicts the contour plot of our 3D PES averaged over the ground vibrational state $v=0$ as a function of $R$ and $\theta$ (hereafter refereed as 3D-ave PES). For this weakly-bound system the global minimum in the interaction energy was found to be -104.138 cm*-1* ($R=6.7$ bohr, $\theta=69^{\circ}$).

Our study is in good agreement with the NH–Ar PES previously published Kendall et al. (1998). Kendall et al. 1998 carried out calculations for the NH–Ar interaction with the supermolecular unrestricted Møller-Plesset (UMP) perturbation theory and a combination of different basis sets. The NH intermolecular distance was fixed at 1.95 bohr. According to the authors the best results have been obtained with the aug-cc-pVTZ(ext-b) basis set, augmented with bond functions, and the global minimum is found at $R=6.75$ bohr and $\theta=67^{\circ}$, with a well-depth of -100.3 cm*-1* and an uncertainty within the 5%. These values are very close to our results for $r=1.95$ bohr ($R=6.7$, $\theta=67^{\circ}$, 103.787 cm*-1*). Furthermore, the results of our 3D-ave PES also agree well with those listed above, confirming the high accuracy of our study. The slightly increased deepness of our well-depth is mostly due to the use of CBS extrapolation, since the energy follows a monotonic trend towards negative values, by approaching the infinite basis set limit. Figure 3 depicts the variation of the angle at which occurs the minimum of the interaction potential, for different NH bond distances. While the equilibrium angle changes substantially over increasing $r$, the $R$ distance is always close to $R=6.7$.

II.2 NH–Ar bound states and dissociation energy

Using the highly correlated 3D-ave PES described in the previous section, we have computed the bound states of NH–Ar complex using a coupled-channel approach, as implemented in the BOUND program bou . The bound state calculations were performed for the main 14N and 40Ar isotopes.

As a first step we performed bound state calculations neglecting the NH fine structure (i.e. NH was considered as a closed shell molecule). The calculations were performed with a propagator step size of 0.01 bohr, and the other propagation parameters were taken as the default BOUND values. The rotational basis includes the rotational states with $N_{max}\leq 10$. The bound energy levels of the NH–Ar complex computed with the 3D-ave PES are listed in Table 1. From the present calculations, dissociation energy ($D_{0}$) of the complex is 73.15 cm*-1* which is slightly larger than the previously calculated value of Kendall et al. Kendall et al. (1998) ($D_{0}=71.5$ cm*-1*). The difference (1.65 cm*-1*) can be mainly attributed to the difference between the two NH–Ar PESs used in the calculations. Indeed, the well depth of the 3D-ave PES considering the vibration motion is slightly deeper (by few cm*-1*) than the rigid rotor one of Kendall et al. Kendall et al. (1998) and this difference leads to a larger estimated value of the dissociation energy.

In order to derive the rotational constant of the NH–Ar complex, we have fitted the energies of Table 1 to the rigid rotor expression: $E_{J}=E_{0}+BJ(J+1)-DJ^{2}(J+1)^{2}$ where $J$ corresponds to the total angular momentum of the complex. We have obtained for the rotational and quartic centrifugal distortion constants, $B=0.1087$ cm*-1* and $D=0.000025$ cm*-1*. Such estimates allow generating the energetic structure of the complex and are useful for the interpretation of future experimental spectra. As a comparison, our rotational constant is in good agreement with the value obtained by Jansen et al. 1993, i.e. B=0.1007.

As previously mentioned, the NH molecule exhibits a fine structure because of the coupling between the rotational angular momentum and the electronic spin. The BOUND program was modified to include this fine structure of the NH molecule Ramachandran et al. (2018). Table 2 presents the bound state energies for the first total angular momentum $J$. The predicted bound energy levels indicate that the coupling of the electron spin to the rotational motion of the complex is very weak. As a consequence, energy levels of NH–Ar are very similar to those obtained by neglecting the fine structure, as already found for the NH–He complex Ramachandran et al. (2018). The dissociation energy is thus not significantly impacted by the fine structure.

III Scattering calculations

Rotational transitions in the NH(${}^{3}\Sigma^{-}$) electronic ground state show fine-structure splitting, due to spin-rotation coupling. The rotational wave function of NH for $j\geq 1$ in the intermediate coupling scheme can be written as Gordy and Cook (1984); Lique et al. (2005):

[TABLE]

where $|N,Sjm\rangle$ denotes pure Hund’s case (b) basis functions and the mixing angle $\alpha$ is obtained by diagonalisation of the molecular Hamiltonian. In this relation corresponding to the Hund’s case (b), the total molecular angular momentum $j$ is defined by:

[TABLE]

where $\mathbf{N}$ and $\mathbf{S}$ are the nuclear rotational and the electronic spin angular momenta. In the pure case (b) limit, $\alpha\to 0$, the $F_{1}$ level corresponds to $N=j-1$ and the $F_{3}$ level to $N=j+1$. The levels in the spin multiplets are usually labeled by the nuclear rotational quantum number $N$ and the spectroscopic index $F_{i}$. This notation will be used hereafter.

Using the new 3D-ave PES, we have studied the collisional excitation of NH by Ar. The scattering calculations were performed for the main 14N and 40Ar isotopes. The detailed description of the Close-Coupling (CC) calculations that consider the fine structure levels of the colliders is given in Ref.32. The quantal coupled equations have been solved in the intermediate coupling scheme using the MOLSCAT codeHutson and Green (1994) modified to take into account the fine structure of the rotational energy levels.

We used a total energy grid with a variable steps. For the energies below 500 cm*-1* the step was equal to 1 cm*-1*, then, between 500 and 1000 cm*-1* it was increased to 2 cm*-1*, and to 20 for the interval 1000-3000 cm*-1*. Using this energy grid, the resonances (shape and Feshbach) that usually appear in the cross sections at low energies were correctly represented.

In order to ensure convergence of the inelastic cross sections, it is necessary to include in the calculations several energetically inaccessible (closed) levels. At the largest energies considered in this work, the NH rotational basis was extended to $N=12$ to ensure convergence of the rotational cross sections between levels with $N<8$. One also needs to converge inelastic cross sections with respect to partial waves. The total angular momentum quantum number $J$ needed for the convergence was set up to 238 for the inelastic cross sections.

Moreover, in MOLSCAT, it is necessary to adjust the propagator’s parameters in order to ensure convergence of cross sections calculations. For all the energies, the minimum and maximum integration distances were $R_{min}=3.0~{}\text{bohr}$ and $R_{max}=50~{}\text{bohr}$, respectively. The STEPS parameter was adjusted for each value of energy in order to obtain a step length of the integrator sufficient to achieve the convergence. In our work, the value of the STEPS parameter decreases with increasing energy, going from 50 to 7 and, therefore, constraining the R spacing below 0.1-0.2 bohr at all energies. The reduced mass of the NH–Ar system is $\mu=10.912$ u and the NH(${}^{3}\Sigma^{-}$) rotational and centrifugal distortion constants have been taken from Ref. 34.

Figure 4 presents the energy variation of the integral cross sections for transitions from the initial rotational level $N=0,F_{1}$ of NH.

The resonances shown at low collisional energies are related to the presence of a $\sim$104 cm*-1* deep attractive potential well. As a consequence, the Ar atom can be temporarily trapped there forming quasi-bond states before dissociation of the complex Smith et al. (1979); Christoffel and Bowman (1983). However, excitation cross sections are less affected and therefore show few resonances. Indeed, the energy spacing between rotational levels is generally larger than the well depth of the PES.

The magnitude of the cross sections shown in figure 4 seems to present the following propensity rules:

(1) Overall decreasing of the cross sections with increasing $\Delta N$, according to the usual trend for rotational excitation. In addition, even $\Delta N$ transitions are favored over odd $\Delta N$ transitions as a consequence of near-homonuclearity of the PES.

(2) Fine-structure conserving transitions are favored, i.e. $\Delta j=\Delta N$ in the case of pure Hund’s case (b).

The same propensity rules are shown in similar systems, such as NH-He and NH-Ne collisions Toboła et al. (2011); Bouhafs and Lique (2015); Ramachandran et al. (2018). In particular, the latter applies in general to molecules in the ${}^{3}\Sigma^{-}$ electronic state. Both porpensity rules have been predicted theoretically (Alexander and Dagdigian, 1983) and also observed for the O2-He Orlikowski (1985); Lique (2010) or SO(X${}^{3}\Sigma^{-}$)–HeLique et al. (2005) collisions.

The thermal rate coefficients, $k_{F_{i}j\to F_{i}^{\prime}j^{\prime}}(T)$, for excitation and de-excitation transitions between fine-structure levels of NH can be calculated by averaging CC excitation cross sections, $\sigma_{F_{i}j\to F_{i}^{\prime}j}$, over a Maxwellian distribution of collision velocities, as follows:

[TABLE]

where $k_{B}$ is the Boltzmann constant, $\mu$ is the reduced mass of the system and $E_{k}$ is the translational energy.

The thermal dependence of these state-to-state rate coefficients for temperatures up to 350 K is shown in Fig. 5 for transitions out of the $N=0,j=1,F_{1}$ level.

The rate coefficients display the same propensity rules as seen for the integral cross sections. In particular, the rate coefficients for $F$-conserving transitions are generally larger than those for $F$-changing transitions.

IV Comparison with experiments

Our new calculated cross sections can be compared with the previous experimental results, obtained for a collisional energy of 410 cm*-1* and for rotational levels up to $N$=4,$F_{1}$ (Ref. 14).

Tab. 3 shows experimental and theoretical values normalized with respect to the $N=0,F_{1}\rightarrow N^{\prime}=1,F_{1}$ cross section. The $F$-conserving propensity rule is overall fullfilled in both the experimental and calculated values. In addition, the $F$-conserving cross sections follow the simple scaling relation observed by Dagdigian 1989, as shown in Fig. 6.

The main discrepancy is the trend of the cross sections over increasing $\Delta N$ and over even/odd $\Delta N$, as discussed for the first propensity rule in Sec. III. Furthermore, according to the results of Ref. 14, the largest cross sections are those with $N^{\prime}=1$, whereas this is not the case in our study. In fact, larger values are related to the transitions involving $N^{\prime}=2$, as also shown in figure 4.

It is likely that these discrepancies are due to a particular feature of the experiment. In fact, as declared by the author, the NH beam was not entirely pure, with 68% of the population in the rotational ground state $N=0,F_{1}$, and approximately 16% and 9% in the $N=1,F_{1}$ and $N=1,F_{2}$ levels, respectively. By taking into account this NH beam population composition, the propensity rules observed in the experiment can be reproduced making a convolution of the various cross sections involved. This is shown in Tab. 4, which gathers values computed using 68% contribution from inelastic cross section for transitions out of the $N=0,F_{1}$, 16% from cross sections involving the $N=1,F_{1}$, and 9% from those involving the $N=1,F_{2}$.

It should be pointed out that there is a 7% population with unknown distribution and thus the theoretical results obtained through convolution are still different in magnitude from the experimental ones.

V Conclusion

We have computed a new highly accurate 3D PES for the NH–Ar collisional system by taking into account the stretching of the NH bond. We carried out these ab initio calculations at the RCCSD(T) level and a complete basis set extrapolation. The results are in good agreement with the most recent PES available Kendall et al. (1998).

Employing our new 3D-ave PES we have calculated the dissociation energy of the NH–Ar van der Waals complex and the corresponding rotational and centrifugal distortion constants. We have also performed scattering calculations at the close-coupling level, obtaining collisional cross sections for energies up to 3000 cm*-1*. We have then determined rate coefficients for temperatures up to 350 K. The resulting values follows the same propensity rules seen in other similar systems (Bouhafs and Lique, 2015; Ramachandran et al., 2018), i.e. overall decreasing with increasing $\Delta N$, even $\Delta N$ favored over odd $\Delta N$ and larger values for $F$-conserving transitions.

Our theoretical results have been compared to a previous experimental studyDagdigian (1989). The discrepancy concerning the $\Delta N$ propensity rules can be explained with the impurity of the NH population of the experimental molecular beam, since we have been able to reproduce the results of the experiment through convolution of various cross sections, as discussed in Sec. IV.

We hope that our results will encourage new experimental studies concerning collisional excitation of NH(${}^{3}\Sigma^{-}$) by Ar. In particular it would be interesting to fill the gap of missing data regarding Ar as a collisional partner, with respect to systems involving He or Ne, more widely studied. Furthermore, a complete overview of these systems could also encourage studies with ortho- and para-H2, highly important for astrophysical environments.

Supplementary Material

The supplementary material provides the analytic form of the NH–Ar potential energy surface and the NH–Ar collisional rate coefficients.

Acknowledgements

F. L. acknowledges the Institut Universitaire de France. We would also thank P. Dagdigian for the interesting discussion about the comparison between our theoretical and his experimental results.