ORIGINAL RESEARCH article

Front. Astron. Space Sci., 02 February 2022

Sec. Astrochemistry

Volume 8 - 2021 | https://doi.org/10.3389/fspas.2021.805162

Rotational Rest Frequencies and First Astronomical Search of Protonated Methylamine

  • 1. I. Physikalisches Institut, Universität zu Köln, Köln, Germany

  • 2. Max Planck Institute for Extraterrestrial Physics, Garching bei München, Germany

Article metrics

View details

1

Citations

3,3k

Views

1k

Downloads

Abstract

We report first laboratory rest frequencies for rotational transitions of protonated methylamine, , measured in a cryogenic 22-pole ion trap machine and employing an action spectroscopy scheme. For this prolate symmetric top molecule thirteen transitions between 80 and 240 GHz were detected in the ground vibrational state, covering JK = 2K − 1K up to JK = 6K − 5K with K = 0, 1, 2. Some transitions exhibit noticeable structure that is attributed to internal rotation splitting. As the CN radical and several of its hydrogenated and protonated forms up to methylamine, CH3NH2, are well known entities in the laboratory and in space, protonated methylamine, , is a promising candidate for future radio astronomical detection.

1 Introduction

The cyano radical, CN, was one of the first molecules detected in the interstellar medium (McKellar, 1940). It is astrochemically linked to its hydrogenated and protonated forms, most of which have been detected in the interstellar medium, typically in the galactic center source Sgr B2, such as HCN (Snyder and Buhl, 1971), HNC (Snyder and Buhl, 1972; Zuckerman et al., 1972), HCNH+ (Ziurys and Turner, 1986), and H2C = NH (Godfrey et al., 1973). The terminal product of this hydrogenation series is methylamine, CH3NH2, which has been well characterized by laboratory spectroscopy (Ohashi et al., 1987; Ilyushin and Lovas, 2007; Motiyenko et al., 2014) and was first detetcted in Sgr B2 and Orion A in 1974 by Kaifu et al. (1974) and Fourikis et al. (1974).

In interstellar environments, molecules may also occur in their protonated forms, generated by a proton transfer from a proton donor like , e.g. CO + + H2. Many of the simpler protonated species have been detected in the interstellar medium, such as HCO+ (Buhl and Snyder, 1970), N2H+ (Turner, 1974), HCNH+ (Ziurys and Turner, 1986) and HOCO+ (Thaddeus et al., 1981), whereas many of the more complex ones so far are only suspected to be present. Examples are , as a product of proton transfer to methyl amine (see, e.g., the KIDA database; Wakelam et al., 2015) or from radiative association of NH3 and (Herbst, 1985), but also protonated methanol, (Jusko et al., 2019), and protonated methane, (Asvany et al., 2012; Asvany et al., 2015). In particular the latter two molecular ions have not yet been searched for in space because their laboratory microwave spectra are predicted to be quite irregular and still not known.

Following up on our recent work on the high-resolution rotational spectra of CN+ (Thorwirth et al., 2019a) and (Markus et al., 2019), in this study, we finally focus on the very last member of the CN hydrogenation/protonation chain, protonated methylamine, . Upon protonation, methylamine, CH3NH2, an asymmetric top molecule featuring two internal large-amplitude motions (internal rotation and inversion), is converted into a much simpler symmetric top of C3v point group symmetry, in which the CH3 and the NH3 groups assume a staggered configuration at the global energy minimum (see inset in Figure 1). Literature on the spectroscopy of is extremely sparse, with only one experimental paper presenting low-resolution infrared (IR) features of the Ar-tagged species (Michi et al., 2003), and one work reporting ab initio values for its IR vibrational frequencies (Zeroka and Jensen, 1998). Therefore, the spectroscopic study in the present work was complemented with new high-level quantum-chemical predictions of the molecular structure and force field (see Section 3). The ion trap experiment is briefly described in Section 2. The experimentally derived rotational transitions of as well as the ground state spectroscopic parameters are summarized in Section 4. A first astronomical search of towards Sgr B2(N) and Sgr B2(M) is finally presented in Section 5.

FIGURE 1

2 Experimental Methods

The rotational transitions of have been measured using an action spectroscopic method which exploits the rotational state dependence of the ternary attachment of He atoms to cations at low temperature (Brünken et al., 2014; Brünken et al., 2017; Doménech et al., 2017; Jusko et al., 2017; Doménech et al., 2018a; Doménech et al., 2018b; Thorwirth et al., 2019a; Asvany et al., 2021). The experiment was performed in the 4 K ion trapping machine COLTRAP described by Asvany et al. (2010, 2014). The ions were generated in a storage ion source by electron impact ionization (Ee ≈ 26–30 eV) of the precursor gas mixture. This mixture consisted of methylamine, CH3NH2 (Aldrich Chem. Corporation, CAS 74-89-5, 98%), and helium (Linde 5.0) which were admitted to the ion source via two separate leakage valves. is generated by a reaction of the type + CH3NH2 + CNH4. A pulse of several ten thousand mass-selected parent ions (m = 32 u) was injected into the 22-pole ion trap filled with about 1014 cm−3 He at 4 K. At the beginning of the trapping time lasting 800 ms, -He complexes formed by three-body collisions with He. The resonant absorption of the cw millimeter-wave radiation by the trapped cold cations is detected by observing the decrease of the number of -He complexes. A rotational line is thus recorded by repeating these trapping cycles (1 Hz) and counting the mass-selected -He complexes (m = 36 u) as a function of the millimeter-wave frequency. This millimeter-wave radiation was supplied by a multiplier chain source (Virginia Diodes, Inc.), covering the ranges 80–125 and 170–1100 GHz. This source was driven by a synthesizer (Rohde&Schwarz SMF100A), which was referenced to a rubidium atomic clock. The beam of the millimeter-wave source has been directed toward the ion trap via an elliptical mirror and a thin diamond vacuum window. The 160 GHz line of was measured with a different multiplier chain (Radiometer Physics GmbH), covering the range of 110–170 GHz.

3 Computational Methods

As no previous high-resolution experimental data were available for , spectroscopic searches and analysis were based on high-level quantum-chemical calculations performed here at the CCSD(T) level of theory (Raghavachari et al., 1989). Equilibrium geometries were calculated using analytic gradient techniques (Watts et al., 1992) and Dunning’s correlation-consistent basis sets as large as cc-pwCVQZ (Peterson and Dunning, 2002). Anharmonic force fields to evaluate the zero-point vibrational contributions (= ΔA0, ΔB0) to the equilibrium rotational constants were calculated using analytic second-derivative techniques (Gauss and Stanton, 1997; Stanton and Gauss, 2000) followed by additional numerical differentiation to calculate the third and fourth derivatives needed for the anharmonic force field (Stanton et al., 1998; Stanton and Gauss, 2000). These calculations were carried out using the frozen core (fc) approximation in combination with the cc-pVTZ basis set (Dunning, 1989).

All calculations were performed using the CFOUR program (Matthews et al., 2020) and strategies summarized elsewhere (Puzzarini et al., 2010). The computed spectroscopic parameters are listed in Table 1 that also provides scaled (best-estimate) parameters obtained from a comparison of experimental and calculated values of isoelectronic ethane, C2H6 (see, e.g., Martinez et al. (2013) for a similar procedure used in the vinyl acetylene/protonated vinyl cyanide family of isoelectronic species). Additional results from the calculations are given in the Supplementary Material. An estimate of the torsional barrier height was obtained using the energy difference between the staggered and the eclipsed forms of calculated at the CCSD(T)/cc-pwCVQZ level and under consideration of harmonic zero-point vibrational contributions calculated at the CCSD(T)/cc-pVTZ level. This procedure results in a barrier of V3 = 1.98 kcal/mol. A similar calculation of ethane, C2H6, yields 2.45 kcal/mol to be compared against an experimental value of 3.00 kcal/mol (Borvayeh et al., 2008). Using the ethane exp/calc-ratio of V3 for the purpose of scaling, a best estimate value of 2.42 kcal/mol (846 cm−1) for the torsional barrier in is obtained.

TABLE 1

ParameterC2H6
exp.acalc.bcalc.bscaledcexp. fit Iexp. fit IId
Ae81 080.84285 189.126
Be20 126.24820 246.782
  ΔA01065.2081222.937
  ΔB0265.982341.653
A080 391.1(38)81 080.84283 966.18983 251.983 251.9e
B019 913.684(24)19 860.26619 905.12919 958.719 920.6224(12)19 941.65(65)
DJ × 10331.026(21)30.5434.8135.435.378(27)35.292(22)
DJK × 10376.37(28)77.8284.4182.879.48(38)87.89(34)
+0.158
μe/D2.15
ρ0.50.540.54d
V3/kcal/mol3.002.451.982.422.42d
F × 10–380.3911f335.152 576f335.152 576f
F3J−397.8(20)−243.0(75)
wrmsg1.440.96

Calculated and experimental spectroscopic parameters (in MHz, unless noted otherwise) of . Two experimental fits to the measured data are presented, one (fit I) ignoring the torsional motion, and one (fit II) including torsional motion. Numbers in parentheses are one standard deviation in units of the last digit.

b This study; Ae, Be, eQq(14N) and dipole moment μe calculated at the CCSD(T)/cc-pwCVQZ, level, zero-point vibrational corrections ΔA0, ΔB0 and centrifugal distortion terms calculated at the CCSD(T)/cc-pVTZ, level. Torsional barriers V3 evaluated from the CCSD(T)/cc-pwCVQZ energy difference of the staggered and eclipsed forms complemented with harmonic vibrational corrections evaluated at the CCSD(T)/cc-pVTZ level.

c Calculated value further scaled with the ratio Xexp/Xcalc of the corresponding parameter of C2H6.

d Fit performed using an in-house program assuming a uniform uncertainty of 15 kHz.

e Kept fixed.

f C2H6: Kept fixed at value of A rotational constant; : Kept fixed at the customary value of A/(ρ(1 − ρ)).

g Weighted rms, dimensionless.

4 Laboratory Results and Spectroscopic Parameters

Using frequency predictions based on the new high-level quantum-chemical calculations, the lines shown in Figure 1 were found subsequently during targeted spectroscopic survey scans. The spectra were recorded in individual measurements in which the frequency was stepped in an up-and-down manner several times. The frequency steps were typically 2 kHz, except for the JK = 4K ← 3K measurement at 159 GHz, for which 10 kHz steps were applied. Such individual measurements were repeated typically ten times. In the depiction of Figure 1, all available spectra were accumulated and rebinned to a stepwidth of 10 kHz. This deep integration was necessary due to the comparably small signal strength (maximum depletion upon photon absorption is on the order of 2%), and due to somewhat noisy signal counts, which had its origin in ion source instabilities caused by the sticky consistency of the methylamine precursor. Also, during the measurements, care was taken to avoid power broadening. Thus, the linewidths were expected to be dominated by the Doppler broadening due to the kinetic temperature of the ions in the trap (nominal temperature T = 4 K with some residual heating to typically 8 K), and a possible contribution due to non-resolved hyperfine splitting. The line patterns detected for a given rotational transition were in line with those expected for a prolate symmetric top molecule and at first sight, no further spectroscopic complexity, as could be expected for resolved nitrogen quadrupole hyperfine structure or torsional motion between the CH3 and NH3 subunits, was detected in individual measurements at our experimental conditions. For all individual measurements, the detected lines were fitted to Gaussian functions, from which line centers and their uncertainties were determined, and the first results were very similar to the frequencies quoted in Table 2 (Model I).

TABLE 2

FrequencyModel IModel II
ococ
79681.065(10)21110.02621,011,00.019
79681.364(10)20100.00720,010,00.010
119517.996(10)3222−0.01032,022,00.006
119519.449(10)31210.01331,021,00.000
119519.449(10)31,121,10.017
119519.621(10)31, −121, −10.008
119519.908(10)3020−0.00530,020,0−0.002
119520.077(10)30, ±120, ±1−0.006
159353.362(30)4232−0.01742,032,00.000
159355.297(10)41310.01141,031,0−0.010
159355.297(10)41,131,10.012
159355.537(10)41, −131, −10.011
159355.899(10)4030−0.02340,030,0−0.023
159356.151(10)40, ±130, ±1−0.002
199185.325(30)5242−0.03152,042,0−0.015
199187.748(10)51410.00851,041,0−0.024
199187.748(10)51,141,10.004
199188.039(10)51, −141, −1−0.007
199188.518(10)5040−0.01750,040,0−0.022
199188.820(10)50, ±140, ±1−0.009
239016.015(30)61510.06761,051,00.019
239016.931(30)60500.02960,050,00.013

Frequencies of pure rotational transitions (in MHz) of , obtained by fitting multiple Gaussians to the spectra in Figure 1. The uncertainties (1 σ) are given in parentheses in units of the last significant digits. The assignment of Model I assumes a standard symmetric top Hamiltonian whereas that of Model II assumes a symmetric top featuring torsional splitting, see text for details.

The measured frequencies of the pure rotational lines collected in Table 2 (Model I) were first fit with a standard symmetric rotor Hamiltonian using the PGOPHER program (Western, 2017), yielding the set of parameters for the ground state given in Table 1 (fit I). Since the transitions of a symmetric rotor obey the ΔK = 0 selection rule, A and DK cannot be determined experimentally, and only B0, DJ and DJK were derived in the least-squares fitting procedure and are presented in Table 1. Computed parameters are also included in Table 1 for comparison. As can be seen, the overall agreement between the calculated and scaled best-estimate values and the experimental values is very good. A simulated stick spectrum and its convolution, based on fit I of Table 1 and accounting for the hyperfine structure due to the quadrupole moment of the nucleus (with spin I = 1), are included in Figure 1. As the computed value = +158 kHz is very small (see also Table 1), it is not surprising that the hyperfine structure is not resolved in our experiment.

Closer inspection after co-addition of all spectra revealed some noticeable blue-shifted shoulders, as discernible in the accumulation of Figure 1. These are particularly evident in the spectra at 119 (J = 3K − 2K) and 159 GHz (J = 4K − 3K), but are also discernible in the weaker spectrum at 199 GHz (J = 5K − 4K). In the spectrum at 119 GHz (second panel in Figure 1 and close-up in Figure 2) both blue shoulders have an offset of about +170 kHz relative to the main peaks (K = 0 and 1), and this offset is found to increase for the higher frequency transitions at 159 (230 kHz) and 199 GHz (290 kHz). While a detailed analysis is hampered by the overall poor signal-to-noise ratio of the spectra, a first analysis was performed here under the assumption that the spectra are affected by internal rotation, an effect that has rarely been observed previously for symmetric top molecules in their ground vibrational states (see, e.g., Ozier and Moazzen-Ahmadi, 2007, and references therein). As indicated in Figure 2, the structure of the J = 3K − 2K spectrum may be decomposed into four features, two strong ones (the putative K = 0 and 1 components used in the first fitting procedure in Table 1, fit I) each of which is accompanied by a weaker satellite that is found blue-shifted relative to the main component. Similar patterns are observed for the J = 4K − 3K and J = 5K − 4K transitions. Ideally, a suitable model Hamiltonian should be able to reproduce both the magnitude of the torsional splitting as well as the intensities of the individual spectroscopic components. First, the magnitude of the torsional splitting to be expected in was estimated based on the so-called hybrid approach described elsewhere (Wang et al., 2001) using the dominant parameters (V3, ρ = Iα/Ic, F = A/(ρ × (1 − ρ)), F3J, F3K) as calculated here and complemented with parameters taken from isoelectronic ethane (Borvayeh et al., 2008). Using this approach, the splitting was predicted in very good agreement with the experimental values. In analogy to other closely related symmetric top molecules such as methyl silane, CH3SiH3, each rotational transition of with K = 1, 2 is expected to split in up to three components (σ = 0, +1, −1) whereas transitions with K = 0 only split into two, σ = 0, ±1 (Pelz et al., 1992; Ozier and Moazzen-Ahmadi, 2007). The symmetries and nuclear spin weights under consideration of torsional splitting that are needed for intensity estimates have been given, for example, in Pelz et al. (1992). From this, it is concluded that in the -spectra two strong torsional components of the K = 1-transitions (σ = 0, 1) are too close in frequency to be spectroscopically resolved and thus only one strong (superposition of σ = 0 and 1) and one weaker torsional component (σ = − 1) are observed. The K = 0 splitting is also (partly) resolved and both components assigned, a stronger (σ = 0) and a weaker one (σ = ±1), respectively. For the K = 2 transitions, only one component (σ = 0) can be assigned in the spectra obtained here with some confidence. Using this spectroscopic knowledge, the accumulated spectra as depicted in Figure 1 were refitted using multiple Gaussian components, and the final frequencies and their assignments are listed in Table 2, Model II. With this improved spectroscopic assignment, a second fit was performed (using an in-house program and assuming a uniform frequency uncertainty of 15 kHz) leading to an alternative set of molecular parameters given in Table 1 (fit II). In this approach, owing to the limited amount of spectroscopic information available to constrain the internal rotation problem more rigorously, several parameters required in the theoretical description were kept fixed at values calculated quantum-chemically (ρ, F, V3). In analogy to fit I, B0, DJ, and DJK were varied in the least squares adjustment and additionally F3J was released, resulting in a good fit result and agreement with the fit I and C2H6 parameter sets. If released also, V3 cannot be determined statistically from the present data set and changing its value from the scaled calculated value of 2.42 kcal/mol to the unscaled value of 1.98 kcal/mol has little impact on the overall fit quality but decreases the value of F3J significantly to −87(3) MHz. Based on this new set of parameters and taking into consideration the proper spin statistical effects (e.g., Pelz et al., 1992), new simulations of spectra have been obtained that are shown in Figure 3. The agreement between the experimental spectra and the simulations is rather compelling and hence speaking very much in favor of internal rotation as to the cause of the peculiar lineshapes and splittings observed.

FIGURE 2

FIGURE 3

5 Interstellar Search

The spectral line transitions reported in Table 2 have frequencies GHz and constitute good targets to be searched for with facilities like ALMA (Atacama Large Millimeter/sub-millimeter Array; ALMA Partnership et al., 2015). While the low-energy transitions are not covered by current ALMA bands, the J = 4K → 3K and JK = 5K → 4K lines can be observed with the new ALMA band 5 receiver (covering a frequency range from 159 to 211 GHz), and the JK = 6K → 5K transitions are observable with the commonly used band 6 receiver (211–276 GHz).

The detection of species closely related to protonated methylamine such as HCN, HCNH+ or CH3NH2 in the star-forming region Sgr B2 (e.g., Schilke et al., 1991; Belloche et al., 2013) motivates the search for protonated methylamine in this region. We have made use of the ALMA spectral line survey published in Sánchez-Monge et al. (2017, see also Schwörer et al., 2019). The observations target the star-forming objects Sgr B2(N) and Sgr B2(M) and cover the whole ALMA band 6 with a spectral resolution of 0.7 km s−1 and an angular resolution of 0.′′4 (corresponding to 3300 au at the distance of the source). This high angular resolution allows to resolve the regions in more than 40 different dense cores (see Sánchez-Monge et al., 2017). Figure 4 presents the spectra extracted towards different selected cores in both regions (see source coordinates in Tables 1, 2 of Sánchez-Monge et al., 2017). The spectra are obtained after averaging the emission inside the 3σ polygon that defines the source size (between 0.′′5 and 0.′′7). We mark with vertical lines the location of the transitions.

FIGURE 4

We have selected objects in different evolutionary stages and with different physical properties. Cores AN02, AN17 and AM03 correspond to dense cores dominated by dust but with a rich chemistry. In particular, core AN02 has been extensively studied in the literature in the search of new chemical species (e.g., Belloche et al., 2014). We have also included objects in which an embedded Hii region has been found (e.g., AN10, AM06, AM15). The presence of the UV radiation from embedded massive stars ionizing the gas and resulting in Hii regions can enhance the protonation of methylamine in the photon-dominated region around the Hii region. As shown in Figure 4, no obvious features are detected at the frequencies of the tabulated transitions, suggesting a low abundance of protonated methylamine in Sgr B2. We have determined an upper limit of cm−2 to the column density of protonated methylamine. This is estimated based on a 3-sigma noise upper limit and assuming a gas temperature of 200–300 K (as derived from other molecular species in Sgr B2, Schwörer et al., 2019). This translates into an upper limit of the fractional abundance of ∼10–10. The high densities of the Sgr B2 region (∼ 105–108 cm−3; see also Schmiedeke et al., 2016) can rapidly attenuate the radiation field and result in a low production of heavy ions like . Moreover, the J = 6K → 5K transitions are located close to bright CH3CN and CN (JK = 13K → 12K) transitions and a strong SO2 feature which may hinder the detection of the weaker in chemically rich, embedded objects. Observations at other frequency ranges, as well as in other astronomical sources associated with highly ionized and less dense gas, may help in the detection of this species in space.

6 Conclusion and Outlook

In the present study, the pure rotational spectrum of protonated methyl amine, , was observed for the first time. It is worthwhile to recall that upon protonation of CH3NH2, forming , the complex dynamical behaviour of the former molecular system is largely simplified, because the inversion at the NH2 group in CH3NH2 is eliminated and only the torsional motion between CH3 and NH3 has to be potentially considered for the pure rotational spectrum of in its ground vibrational state. In this study, a total of five pure rotational transitions from J = 2K − 1K to J = 6K − 5K and K up to 2 were identified. As indicated through complementary quantum-chemical calculations, explicit treatment of nuclear quadrupole hyperfine structure from the presence of the nitrogen nucleus was found to be negligible. While the strongest spectroscopic features in the spectrum can be assigned and fitted reasonably well in a straightforward fashion using a standard symmetric-top Hamiltonian, peculiar lineshapes and weak substructure identified in the spectra upon close inspection indeed required extension of the theoretical treatment to account for internal rotation. This refined model description permits very convincing reproduction of the experimental line profiles. A more comprehensive treatment of the internal rotational problem in would certainly benefit from extension of the present work to a higher degree of rotational excitation but also from detection of vibrational satellites of in torsionally excited states.

Future work towards the spectroscopic characterization of protonated amines offers many possibilities, not only in the millimeter-wave regime. For example, the low-resolution vibrational spectra of and also protonated ethyl amine, , were observed recently in the range from 700 to 1750 cm−1 (Thorwirth et al., 2019b) using the FELion ion trap apparatus (Jusko et al., 2019). For the purpose of millimeter-wave radio astronomical searches of the data presented here should already suffice. Although the column density of in the interstellar medium may be some orders of magnitude lower than that of CH3NH2, the more favorable partition function of will support radio astronomical detectability. In the search towards SgrB2 presented in this work was unfortunately not found.

Statements

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

PS: measurement, manuscript writing MT: measurement ST: quantum chemical computations, data evaluation, manuscript writing CE: data evaluation, manuscript writing AS-M: data evaluation, manuscript writing AS: data evaluation, manuscript writing PS: supervision SS: supervision, manuscript writing OA: measurement, data evaluation, manuscript writing, supervision.

Funding

This work has been supported via Collaborative Research Centre 956, sub-projects A6 and B2, funded by the Deutsche Forschungsgemeinschaft (DFG, project ID 184 018 867), as well as DFG SCHL 341/15-1 (Gerätezentrum “Cologne Center for Terahertz Spectroscopy”).

Acknowledgments

The authors gratefully acknowledge the work done over the last years by the electrical and mechanical workshops of the I. Physikalisches Institut. We thank Frank Lewen for assistance using the mm-wave setup. CE is grateful for support from the Max Planck Society.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2021.805162/full#supplementary-material

References

  • 1

    AsvanyO.BielauF.MoratschkeD.KrauseJ.SchlemmerS. (2010). Note: New Design of a Cryogenic Linear Radio Frequency Multipole Trap. Rev. Sci. Instrum.81, 076102. 10.1063/1.3460265

  • 2

    AsvanyO.BrünkenS.KlugeL.SchlemmerS. (2014). COLTRAP: A 22-pole Ion Trapping Machine for Spectroscopy at 4 K. Appl. Phys. B114, 203211. 10.1007/s00340-013-5684-y

  • 3

    AsvanyO.KriegJ.SchlemmerS. (2012). Frequency Comb Assisted Mid-Infrared Spectroscopy of Cold Molecular Ions. Rev. Sci. Instrum.83, 093110. 10.1063/1.4753930

  • 4

    AsvanyO.MarkusC. R.RoucouA.SchlemmerS.ThorwirthS.LauzinC. (2021). The Fundamental Rotational Transition of NO+. J. Mol. Spectrosc.378, 111447. 10.1016/j.jms.2021.111447

  • 5

    AsvanyO.YamadaK. M. T.BrünkenS.PotapovA.SchlemmerS. (2015). Experimental Ground-State Combination Differences of CH5+. Science347, 13461349. 10.1126/science.aaa3304

  • 6

    BellocheA.GarrodR. T.MüllerH. S. P.MentenK. M. (2014). Detection of a Branched Alkyl Molecule in the Interstellar Medium: Iso-Propyl Cyanide. Science345, 15841587. 10.1126/science.1256678

  • 7

    BellocheA.MüllerH. S. P.MentenK. M.SchilkeP.ComitoC. (2013). Complex Organic Molecules in the Interstellar Medium: IRAM 30 m Line Survey of Sagittarius B2(N) and (M). Astron. Astrophys.559, A47. 10.1051/0004-6361/201321096

  • 8

    BorvayehL.Moazzen-AhmadiN.HornemanV.-M. (2008). The ν12ν9 Band of Ethane: A Global Frequency Analysis of Data from the Four Lowest Vibrational States. J. Mol. Spectrosc.250, 5156. 10.1016/j.jms.2008.04.009

  • 9

    BrünkenS.KlugeL.StoffelsA.AsvanyO.SchlemmerS. (2014). Laboratory Rotational Spectrum of l-C3H+ and Confirmation of its Astronomical Detection. Astrophys. J. Lett.783, L4. 10.1088/2041-8205/783/1/l4

  • 10

    BrünkenS.KlugeL.StoffelsA.Pérez-RíosJ.SchlemmerS. (2017). Rotational State-dependent Attachment of He Atoms to Cold Molecular Ions: An Action Spectroscopic Scheme for Rotational Spectroscopy. J. Mol. Spectrosc.332, 6778. 10.1016/j.jms.2016.10.018

  • 11

    BuhlD.SnyderL. E. (1970). Unidentified Interstellar Microwave Line. Nature228, 267269. 10.1038/228267a0

  • 12

    DoménechJ. L.JuskoP.SchlemmerS.AsvanyO. (2018a). The First Laboratory Detection of Vibration-Rotation Transitions of 12CH+ and 13CH+ and Improved Measurement of Their Rotational Transition Frequencies. Astrophys. J.857, 61. 10.3847/1538-4357/aab36a

  • 13

    DoménechJ. L.SchlemmerS.AsvanyO. (2017). Accurate Frequency Determination of Vibration-Rotation and Rotational Transitions of SiH+. Astrophys. J.849, 60. 10.3847/1538-4357/aa8fca

  • 14

    DoménechJ. L.SchlemmerS.AsvanyO. (2018b). Accurate Rotational Rest Frequencies for Ammonium Ion Isotopologues. Astrophys. J.866, 158. 10.3847/1538-4357/aadf83

  • 15

    DunningT. H. (1989). Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen. J. Chem. Phys.90, 10071023. 10.1063/1.456153

  • 16

    ALMA PartnershipFomalontE. B.VlahakisC.CorderS.RemijanA.BarkatsD.et al (2015). The 2014 ALMA Long Baseline Campaign: An Overview. Astrophys. J. Lett.808, L1. 10.1088/2041-8205/808/1/L1

  • 17

    FourikisN.TakagiK.MorimotoM. (1974). Detection of Interstellar Methylamine by its 202 → 110Aa-State Transition. Astrophys. J. Lett.191, L139. 10.1086/181570

  • 18

    GaussJ.StantonJ. F. (1997). Analytic CCSD(T) Second Derivatives. Chem. Phys. Lett.276, 7077. 10.1016/s0009-2614(97)88036-0

  • 19

    GodfreyP. D.BrownR. D.RobinsonB. J.SinclairM. W. (1973). Discovery of Interstellar Methanimine (Formaldimine). Astrophys. Lett.13, 119.

  • 20

    HerbstE. (1985). The Rate of the Radiative Association Reaction between CH3+ and NH3 and its Implications for Interstellar Chemistry. Astrophys. J.292, 484486. 10.1086/163179

  • 21

    IlyushinV.LovasF. J. (2007). Microwave Spectra of Molecules of Astrophysical Interest. Xxv. Methylamine. J. Phys. Chem. Ref. Data36, 11411276. 10.1063/1.2769382

  • 22

    JuskoP.BrünkenS.AsvanyO.ThorwirthS.StoffelsA.van der MeerL.et al (2019). The FELion Cryogenic Ion Trap Beam Line at the FELIX Free-Electron Laser Laboratory: Infrared Signatures of Primary Alcohol Cations. Faraday Discuss.217, 172202. 10.1039/c8fd00225h

  • 23

    JuskoP.StoffelsA.ThorwirthS.BrünkenS.SchlemmerS.AsvanyO. (2017). High-resolution Vibrational and Rotational Spectroscopy of CD2H+ in a Cryogenic Ion Trap. J. Mol. Spectrosc.332, 5966. 10.1016/j.jms.2016.10.017

  • 24

    KaifuN.MorimotoM.NaganeK.AkabaneK.IguchiT.TakagiK. (1974). Detection of Interstellar Methylamine. Astrophys. J. Lett.191, L135L137. 10.1086/181569

  • 25

    MarkusC. R.ThorwirthS.AsvanyO.SchlemmerS. (2019). High-Resolution Double Resonance Action Spectroscopy in Ion Traps: Vibrational and Rotational Fingerprints of CH2NH2+. Phys. Chem. Chem. Phys.21, 2640626412. 10.1039/c9cp05487a

  • 26

    MartinezO.Jr.LattanziV.ThorwirthS.McCarthyM. C. (2013). Detection of Protonated Vinyl Cyanide, CH2CHCNH+, a Prototypical Branched Nitrile Cation. J. Chem. Phys.138, 094316. 10.1063/1.4793316

  • 27

    MatthewsD. A.ChengL.HardingM. E.LippariniF.StopkowiczS.JagauT.-C.et al (2020). Coupled-Cluster Techniques for Computational Chemistry: The CFOUR Program Package. J. Chem. Phys.152, 214108. 10.1063/5.0004837

  • 28

    McKellarA. (1940). Evidence for the Molecular Origin of Some Hitherto Unidentified Interstellar Lines. Publ. Astron. Soc. Pac.52, 187. 10.1086/125159

  • 29

    MichiT.OhashiK.InokuchiY.NishiN.SekiyaH. (2003). Infrared Spectra and Structures of (CH3NH2)nH+ (n = 1-4). Binding Features of an Excess Proton. Chem. Phys. Lett.371, 111117. 10.1016/s0009-2614(03)00217-3

  • 30

    MotiyenkoR. A.IlyushinV. V.DrouinB. J.YuS.MargulèsL. (2014). Rotational Spectroscopy of Methylamine up to 2.6 THz. Astron. Astrophys.563, A137. 10.1051/0004-6361/201323190

  • 31

    OhashiN.TakagiK.HougenJ. T.OlsonW. B.LaffertyW. J. (1987). Far-Infrared Spectrum and Ground State Constants of Methyl Amine. J. Mol. Spectrosc.126, 443459. 10.1016/0022-2852(87)90249-9

  • 32

    OzierI.Moazzen-AhmadiN. (2007). “Internal Rotation in Symmetric Tops,” in Advances in Atomic, Molecular, and Optical Physics. Editors BermanP. R.LinC. C.ArimondoE. (Amsterdam: Elsevier), 54, 423509. 10.1016/s1049-250x(06)54007-8

  • 33

    PelzG.MittlerP.YamadaK. M. T.WinnewisserG. (1992). Millimeter-Wave Spectra and Revised Spectroscopic Parameters of Methylsilane. J. Mol. Spectrosc.156, 390402. 10.1016/0022-2852(92)90240-o

  • 34

    PetersonK. A.DunningT. H. (2002). Accurate Correlation Consistent Basis Sets for Molecular Core-Valence Correlation Effects: The Second Row Atoms Al-Ar, and the First Row Atoms B-Ne Revisited. J. Chem. Phys.117, 1054810560. 10.1063/1.1520138

  • 35

    PuzzariniC.StantonJ. F.GaussJ. (2010). Quantum-Chemical Calculation of Spectroscopic Parameters for Rotational Spectroscopy. Int. Rev. Phys. Chem.29, 273367. 10.1080/01442351003643401

  • 36

    RaghavachariK.TrucksG. W.PopleJ. A.Head-GordonM. (1989). A Fifth-Order Perturbation Comparison of Electron Correlation Theories. Chem. Phys. Lett.157, 479483. 10.1016/s0009-2614(89)87395-6

  • 37

    Sánchez-MongeÁ.SchilkeP.SchmiedekeA.GinsburgA.CesaroniR.LisD. C.et al (2017). The Physical and Chemical Structure of Sagittarius B2. II. Continuum Millimeter Emission of Sgr B2(M) and Sgr B2(N) with ALMA. Astron. Astrophys.604, A6. 10.1051/0004-6361/201730426

  • 38

    SchilkeP.WalmsleyC. M.MillarT. J.HenkelC. (1991). Protonated HCN in Molecular Clouds. Astron. Astrophys.247, 487.

  • 39

    SchmiedekeA.SchilkeP.MöllerT.Sánchez-MongeÁ.BerginE.ComitoC.et al (2016). The Physical and Chemical Structure of Sagittarius B2. I. Three-Dimensional Thermal Dust and Free-Free Continuum Modeling on 100 au to 45 pc Scales. Astron. Astrophys.588, A143. 10.1051/0004-6361/201527311

  • 40

    SchwörerA.Sánchez-MongeÁ.SchilkeP.MöllerT.GinsburgA.MengF.et al (2019). The Physical and Chemical Structure of Sagittarius B2. IV. Converging Filaments in the High-Mass Cluster Forming Region Sgr B2(N). Astron. Astrophys.628, A6. 10.1051/0004-6361/201935200

  • 41

    SnyderL. E.BuhlD. (1972). Detection of Several New Interstellar Molecules. Ann. N. Y. Acad. Sci.194, 1724. 10.1111/j.1749-6632.1972.tb12687.x

  • 42

    SnyderL. E.BuhlD. (1971). Observations of Radio Emission from Interstellar Hydrogen Cyanide. Astrophys. J. Lett.163, L47. 10.1086/180664

  • 43

    StantonJ. F.GaussJ. (2000). Analytic Second Derivatives in High-Order many-body Perturbation and Coupled-Cluster Theories: Computational Considerations and Applications. Int. Rev. Phys. Chem.19, 6195. 10.1080/014423500229864

  • 44

    StantonJ. F.LopreoreC. L.GaussJ. (1998). The Equilibrium Structure and Fundamental Vibrational Frequencies of Dioxirane. J. Chem. Phys.108, 71907196. 10.1063/1.476136

  • 45

    ThaddeusP.GuélinM.LinkeR. A. (1981). Three New 'nonterrestrial' Molecules. Astrophys. J.246, L41L45. 10.1086/183549

  • 46

    ThorwirthS.SchreierP.SalomonT.SchlemmerS.AsvanyO. (2019a). Pure Rotational Spectrum of CN+. Astrophys. J. Lett.882, L6. 10.3847/2041-8213/ab3927

  • 47

    ThorwirthS.SchmidP. C.TöpferM.BrünkenS.AsvanyO.SchlemmerS. (2019b). “Spectroscopic Studies of Protonated Amines: CH3NH3+ and C2H5NH3+,” in International Symposium on Molecular Spectroscopy, 74th meeting, Champaign-Urbana, IL, U.S.A., June 17-21, Talk WF04.

  • 48

    TurnerB. E. (1974). U93.174 - A New Interstellar Line with Quadrupole Hyperfine Splitting. Astrophys. J. Lett.193, L83. 10.1086/181638

  • 49

    WakelamV.LoisonJ. C.HerbstE.PavoneB.BergeatA.BeroffK.et al (2015). The 2014 KIDA Network for Interstellar Chemistry. Astrophys. J. Suppl. Ser.217, 20. 10.1088/0067-0049/217/2/20

  • 50

    WangS.-X.SchroderusJ.OzierI.Moazzen-AhmadiN.McKellarA. R. W.IlyushynV. V.et al (2001). Infrared and Millimeter-Wave Study of the Four Lowest Torsional States of CH3CF3. J. Mol. Spectrosc.205, 146163. 10.1006/jmsp.2000.8235

  • 51

    WattsJ. D.GaussJ.BartlettR. J. (1992). Open-Shell Analytical Energy Gradients for Triple Excitation many-body, Coupled-Cluster Methods: MBPT(4), CCSD+T(CCSD), CCSD(T), and QCISD(T). Chem. Phys. Lett.200, 17. 10.1016/0009-2614(92)87036-o

  • 52

    WesternC. M. (2017). Pgopher: A Program for Simulating Rotational, Vibrational and Electronic Spectra. J. Quant. Spectrosc. Radiat. Transf.186, 221242. 10.1016/j.jqsrt.2016.04.010

  • 53

    ZerokaD.JensenJ. O. (1998). Infrared Spectra of Some Isotopomers of Methylamine and the Methylammonium Ion: A Theoretical Study. J. Mol. Struct. THEOCHEM425, 181192. 10.1016/s0166-1280(97)00109-7

  • 54

    ZiurysL. M.TurnerB. E. (1986). HCNH(+) - A New Interstellar Molecular Ion. Astrophys. J.302, L31L36. 10.1086/184631

  • 55

    ZuckermanB.MorrisM.PalmerP.TurnerB. E. (1972). Observations of Cs, HCN, U89.2, and U90.7 in NGC 2264. Astrophys. J. Lett.173, L125. 10.1086/180931

Summary

Keywords

rotational spectroscopy, protonated methylamine, ion trap, astrochemistry, symmetric top molecule

Citation

Schmid PC, Thorwirth S, Endres CP, Töpfer M, Sánchez-Monge Á, Schwörer A, Schilke P, Schlemmer S and Asvany O (2022) Rotational Rest Frequencies and First Astronomical Search of Protonated Methylamine. Front. Astron. Space Sci. 8:805162. doi: 10.3389/fspas.2021.805162

Received

29 October 2021

Accepted

22 December 2021

Published

02 February 2022

Volume

8 - 2021

Edited by

André Canosa, UMR6251 Institut de Physique de Rennes (IPR), France

Reviewed by

Roman Motiyenko, Université de Lille, France

Anthony Remijan, National Radio Astronomy Observatory, United States

Updates

Copyright

*Correspondence: Oskar Asvany,

† Present Address: Matthias Töpfer, Vinnolit GmbH & Co. KG, Hürth, Germany

This article was submitted to Astrochemistry, a section of the journal Frontiers in Astronomy and Space Sciences

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics