arXiv:1809.02083v1 [astro-ph.GA] 6 Sep 2018 [email protected] orsodn uhr sa Asvany Oskar author: Corresponding oetfeunytasto fNH of transition frequency lowest rniin r eetd epciey n v rniin r d are transitions five and respectively, detected, are transitions ai srnmcldtcino h w te isotopologues. other two the of detection astronomical radio n NHD and 1 2 nttt eEtutr el aei IMCI) ern 1 Serrano (IEM-CSIC), Materia la de Estructura de Instituto .Pyiaice nttt nvri¨tz K¨oln, Z¨ulpi Universit¨at zu Institut, Physikalisches I. erpr etfeunisfrrttoa rniin ftedeu the of transitions rotational for frequencies rest report We yee sn L using 2018 Typeset 7, September version Draft CUAERTTOA ETFEUNISFRAMNU O ISOT ION AMMONIUM FOR FREQUENCIES REST ROTATIONAL ACCURATE Keywords: 3 + esrdi roei o rpmcie o h ymti osN tops symmetric the For machine. trap ion cryogenic a in measured , A T E X S:mlcls—mtos aoaoy oeua oeua aa—t — data molecular — molecular laboratory: methods: — molecules ISM: troscopic twocolumn Jos .Dom L. e ´ tl nAASTeX61 in style 3 D + umte oTeAtohsclJournal Astrophysical The to Submitted enech ´ a led nw ntelbrtr n pc,ti okealsth enables work this space, and laboratory the in known already was Acpe etme ,2018) 4, September (Accepted , 1 hrSr 7 03 K¨oln, Germany 50937 77, Str. cher tpa Schlemmer Stephan ABSTRACT 3 206Mdi,Spain Madrid, E28006 23, tce o h smerctpNH top asymmetric the for etected eae moimioooousNH isotopologues ammonium terated , 2 and sa Asvany Oskar H 3 D + 2 n NHD and cnqe:spec- echniques: OPOLOGUES 2 3 + D n n three and one 2 + 3 D hl the While . + NH , future e 2 D 2 + 2 Domenech,´ Schlemmer and Asvany

1. INTRODUCTION wards the D atoms), thus making their detection feasi- Nitrogen is one of the most abundant elements in ble. Indeed, the detection of ammonium in space was the local universe, and has a notably rich chemistry, claimed through the assignment of an emission line cen- with more than seventy nitrogen-containing molecules tered at 262817 GHz (observed both in Orion IRc2 and − + identified in space to date (CDMS 2018). Two of the in Barnard B1-bS) to the 10 00 transition of NH3D most abundant nitrogen-bearing molecules in the inter- (Cernicharo et al. 2013). The laboratory rest frequency was derived from the analysis of the high resolution stellar medium (ISM) are N2 and NH3, which are pre- ν dicted to be present in many different media, from cold infrared spectrum of the 4 band, originally made by dark clouds (see, e.g., Nejad et al. (1990)) or protostel- Nakanaga & Amano (1986) and considerably improved lar cores (see, e.g., Aikawa et al. (2008)) to active galac- by Dom´enech et al. (2013). The rest frequency was later confirmed, and its uncertainty substantially decreased, tic nuclei (see, e.g., Harada et al. (2010)). N2, being a homonuclear diatomic molecule, can not be observed by the direct measurement of the rotational transition by Stoffels et al. (2016). through its (electric-dipole) rotation or vibration tran- + sitions. It can be observed, however, through electronic Besides their importance as tracers of NH4 , the transitions in the far-ultraviolet, with only one direct ob- deuterated isotopologues are also relevant regarding the servation reported by Knauth et al. (2004), or by its fre- subject of deuterium fractionation. In cold dark clouds, + when CO freezes onto grains and its abundance falls quently used proxy, the diazenylium ion (N2H ). Am- monia, in turn, is the first polyatomic molecule identified below that of HD, deuteration of other species becomes ∼ in space (Cheung et al. 1968), and since then has been efficient. It can be explained by the exothermicity ( + → + observed in many environments. To date, all deuterated 230 K) of the reaction H3 + HD H2D + H2 and sub- sequent deuteron transfer to other molecules. Similarly, isotopologues of ammonia, NH2D, NHD2 and ND3 have + + the exothermic reactions of H2D and HD2 with HD ef- been observed (see e.g. Harju et al. (2017) for a recent + + 15 15 ficiently produce HD and D (Hugo et al.2009) in case account), as well as NH3 and NH2D (Gerin et al. 2 3 of extreme depletion, and subsequent deuteron transfer 2009). + + to ammonia, e.g. NH3 + H2D → NH3D + H2, lead The astrophysical relevance of the ammonium ion, + + to the deuterated forms of NH4 . Furthermore, the NH4 , stems from its role as gas phase forebear of the ammonia molecule through its dissociative recombina- zero point energy differences of other reactions involv- tion. The full family of nitrogen hydrides is initiated by ing H or D atoms at temperatures of cold dark clouds + + (T ∼ 10 K) also enhance deuteration (Millar 2003). For the reaction N +H2 → NH + H (Zymak et al. 2013). the case of nitrogen, the reaction N+ + HD → ND+ + H Subsequent exothermic H2 abstraction reactions lead to + + + has a lower endothermicity than the corresponding the formation of NH2 , NH3 and NH4 , and their re- combination with electrons form the neutral hydrides reaction with H2 mentioned above (Grozdanov et al. 2016), leading to an enrichment in deuterated ammo- NH, NH2 and NH3. Ammonia is supposed to be mainly depleted onto grains at the temperatures of cold dark nium and ammonia isotopologues. The detected abun- clouds (T ∼ 10 K). On the other hand, the grains can dances of deuterated variants of ammonia are orders of magnitude higher compared to what is expected free NH3 molecules upon heating by shocks or in irradi- × −5 ated regions, increasing its gas-phase concentration, and based on the [D]/[H] abundance ratio (2.35 10 in + the local universe (Linsky et al. 2006)). Detection leading to appreciable quantities of NH4 , formed by the + of other deuterated ammonium isotopologues besides proton transfer from H3 to NH3. Actually, NH3 has one + of the highest proton affinities (PA) of simple interstellar NH3D will clearly help to constrain and understand the conditions of formation and distribution of ammo- molecules (PA[NH3]=8.85 eV) including H2 (which has + nia molecules, and, possibly, other nitrogen-containing PA[H2]=4.39 eV), so NH4 , once formed, remains stable prebiotic species. against further collisions with H2. + The only previous spectroscopic laboratory data on NH4 has tetrahedral symmetry with no permanent + the doubly and triply deuterated species (NH2D2 and electric dipole moment. Therefore, it is untraceable by + NHD3 ) that we are aware of are the recent publication radio astronomy through its rotational transitions. On + + of the analysis of the ν1 and ν6 infrared bands of NH2D the other hand, the deuterated isotopologues NH3D , 2 + + (Chang & Nesbitt 2018), and the communication of pre- NH2D and NHD have sizeable permanent electric 2 3 ν + dipole moments (approximately 0.26 D, 0.29 D and 0.24 liminary results on the analysis of the 1 band of NHD3 D, respectively) due to the separation between the cen- (Chang & Nesbitt 2013), both recorded at high resolu- ter of charge (which remains centered on the central N tion in a supersonic slit-jet discharge with an infrared atom) and the center of mass (which is displaced to- difference-frequency laser spectrometer. In this work, we Accurate rotational rest frequencies for ammonium ion isotopologues 3 present direct and accurate laboratory measurements of Table 1. + + the lowest frequency rotational transitions of NH3D , For NH2D2 , we based our search on the predictions + + NH2D2 and NHD3 . The fundamental rotational fre- from the high resolution infrared vibration-rotation + quency for NH3D has already been published previ- spectra of the ν1 and ν6 bands recently published by ously, and we only confirm its value, while our reported Chang & Nesbitt (2018). In total, we searched for eight + values for the other two isotopologues represent the first rotational lines of NH2D2 of which five were detected. available sub-mm wave data. All measured and predicted transitions are included in Table 1. Figure 1 shows two example measurements + 2. EXPERIMENTAL METHODS for NH2D2 . All lines have been measured at least 7 The rotational transitions of the deuterated ammo- times and fitted to Gaussian functions, from which line nium isotopologues have been measured in the K¨oln centers and widths were determined. The values quoted laboratories exploiting the rotational state dependence in Table 1 represent the combined mean and standard of the attachment of He atoms to cations (Br¨unken et al. deviation of all measurements. During the measure- (2014, 2017); Jusko et al. (2017); Dom´enech et al. ments, care was taken to avoid power broadening. By (2017, 2018)). The experiment was performed in the 4 K this, the linewidths are determined only by the Doppler trapping machine COLTRAP described by Asvany et al. broadening due to the kinetic temperature of the ions T (2010, 2014). The ions were generated in a storage ion in the trap (nominal temperature = 4 K), and also source by bombarding the precursor gas mixture with by a small contribution from the unresolved hyperfine 14 + splitting due to the quadrupole moments of the N and electrons (energy 30-40 eV). For generating NH3D and + 2H nuclei, both with spin I=1. Effectively, we measure NH2D2 , NH3 (Messer, 99.8 % purity) and a 1:5 mix- + T − 2 ture of D2 (Linde, 99.8 %) and He (Linde 99.999 %) linewidths corresponding to = 10 20 K for NH D2 , − were admitted to the ion source via two separate leak- depending on the observed line. The 10 00 line of + age valves. The approximate proportions were 1:1:5 NH3D seems to be more affected by hyperfine split- + + + ting, with an effective temperature T ≈ 30 K. for NH3D and 1:2:10 for NH2D2 . For NHD3 we applied a similar mixture using ND3 (Campro Scien- + tific, 99 % atom D) and H2 (Linde, 99.9999%). A For NHD3 , our search was guided by mass-scaling pulse of several ten thousand mass-selected parent ions the rotational constants from the other isotopologues, was injected into the 4 K cold 22-pole ion trap filled as well as by an estimation based on the available with about 1014 cm−3 He. During the trapping time IR-derived rotational constants from the Nesbitt group + of 750 ms, cation-helium complexes formed by three- (Chang & Nesbitt (2013)). In contrast to NH3D , we − body collisions. The detection of the resonant absorp- were able to measure not only the fundamental 10 00 + − − tion of the admitted cw (sub)millimeter radiation by for NHD3 , but also detected the 20 10 and 21 11 transitions. The reason is a more advantageous partition the trapped cold parent cations was achieved by ob- + serving the decrease of the number of cation-helium function for NHD3 , as well as more available millimeter- wave power at the lower transition frequencies of the complexes formed. For example, for recording a ro- + + + heavier NHD3 . tational transition of NH2D2 , the NH2D2 −He com- plexes are counted (typical counts are on the order of 2000) as a function of the (sub)millimeter-wave fre- 4. SPECTROSCOPIC PARAMETERS quency. The (sub)millimeter-wave radiation was sup- As accurate spectroscopic parameters and predictions + plied by a rubidium atomic clock-referenced synthesizer for the ground state of NH3D were already given by (Rohde&Schwarz SMF100A) driving a multiplier chain Stoffels et al. (2016) and Dom´enech et al. (2013), we source (Virginia Diodes, Inc.), covering the range 80- + + treat here only the isotopologues NH2D2 and NHD3 . 1100 GHz. + For NH2D2 , the measured frequencies of the pure rota- tional lines collected in Table 1 were fit together with 3. RESULTS the ground state combination differences derived from We started our measurements by optimizing the ex- Chang & Nesbitt (2018)(two duplicated entries had to perimental conditions using the known 10 − 00 line of be trimmed) using the program PGOPHER (Western + NH3D , which has been measured by Stoffels et al. 2017), rendering the final set of parameters for the (2016) with the same ion trap technique. We confirm ground state given in Table 2. Both sets of data their value for 10 − 00, as well as the inability of the ap- were weighted according to their different uncertain- plied setups to detect the 20 −10 and 21 −11 transitions. ties, in the kHz range for the present measurements, All experimental values of this work are summarized in and ∼ 7 MHz for the IR measurements, resulting in a 4 Domenech,´ Schlemmer and Asvany

fit. It may be also noted that the δk parameter is not significantly different from zero in this fit. + For NHD3 , we only have the three presently measured 1.05 transition frequencies available therefore we cannot pro- vide a least squares fit. Since, furthermore, the transi-

1.00 He counts He - tions obey ∆K = 0, we cannot determine C and DK . + 2 D 2 With the three experimental frequencies we have de-

0.95 termined B, DJ and DJK , shown in Table 2 together with the values quoted by Chang & Nesbitt (2013). Our

2 -1 normalized NH normalized

12 01 1 -0 01 00 quoted errors are obtained propagating the experimental

0.90 uncertainties through the closed linear relations exist- 455.5345 455.5350 248.9180 248.9185 ing between the three observed frequencies and the three frequency / GHz retrieved constants. We have assumed that there is no Figure 1. Example measurements of the 101 −000 and 212 − correlation between the experimental frequencies, thus + 101 rotational transitions of NH2D2 , recorded as depletion the uncertainties are just added in quadrature with the +− signal of the normalized NH2D2 He counts. Grey dots are appropriate coefficients. single measurements, the blue line is the average binned in 25 kHz steps, and the black trace is a Gaussian fit. Table 2. Derived spectroscopic parameters (in MHz) from a fit to all measured transitions. Numbers in parentheses are one standard deviation in units of the last digit. Table 1. Frequencies of pure rotational transitions (in MHz) + a of deuterated ammonium isotopologues. The final error is NH2D2 Parameter this work former work given in parentheses in units of the last digit. A 145600.69(33) 145601.7(12) B 118963.31(77) 118966.3(12) + J ′ ← J ′′ a NH3D K K this work former work C 103324.30(63) 103328.6(18) 10 ← 00 262816.8864(8) 262816.904(15) ∆J 1.698(46) 2.10(9) + ′ ′′ b NH2D J ′ ′ ← J ′′ ′′ this work former prediction 2 KaKc Ka Kc ∆JK -0.85(24) -1.9(4) ← c 110 101 42273.4715(47) 42270(6) ∆K 3.08(24) 3.87(3) ← 111 000 248918.2760(8) 248924(6) δJ 0.371(15) 0.18(6) ← 202 111 412755.5374(27) δK 0.01(20) 0.9(3) 212 ← 101 455534.5852(16) 455540(12) + b NHD3 Parameter this work former work ← 221 110 540046.0112(75) B 111117.9794(6) 111120.9(8) 220 ← 111 560809.9582(37) DJ 1.7539(2) 2.01(7) c 303 ← 212 632874.2(17) DJK -1.3763(11) -2.0(2) c 313 ← 202 656593.10(82) a + ′ ′′ Chang & Nesbitt (2018) NHD JK ← JK this work b 3 Chang & Nesbitt (2013) 10 ← 00 222228.9432(7)

20 ← 10 444415.7929(42) 5. CONCLUSION AND OUTLOOK 21 ← 11 444421.2980(18) We have measured, with kHz-level accuracies, rota- a Stoffels et al. (2016) tional transitions originating from the lowest energy lev- b Chang & Nesbitt (2018) + c els of the deuterated ammonium isotopologues NH3D , our prediction + + + NH2D2 and NHD3 . Prior to this work, NH3D was already identified in Orion IRc2 and Barnard B1b + + weighted standard deviation of the fit σw = 1.02. The (Cernicharo et al. 2013). For the NH2D2 and NHD3 accuracy of the spectroscopic parameters of this work is isotopologues these rest frequencies will facilitate the marginally improved (a factor of ∼ 2 − 3) with respect search for these ions in the interstellar medium. Cold to the former work, since only five rotational lines could dark clouds, highly depleted in C and O are the kind of be recorded in high resolution but eight parameters are sources where high deuterium fractionation is expected. Accurate rotational rest frequencies for ammonium ion isotopologues 5

+ The cold protostellar core B1b is a source where the higher symmetry (it is a symmetric top, like NH3D ) abundance of deuterated isotopologues of ammonia and and, therefore, a more favorable partition function, so + the diazenylium ion (N2H ) is quite remarkable, as it may also be soon within the reach of observational is the presence of other multiply deuterated species, radio astronomy. In any case the frequencies provided like D2CS (Marcelino et al. 2005) or D2S (Vastel et al. in this work will enable those detections. 2003). The chemical model used in the interpretation + of NH3D observations in B1b predicted that the abun- + This work (including a research stay of J.L.D. in dance of NH2D2 could be just a factor of two lower + K¨oln) has been supported by the Deutsche Forschungs- than that of NH3D . Therefore, one could expect the gemeinschaft (DFG) via SFB 956 project B2 and the detection of more deuterated species of the ammonium Ger¨atezentrum ”Cologne Center for Terahertz Spec- ion in this or other cold dark clouds. We note, how- + troscopy”. J.L.D. acknowledges partial financial support ever, that the partition function of NH2D , being an 2 from the Spanish MINECO through grant FIS2016- asymmetric top, is more unfavorable, and, even at the 77726-C3-1-P and from the European Research Coun- low temperature of B1b (T ∼ 12 K), spectral dilution cil through grant agreement ERC-2013-SyG-610256- is to be expected. The slightly higher electric dipole NANOCOSMOS. The authors gratefully acknowledge moment (0.29 D vs 0.26 D) just marginally improves the work done over the last years by the electrical and this situation. Although the necessary sensitivity for mechanical workshops of the I. Physikalisches Institut. this detection seems demanding, the ever increasing We thank Marius Hermanns, Frank Lewen and Matthias capabilities of receivers and observatories will likely + T¨opfer for assistance in the measurements. make it possible in a not distant future. NHD3 has a smaller abundance and dipole moment, but it has Software: PGOPHER (Western 2017)

REFERENCES

Aikawa, Y., Wakelam, V., Garrod, R. T., & Herbst, E. Dom´enech, J. L., Cueto, M., Herrero, V. J., et al. 2013, 2008, Astrophys. J., 674, 984 Astrophys. J., 771, L11 Asvany, O., Bielau, F., Moratschke, D., Krause, J., & Dom´enech, J. L., Jusko, P., Schlemmer, S., & Asvany, O. Schlemmer, S. 2010, Rev. Sci. Instr., 81, 076102 2018, Astrophys. J., 857, 61 Asvany, O., Br¨unken, S., Kluge, L., & Schlemmer, S. 2014, Dom´enech, J. L., Schlemmer, S., & Asvany, O. 2017, Appl. Phys. B, 114, 203 Astrophys. J., 849, 60 Br¨unken, S., Kluge, L., Stoffels, A., Asvany, O., & Gerin, M., Marcelino, N., Biver, N., et al. 2009, Astron. Schlemmer, S. 2014, Astrophys. J. Lett., 783, L4 Astrophys., 498, L9 Br¨unken, S., Kluge, L., Stoffels, A., P´erez-R´ıos, J., & Grozdanov, T. P., McCarroll, R., & Roueff, E. 2016, Schlemmer, S. 2017, Journal of Molecular Spectroscopy, Astron. Astrophys., 589, A105 332, 67 Harada, N., Herbst, E., & Wakelam, V. 2010, Astrophys. J., CDMS. 2018, Cologne Database for Molecular 721, 1570 Spectroscopy. Molecules in the Interstellar Medium or Harju, J., Daniel, F., Sipil¨a, O., et al. 2017, Astron. Circumstellar Shells (as of 05/2018), University of Astrophys., 600, A61 Cologne, Hugo, E., Asvany, O., & Schlemmer, S. 2009, J. Chem. https://www.astro.uni-koeln.de/cdms/molecules Phys., 130, 164302 Cernicharo, J., Tercero, B., Fuente, A., et al. 2013, Jusko, P., Stoffels, A., Thorwirth, S., et al. 2017, J. Mol. Astrophys. J., 771, L10 Spectrosc., 332, 59 Chang, C.-H., & Nesbitt, D. J. 2013, in 68th OSU Knauth, D. C., Andersson, B.-G., McCandliss, S. R., & International Symposium on Molecular Spectroscopy, Warren Moos, H. 2004, Nature, 429, 636 Contribution MG03. Linsky, J. L., Draine, B. T., Moos, H. W., et al. 2006, https://kb.osu.edu/handle/1811/55596 Astrophys. J., 647, 1106 Chang, C.-H., & Nesbitt, D. J. 2018, J. Chem. Phys., 148, Marcelino, N., Cernicharo, J., Roueff, E., Gerin, M., & 014304 Mauersberger, R. 2005, Astrophys. J., 620, 308 Cheung, A. C., Rank, D. M., Townes, C. H., Thornton, Millar, T. J. 2003, Space Sci. Rev., 106, 73 D. D., & Welch, W. J. 1968, Phys. Rev. Lett., 21, 1701 Nakanaga, T., & Amano, T. 1986, Can. J. Phys., 64, 1356 6 Domenech,´ Schlemmer and Asvany

Nejad, L., Williams, D. A., & Charnley, S. B. 1990, Mon. Western, C. M. 2017, Journal of Quantitative Spectroscopy Not. R. Astron. Soc., 246, 186 and Radiative Transfer, 186, 221 Stoffels, A., Kluge, L., Schlemmer, S., & Br¨unken, S. 2016, Zymak, I., Hejduk, M., Mulin, D., et al. 2013, Astrophys. Astron. Astrophys., 593, A56 Vastel, C., Phillips, T. G., Ceccarelli, C., & Pearson, J. J., 768, 86 2003, Astrophys. J., 593, L97