Issue 
A&A
Volume 683, March 2024



Article Number  A155  
Number of page(s)  12  
Section  Atomic, molecular, and nuclear data  
DOI  https://doi.org/10.1051/00046361/202348865  
Published online  15 March 2024 
Hyperfine excitation of NH and ND by molecular hydrogen
Rate coefficients and astrophysical modeling
^{1}
Univ. Rennes, CNRS, IPR (Institut de Physique de Rennes) – UMR 6251,
35000
Rennes,
France
email: paul.pirlot@univrennes.fr; francois.lique@univrennes.fr
^{2}
Instituto de Física Fundamental (CSIC),
Calle Serrano 121123,
28006
Madrid,
Spain
Received:
6
December
2023
Accepted:
8
January
2024
The NH and ND radicals are of key importance in the comprehension of nitrogen chemistry and the enhancement of deuterated molecules in the interstellar medium. Observations by space telescopes yield spectra that can resolve the fine and hyperfine structure of these radicals, a consequence of the electronic and magnetic interactions of nitrogen, hydrogen, and deuterium nuclei. Accurate rate coefficients, induced by collisions with H_{2}, are required to interpret spectra of these radicals. We report the first rate coefficients for fine and hyperfine transitions of NH and ND in collision with both ortho and paraH_{2}. Based on a recent fourdimensional potential energy surface, finestructure resolved cross sections and rate coefficients are computed with the timeindependent closecoupling method over a temperature range of 5–300 K. Our calculations include the first 25 energy levels of NH and ND. Hyperfine resolved cross sections and rate coefficients are determined using the infiniteorder sudden (IOS) approximation between 5 and 200 K for NH and 100 K for ND. We consider the first 71 and 105 energy levels of NH and ND, respectively. General propensity rules are discussed. We found a significant isotopic substitution effect in the rate coefficients. In addition, the rate coefficients for collisions with H_{2} are larger than those with He by a factor of up to 5, leading to lower critical densities for collisional excitation with H_{2} than He. The impact of the new set of collisional data has been investigated in simple radiative transfer models of the NH emission seen toward the Orion Bar and the ejecta of the η Carinae binary star. We observed significant differences by a factor of 5 between the presently determined column densities for NH compared to those from the literature using He as a collider.
Key words: astrochemistry / molecular data / molecular processes / radiative transfer / scattering
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1 Introduction
Nitrogen hydrides are key species in the formation of complex nitrogenbearing molecules observed in the interstellar medium (ISM; Gerin et al. 2016). An accurate determination of the abundances of these species is essential for understanding nitrogen chemistry, which remains a subject of ongoing exploration (Bacmann et al. 2016). Among nitrogen hydrides, the NH radical is of particular interest since it is an important intermediate in the synthesis of heavier molecules (Bacmann et al. 2010; Gerin et al. 2016). Observations of NH were first reported in comets (Swings et al. 1941) and stellar atmospheres (Schmitt 1969; Lambert & Beer 1972). Meyer & Roth (1991) detected NH in the ISM in ζ Per and HD 27778 diffuse clouds through the A^{3}Π − X^{3}Σ UV absorption band. The high rotational constant of NH leads to rotational lines in the submillimeter frequency range (≳ 1 THz), making groundbased observations not suitable for detecting rotational transitions. Thanks to the Infrared Space Observatory (ISO), Cernicharo et al. (2000) and Goicoechea et al. (2004) reported the first detection of farinfrared finestructure components 2_{3} → 1_{2} and 2_{2} → 1_{1} toward the Sgr B2 dense molecular cloud.
From the first detection of NH, large divergences between the determined abundances from observations and those predicted by chemical models were pointed out (Wagenblast et al. 1993). The much higher spectral resolution observations with the heterodyne receiver HIFI on board the Herschel space telescope led to the suggestion that the formation of NH in cold molecular clouds and prestellar cores may occur through dissociative recombination of N_{2}H^{+} (HilyBlant et al. 2010; Dislaire et al. 2012; Le Gal et al. 2014) in addition to the well known N^{+} + H_{2} reaction followed by dissociative recombination. Various models confirmed the poor understanding of the NH abundance and other nitrogen hydrides despite the many observations in the solarmass Class 0 protostar IRAS 162932422 (Bacmann et al. 2010; HilyBlant et al. 2010), in the H II region G10.60.4 (W31; Persson et al. 2012), the starforming region W49N (Persson et al. 2012), the binary star η Carinae, and the Homunculus nebula (Gull et al. 2020). Goicoechea & Roncero (2022) recently detected NH emission lines toward the Orion Bar photodissociation region (PDR), where their modeling based on statespecific rate coefficients calculations for the N + H_{2} (υ) → NH + H reaction lead to an enhancement of the NH column density by a factor ~25 compared to chemical models using thermal rate coefficients and showing better agreement with observations. The study revealed the importance of UVpumped vibrationally excited H_{2} in the formation of NH in strongly UVirradiated molecular clouds.
The deuterium fractionation is also a fundamental parameter of the study of the ISM. It is well known that the abundance ratio between a deuterated molecule and its hydrogenated counterpart can deviate by several orders of magnitude from the elemental ratio [D]/[H] ~ 2 × 10^{−5} (Linsky et al. 2006) depending on the molecule and on the source. High degrees of molecular deuteration have been observed in various cold environments (Ceccarelli et al. 1998; Caselli et al. 2003; Bacmann et al. 2003).
In this context, the study of the ND radical is of interest to infer the deuterium enhancement in the chemistry of nitrogen hydrides. Bacmann et al. (2010) reported the first detection of ND in IRAS 162932422 resolving the hyperfine multiplet of the 1 → 0 rotational line and deriving an abundance ratio [ND]/[NH] ≥ 30%. Bacmann et al. (2016) also detected ND toward the 16293E prestellar core and confirmed the high degree of deuteration of the NH radical.
At low densities, the determination of molecular abundances cannot be achieved assuming local thermodynamic equilibrium (LTE). The computation of reliable statetostate rate coefficients of NH and ND induced by collisions with the main collider H_{2} in molecular clouds is required. Interactions of the molecular electronic spin and nuclear spin of the nitrogen and the hydrogen or deuterium atoms to the molecular rotation leads to a large number of energy levels to consider in scattering calculations. In addition, the inclusion of the internal structure of H_{2} requires a large number of channels to consider in quantum dynamical calculations that are not achievable with current computational resources. There are presently no collisional data reported about fine and hyperfine resolved rate coefficients for NH and ND in collision with both ortho and paraH_{2}. Because of this lack of data, Bacmann et al. (2016) used scaled NH/ND– He hyperfine rate coefficients computed by Dumouchel et al. (2012) to model molecular spectra. In addition, Goicoechea & Roncero (2022) scaled finestructure NHHe rate coefficients from Tobola et al. (2011). These estimations of NH–H_{2} rate coefficients may not be accurate enough, especially in the case of light hydrides. Hyperfine lineoverlap and opacity effects may also lead to subtle radiative effects and anomalous hyperfine line emission (Goicoechea et al. 2022).
To date, only a few investigations have reported hyperfine resolved rate coefficients for collisions involving a molecule possessing two nonzero nuclear spins such as N_{2}H^{+}–He (Daniel et al. 2005), NH/ND–He (Dumouchel et al. 2012), or ^{13}CCH/C^{13}CH–paraH_{2} (Pirlot Jankowiak et al. 2023). Because NH is detected in both cold and warm astrophysical environments, it is essential to take into account explicitly the internal structure of H_{2} in the scattering calculations. This work aims to fill this objective by providing fine and hyperfine resolved rate coefficients for NH and ND in collisions with molecular hydrogen (both ortho and paraH_{2}). Finestructure calculations are achieved with the closecoupling (CC) method, while the infiniteorder sudden (IOS) approximation (Faure & Lique 2012; Lanza & Lique 2014) is employed to overcome the computational challenges posed by the inclusion of the hyperfine structure in the scattering calculations. The paper is organized as follows. Section 2 presents the potential energy surface (PES) used in this work. Details about the treatment for the ND isotopolog is also highlighted. Section 3 describes the scattering formalism for the obtention of both fine and hyperfinestructure resolved rate coefficients. Finestructure results are presented and discussed. Then hyperfine resolved rate coefficients are presented. Section 4 presents an illustration of the impact of the new computed rate coefficients in a simple radiative transfer modeling under nonLTE conditions. A summary of the work and our conclusions are presented in Sect. 5.
2 Potential energy surface
In this work we use the fourdimensional PES computed by Pirlot Jankowiak et al. (2021, hereafter Paper I). Ab initio calculations were done for 33 516 geometries with the explicitly correlated coupled cluster single, double, and a perturbative treatment of triple excitations [RCCSD(T)F12a] (Knizia et al. 2009) using the MOLPRO software (Werner et al. 2020). The augmented correlationconsistent polarized valence triple zeta (augccpVTZ) basis set (Dunning 1989) was used and the counterpoise scheme of Boys & Bernardi (1970) was considered to take into account the basis set superposition error.
This PES is provided as a function of the Jacobi coordinate system. The origin of this coordinate system lies at the center of mass of NH, with the intermolecular separation between NH and H_{2} centers of mass denoted as the vector R, aligned along the 𝓏axis. The rotational motion of NH around the R axis is described by the angle θ_{1}, while the angles (θ_{2}, ϕ) specify the orientation of H_{2} (see Fig. 1).
In this PES, both NH and H_{2} are treated as rigid rotors. We consider the internuclear distances taken in the average of the vibrational ground state for both species, with r_{NH} = 1.958a_{0} and ^{r}HH = 1.449a_{0} (Huber & Herzberg 1979).
To make this PES suitable for timeindependent quantum scattering calculations, the analytical representation of the potential was done in terms of bispherical harmonics, given by (1) (2)
Here P_{lm} denotes associated Legendre polynomials, and (…) represents 3jWigner symbols; l_{1}, l_{2}, and l are defined such that l_{1} − l_{2}1 < l < h + l_{2}. Then, the expansion coefficients are determined through a leastsquares fitting procedure, setting and and resulting in a total of 86 terms for subsequent scattering calculations.
The global minimum was found to be D_{e} = 149.10 cm^{−1} corresponding to a linear geometry characterized by R = 6.30a_{0}, θ_{1} = 180°, θ_{2} = 0°, and ϕ = 0°, where one hydrogen atom is oriented toward the nitrogen end. Additionally, a secondary minimum was observed for a Tshaped geometry with R = 6.77a_{0}, θ_{1} = 0°, θ_{2} = 90°, and ϕ = 0°, exhibiting a well depth D_{e} = 109.52 cm^{−1}. This PES has undergone validation, displaying good agreement between bound state calculations and experimental measurements from Fawzy et al. (2005).
The deuteration of NH yields the ND molecule by a displacement of the center of mass of δr = −0.1149a_{0} toward the deuterium end. To assess the interaction potential for the ND– H_{2} complex, we introduce new coordinates and transform the primary isotopolog coordinates as follows: (3) (4)
Within the rigid rotor approximation, we assume that the bond lengths for NH and ND remain identical. Notably, we found that the transformation of the (θ_{2}, ϕ) angles have no impact on the dynamics calculations and was not necessary. In order to obtain an analytical representation suitable for scattering calculations, a GaussLegendre quadrature procedure was performed for 588 geometries from the transformed NH–H_{2} PES with and .
Figure 2 displays the first radial coefficients for both NH–H_{2} and ND–H_{2} complexes. The most significant difference arises from the υ_{101} term. This term exhibits a lower magnitude for the ND–H_{2} PES than the NH–H_{2} PES, indicating a more pronounced odd anisotropy for the NH–H_{2} complex. This effect was already seen for the NH/ND–He collisional systems in the study of Dumouchel et al. (2012). We anticipate that the incorporation of these radial coefficients in scattering calculations will lead to differences in cross sections between the two isotopologs.
Fig. 1 Description of the NH–H_{2} and ND–H_{2} complexes in Jacobi coordinates. Unprimed and primed coordinates relate to NH and ND, respectively. 
Fig. 2 First expansion coefficients for the NH–H_{2} (solid) and ND–H_{2} (dashed) collisional systems. 
3 Scattering calculations
3.1 Finestructure excitation
The NH and ND radicals are both openshell molecules, and possess a nonzero electronic spin S = 1. This is coupled to the spinfree rotational angular momentum n_{1} such that
where j_{1} is the total rotational angular momentum of the target. Then, in the ^{3}Σ^{−} ground electronic state, the NH and ND rotational levels are split by spinrotation and spinspin interaction. In the intermediate coupling scheme, the wave function can be defined for j_{1} ≥ 1 as (Gordy & Cook 1984) (5)
where n_{1;} j_{1}m〉 denotes the wave function in a pure Hund’s case (b), and α is the angle arising from the diagonalization of the molecular Hamiltonian. In the pure Hund’s case (b) limit (α → 0), the F_{i} labels are associated with n_{1} = j_{1} − S for F_{1} and n_{1} = j_{1} + S for F_{3}. For simplicity, finestructure levels are labeled as corresponding to the pure Hund’s case (b) limit, although all calculations were done in the intermediate coupling scheme. In the CC approach, the NH or ND angular momentum couples with the rotational angular momentum j_{2} of the H_{2} collider to form the resultant sum j_{12}. This couples to the orbital angular momentum L to yield the total angular momentum J:
This is the general approach used to compute finestructure resolved cross sections with the CC procedure.
By the resolution of the CC equations, it is possible to determine finestructure statetostate cross sections defined as (6)
where is the wave vector of the initial state, [X] ≡ (2X + 1) and are the Tmatrix elements (i.e., the solutions of the CC equations). The rate coefficient from an initial state (n_{1} j_{1} j_{2}) to a final state can be obtained by integrating the CC cross section in Eq. (6) over a MaxwellBoltzmann distribution of the collisional energy E_{c}: (7)
with µ the reduced mass of the collisional system and k_{B} the Boltzmann constant.
Details of the scattering calculations can be found in Appendix A. The spectroscopic parameters used in finestructure calculations for NH and ND are summarized in Table 1. The scattering calculations were done using the HIBRIDON package of programs (Alexander et al. 2023).
Spectroscopic constants of NH and ND radicals in their ^{3}Σ^{−} electronic state.
3.2 Results
Figures 3 and 4 illustrate the temperature dependence of the finestructure resolved rate coefficients for NH and ND in collision with ortho and paraH_{2}. The magnitude of the rate coefficients can vary up to a factor of 3 between low and high temperatures. In addition, most of the rate coefficients for transitions induced by orthoH_{2} are larger than those induced by paraH_{2} by a factor of 3–10. This behavior arises due to the contributions of radial coefficients to the coupling between energy levels in scattering calculations. Specifically, orthoH_{2} collisions involve coefficients with l_{2} = 0, 2, whereas paraH_{2} collisions only consider coefficients with l_{2} = 0. The orthoH_{2} collisions are then characterized by enhanced anisotropy for j_{2} > 0 in the PES and tend to yield larger rate coefficients than for paraH_{2} collisions, a wellknown trend in neutral collisional systems involving H_{2} (Kalugina & Lique 2015; Desrousseaux et al. 2021; Demes et al. 2023).
As previously reported in Paper I, NH–orthoH_{2} collisions exhibit larger rate coefficients for odd ∆n_{1} transitions compared to NHparaH_{2} collisions. This is related to the magnitude of the υ_{121} radial coefficient, which has a greater intensity compared to those with l_{2} = 0 and provides a larger contribution in the coupling terms of the CC equations. A similar behavior is also observed for collisions involving ND.
Furthermore, Figs. 3 and 4 display that finestructure conserving transitions (∆n_{1} = ∆j_{1}) dominate over finestructure changing transitions (∆n_{1} ≠ ∆j_{1}). Alexander & Dagdigian (1983) already investigated this propensity rule. This rule implies that the electronic spin acts as a spectator during collisions, conserving its projection, due to the absence of electronic spin dependence in the potential energy surface. The rule can be seen clearly in Fig. 5.
It is interesting to have a look at the impact of the isotopic substitution on the rate coefficients. Figure 6 compares rate coefficients for NH and ND in collisions with ortho and paraH_{2}. The ND transitions exhibit rate coefficients that are nearly three times larger than NH transitions for orthoH_{2} collisions, a disparity that jumps to a factor of 5 for paraH_{2} collisions. This discrepancy arises from the displacement of the center of mass between NH–H_{2} and ND–H_{2} PESs and differences in NH and ND spectroscopic parameters (see Table 1). Similar isotopic differences have been observed in various collisional systems, particularly with light molecules such as NH/ND–He (Dumouchel et al. 2012) and OH/OD–H_{2} (Dagdigian 2021). These differences become less pronounced for heavier molecules, as seen in HCO^{+}/DCO^{+}–H_{2} (DenisAlpizar et al. 2020).
When collisional data involving H_{2} are not available, it is customary to use He as a proxy to estimate rate coefficients for the paraH_{2}(j_{2} = 0) collider, based on the assumption that their cross sections are similar. Both He and paraH_{2}(j_{2} = 0) possess two valence electrons and have a spherical shape. For NH or ND collisions, a massscaling relation can be used to model paraH_{2} with He (8)
where X is either NH or ND.
Figure 7 displays the discrepancies between finestructure rate coefficients computed in this work with those available from Dumouchel et al. (2012).^{1} Specifically, we compare NH and ND in collision with paraH_{2} and scaledHe. In the case of NH, substantial differences are observed. The paraH_{2} rate coefficients are larger than those of He by more than a factor of 5, particularly at low temperatures (T ≤ 50 K). Some transitions can exhibit larger deviations even at higher temperatures, showing a strong temperaturedependence between the datasets. A similar behavior is observed in the case of ND. This could be explained by the contrasts in the interaction potentials. The PES calculated for NH–He by Cybulski et al. (2005) presents a well depth smaller by a factor ~7 compared to the NH–H_{2} PES computed in Paper I. This difference leads to larger cross sections for paraH_{2} collisions. As a conclusion, He is an unsuitable template for paraH_{2}(j_{2} = 0) collisions, necessitating the explicit computation of these rate coefficients.
Fig. 3 Temperature dependence of finestructure resolved deexcitation rate coefficients for NH in collision with paraH_{2} (solid) and orthoH_{2} (dashed). The finestructure conserving transitions (∆n_{1} = ∆j_{1}) are presented in the top panel; the finestructure changing transitions (∆n_{1} ≠ ∆j_{1}) are shown in the bottom panel. 
Fig. 4 Temperature dependence of finestructure resolved deexcitation rate coefficients for ND in collision with paraH_{2} (solid) and orthoH_{2} (dashed). The finestructure conserving transitions (∆n_{1} = ∆j_{1}) are presented in the top panel; the finestructure changing transitions (∆n_{1} ≠ ∆j_{1}) are shown in the bottom panel. 
Fig. 5 Finestructure resolved rate coefficients at 50 K for NH (top) and ND (bottom) deexcitations out of the rotational level. 
Fig. 6 Systematic comparison between NH and ND rate coefficients at 10 and 150 K. The top panel shows orthoH_{2} collisions, where the green dashed lines quantify deviations by a factor of 3. The bottom panel is related to paraH_{2} collisions and the green dashed lines present differences of a factor of 5. 
Fig. 7 Comparison between the present finestructure resolved rate coefficients for NH and ND in collision with paraH_{2} and those in collision with He done by Dumouchel et al. (2012). The dashed lines represent deviations of the data of a factor of 5. 
3.3 Hyperfinestructure excitation
Additional energy splittings occur for NH and ND due to the presence of the nuclear spins of the hydrogen, deuterium, and nitrogen atoms, I_{H} = 1/2, I_{D} = 1, and I_{N} = 1, respectively. These spins couple to the rotational quantum number, as
with F_{1} and F being the rotational quantum numbers including nuclear spins of the targeted isotopolog. Such splittings are generally small (10^{−3}–10^{−4} cm^{−1}) due to the magnitude of the possible electric quadrupole and magnetic moments. In the case of NH, each rotational level (n_{1} > 2) possesses 18 hyperfine components except for n_{1} = 0,1,2, which have respectively 5, 13, and 17 hyperfine energy levels. For ND these splittings lead to 27 levels for n_{1} > 2 and 7, 19, and 26 levels for n_{1} = 0,1,2.
Collisional excitation between two molecules is a computational challenge, particularly in calculations involving openshell molecules with two nonzero nuclear spins. The computational cost of scattering calculations becomes prohibitive with our current computational resources, due to the extensive number of channels that must be considered. For NH and ND, our objective is to provide hyperfine rate coefficients up to 200 K, taking into account energy levels up to 326 cm^{−1} for NH and up to 100 K with energy levels up to 176 cm^{−1} for ND. Achieving this goal entails accounting for 71 energy levels in the case of NH and 105 for ND. To account both the nuclear spins of the target and the internal structure of the collider, the IOS limit is tested and used.
3.3.1 Infiniteorder sudden limit approach
It is possible to estimate statetostate hyperfine resolved rate coefficients only by considering the finestructure rate coefficient out of the fundamental energy level. Based on several works (Alexander 1982; Orlikowski 1985; Lanza & Lique 2014), the hyperfine resolved rate coefficient for a ^{3}σ electronic state molecule in collision with ortho and paraH_{2} can be written for an initial state i = (n_{1} j_{1} F_{1} F j_{2}) to a final one as (see details in Appendix B) (9) (10)
where [abc…] ≡ (2a + 1)(2b + 1)(2c + 1)…, λ_{1}, λ_{2} are chosen so that and or D. Moreover, IOSp stands for collisions with paraH_{2}, whereas IOSo stands for orthoH_{2} collisions. Equations (9) and (10) are used considering NH and ND in a pure Hund’s case (b) limit. It should be noted that the rate coefficient can be replaced by the cross section for IOS calculations. Equations (9) and (10) can also be used by considering deexcitations rate coefficients through the detailed balance: (11) (12)
It is worth noting that the gap between hyperfine energy levels remains relatively small compared to the collisional energies involved. Consequently, it is reasonable to anticipate that the IOS limit can provide reliable predictions for hyperfine resolved transitions. Then, the scaling relation can be used, first introduced by Neufeld & Green (1994, hereafter NG) as (13)
This relation is constructed so that the summation of the rate coefficients over the final hyperfine labels for a given transition gives the CC finestructure rate coefficient: (14)
One can see that using Eq. (13) implies that the calculation of quasielastic transitions (i.e., and with and F ≠ F’) requires the determination of finestructure elastic transitions. However, these transitions are usually difficult to converge. In the following we provide quasielastic transitions without the NG correction and only using the IOS limit through Eqs. (9) and (10). In addition, the presence of 3j and 6j Wigner coefficients involves strict propensity rules. In particular, the 3j coefficient sets all hyperfine rate coefficients to zero for the case where . The corresponding rate coefficients are only nonvanishing for hyperfine elastic transitions, which are beyond the scope of this work.
3.3.2 Validation of the IOS approach
The IOS approximation described previously can be tested and compared with the recoupling approach. This reference method, for hyperfine treatment can be found in detail in Alexander & Dagdigian (1985), among others. This comparison is applied to the NH–H_{2} and ND–H_{2} collisional systems, but taking into account only the nuclear spin of the hydrogen (or the deuterium). The IOS approach was tested for several cross sections at two different collisional energies.
Figure 8 displays a comparison of cross sections computed using the NG correction and the recoupling approach for NH and ND in collision with both ortho and paraH_{2}. Quasielastic transitions, calculated with the IOS approximation, are also included in these plots. Of all the tested collisional systems, the figure shows that most of the largest transitions match almost perfectly those obtained with the recoupling approach. It is possible to quantify the discrepancies between methods by introducing the weighted mean error factor (WMEF; Loreau et al. 2018), (15)
with being the hyperfine cross section obtained with the recoupling approach, . Larger discrepancies are observed for low collisional energies, while most transitions exhibit better agreement at higher energies. For paraH_{2} collisions, the WMEFs remain of the same order of magnitude for NH and ND, regardless of the energy. The inclusion of quasielastic transitions has minimal impact in this case. Among all the tested collisional systems, orthoH_{2} collisions display the most noticeable deviations, mostly due to the inclusion of quasielastic transitions, which are among the most divergent transitions in terms of intensity. However, it is interesting to note that, with the exception of these transitions, the NG correction performs surprisingly well in the case of orthoH_{2} collisions.
Fig. 8 Comparison between the NG correction and the recoupling approach for NH (top) and ND (bottom) cross sections for collisions with ortho(left) and paraH_{2} (right). Cross sections are displayed for collisional energies of 50 and 200 cm^{−1}. Quasielastic transitions computed in the IOS limit are also present. The dashed lines represent deviations of the data of a factor of 2. 
3.3.3 NH–H_{2} and ND–H_{2} rate coefficients
Figure 9 displays the propensity rules discussed in Sect. 3.1 and applied here to the hyperfine rate coefficients computed through the IOS approach using Eqs. (9), (10), and (13). Among all collisional systems, the strong propensity rule governing Δn_{1} = Δj_{1} = ΔF_{1} = ΔF transitions is found. This rule can be explained by vector coupling arguments. Given the independence of nuclear spins from the potential, they remain spectators during the collision, and they cannot be reoriented. As in Sect. 3.1, orthoH_{2} collisions are governed by the odd Δn_{1} propensity rule, except for NH for ΔF_{1} ≠ ΔF transitions, where this rule stands for paraH_{2} collisions. These trends are well known for various types of systems having one or two nuclear spins (Daniel et al. 2004; Dumouchel et al. 2012; Kalugina & Lique 2015; Lique et al. 2016; Klos et al. 2020; LaraMoreno et al. 2021; Ndaw et al. 2021; Godard Palluet & Lique 2024).
4 Rotational excitation of NH in the ISM
Detections of several NH finestructure lines in the ejecta of the massive binary star η Carinae have been reported by Gull et al. (2020). Additionally, Goicoechea & Roncero (2022) detected the submillimeter 1_{2} → 0_{1} finestructure line in the Orion Bar PDR. To our knowledge, these are the only two sources where NH rotational lines are observed in emission, while all other previous detections report absorption lines. In both studies, modeling based on accurate NH–H_{2} collisional data was not possible because of the lack of these data. In the absence of rate coefficients for NH–H_{2} collisions, astrophysical modeling relied on He collisional data using the scaling relation given by Eq. (8) or based on crude estimation when even these data are missing.
In the following discussions, we focus on the three finestructure lines 1_{0} → 0_{1} (946 GHz), 1_{2} → 0_{1} (974 GHz) and 1_{1} → 0_{1} (1000 GHz) detected by the Herschel/SPIRE survey in the study of Gull et al. (2020) and later toward the Orion Bar PDR with the Herschel/HIFI spectrometer by Goicoechea & Roncero (2022). We carried out simple radiative transfer modeling to assess the impact of the new NH–H_{2} rate coefficients in comparison to scaled NH–He data. We explored the impact of the new NH–H_{2} rate coefficients in modeling for typical conditions of these two environments with the RADEX nonLTE radiative transfer model using the escape probability approximation (Van der Tak et al. 2007).
Fig. 9 Hyperfine rate coefficients of NH (top) out of the (n_{1} = 4, j_{1} = 5, F_{1} = 5.5, F = 6.5) level and ND (bottom) out of the (n_{1} = 4, j_{1} = 5, F_{1} = 5, F = 6) level at 50 K as a function of the variation in the rotational angular momentum n_{1}. The solid lines stand for paraH_{2} collisions and the dashed lines for orthoH_{2} collisions. The ΔF_{1} = ΔF − 2 and Δn_{1} = Δ j_{1} − 2 transitions are represented by triangles, whereas the ΔF_{1} = ΔF − 1 and Δn_{1} = Δ j_{1} − 1 transitions are represented by squares. 
Fig. 10 Evolution of the critical density of NH with the temperature for the finestructure components of the 1 → 0 line. The solid lines stand for NH–H_{2} rate coefficients, whereas the dashed lines stand for NH–He coefficients. 
4.1 General features: Critical densities and background continuum
A useful parameter for qualitatively understanding the excitation behavior of molecular lines is the critical density from an upper level u to a lower one l, describing the gas density where radiative and collisional processes are in competition. Using the new finestructure NH–H_{2} rate coefficients, we estimated these critical densities to be (16)
where OPR is the orthotoparaH_{2} ratio assuming a thermalized distribution of the rotational population of H_{2}. We note, however, that the OPR can be out of equilibrium in various astrophysical environments, and can reach ~10^{−3} at low temperatures (Faure et al. 2019). The A_{ul} are the finestructure computed Einstein coefficients for NH radiative transitions (see Appendix C). These critical densities are compared to those computed with NH–He rate coefficients, also using Eq. (16), excluding the orthoH_{2} term and the OPR.
Figure 10 displays the critical density n_{cr} for each finestructure component of the 1 → 0 rotational line. Due to the high rotational constant of NH, Einstein coefficients corresponding to these lines are large, of the order of A_{ul} ~ 10^{−3}–10^{−2} s^{−1}. This is why the computed critical densities are high compared to the density of molecular clouds (n(H_{2}) ~ 10^{3}–10^{6} cm^{−3}). At very low temperature, a very high gas density is needed to reach the LTE regime. In this case, most of the excitation is due to radiative processes (i.e., when inelastic collisions do not dominate) and NH can be easily detected in absorption toward strong submillimeter continuum sources. It can be seen from Eq. (16) that the evolution of the critical density strongly depends on the OPR, reflecting the importance of collisions induced by orthoH_{2} at high temperatures because of their different excitation efficiency. Discrepancies between H_{2} and He critical densities are significant, exceeding an order of magnitude especially for the 1_{0} → 0_{1} line. These variations persist, even at low temperatures when paraH_{2} collisions are dominant, reflecting the previously discussed differences between paraH_{2} and He rate coefficients shown in Fig. 7.
It is interesting to assess the impact of rate coefficients on nonLTE modeling using He or H_{2} sets of collisional data, particularly in environments influenced by a strong farIR and submillimeter continuum. As an illustrative example, in these first models we adopt the physical parameters of Gull et al. (2020) for η Carinae; that is, we include a fixed column density N(NH) = 5 × 10^{15} cm^{−2}. We adopt a line width of 188 km/s in the ejecta (see Gull et al. 2020). These broad line widths ensure NH optically thin line emission. The kinetic temperature was chosen to be T_{k} = 150 K.
Figure 11 presents the variations in the excitation temperatures derived from models using NH–H_{2} and NH–He rate coefficients. We considered two distinct continua: the cosmic microwave background (CMB) described by a single blackbody at T_{bg} = 2.73 K and the modified source continuum that Gull et al. (2020) adopted. We found that excitation temperatures for individual finestructure lines did not exhibit significant relative differences. Then we computed an average excitation temperature over these lines. Radiative processes dominate for very low gas densities, with 〈T_{ex}〉 ≃ T_{bg} (i.e., radiative thermalization). Conversely, 〈T_{ex}〉 thermalizes to T_{k} when collisional processes are dominant for high densities. Given the typical conditions in η Carinae, a nonLTE behavior is obvious for both continua and sets of collisional data. The influence of the background continuum on excitation temperature is especially pronounced for low densities. When we adopt the specific source continuum, the averaged excitation temperature computed with NH–H_{2} rate coefficients depart from a radiative regime at n(H_{2}) ~ 10^{7} cm^{−3}, whereas when the continuum is only the CMB, the NH excitation is in the nonLTE regime for a gas density of ~10^{4}–10^{5} cm^{−3}. However, the impact of the continuum on the excitation temperature using H_{2} rate coefficients is negligible for n(H_{2}) ≥ 10^{8} cm^{−3} and excitation temperatures become identical between the two continua. On the other hand, for n(H_{2}) ~ 4 × 10^{8}−10^{10} cm^{−3}, models computed using the CMB and NH–He rate coefficients predict a suprathermal excitation (〈T_{ex}〉 > T_{k}·) which becomes moderate when we adopt the source continuum. This difference in behavior with the H_{2} collisional data is related to the lower magnitude of He rate coefficients. Consequently, excitation temperatures computed from He rate coefficients are more influenced by the strength of the source continuum than those computed using H_{2} rate coefficients.
Since the same excitation trends apply for both NH and ND, these results imply that the ND column densities inferred by Bacmann et al. (2016) in the prestellar core 16293E may be revised. Since ND is observed there in emission, a decrease in its calculated abundance can be anticipated. In addition, due to the larger values of the NHH_{2} rate coefficients compared to the NH–He values used in their work, models using the new collisional data may not be as sensitive to the modification of the background source continuum as they suggested.
Fig. 11 Dependence of the excitation temperature of the 1 → 0 rotational line on the gas density. The solid lines stand for NH–H_{2} rate coefficients, whereas the dashed lines stand for NH–He coefficients. 
Fig. 12 Models of the NH 1 → 0 rotational line toward η Carinae. The solid curve represents the observations of Gull et al. (2020), while the dashed curves refer to modeling using finestructure NH–H_{2} rate coefficients. 
4.2 Applications to the binary star η Carinae and the Orion Bar PDR
An essential parameter is the column density of the molecules, which allows us to infer the fractional abundance of a given molecule. For the η Carinae ejecta, Gull et al. (2020) determined a column density of N(NH) = 5 × 10^{15} cm^{−2} for a gas density of n(H_{2}) = 10^{8} cm^{−3} and a kinetic temperature T_{k} = 200 K. Using the new set of finestructure rate coefficients, we investigate the differences between the column density determined in their work and the value estimated in this study.
Figure 12 displays a simple model of the 1 → 0 rotational line compared to the model determined by Gull et al. (2020). Line fluxes computed using the previously determined column density and the new set of finestructure rate coefficients show a significant deviation from the observations. These differences are likely to be produced by an overestimation of the inelastic rate coefficients by these authors. By fixing the kinetic temperature to T_{k} = 200 K and a gas density to n(H_{2}) = 10^{8} cm^{−3}, it is possible to set the column density by minimizing the χ^{2} parameter. Based on our new collisional data, the minimization of this parameter must give an estimation of the integrated lines that best reproduce the observations. We found a column density of N(NH) = 2.43 × 10^{16} cm^{−2}, increasing the previous estimation by a factor of ~5. It looks like the two most intense lines are actually able to reproduce much better the observations than the 1_{0} → 0_{1} line. This may be due to the larger intensity uncertainty of this weak line or due to anomalous hyperfine emission produced by hyperfinestructure line overlap effects not treated by RADEX.
A 3σ detection of the 1_{2} → 0_{1} (974 GHz) line in emission toward the Orion Bar was recently reported by Goicoechea & Roncero (2022). These authors estimated the NH excitation conditions and NH column density using a Monte Carlo nonLTE radiative transfer model. For these calculations, the authors used the available finestructure NH–He collisional data from Toboła et al. (2011) to estimate the NH–H_{2} rate coefficients through Eq. (8). These rate coefficients were also extrapolated to high temperatures. They estimated an excitation temperature of 10–15 K for the 974 GHz line and a column density of 1.3 × 10^{13} cm^{−2}.
Figure 13 shows an illustration of the brightness temperatures for hyperfine transitions corresponding to the 1_{2} → 0_{1} line. These models were carried out using the new sets of hyperfine resolved rate coefficients, assuming a kinetic temperature of T_{k} = 150 K, a gas density n(H_{2}) = 10^{5} cm^{−3}, and a typical column density of N(NH) = 10^{13} cm^{−2}. The CMB was used as the background continuum and the hyperfine Einstein coefficients were extracted from the CDMS database (Endres et al. 2016). For the different sets of rate coefficients, we found quite important disparities in brightness temperatures. In particular, all the hyperfinestructure line peaks are brighter by a factor of 5 when using H_{2} collisional data compared to those obtained with He rate coefficients. It was found that matching brightness temperatures for all components between the two sets of collisional data implied a decrease in the column density by a factor of ~5 when using NH–H_{2} rate coefficients, leading to N(NH) ~2 × 10^{12} cm^{−2}. This discrepancy in column densities is due to the high contribution of orthoH_{2} rate coefficients. At T_{k} = 150 K, the thermalized OPR is ~2.51. Then each rate coefficient contributes to the radiative transfer. As mentioned in the previous sections, the NH–orthoH_{2} rate coefficients are larger than those for NH–paraH_{2} by a factor of more than 3, while the NH–paraH_{2} rates, being larger than NH–scaledHe rates, are larger by up to a factor of 5.
Fig. 13 Model of the NH hyperfine components of the 1_{2} → 0_{1} line. A line width value of ∆v = 0.5 km s^{−1} is used to avoid overlap of the hyperfine lines displayed in the spectrum. The column density is set to N(NH) = 10^{13} cm^{−2}. 
5 Conclusion
We computed the first fine and hyperfine resolved rate coefficients for NH and ND in collision with both ortho and paraH_{2}. These data are based on the 4D PES reported in Paper I. Finestructure rate coefficients were determined with the CC method up to 300 K, taking into account NH and ND energy levels up to n_{1} = 8. Both NH and ND displayed significantly larger finestructure transitions when interacting with orthoH_{2} than with paraH_{2}. This trend is standard for light neutral collisional systems. Moreover, these transitions display the usual properties of ^{3}Σ systems where the magnitude of the transitions varies with the transfered angular momentum as (∆n_{1} = ∆j_{1}) > (∆n_{1} ≠ ∆j_{1})
ND collisions with H_{2} yield rate coefficients that differ by a factor of 3–5 with respect to those of NH. These large differences could have been anticipated for such light collisional systems due to the large differences in spectroscopic parameters and the large displacement of the center of mass. Consequently, these findings highlight the importance of computing rate coefficients explicitly for deuterated isotopologs. Additionally, a comparison between H_{2} and scaledHe colliders has shown deviations of a factor >5, emphasizing the inadvisability of using of He as a proxy for rate coefficients involving H_{2}.
The treatment of the hyperfine structure was done using the IOS approach (Faure & Lique 2012) with the NG correction (Neufeld & Green 1994), adapted for the first time to collisions between an openshell target in a ^{3}Σ electronic state and a linear molecule. This allowed us to provide hyperfine collisional data up to 200 K for NH collisions and 100 K for ND collisions, including energy levels up to n_{1} = 4. The computed hyperfine rate coefficients follow a strong propensity for ∆n_{1} = ∆j_{1} = ∆F_{1} = ∆F transitions. A larger gap is observed between ∆n_{1} ≠ ∆ j_{1} = ∆F_{1} = ∆F and ∆n_{1} = ∆ j_{1} = ∆F_{1} ≠ ∆F transitions in the case of NH as opposed to ND collisions.
The new sets of NHH_{2} fine and hyperfine resolved rate coefficients were included in a radiative transfer modeling for typical conditions of the ejecta of the η Carinae binary star and the Orion Bar PDR where NH submillimeter lines have been detected in emission. We found important differences between the present column densities and those derived in the literature (up to a factor of ~5). These differences are mainly due to the use of the new NHH_{2} rate coefficients. These results suggest a downward revision of the calculation of the abundance of NH (and possibly ND) in these environments using more robust radiative transfer models. These data are expected to be useful for interpreting observations of these species detected in both cold and warm environments where both ortho and paraH_{2} coexist.
The computed collisional data for the fine and hyperfine structure of NH and ND in collision with both ortho and paraH_{2} will be available on the following data bases: EMAA^{2}, BASECOL (Dubernet et al. 2024), and LAMDA (Van der Tak et al. 2020).
Acknowledgements
We acknowledge financial support from the European Research Council (Consolidator Grant COLLEXISM, Grant Agreement No. 811363). We wish to acknowledge the support from the CEA/GENCI (Grand Equipement National de Calcul Intensif) for awarding us access to the TGCC (Très Grand Centre de Calcul) Joliot Curie/IRENE supercomputer within the A0110413001 project. J.R.G. thanks the Spanish MICINN for funding support under grant PID2019106110GBI00.
Appendix A Scattering parameters used for dynamical calculations
The resolution of the CC equations requires optimizing several parameters in order to converge finestructure cross sections for different total energies E_{tot}. For NH and ND in collision with ortho and paraH_{2}, the total energy grid is chosen with an increasing step ∆E with the energy. The step is small at low energies to be able to describe resonances and becomes larger at higher energies to save computational time since the energy variation of the cross sections become smoother. The converged rotational basis n_{max} of the target considered and the total angular momentum J_{tot} are presented in the Tables A.1–A.2. These parameters were chosen in order to converge integral cross sections to 1% accuracy per parameter. The wavefunctions were propagated for all collisional systems from 3.75a_{0} to 60a_{0}. The inclusion of the rotational basis for paraH_{2} and for orthoH_{2} were found to be enough for converged scattering calculations. As reported in Paper I, the impact of higher basis for both colliders on cross sections is less than 5% around E_{tot} = 500 cm^{−1} and can be neglected. The reduced mass used for the NH–H_{2} and ND–H_{2} collisional systems are respectively 1.777 and 1.790 amu.
Converged parameters used in closecoupling calculations for the NH–H_{2} collisional system.
Converged parameters used in closecoupling calculations for the ND–H_{2} collisional system.
Appendix B Details about IOS scaling relations
We describe here a more detailed derivation of the IOS scaling relations. From Eq. (6), it is possible to use angular momentum algebra in order to write the collisional integral cross section between a ^{3}Σ target and a linear collider so that (B.1)
where is a quantity defined as (B.2)
and represents the IOS Tmatrix coefficients excluding the angular dependence, similar to the coefficients defined by Corey & McCourt (1983) applied for atommolecule collisions.
Within the IOS approximation, we can express the integral cross section for a transition out of the finestructure level. For collisions induced by H_{2}, we detail two cases. These are treated in a pure Hund’s case (b) even if the intermediate coupling is used for the calculations of the CC cross sections.
B.1 paraH_{2}
Setting n_{1} = 1, j_{1} = 0 and j_{2} = 0 in (B.1) gives (B.3)
It is possible to see that by setting , the 3j and 6j coefficients are nonzero for and so that (B.4)
It should be noted that an equivalent set of equations like (B.4) can be found for . In the following and for the determination of the hyperfine rate coefficients, we chose to use the case where . The evaluation of the 3j and 6j coefficients in (B.4) leads to the expression for : (B.5)
This relation stands for S = 1. By reintroducing (B.5) into (B.1), we obtain the IOS cross section scaling relation out of the in the case of a collision with paraH_{2}: (B.6)
where the IOS cross section is replaced by the CC cross section. It is also possible to replace the cross section by the rate coefficient. The extension to hyperfine scattering calculation gives directly Eq. (9).
B.2 orthoH_{2}
Similarly, by applying n_{1} = 1, j_{1} = 0 and j_{2} = 1 in (B.1), we find (B.7)
In this case, for , the 3j and 6j coefficients give a contribution for and . Replacing in (B.7), we find (B.8)
The 3j and 6j coefficients in (B.8) are analyzed so that the expression for becomes (B.9)
We can reintroduce (B.9) into (B.1) and obtain the IOS cross section scaling relation out of the in the case of a collision with orthoH_{2}: (B.10)
Similarly to the case of paraH_{2} the extension to hyperfine scattering calculation is given by Eq. (10).
Appendix C Details about finestructure Einstein coefficients
For an openshell molecule following the coupling scheme given in Sect. 3.1, the general form of the Einstein coefficient is given by (see, e.g., Sobelman 1979) (C.1)
where v is the frequency of the transition, h is the Planck constant, c is the speed of light, and , and are the line strength factors for a finestructure and a rotational transition, respectively. The term depends on the electric dipole moment matrix elements, which allow transitions for ∆n_{1} = ± 1 (Axner et al. 2004). It turns out that (C.2)
Reintroduce (C.2) into (C.1) simplifies to (C.3)
In the case of NH and ND, S = 1 andµ_{D} = 1.39 D (Keun Park & Sun 1993).
References
 Alexander, M. H. 1982, J. Chem. Phys., 76, 429 [NASA ADS] [CrossRef] [Google Scholar]
 Alexander, M. H., & Dagdigian, P. J. 1983, J. Chem. Phys., 79 [Google Scholar]
 Alexander, M. H., & Dagdigian, P. J. 1985, J. Chem. Phys., 83, 2191 [NASA ADS] [CrossRef] [Google Scholar]
 Alexander, M. H., Dagdigian, P. J., Werner, H. J., et al. 2023, Comp. Phys. Commun., 289, 108761 [NASA ADS] [CrossRef] [Google Scholar]
 Axner, O., Gustafsson, J., Omenetto, N., & Winefordner, J. D. 2004, Spectrochim. Acta, 59, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Bacmann, A., Lefloch, B., Ceccarelli, C., et al. 2003, ApJ, 585, L55 [NASA ADS] [CrossRef] [Google Scholar]
 Bacmann, A., Caux, E., HilyBlant, P., Parise, B., & Pagani, L. 2010, A&A, 521, L42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bacmann, A., Daniel, F., Ceccarelli, C., et al. 2016, A&A, 587, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Boys, S. F., & Bernardi, S. 1970, Mol. Phys., 19, 553 [NASA ADS] [CrossRef] [Google Scholar]
 Caselli, P., Van der Tak, F. F. S., Ceccarelli, C., & Bacmann, A. 2003, A&A, 403, L37 [CrossRef] [EDP Sciences] [Google Scholar]
 Ceccarelli, C., Castets, A., Loinard, L., Caux, E., & Tielens, A. G. G. M. 1998, A&A, 338, L43 [NASA ADS] [Google Scholar]
 Cernicharo, J., Goicoechea, J. R., & Caux, E. 2000, ApJ, 534, L199 [NASA ADS] [CrossRef] [Google Scholar]
 Corey, G. C., & McCourt, F. R. 1983, J. Chem. Phys., 87, 2723 [Google Scholar]
 Cybulski, H., Krems, R. V., Sadeghpour, H. R., et al. 2005, J. Chem. Phys., 122 [Google Scholar]
 Dagdigian, P. J. 2021, MNRAS, 505, 1987 [CrossRef] [Google Scholar]
 Daniel, F., Dubernet, M. L., & Meuwly, M. 2004, J. Chem. Phys., 121, 4540 [NASA ADS] [CrossRef] [Google Scholar]
 Daniel, F., Dubernet, M. L., Meuwly, M., Cernicharo, J., & Pagani, L. 2005, MNRAS, 363, 1083 [NASA ADS] [CrossRef] [Google Scholar]
 Demes, S., Lique, F., Loreau, J., & Faure, A. 2023, MNRAS, 524, 2368 [NASA ADS] [CrossRef] [Google Scholar]
 DenisAlpizar, O., Stoecklin, T., Dutrey, A., & Guilloteau, S. 2020, MNRAS, 497, 4276 [Google Scholar]
 Desrousseaux, B., QuintasSánchez, E. L., Dawes, R., Marinakis, S., & Lique, F. 2021, J. Chem. Phys., 154, 034304 [NASA ADS] [CrossRef] [Google Scholar]
 Dislaire, V., HilyBlant, P., Faure, A., et al. 2012, A&A, 537, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dubernet, M. L., Boursier, C., DenisAlpizar, O., et al. 2024, A&A, 683, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dumouchel, F., Klos, J., Tobola, R., et al. 2012, J. Chem. Phys., 137, 114306 [NASA ADS] [CrossRef] [Google Scholar]
 Dunning, T. H. 1989, J. Chem. Phys., 90, 1007 [Google Scholar]
 Endres, C. P., Schlemmer, S., Schilke, P., Stutzki, J., & Müller, H. S. P. 2016, J. Mol. Spectrosc., 327, 95 [NASA ADS] [CrossRef] [Google Scholar]
 Faure, A., & Lique, F. 2012, MNRAS, 425, 740 [NASA ADS] [CrossRef] [Google Scholar]
 Faure, A., HilyBlant, P., Rist, C., et al. 2019, MNRAS, 487, 3392 [Google Scholar]
 Fawzy, W. M., Kerenskaya, G., & Heaven, M. C. 2005, J. Chem. Phys., 122, 144318 [NASA ADS] [CrossRef] [Google Scholar]
 FloresMijangos, J., Brown, J. M., Matsushima, F., et al. 2004, J. Mol. Spectrosc., 225, 189 [NASA ADS] [CrossRef] [Google Scholar]
 Gerin, M., Neufeld, D. A., & Goicoechea, J. R. 2016, ARA&A, 54, 181 [NASA ADS] [CrossRef] [Google Scholar]
 Godard Palluet, A., & Lique, F. 2024, MNRAS, 527, 6702 [Google Scholar]
 Goicoechea, J. R., & Roncero, O. 2022, A&A, 664, A190 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Goicoechea, J. R., RodriguezFernandez, N. J., & Cernicharo, J. 2004, ApJ, 600, 214 [CrossRef] [Google Scholar]
 Goicoechea, J. R., Lique, F., & SantaMaria, M. G. 2022, A&A, 658, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gordy, W., & Cook, R. L. 1984, Microwave Molecular Spectra (New Jersy: John Wiley & Sons, Inc) [Google Scholar]
 Gull, T. R., Morris, P. W., Black, J. H., et al. 2020, MNRAS, 499, 5269 [NASA ADS] [CrossRef] [Google Scholar]
 HilyBlant, P., Maret, S., Bacmann, A., et al. 2010, A&A, 521, L52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Huber, K. P., & Herzberg, G. 1979, Molecular Spectra And Molecular Structure, IV. Constants Of Diatomic Molecules (New York: Van Nostrand Reinhold) [CrossRef] [Google Scholar]
 Kalugina, Y., & Lique, F. 2015, MNRAS, 446, 21 [Google Scholar]
 Keun Park, J., & Sun, H. 1993, Chem. Phys. Lett., 211, 6 [Google Scholar]
 Klos, J., Dagdigian, P. J., Alexander, M. H., Faure, A., & Lique, F. 2020, MNRAS, 493, 3491 [CrossRef] [Google Scholar]
 Knizia, G., Adler, T. B., & Werner, H. J. 2009, J. Chem. Phys., 130, 054104 [NASA ADS] [CrossRef] [Google Scholar]
 Lambert, D. L., & Beer, R. 1972, ApJ, 177, 514 [Google Scholar]
 Lanza, M., & Lique, F. 2014, J. Chem. Phys., 141, 164321 [NASA ADS] [CrossRef] [Google Scholar]
 LaraMoreno, M., Stoecklin, T., & Halvick, P. 2021, MNRAS, 507, 4086 [NASA ADS] [CrossRef] [Google Scholar]
 Le Gal, R., HilyBlant, P., Faure, A., et al. 2014, A&A, 562, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Linsky, J. L., Draine, B. T., Moos, H. W., et al. 2006, ApJ, 647, 1106 [Google Scholar]
 Lique, F., Bulut, N., & Roncero, O. 2016, MNRAS, 461, 4477 [CrossRef] [Google Scholar]
 Loreau, J., Lique, F., & Faure, A. 2018, ApJ, 853, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Meyer, D. M., & Roth, K. C. 1991, ApJ, 376, L49 [NASA ADS] [CrossRef] [Google Scholar]
 Ndaw, D., Bop, C. T., Dieye, G., Boye Faye, N. A., & Lique, F. 2021, MNRAS, 503, 5976 [NASA ADS] [CrossRef] [Google Scholar]
 Neufeld, D. A., & Green, S. 1994, ApJ, 432, 158 [NASA ADS] [CrossRef] [Google Scholar]
 Orlikowski, T. 1985, Mol. Phys., 56, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Persson, C. M., De Luca, M., Mookerjea, B., Olofsson, A. O. H., & Black, J. H. 2012, A&A, 543, A145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pirlot Jankowiak, P., Kalugina, Y., Ramachandran, R., et al. 2021, J. Chem. Phys., 155, 134303 [NASA ADS] [CrossRef] [Google Scholar]
 Pirlot Jankowiak, P., Lique, F., & Dagdigian, P. J. 2023, MNRAS, 523, 3732 [NASA ADS] [CrossRef] [Google Scholar]
 Ramachandran, R., Klos, J., & Lique, F. 2018, J. Chem. Phys., 148, 084311 [NASA ADS] [CrossRef] [Google Scholar]
 Schmitt, J. L. 1969, PASP, 81, 657 [NASA ADS] [CrossRef] [Google Scholar]
 Sobelman, I. I. 1979, Atomic Spectra and Radiative Transitions (New York: SpringerVerlag Berlin Heidelberg New York), 1 [Google Scholar]
 Swings, P., Elvey, C. T., & Babcock, H. W. 1941, ApJ, 94, 320 [NASA ADS] [CrossRef] [Google Scholar]
 Takano, S., Klaus, T., & Winnewisser, G. 1998, J. Mol. Spectrosc., 192, 309 [NASA ADS] [CrossRef] [Google Scholar]
 Tobola, R., Dumouchel, F., Klos, J., & Lique, F. 2011, J. Chem. Phys., 134, 024305 [CrossRef] [Google Scholar]
 Van der Tak, F. F. S., Black, J. H., Schöier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, MNRAS, 468, 627 [Google Scholar]
 Van der Tak, F. F. S., Lique, F., Faure, A., Black, J. H., & van Dishoeck, E. F. 2020, Atoms, 8, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Wagenblast, R., Williams, D. A., Millar, T. J., & Nejad, L. A. M. 1993, MNRAS, 260, 420 [NASA ADS] [Google Scholar]
 Werner, H. J., Knowles, P. J., Manby, F. R., et al., F. R. M. 2020, J. Chem. Phys., 152, 144107 [NASA ADS] [CrossRef] [Google Scholar]
Ramachandran et al. (2018) published more recent data about the NH–He complex that are consistent with those of Dumouchel et al. (2012), but did not provide finestructure rate coefficients for the ND–He complex.
All Tables
Spectroscopic constants of NH and ND radicals in their ^{3}Σ^{−} electronic state.
Converged parameters used in closecoupling calculations for the NH–H_{2} collisional system.
Converged parameters used in closecoupling calculations for the ND–H_{2} collisional system.
All Figures
Fig. 1 Description of the NH–H_{2} and ND–H_{2} complexes in Jacobi coordinates. Unprimed and primed coordinates relate to NH and ND, respectively. 

In the text 
Fig. 2 First expansion coefficients for the NH–H_{2} (solid) and ND–H_{2} (dashed) collisional systems. 

In the text 
Fig. 3 Temperature dependence of finestructure resolved deexcitation rate coefficients for NH in collision with paraH_{2} (solid) and orthoH_{2} (dashed). The finestructure conserving transitions (∆n_{1} = ∆j_{1}) are presented in the top panel; the finestructure changing transitions (∆n_{1} ≠ ∆j_{1}) are shown in the bottom panel. 

In the text 
Fig. 4 Temperature dependence of finestructure resolved deexcitation rate coefficients for ND in collision with paraH_{2} (solid) and orthoH_{2} (dashed). The finestructure conserving transitions (∆n_{1} = ∆j_{1}) are presented in the top panel; the finestructure changing transitions (∆n_{1} ≠ ∆j_{1}) are shown in the bottom panel. 

In the text 
Fig. 5 Finestructure resolved rate coefficients at 50 K for NH (top) and ND (bottom) deexcitations out of the rotational level. 

In the text 
Fig. 6 Systematic comparison between NH and ND rate coefficients at 10 and 150 K. The top panel shows orthoH_{2} collisions, where the green dashed lines quantify deviations by a factor of 3. The bottom panel is related to paraH_{2} collisions and the green dashed lines present differences of a factor of 5. 

In the text 
Fig. 7 Comparison between the present finestructure resolved rate coefficients for NH and ND in collision with paraH_{2} and those in collision with He done by Dumouchel et al. (2012). The dashed lines represent deviations of the data of a factor of 5. 

In the text 
Fig. 8 Comparison between the NG correction and the recoupling approach for NH (top) and ND (bottom) cross sections for collisions with ortho(left) and paraH_{2} (right). Cross sections are displayed for collisional energies of 50 and 200 cm^{−1}. Quasielastic transitions computed in the IOS limit are also present. The dashed lines represent deviations of the data of a factor of 2. 

In the text 
Fig. 9 Hyperfine rate coefficients of NH (top) out of the (n_{1} = 4, j_{1} = 5, F_{1} = 5.5, F = 6.5) level and ND (bottom) out of the (n_{1} = 4, j_{1} = 5, F_{1} = 5, F = 6) level at 50 K as a function of the variation in the rotational angular momentum n_{1}. The solid lines stand for paraH_{2} collisions and the dashed lines for orthoH_{2} collisions. The ΔF_{1} = ΔF − 2 and Δn_{1} = Δ j_{1} − 2 transitions are represented by triangles, whereas the ΔF_{1} = ΔF − 1 and Δn_{1} = Δ j_{1} − 1 transitions are represented by squares. 

In the text 
Fig. 10 Evolution of the critical density of NH with the temperature for the finestructure components of the 1 → 0 line. The solid lines stand for NH–H_{2} rate coefficients, whereas the dashed lines stand for NH–He coefficients. 

In the text 
Fig. 11 Dependence of the excitation temperature of the 1 → 0 rotational line on the gas density. The solid lines stand for NH–H_{2} rate coefficients, whereas the dashed lines stand for NH–He coefficients. 

In the text 
Fig. 12 Models of the NH 1 → 0 rotational line toward η Carinae. The solid curve represents the observations of Gull et al. (2020), while the dashed curves refer to modeling using finestructure NH–H_{2} rate coefficients. 

In the text 
Fig. 13 Model of the NH hyperfine components of the 1_{2} → 0_{1} line. A line width value of ∆v = 0.5 km s^{−1} is used to avoid overlap of the hyperfine lines displayed in the spectrum. The column density is set to N(NH) = 10^{13} cm^{−2}. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.