Experimental and Computational Studies on the Reactivity of Methanimine Radical Cation (H2CNH+•) and its Isomer Aminomethylene (HCNH2+•) With C2H2

1 Introduction

The reaction of ions is thought to play a crucial role in the synthesis of complex organic species in both the interstellar medium (ISM) and the atmospheres of planets and their satellites Larsson et al. (2012). As increasingly complex species are identified in such environments, the potential impact of isomers increases as well. The existence of isomeric ions in the ISM has been known for a long time and in some cases, such as HCO⁺ and HOC⁺, a determination has been made of their relative abundances Woods et al. (1983). As isomers have different spectroscopic and chemical properties, and, in many cases, isomerization barriers are too high to be overcome in interstellar conditions, they must be treated separately.

Notably, the relative abundance of isomers in a given environment cannot usually be determined by their thermodynamic properties, with the feasibility of formation and destruction mechanisms often more important than relative energies. Furthermore, while isomers often have quite different dipole moments and other spectroscopic properties and therefore have distinct radioastronomic profiles Woods et al. (1983), in situ measurements of ion abundances in planetary and satellite atmospheres often employ mass spectrometers which bring extreme sensitivity and time/localisation specific information, but cannot distinguish between isomers. The novel Atacama Large Interferometer Array (ALMA) has enabled data collection on the density as a function of altitude for HCN and HNC Cordiner et al. (2019) and the unprecedented resolution of this device should also allow for the retrieval of data about the abundance and distribution of isomeric ions in different astronomic objects. However, at present, the relative abundance of isomers in such environments has proven difficult to assess. So, in order to develop accurate models to predict these abundances, it is important to have reliable data on the reactivity of ionic isomers with common interstellar and atmospheric molecules.

In recent years, efforts have been made to selectively generate relevant isomers Fathi et al. (2016), with VUV photoionization in particular proving an effective tool Polášek et al. (2016) while noble gas tagging has been successful in characterizing the presence and the relative abundance of ions produced by electron impact ionization Brünken et al. (2019).

Like Earth, Titan has a dense nitrogen-dominated atmosphere Hörst (2017), which is one of the most chemically complex atmospheres in the solar system, as demonstrated by the data from the ion and neutral mass-spectrometer (INMS) of the Cassini mission. Titan’s atmospheric chemistry is driven by the ionization and/or dissociation of its two primary components, N₂ and CH₄, through a combination of extreme ultraviolet (EUV) radiation Imanaka and Smith (2007) and magnetospheric electrons. The resultant ions and radicals can then react to form more complex species Vigren et al. (2012), Westlake et al. (2014), including large hydrocarbons and nitrogen-bearing compounds Vinatier et al. (2007), Nixon et al. (2018), with cations of masses up to m/z 99 having been identified at altitudes of 950 km above the surface while cations of masses up to m/z 350 were detected by the Cassini plasma spectrometer ion beam sensor (CAPS-IBS) Vuitton et al. (2007).

Reactions of nitrogen-containing ions with hydrocarbons, especially unsaturated hydrocarbons such as C₂H₂ and C₂H₄, have been suggested as the basis for the formation of some of the complex ions that have been detected in Titan’s atmosphere Westlake et al. (2014), Vuitton et al. (2019). However, protonated nitriles, a major part of the nitrogen-containing ions in Titan’s atmosphere, are less reactive with their major reaction pathway in Titan’s atmosphere being dissociative recombination Vuitton et al. (2007). However, this is not the case for the methylenimine cation H₂CNH^+• and its isomer the aminomethylene ion HCNH₂^+•, which are both reactive radical cations and could form the basis of chain elongation reactions with unsaturated and saturated hydrocarbons that produce larger ions, from which neutral species can then be generated by dissociative recombination. Though the m/z 29 signal recorded by INMS is regarded as mostly arising from the presence of C₂H₅⁺, models predict a density of HCNH₂^+• and its isomers amounting to 1.1 × 10⁻² cm⁻³ in Titan’s ionosphere Vuitton et al. (2007).

Although the H₂CNH^+• and HCNH₂^+• ions have not yet been identified in the ISM, neutral methanimine (i.e. methylene imine, H₂CNH) is ubiquitous, having been detected in several objects including giant molecular clouds as well as both high mass and solar-type protostellar systems Dickens et al. (1997), Suzuki et al. (2016), Widicus Weaver et al. (2017), Ligterink et al. (2018), Bogelund et al. (2019). Methanimine, which could form either on the surface of dust grains Bernstein et al. (1995), Theule et al. (2011) or via gas phase reactions Suzuki et al. (2016), could serve as the basis for synthetic routes to form more complex nitrogen-containing molecules, including biomolecule precursors such as glycine and its corresponding α-aminonitrile Basiuk and Bogillo (2002), Aponte et al. (2017).

Due to its low ionization energy of 9.99 eV Holzmeier et al. (2013), methanimine can be ionized by VUV photons in the ISM. A possible pathway for the synthesis of larger nitrogen-containing species could therefore be the ionization of H₂CNH followed by the reaction of the resulting radical cation with an unsaturated hydrocarbon molecule, such as ethene or ethyne. The resultant ion could then either undergo dissociative recombination to form a complex neutral molecule through elimination of a hydrogen atom, or react further with other interstellar compounds to generate an even larger ionic species. Since methanimine has been detected also in ion-rich regions [e.g., at the “radical-ion peak” in the Orion molecular cloud Dickens et al. (1997)], such a course of reaction does not seem unlikely.

This paper presents a reactivity study of both H₂CNH^+• and HCNH₂^+• with ethyne using synchrotron radiation and apt neutral precursors to selectively generate the charged species. It aims to establish the occurrence of product channels leading to the formation of new C-C or C-N bonds through chain elongation reactions with ethyne, thereby leading to an increase in chemical complexity.

In Section 2 and Section 3 we present the experimental and theoretical methodologies. In Section 4 we provide a detailed description of the experimental results while in Section 5 the computational results are presented alongside proposed reaction mechanisms. Experimental results are discussed and interpreted in light of calculations of the most probable reaction mechanisms in Section 6. Finally, the conclusions are summarized in Section 7, where implications of present results on the formation of complex N-containing ions in Titan’s atmosphere and in the ISM are briefly addressed.

2 Experimental Methodology

Experiments were performed using the CERISES apparatus, an associated experiment of the DESIRS beamline Nahon et al. (2012) of the SOLEIL synchrotron radiation facility. Since the experimental set-up has been described in detail previously Alcaraz et al. (2004), Cunha de Miranda et al. (2015), only the most relevant details will be given here. CERISES is a guided ion beam tandem mass spectrometer composed of two octopoles located between two quadrupole mass filters in a Q1-O1-O2-Q2 configuration that permits investigation of bimolecular reactions of mass-selected ions occurring in a 4 cm long reaction cell located after the first octopole. Absolute reaction cross sections (CSs) and branching ratios (BRs) as a function of both photon (E_phot) and collision (E_CM) energies are derived by measuring the yields of parent- and product-ions.

The HCNH₂^+• and H₂CNH^+• isomers were produced by dissociative photoionization of the gaseous precursors cyclopropylamine (c-C₃H₅NH₂) and azetidine (c-CH₂CH₂CH₂NH) respectively, according to the following reactions:

$\text{c}-{\text{C}}_3{\text{H}}_5{\text{NH}}_2+h\nu \to {{\text{HCNH}}_2}^{+•}+{\text{C}}_2{\text{H}}_4(1)$

$\text{c}-{\text{CH}}_2{\text{CH}}_2{\text{CH}}_2\text{NH}+h\nu \to {\text{H}}_2{\text{CNH}}^{+•}+{\text{C}}_2{\text{H}}_4(2)$

These commercially available gaseous precursors were introduced into the ion source at roughly 10⁻⁶ mbar [see Richardson et al. (2021), Sundelin et al. (2021) for a complete description of the generation of the ions]. In addition, photoionization of methanimine (H₂CNH) was used to selectively generate the H₂CNH^+• isomer, with the details of the ion generation processes presented in Section 4.

The vacuum ultraviolet (VUV) radiation for photoionization is provided by the undulator-based DESIRS beamline Nahon et al. (2012), operating between 5 and 40 eV. Photons at the desired energies are selected and scanned simultaneously with the undulator peak energy by a normal incidence monochromator equipped with a low dispersion uncoated SiC grating (200 grooves/mm) optimized to provide photon flux up to 10¹³ photons/s. In the present experiments, the photon energies required to produce the [CH₃N]^+• isomer ions are in the region 9.5–14.0 eV. The monochromator slits were set in the range 300–600 μm, corresponding to a resolution of 20–40 meV, in order to balance the maximization of signal with avoiding detector saturation. The monochromator was operated at its first order of diffraction. Higher orders were removed from the incident beam using a gas filter Mercier et al. (2000) installed on the beam line and filled with Argon (at a pressure of 0.2 mbar) to efficiently remove any photons with E_phot larger than 15.7 eV. The absolute scale of the photon energy was checked using sharp absorption lines of atomic Argon around 11.828 and 14.304 eV Minnhagen (1973), Kramida et al. (2020) that were observed with systematic shifts of about 10–20 meV above their tabulated values.

The resolution of the mass filters Q1 and Q2 were adjusted in order to provide a suitable compromise between signal and separation of adjacent masses. We estimate that in this way we are able to remove 99.99% of neighboring ions for Q1 and 99.90% of neighboring ions for Q2 at the lowest collision energy, with a loss of performance reaching 98.80% of neighboring ions at the highest acceleration potentials of 20 V.

The collision energy in the lab is dependent on the reagent ion charge (+1 in our case) and on the difference between the ion source and the reaction cell potentials. The retarding potential method was applied to the reagent ion beam Teloy and Gerlich (1974) to determine the maximum of the first derivative of the reagent ion yield, which defines the zero of the kinetic energy. In this way we have estimated an average reagent ion beam FWHM of ∼0.4 eV, equivalent to ∼0.19 eV in the center-of-mass frame. By changing the potentials of the reaction cell and all subsequent elements, we were able to scan a collision energy range from ∼0.06 to ∼10 eV in the center of mass frame (E_CM).

C₂H₂ was stored in a dilute mixture with acetone, which was removed using a coil cooled with dry ice before insertion into the instrument. C₂H₂ was introduced in the reaction cell at a dynamic pressure of up to 2.0 × 10⁻⁷ bar, which guarantees operation close to a single collision regime, keeps secondary reactions at a reduced level and limits the attenuation of the reagent ion to less than 10%, while still providing an accurate pressure measurement and a dynamic range on the absolute cross sections determination in the range 0.01–100 Å². The absolute pressure in the reaction cell was measured using an MKS 398H differential manometer instrument. All data were collected with either reaction gas in the cell or in the surrounding chamber, therefore allowing the removal of any contribution resulting from collisions occurring outside the reaction cell, corresponding to no more than 5% of product ion counts. The absolute cross section acquisitions were performed in the so-called multi-scan mode, which involves recording signals for all ionic species at a single point before moving to the next point. This method has the benefit of drastically reducing the effects of potential drifts of source pressure, reaction cell pressure and photon flux.

3 Theoretical Methodology

The mechanisms for the reactions of H₂CNH^+• and HCNH₂^+• were studied using GAUSSIAN 09, Revision D.01 Frish et al. (2013). Geometries for intermediate structure, i.e. minima and transition states (TSs), were calculated at the MP2/6-31G(d) and MP2/6-311++G(d,p) levels. The identity of TS and minima were checked by frequency calculations and zero-point energy corrections were applied to the obtained energies. IRC calculations were performed at the MP2/6-31G(d) level to ensure that the TS connect the correct minima. Single point energy calculations were carried out for all stationary points at the CCSD(T)/6-311++G(d,p) level with zero-point energy corrections taken from the MP2/6-311++G(d,p) level. Reactants and product geometries and energies were also calculated at the same level (including zero-point energy corrections) and the reactants’ energies subtracted from the optimizations of the TS and minima, resulting in the relative energies E_rel. Cartesian coordinates, structures and electronic energies of reagents, products, minima and TS are reported in the Supplementary Material.

4 Experimental Results

The characterization of the fragmentation pathways for generation of the two isomers via dissociative photoionization of cyclopropylamine (c-C₃H₅NH₂) and azetidine (c-CH₂CH₂CH₂NH) has been described in previous publications Richardson et al. (2021), Sundelin et al. (2021) to which the interested reader is referred to for details. In the following subsection we will focus on the formation of the H₂CNH^+• isomer through direct ionization of methanimine.

4.1 Generation of H₂CNH^+• From Photoionization of Methanimine

H₂CNH was synthesized by dehydrocyanation of aminoacetonitrile performed under vacuum, following an improved version of the original method Guillemin and Denis (1988) and the interested reader is referred to Gans et al. (2019) for details of the chemical procedure. After synthesis, methanimine was stored under liquid nitrogen, and transferred to dry ice for introduction in the ion source of the CERISES set-up using its vapor pressure at dry ice temperature (T ∼ −78°C).

The photoionization efficiency (PIE) curves (i.e. ion yields as a function of photon energy) for the parent radical cation H₂CNH^+• (m/z 29) and the photo-dissociation product HCNH⁺ (m/z 28) are reported in Figure 1, where the observed Appearance Energies (AEs) are given. For H₂CNH^+•, an Ionization Energy (IE) of 10.01 ± 0.08 eV was obtained using the method of linear threshold extrapolation of the PIE curve at m/z 29 Chupka (1971), Traeger and McLoughlin (1981), Castrovilli et al. (2014), a value that is fully consistent with the adiabatic ionization energy of methanimine of 9.99 eV as determined via mass-selected threshold photoelectron spectroscopy Holzmeier et al. (2013). For the HCNH⁺ photo-fragment, an experimental AE = 11.02 ± 0.08 eV was estimated taking the first point above noise level of the PIE curve at m/z 28, after signal correction due to “spillover” of the more abundant peak at m/z 29. This is also consistent with the onset of dissociative photoionization of methanimine via H^• loss, observed at 11.13 eV Holzmeier et al. (2013).

FIGURE 1

FIGURE 1. PIE curves for the photoionization of methanimine in the 9.8–14 eV photon energy region: open black circles are for H₂CNH^+• radical cations (m/z 29) and open red triangles for HCNH⁺ dissociation products (m/z 28). Insets (A) and (B) show enlarged plots of the threshold regions for the corresponding m/z values. The dashed line indicates the linear extrapolation used to extract the experimental IE.

We emphasize that the alternative generation method starting from neutral methanimine is explored as a way to produce the methanimine radical cation isomer free from interference from isobaric contaminants (C₂H₅⁺, H¹³CNH⁺ and ¹³CCH₄^+•), as is detailed in Section 4.3. Conversely, since the methanimine sample is free from heavier mass contaminants, the only expected interference at m/z 29 is limited to the ¹³C contribution of the HCNH⁺ photo-fragment at m/z 28. However, from Figure 1 it can be clearly seen that this fragment is less intense than the m/z 29 cation at least up to a photon energy of ∼13.1 eV. As a consequence, any contamination due to H¹³CNH⁺ at m/z 29 should account for less than 1% of the m/z 29 signal below E_phot = 13.1 eV.

4.2 Reactions of H₂CNH^+• and HCNH₂^+• With C₂H₂

The reaction of both isomers yields products at m/z 54, 39 and 28, while for the H₂CNH^+• isomer an additional minor product was observed at m/z 27. The most probable product channels are indicated as Eqs 3–11, and the reaction enthalpies estimated from literature values are summarized and compared with results from our theoretical calculations in Table 1.

$\begin{array}{ll}\hfill {\text{H}}_2{\text{CNH}}^{+•}/{\text{HCNH}}_2^{+•}+{\text{C}}_2{\text{H}}_2& \to {\text{C}}_2{\text{H}}_3^{+}(m/z27)+{\text{H}}_2{\text{CN}}^{•}\hfill \end{array}(3)$

$\begin{array}{ll}\hfill & \to {\text{C}}_2{\text{H}}_3^{+}(m/z27)+{\text{HCNH}}^{•}\hfill \end{array}(4)$

$\begin{array}{ll}\hfill & \to {\text{HCNH}}^{+}(m/z28)+{\text{C}}_2{\text{H}}_3^{•}\hfill \end{array}(5)$

$\begin{array}{ll}\hfill & \to {\text{C}}_2{\text{H}}_4^{+•}(m/z28)+\text{HCN}\hfill \end{array}(6)$

$\begin{array}{ll}\hfill & \to \text{c}-{\text{C}}_3{\text{H}}_3^{+}(m/z39)+{\text{NH}}_2^{•}\hfill \end{array}(7)$

$\begin{array}{ll}\hfill & \to {\text{H}}_2{\text{CCHCNH}}^{+}(m/z54)+{\text{H}}^{•}\hfill \end{array}(8)$

$\begin{array}{ll}\hfill & \to \text{c}-\text{CHCHC}{({\text{NH}}_2)}^{+}(m/z54)+{\text{H}}^{•}\hfill \end{array}(9)$

$\begin{array}{ll}\hfill & \to {{\text{HCCCHNH}}_2}^{+}(m/z54)+{\text{H}}^{•}\hfill \end{array}(10)$

$\begin{array}{ll}\hfill & \to {\text{H}}_2{\text{CCHNCH}}^{+}(m/z54)+{\text{H}}^{•}\hfill \end{array}(11)$

$\begin{array}{ll}\hfill & \to {\text{H}}_2{{\text{CCCNH}}_2}^{+}(m/z54)+{\text{H}}^{•}\hfill \end{array}(12)$

$\begin{array}{ll}\hfill & \to {\text{CH}}_2{{\text{CNCH}}_2}^{+}(m/z54)+{\text{H}}^{•}\hfill \end{array}(13)$

TABLE 1

TABLE 1. Reaction enthalpies for the reaction: HCNH₂^+•/H₂CNH^+• + C₂H₂ → products.

Although the m/z 39 product could alternatively correspond to the formation of HCCN^+• or HCNC^+• ions formed together with a CH₄ molecule, these processes are not only endothermic by about 2 eV [on the basis of the calculated energies of HCCN^+• and HCNC^+• Harland and McIntosh (1985), Goldberg et al. (1995), Holmes et al. (2006)] but also very difficult to account for from a mechanistic point of view. Hence this possibility is not considered any further and the corresponding products are not included in Table 1.

4.3 Considerations About Impurities in the Reagent Ions

When selecting m/z 29 from the neutral precursors (cyclopropylamine and azetidine), there is no way to discriminate against the presence of isobaric contaminants, namely C₂H₅⁺, H¹³CNH⁺ and ¹³CCH₄^+•. Calculations and experimental evidence from our previous works Richardson et al. (2021), Sundelin et al. (2021) show that dissociative photoionization of azetidine to give H₂CNH^+• plus C₂H₄ has an experimental AE (10.2 ± 0.1 eV) very close to the fragmentation channel giving HCNH⁺ plus C₂H₅^+•. Additionally, the azetidine radical cation can fragment into C₂H₄^+• plus H₂CNH at photon energies above ∼11.3 eV. Hence both H¹³CNH⁺ and ¹³CCH₄^+• might be present as impurities when mass selecting reagent ions at m/z 29 from azetidine. For cyclopropylamine, dissociative photoionization to HCNH₂^+• (plus C₂H₄) has a similar AE (10.2 ± 0.1 eV) also very close to fragmentation into HCNH⁺ plus C₂H₅^•. Additionally, the formation of the isobaric cation C₂H₅⁺ is possible at photon energies higher than ∼12.4 eV. Since fragmentation into C₂H₄^+• plus HCNH₂ opens up at very high photon energies, any contamination due to ¹³CCH₄^+• is not relevant in this case. Hence, for cyclopropylamine the relevant contaminants at m/z 29 are H¹³CNH⁺ and C₂H₅⁺.

In the following the reactivity of such isobaric impurities with C₂H₂ will be addressed.

4.3.1 Contamination due to C₂H₅⁺

Literature data Kim et al. (1977), Anicich (1993), Anicich (2003) indicate that C₂H₅⁺ reacts with C₂H₂ with an overall rate k = 1.9 × 10⁻¹⁰ cm³⋅molec⁻¹⋅s⁻¹ to give the following products (with BRs indicated in parenthesis):

$\begin{array}{ll}\hfill {\text{C}}_2{\text{H}}_5^{+}+{\text{C}}_2{\text{H}}_2& \to {\text{C}}_3{\text{H}}_3^{+}(m/z39)+{\text{CH}}_4(\text{BR}=0.36)\hfill \end{array}(14)$

$\begin{array}{ll}\hfill & \to {\text{C}}_4{\text{H}}_5^{+}(m/z53)+{\text{H}}_2(\text{BR}=0.64)\hfill \end{array}(15)$

A recent re-evaluation Richardson et al. (n.d.) of this reaction using the same set-up employed in this work where C₂H₅⁺ was generated via VUV dissociative photoionization of C₂H₅Br at a photon energy of 11.3 eV and at the lowest achievable collision energy (E_CM ∼ 0.08 eV) gives quite different product branching ratios. Most notably, Eq. 14 is found to be the major channel (BR = 0.71), while Eq. 15 is the second most abundant channel (BR = 0.22), with minor amounts of the other products (C₄H₇⁺ and C₂H₃⁺) accounting for the remaining product flux.

Since Eq. 14 from C₂H₅⁺ leads to a product of the same mass as Eq. 7, contamination at m/z 39 will emerge as the photon energy is increased when working with HCNH₂^+• from cyclopropylamine.

4.3.2 Contamination due to H¹³CNH⁺ and ¹³CCH₄^+•

As already mentioned, a contribution from H¹³CNH⁺ at all photon energies is expected for both azetidine and cyclopropylamine experiments. However, literature reports indicate that HCNH⁺ is unreactive with C₂H₂ [k < 1.0 × 10⁻¹³ cm³⋅molec⁻¹⋅s⁻¹ at room temperature Anicich et al. (2000), Anicich and McEwan (2001), Milligan et al. (2001)], as expected from thermochemical considerations, since the most feasible channel (i.e. proton transfer to give C₂H₃⁺ plus HCN) is endothermic by approximately 0.75 eV. In our experiment, this channel could potentially open but only at high photon energies and should be accounted for when looking at data collected at m/z 27 under such conditions. Contamination due to ¹³CCH₄^+• is expected only when working with azetidine above ∼11 eV photon energy. Literature data Anicich (1993) indicate that the ethylene radical cation C₂H₄^•+ reacts with C₂H₂ with an overall rate k = 8.4 × 10^–10 cm³⋅molec⁻¹⋅s⁻¹ to give the following products:

$\begin{array}{ll}\hfill {\text{C}}_2{\text{H}}_4^{+•}+{\text{C}}_2{\text{H}}_2& \to {\text{C}}_3{\text{H}}_3^{+}(m/z39)+{\text{CH}}_3^{•}(\text{BR}=0.77)\hfill \end{array}(16)$

$\begin{array}{ll}\hfill & \to {\text{C}}_4{{\text{H}}_5}^{+}(m/z53)+{\text{H}}^{•}(\text{BR}=0.23)\hfill \end{array}(17)$

The reaction of ¹³CCH4^+• should therefore lead to products at m/z 39 and 40 (i.e. C₃H₃⁺ and ¹³CC₂H₃⁺ from Eq. 16) plus a minor amount of m/z 54 (¹³CC₃H₅⁺) via Eq. 17.

To further check for the reactivity of ¹³CCH₄^+• contaminant, control experiments have been performed by selecting m/z 28 ions (a combination of HCNH⁺ and C₂H₄^•+) from dissociative ionization of azetidine and reacting them with C₂H₂. Products observed at m/z 39 and 53 are indicative of the occurrence of Eqs 16, 17 and the BR is consistent with the literature value reported above. Hence contamination due to ¹³CCH₄^+• might overlap with products at m/z 39 and 54 in the reaction of H₂CNH^+• generated from azetidine. As already discussed in Sundelin et al. (2021), since the ratio of the photodissociation yields for m/z 28 and 29 at the explored photon energies is in the range 0.3–1.0, only a small contamination of ¹³CCH₄^+• should be present in the reagent beam.

4.4 Experimental Results: Data as a Function of the Photon Energy (E_phot)

Product cross sections (CSs) as a function of E_phot for the reactions of HCNH₂^+• and H₂CNH^+• isomers generated through dissociative photoionization of cyclopropylamine and azetidine are given in Figure 2, while in Figure 3 analogous data are presented for H₂CNH^+• generated via photoionization of methanimine. In both figures the large uncertainties at the lowest photon energies are due to the low parent ion intensity near threshold Richardson et al. (2021). Trends as a function of E_phot give qualitative indications on how the reactivity is affected by changes in the internal energy content of the reagent ions.

FIGURE 2

FIGURE 2. Reactive cross sections as a function of E_phot for the reaction of HCNH₂^+• (top) and H₂CNH^+• (bottom) with C₂H₂. The experimental AEs of parent ions at m/z 29 from both precursors are at E_phot = 10.2 ± 0.1 eV (outside the x axis range). Collision energies are fixed at E_CM = 0.16 eV (top) and 0.12 eV (bottom). The vertical dashed lines indicate the photon energies at which data have been collected as a function of E_CM (see Figure 4).

FIGURE 3

FIGURE 3. Reactive cross sections as a function of E_phot for the reaction of H₂CNH^+• with C₂H₂. The collision energy is fixed at E_CM ∼ 0.1 eV. The x axis scale at the top is the difference between E_phot (bottom x scale) and the AE for production of the parent ion from photoionization of methanimine, and equal to 10.01 ± 0.08 eV, as detailed in the text. The vertical dashed lines indicate the photon energies at which data have been collected as a function of E_CM (see Figure 5).

The reaction of HCNH₂^+• (top part of Figure 2) gives a major product at m/z 54 with no significant dependence on E_phot and, by extension, on the internal energy of the ion, at least up to ∼12.5 eV. The decrease in CS at E_phot > 12.5 eV is an artefact due to the increasing contamination of C₂H₅⁺ in the parent beam (as discussed in Section 4.3), which reduces the relative amount of HCNH₂^+• ion available for reaction. The other products at m/z 28 and 39 are minor, but they also do not show any significant change as a function of the photon energy, with the rise in the m/z 39 product above E_phot > 12.5 being due to the emergence of the C₂H₅⁺ contaminant.

The reaction with H₂CNH^+• from azetidine (bottom part of Figure 2) also gives a major product at m/z 54, but differently from the previous isomer, its CS shows a marked decrease with increasing E_phot before plateauing, thus pointing to a barrierless pathway somehow inhibited by the increase in internal energy of the reactant ion. The other main channel is at m/z 28, with the CS showing no change with increasing E_phot, making it the major product at E_phot above ∼11.7 eV. The m/z 39 channel shows a very small cross section up to E_phot = 11.4–11.5 eV, though above this energy a slight increase is observed which might be the result of the minor ¹³CCH₄^+• impurity (see discussion in Section 4.3). Being very minor, this channel is discussed here only due to its relevance at higher collision energies. The reaction with H₂CNH^+• shows an additional minor product at m/z 27, whose CS increases with increasing photon energy.

CSs for the reaction of H₂CNH^+• from methanimine (Figure 3) show the same four products, with channels at m/z 54 and 28 being the major ones at all photon energies, and those at m/z 39 and 27 being minor throughout the whole photon energy range, in good agreement with the data collected using azetidine. As already mentioned, products observed using methanimine as precursor do not suffer from the contamination issues mentioned in Section 4.3, however a quantitative comparison as a function of E_phot should be performed with care, since the photon energy scales are not directly related. In fact, while in the case of methanimine photoionization the excess energy carried by the photon will contribute to increase the internal energy of the reacting H₂CNH^+•, in the case of dissociative photoionization from azetidine the excess energy might be dissipated also as internal and kinetic energy of the neutral co-fragment. Therefore, the change in internal energy of H₂CNH^+• with increasing E_phot is expected to be smaller in the azetidine experiment than for the methanimine precursor. CS for product at m/z 54 decreases sharply with increasing E_phot up to about 10.8 eV, after which it levels off and is approximately constant. Possible explanations of the observed trend are discussed in Section 6, here we only highlight the good qualitative agreement with the decrease as a function of E_phot observed with the dissociative photoionization of azetidine.

Similarly to the azetidine experiment, the CS for the m/z 28 product shows no notable change as a function of E_phot. The CS for the very minor product at m/z 39 also shows no significant dependence on E_phot. This is at odds with the azetidine data in which a very slight increase was observed at higher photon energies, further supporting the idea that this increase is due to the emergence of a pathway with the ¹³CCH₄^+• impurity. Finally, CS for the other minor channel at m/z 27 shows a similar rise with increasing E_phot as observed with the azetidine precursor.

The BRs for the various product channels, obtained from the absolute cross section measurements, are summarized in Table 2 for both isomers. In the same Table estimates of the total (i.e. summed over all the product channels) rate constants at a fixed average energy k_tot (E_ave) have been given. For details on how cross sections have been converted into rate constants see our previous paper Sundelin et al. (2021). For the HCNH₂^+• isomer, BRs and k_tot values were obtained via dissociative photoionization of cyclopropylamine at E_phot = 12.0 eV, while for the H₂CNH^+• the methanimine data at E_phot = 11.1 eV were used. For both isomers the collision energy was fixed at E_CM = 0.08 ± 0.01 eV, corresponding to an average energy E_ave = 0.10 ± 0.01 eV Ervin and Armentrout (1985), Nicolas et al. (2002). The ratio between k_tot and the Langevin collision rate constant (k_L), also reported in Table 2, is a useful estimate of the overall efficiency of the reactions, which proceed at a rate of ∼10% of the Langevin limit for both isomers.

TABLE 2

TABLE 2. Total rate constants at fixed average energy k_tot (E_ave) and branching ratios (BRs) for the reaction of HCNH₂^+• and H₂CNH^+• with C₂H₂. Results have been obtained at E_CM = 0.08 ± 0.01 eV, corresponding to an average energy E_ave = 0.10 ± 0.01. The HCNH₂^+• isomer was produced via dissociative photoionization of cyclopropylamine at E_phot = 12.0 eV while the H₂CNH^+• isomer was produced by photoionization of methanimine at E_phot = 11.1 eV.

4.5 Experimental Results: Data as a Function of the Collision Energy (E_CM)

The CSs as a function of E_CM for the reaction of both HCNH₂^+• and H₂CNH^+• isomers generated through dissociative photoionization of cyclopropylamine and azetidine are shown in Figure 4, while in Figure 5 analogous data are presented for H₂CNH^+• generated via photoionization of methanimine.

FIGURE 4

FIGURE 4. Reactive cross sections as a function of the collision energy (E_CM) for the reaction of HCNH₂^+• (top) and H₂CNH^+• (bottom) with C₂H₂. The photon energy is fixed at E_phot = 12.0 eV (top) and 11.5 eV (bottom). The vertical dashed lines indicate the collision energy at which data have been collected as a function of E_phot (see Figure 2). In the top panel: the grey solid line is the Langevin CS (rescaled by a factor 0.1 to fit in the figure). The magenta point is the total reactive CS (σ_tot) used to calculate the total rate constant k_tot reported in Table 2.

FIGURE 5

FIGURE 5. Reactive cross sections as a function of E_CM for the reaction of H₂CNH^+• (from methanimine) with C₂H₂. The photon energy is fixed at E_phot = 11.1 eV. The vertical dashed line is at the collision energy at which data have been collected as a function of E_phot (see Figure 3). The grey solid line is the Langevin CSs (rescaled by a factor 0.1 to fit in the figure). The magenta point is the total reactive CS (σ_tot) used to calculate the total rate constant k_tot reported in Table 2.

For the HCNH₂^+• isomer, the m/z 54 product shows a marked decrease with increasing E_CM, indicative of a barrierless pathway most likely proceeding via a complex-mediated mechanism. The m/z 28 and 39 channels show a combination of two trends, with a similar decrease below 1 eV followed by a gradual increase with increasing E_CM, indicating the presence of a barrierless pathway, mostly contributing at low collision energy, superimposed on an endothermic process (or a process involving a barrier) appearing at high E_CM. It should be noted that CS at m/z 28 are affected by collisional fragmentation of the parent which is discussed in more detail below.

For the H₂CNH^+• isomer with the azetidine precursor (bottom of Figure 4), CSs were measured only for products at m/z 54 and 39, with the former decreasing with increasing E_CM and the latter showing a decrease with increasing E_CM, at low collision energies, followed by a ramp at very high collision energies. A more complete picture for the reactivity of the H₂CNH^+• isomer is obtained from the experiments using methanimine (see Figure 5) since CS as a function of E_CM have been measured for all of the four product channels. The trend for m/z 54 product, decreasing with increasing E_CM, confirms the results from the azetidine experiment. Also similar to the azetidine experiment, the m/z 39 product is minor at all collision energies and shows an initial decrease with increasing E_CM followed by an inversion in the E_CM trend at high energies. The m/z 28 channel shows a slight initial decrease followed by a gradual increase with increasing E_CM that resembles the trend observed for the HCNH₂^+• isomer. The m/z 27 product shows a similar increase with increasing E_CM as it does with photon energy, indicating that the barrier to this pathway can be overcome equally effectively with collision energy as internal energy. A detailed discussion with explanations of the trends as a function of collision energy is given in Section 6 on the basis of the results obtained from computations and presented in the next section.

5 Computational Results

Our calculations show that HCNH₂^+• is slightly more stable (by 18.9 kJ/mol) than H₂CNH^+•, consistent with a previous literature value Nguyen et al. (1994) in which the difference in the heat of formations of the two isomers, calculated at the UQCISD(T)/6-311++G(3df, 2p) and RCCSD(T)/cc-pVTZ levels of theory, is equal to 16 kJ/mol. Given the open shell nature and the presence of a charge such differences in the calculations are not surprising. All enthalpies are calculated with the GAUSSIAN suite of programs at the CCSD(T)/6-311++G(d,p) level of theory. Unless explicitly noted otherwise, relative enthalpies are given (in kJ/mol and in parenthesis after each structure) with respect to the sum of the separated reactants having the lowest energy, i.e. HCNH₂^+• plus C₂H₂.

For both isomers, the reaction mechanisms proceed via the formation of covalently bound adducts through a radical-terminal approach, with the barriers for the non-radical terminal attack being endothermic by more than 45 kJ/mol in both cases.

HCNH₂^+•isomer (pathways are described graphically in Figures 6, 7, while reaction enthalpies are compared with literature values in Table 1): for this isomer, an initial van der Waals cluster A1 (−27.9) can proceed to a covalently-bound adduct A2 (−227.6) via the transition state (TS) A1-2 (−27.4). From A2, a series of pathways are possible, all of which are shown diagrammatically in Figure 7 as well as being described here.

FIGURE 6

FIGURE 6. Structures and energies of the van der Waals adducts (and their evolution into covalently bound intermediates) for the entrance channel of the reactions of HCNH₂^+• (panel A) and H₂CNH^+• (panel B) plus C₂H₂. Energies (in kJ/mol) are defined with respect to the separated HCNH₂^+• plus C₂H₂ reactants. Calculations are at the MP2/6-311++G(d,p)//CCSD(T)/6-311++G(d,p) level of theory.

FIGURE 7

FIGURE 7. Energies and reaction pathways for the HCNH₂^+• plus C₂H₂. Energies (in kJ/mol) are defined with respect to the separated HCNH₂^+• plus C₂H₂ reactants. Calculations are at the MP2/6-311++G (d,p)//CCSD(T)/6-311++G(d,p) level of theory. Molecular structures are shown only for intermediates and products. For TSs’ structures see the Supplementary Material.

The most favorable pathway, from the thermochemical point of view, to form the observed major product at m/z 54 proceeds from a rotamer of A2, named A3 (−217.1), with the two structures being separated by the TS A2-3 (−191.0). A3 can then undergo a [1,4] H-shift (from the N atom to the terminal C) to give A4 (−150.1) via the TS A3-4 (−71.7). From A4 the ejection of an H^• via A4-P (−49.7) gives the m/z 54 product in the form of the lowest energy isomer H₂CCHCNH⁺ (protonated vinyl cyanide) in combination with H^• (−107.5).

While producing a higher energy C₃H₄N⁺ isomer, the lowest energy pathway for m/z 54 proceeds through the cyclic intermediate A7 (−154.2) which can be reached via A2-7 (−147.6). A7 can then barrierlessly eject the H^• from the central carbon to give a three-membered ring aminic ion c-CHCHC(NH₂)⁺ (−73.6).

A different C₃H₄N⁺ isomer at m/z 54 can be formed from A3 by the direct ejection of H^• via the TS A3-P (−47.8) to give a linear aminic ion HCCCHNH₂⁺ plus H^• (−65.3). A fourth pathway leading to the formation of a product at m/z 54 proceeds from A2, which can undergo a [1,3] H-shift via the TS A2-9 (−43.1) to give A9 (−222.8), which can then eject H^• via TS A9-P (−32.8) to give the already mentioned H₂CCHCNH⁺ plus H^•. All these pathways are fully submerged, hence they can all contribute to the observed m/z 54 product.

Though not included in Figure 7, our calculations have identified another pathway, higher in energy with respect to the ones above, that leads to the production of a fourth high energy C₃H₄N⁺ isomer, the linear CH₂CCNH₂⁺ (protonated ketenimine). This pathway (shown in Supplementary Figure S30 of the Supplementary Material) starts from A9 and proceeds via successive interconversions into the rotamers A10 (−230.5), via TS A9-10 (−226.7), and then into A11 (−180.1), via TS A10-11 (−188.9), with the lower energy of the TS being due to the zero-point correction. A11 can undergo a hydrogen migration to give A12 (−226.2), via TS A11-12 (−60.5), which can finally eject H^• to give CH₂CCNH₂⁺ plus H^• (−23.6) via TS A12-P (−3.3).

The next major product from the HCNH₂^+• isomer is at m/z 28. A possible pathway has been identified stemming from A4. This intermediate can interconvert to give its rotamer A5 (−138.8) via the TS A4-5 (−157.6). Note that the lower energy of the TS compared to A4 is due to the zero point correction. A5 can then cleave the central C-C bond to give a van der Waals complex of the products, A6 (−49.5), via the TS A5-6 (−30.9). The van der Waals adduct can subsequently fragment, with no energy barrier, into the products HCNH⁺ (protonated hydrogen cyanide) plus C₂H₃^• (−9.0).

Finally, a pathway should be accounted for the extremely minor product detected at m/z 39 and attributed to c-C₃H₃⁺. It starts from the cyclic A7 intermediate, which can cleave the C-N bond to give the van der Waals complex A8 (−46.1) via the TS A7-8 (−37.7). A8 can then fragment barrierlessly into the products c-C₃H₃⁺ and NH₂^+• (+3.9). The small yield of the m/z 39 product is compatible with the calculated endothermicity of the channel.

H₂CNH^+• isomer (pathways are described graphically in Figures 6, 8, while reaction enthalpies are compared with literature values in Table 1): for this isomer the reagents (at +18.9 kJ/mol above the other isomer) come together to form an initial van der Waals complex B1 (−32.6) which can then evolve, by formation of a C-N bond, into the covalently-bound complex B2 (−146.4) via the TS B1-2 (−33.6), with the lower energy of the TS compared to the van der Waals complex arising from the zero-point energy correction. B2 can then interconvert into its rotamer B3 (−154.7) via TS B2-3 (−141.9) from which the various pathways then proceed, as shown diagrammatically in Figure 8 and described in the following.

FIGURE 8

FIGURE 8. Energies and reaction pathways for the H₂CNH^+• plus C₂H₂. Energies (in kJ/mol) are defined with respect to the separated HCNH₂^+• plus C₂H₂ reactants. Calculations are at the MP2/6-311++G(d,p)//CCSD(T)/6-311++G(d,p) level of theory. Molecular structures are shown only for intermediates and products. For TSs’ structures see the Supplementary Material.

The lowest energy pathway to form the m/z 54 product proceeds via a cyclisation to give a four-membered ring intermediate, B5 (−239.3), via B3-5 (−48.3). B5 can open to give the N-terminal structure A4, via B5-A4 (−108.6). The A4 intermediate is a common structure with the potential energy hypersurface of the other isomer, hence from there the reaction can proceed to give H₂CCHCNH⁺ (protonated vinyl cyanide) plus H^• (−107.5), as already described above. Other C₃H₄N⁺ isomers could be produced from B3 via an interconversion into a second rotamer B4 (−154.8) via the TS B3-4 (−151.4). The latter rotamer can undergo a [1,3] H-shift to give B9 (−240.2) via B4-9 (+32.8). B9 can then eject H^• to give H₂CCHNCH⁺ (protonated vinyl isonitrile) plus H^• (−65.2) via the TS B9-P (−24.4). While this pathway involves a barrier higher in energy than the reactants (i.e. B4-9 that is 13.9 kJ/mol above the reactants) it is potentially relevant at higher photon energies where the HCNH₂^+• isomer might be formed with sufficient internal energy to both overcome the barrier and to potentially inhibit the alternative low energy pathway going via the cyclisation to give B5.

There are also two higher-energy pathways leading to another C₃H₄N⁺ isomer at m/z 54 that proceed via B9. For clarity reasons they have not been included in Figure 8, but they are shown in Supplementary Figure S31 and are described briefly here. The first pathway involves a [1,2] H-shift to give B11 (−277.1) via a TS B9-11 (−35.5), with B11 then ejecting H^• to give CH₂CNCH₂⁺ (−24.4) via B11-P (+3.7). The second involves a [1,3] H-shift from B9 to give B12 (-276.1) via TS B9-12 (+10.2), which then also eject H^• to give CH₂CNCH₂⁺, this time via TS B12-P (+3.2).

The m/z 28 product, corresponding to HCNH⁺ formed in association with C₂H₃^•, can proceed from the A4 intermediate, as already described for the HCNH₂^+• isomer, but it can also form from B3 via a [1,4] H-shift to give B6 (−133.1) via the TS B3-6 (−31.7). B6 can then interconvert into its rotamer B7 (−138.7), via the TS B6-7 (−98.5), which can then cleave the central C-N bond to give B8 (−69.8), a van der Waals complex of the products, via B7-8 (6.4). B8 can then dissociate barrierlessly into HCNH⁺ and C₂H₃^•.

As A4 is a common structure, the whole of the HCNH₂^+• PES is also accessible for the H₂CNH^+• isomer, therefore creating a potential pathway to a m/z 39 product through a series of rearrangements via A2, as previously described for the HCNH₂^+• isomer. However, due to the large number of interconversions, any contribution from this pathway is expected to be extremely minor despite being fully submerged with respect to the reactants.

While not included in Figure 8, there is also a pathway to interconvert between B5 and the rotamers B6 and B7 and a diagram is shown in Supplementary Figure S32. B5 can undergo a [1,2] H-shift to give B10 (−231.4) via B5-10 (−71.7). This can then ring-open via C-C bond cleavage to give B6 via B10-6 (−48.3) or B7 via B10-7 (−24.3). For the very minor m/z 27 product, only observed at high photon and high collision energies, two endothermic pathways have been found, one proceeding via intermediate B6 without a barrier to give C₂H₃⁺ plus HCNH (+145.1) and the other one going from B9 to give C₂H₃⁺ plus H₂CN (+104.0), also with no further barrier. Since the energies of both products are out of scale they have not been reported in Figure 8.

6 Discussion

Comparison of the two sources of H₂CNH^+• shows a high level of qualitative agreement with regards to the trends as a function of both E_phot and E_CM, as well as broad quantitative agreement for both CSs and BRs of the various products. Exact quantitative comparison is not possible due to the different scaling of internal energy with photon energy. Similarly, comparison of the data collected with the two different isomers indicates a high level of selectivity for both isomers over the whole photon energy range.

For both isomers the major product is at m/z 54, corresponding to C₃H₄N⁺. Previous computational studies Wang et al. (2015), Heerma et al. (1986) have provided strong evidence for the existence of several C₃H₄N⁺ isomers in the gas phase, the most stable of which being protonated vinyl cyanide H₂CCHCNH⁺. The next most stable is a three-membered ring structure c-(CHCHC)–NH₂⁺ that is 22.2 kJ/mol higher in energy Wang et al. (2015), followed by the linear aminic structure HCCCHNH₂⁺ [37.7 kJ/mol above the most stable isomer Wang et al. (2015)]. At higher energies there is protonated vinyl isonitrile, H₂CCHNCH⁺ [49.8 kJ/mol above the most stable isomer Wang et al. (2015)], protonated ketenimine H₂CCCNH₂⁺ [57.4 kJ/mol above the most stable isomer Wang et al. (2015)] and CH₂CNCH₂⁺ [61.1 kJ/mol above the most stable isomer Wang et al. (2015)]. Although our experimental technique can not distinguish between the different C₃H₄N⁺ isomers, ab-initio calculations can help in differentiating among the pathways leading to them.

For the reaction of HCNH₂^+• the m/z 54 product is approximately constant in cross section as a function of the photon energy while it decreases with increasing collision energy. These findings can be rationalized by the fact that the m/z 54 formation pathways have only submerged barriers and proceed via complex-mediated mechanisms, as confirmed by our calculations. In fact, not just one, but four different pathways with fully submerged barriers have been found (see Section 5) to yield the three lowest energy isomers CH₂CHCNH⁺, c-CHCHC(NH₂)⁺ and HCCCHNH₂⁺. Hence, the observed m/z 54 product is expected to correspond to a combination of the three isomers.

For the H₂CNH^+• isomer, the CS trend with E_phot for the m/z 54 product can be explained through a combination of multiple channels. At low internal energies, a barrierless pathway proceeds via the cyclized adduct B5 but this decreases in relevance as the internal energy increases due to the inhibition of the cyclisation. At higher E_phot, a pathway with a slight barrier via B9 opens, with this pathway giving the plateau that is observed at higher internal energies. Hence, for H₂CNH^+•, in the absence of internal excitation, the observed m/z 54 product is exclusively due to protonated vinyl cyanide H₂CCHCNH⁺.

As a function of the collision energy, the m/z 54 channels for both isomers show an initial decrease in CS section, indicative of barrierless pathways as suggested by our calculations. This is due to the higher collision energy inhibiting the formation of the van der Waals clusters (A1 and B1) in the entrance channels.

For the reaction of HCNH₂^+•, the m/z 28 and 39 products are approximately constant in cross section as a function of E_phot. This is in agreement with the computations that identify fully submerged pathways for the m/z 28 product and for formation of the van der Waals complex of the products for the m/z 39 channel, with the barrierless fragmentation of this complex lying approximately equal in energy with the reactants. The lower cross section of the m/z 28 product for this isomer than for H₂CNH^+• is expected to be due to the presence of the multiple lower-lying channels to give products at m/z 54 which are discussed above.

For the reaction of H₂CNH^+•, the m/z 28 channel corresponds to a combination of collision induced dissociation, especially at higher collision energies, with the two fully submerged pathways (see blue and red lines in Figure 8) leading to a CS trend that is independent of the photon energy at low collision energies. It is not clear what contribution there might be at low collision energies from the fragmentation of the parent at higher photon energies.

The small signal observed for the m/z 39 product for H₂CNH^+• is expected to indicate the very small level of interconversion between the isomers via A4. While this interconversion is barrierless for both ions, the fact that this requires passage over twelve TSs and intermediate structures is expected to lead to a significant level of kinetic inhibition, especially compared to other barrierless pathways present.

In general, products at m/z 28 and 39 from both isomers show an initial decrease in CS section with increasing collision energy, indicative of barrierless pathways as suggested by our calculations. This is due to the higher collision energy inhibiting the formation of the van der Waals clusters (A1 and B1) in the entrance channels. For both isomers, above ∼0.7 eV, the m/z 28 product shows a gradual rise with increasing collision energy which is probably due to the collisional fragmentation of the parent ion into HCNH⁺ and H^•. The very minor m/z 39 channel also shows a higher collision energy feature for both isomers, with appearance energies of ∼1 eV for HCNH₂^+• and ∼4 eV for H₂CNH^+•. This is believed to correspond to a mechanism to give either c-C₃H₃⁺ or its high energy isomer CH₂CCH⁺ via a significant barrier, though no individual channel has been identified.

The m/z 27 channel is only observed with the H₂CNH^+• isomer, where it shows an increase with increasing photon energy, a trend that relates well with the endothermic channels Eqs 3 and 4 identified both in literature and our calculations. As a function of the collision energy, this channel shows a steady increase, again indicative of an endothermic proton transfer process which can proceed directly at higher collision energies.

7 Conclusion

Absolute cross sections for the reaction of methanimine (H₂CNH^+•) and aminomethylene (HCNH₂^+•) radical cations with C₂H₂ have been measured as a function of photon and collision energies under a single collision regime using a guided ion beam tandem mass spectrometer with VUV ionization for ion generation.

Through comparison of H₂CNH^+• formed both via direct and dissociative ionization and HCNH₂^+• formed exclusively through dissociative ionization, we are able to conclude that both isomers are formed selectively with a high level of correspondence between the two H₂CNH^+• sources. Ab-initio calculations of the most relevant stationary points on the potential energy surfaces allow for the rationalization of experimental results.

Our experiments clearly show that the reaction of both isomers with C₂H₂ lead to the formation of C₃H₄N⁺ products (detected at m/z 54), and by extension to an elongation of the main skeleton of the ion. Measured CSs show a negative dependence on the collision energy, indicative of barrierless and exothermic processes that are feasible at low temperatures. Thus, the title reactions can be viable intermediate steps to produce larger ions in the ISM as well as ionospheres of planets and their satellites.

But what about the structures of the C₃H₄N⁺ products? Theoretical calculations show that, for the reaction of HCNH₂^+•, the lowest energy isomer (H₂CCHCNH⁺, protonated vinyl cyanide, a.k.a. acrylonitrile) as well as two higher energy isomers, namely c-CHCHC(NH₂)⁺ and HCCCHNH₂⁺, can be generated via barrierless and exothermic pathways. It should, however, be noted that our experiment cannot yield any information about the branching ratios for the formation of the different C₃H₄N⁺ isomers. Since all of these pathways involve submerged barriers, no preference for the production of individual ionic products could be inferred from our calculations starting from the HCNH₂^+• isomer. Conversely, for the H₂CNH^+• reagent, the only barrierless pathway is the one leading to the production of protonated vinyl cyanide, which somewhat restricts the potential of the methanimine radical cation in the synthesis of interstellar nitrogen compounds.

Protonated vinyl cyanide has been proposed as a likely intermediate in star forming regions and in the atmosphere of Titan, where the title reaction can contribute to its formation. In particular, Cassini flyby missions to Titan’s upper atmosphere have detected the presence of a strong signal at m/z 54 that has been tentatively assigned to protonated vinyl cyanide on the basis of the large proton affinity of the neutral Vuitton et al. (2007). More recently, the spectroscopic detection of vinyl cyanide in Titan’s atmosphere has been confirmed using ALMA data Palmer et al. (2017). In the ISM, while vinyl cyanide is one of the most abundant nitriles detected in several environments such as the circumstellar shell of evolved carbon stars and dark clouds, the presence of its protonated form has been proposed on the basis of laboratory spectra but still not confirmed Martinez et al. (2013). It is therefore noteworthy that our title reaction, starting from the methanimine radical cation (H₂CNH^+•), could lead to protonated vinyl cyanide, from which acrylonitrile can be formed either through dissociative recombination with electrons Vigren et al. (2009) or proton transfer to species such as NH₃ Taquet et al. (2016).

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Datasets shown in Figures 2-5 are made available as ASCII files in the Zenodo repository accessible via DOI: 10.5281/zenodo.5419565. Further inquiries can be directed to the corresponding author.

Author Contributions

J-CG synthesized the reagent methanimine. CA, CR, RT, and JZ developed the experimental set-up and RT developed the data acquisition and analysis suite. All authors were involved in collecting the data. VR performed the data analysis of the experimental data. WG and DS worked on the PESs for both isomers while PT worked on the PESs for the fragmentation of reagent species. DA and VR worked on the interpretation of the results while DA, WG, and VR wrote the manuscript draft. All authors were involved in reviewing the manuscript.

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.

Acknowledgments

We are grateful to the DESIRS beamline team for their assistance during the synchrotron measurements and to the technical staff of SOLEIL for the smooth running of the facility under projects nos. 20180118 and 20190249. This work was supported by the European Union’s Horizon 2020 research and innovation program “Astro-Chemistry Origins” (ACO), Grant No. 811312. WG thanks the Swedish Research Council for a project grant (grant number 2019-04332). MP and JZ acknowledge support from the Ministry of Education, Youth and Sports of the Czech Republic (grant No. LTC20062). VR acknowledges funding for a PhD fellowship from the Dept. Physics, University of Trento. J-CG thanks the Centre National d’Etudes Spatiales (CNES) for a grant. CA, RT and CR acknowledge the synchrotron SOLEIL for the support to the associated CERISES setup since 2008 and subsistence expenses during beamtime periods. The Orsay group acknowledges support from the LabEx PALM (ANR-10-LABX-0039-PALM, Project ERACOP).

Supplementary Material

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

References

Alcaraz, C., Nicolas, C., Thissen, R., Žabka, J., and Dutuit, O. (2004). ¹⁵N⁺ +CD₄ and O⁺+ ¹³CO₂ State-Selected Ion−Molecule Reactions Relevant to the Chemistry of Planetary Ionospheres. J. Phys. Chem. A. 108, 9998–10009. doi:10.1021/jp0477755

CrossRef Full Text | Google Scholar

Anicich, V. G. (2003). An index of the Literature for Bimolecular Gas Phase Cation-Molecule Reaction Kinetics. JPL Publication 03-19. Available at: http://hdl.handle.net/2014/7981.

Google Scholar

Anicich, V. G. (1993). Evaluated Bimolecular Ion‐Molecule Gas Phase Kinetics of Positive Ions for Use in Modeling Planetary Atmospheres, Cometary Comae, and Interstellar Clouds. J. Phys. Chem. Reference Data. 22, 1469–1569. doi:10.1063/1.555940

CrossRef Full Text | Google Scholar

Anicich, V., and McEwan, M. J. (2001). Termolecular Ion-Molecule Reactions in Titan's Atmosphere III. Clustering of Ions. Icarus 154, 522–530. doi:10.1006/icar.2001.6696

CrossRef Full Text | Google Scholar

Anicich, V., Milligan, D. B., Fairley, D. A., and McEwan, M. J. (2000). Termolecular Ion-Molecule Reactions in Titan's Atmosphere, I Principal Ions with Principal Neutrals. Icarus 146, 118–124. doi:10.1006/icar.2000.6353

CrossRef Full Text | Google Scholar

Aponte, J. C., Elsila, J. E., Glavin, D. P., Milam, S. N., Charnley, S. B., and Dworkin, J. P. (2017). Pathways to Meteoritic Glycine and Methylamine. ACS Earth Space Chem. 1 (1), 3–13. doi:10.1021/acsearthspacechem.6b00014

PubMed Abstract | CrossRef Full Text | Google Scholar

Basiuk, V. A., and Bogillo, V. I. (2002). Theoretical Study of Amino Acid Precursor Formation in the Interstellar Medium. 2. Reaction of Methylenimine with CN Radical. Adv. Space Res. 30 (1), 1445–1450. doi:10.1016/S0273-1177(02)00511-2

CrossRef Full Text | Google Scholar

Bernstein, M. P., Sandford, S. A., Allamandola, L. J., Chang, S., and Scharberg, M. A. (1995). Organic Compounds Produced by Photolysis of Realistic Interstellar and Cometary Ice Analogs Containing Methanol. Ap. J. 454 (1), 327–344. doi:10.1086/176485

CrossRef Full Text | Google Scholar

Bøgelund, E. G., McGuire, B. A., Hogerheijde, M. R., van Dishoeck, E. F., and Ligterink, N. F. W. (2019). Methylamine and Other Simple N-Bearing Species in the Hot Cores NGC 6334I MM1-3. A. A 624, A82. doi:10.1051/0004-6361/201833676

CrossRef Full Text | Google Scholar

Brünken, S., Lipparini, F., Stoffels, A., Jusko, P., Redlich, B., Gauss, J., et al. (2019). Gas-Phase Vibrational Spectroscopy of the Hydrocarbon Cations l-C₃H⁺, HC₃H⁺, and c-C₃H₂⁺: Structures, Isomers, and the Influence of Ne-Tagging. J. Phys. Chem. A. 123 (37), 8053–8062. doi:10.1021/acs.jpca.9b06176

CrossRef Full Text | Google Scholar

Castrovilli, M. C., Bolognesi, P., Cartoni, A., Catone, D., O’Keeffe, P., Casavola, A. R., et al. (2014). Photofragmentation of Halogenated Pyrimidine Molecules in the VUV Range. J. Am. Soc. Mass. Spectrom. 25, 351–367. doi:10.1007/s13361-013-0783-x

CrossRef Full Text | Google Scholar

Chupka, W. A. (1971). Effect of thermal Energy on Ionization Efficiency Curves of Fragment Ions. J. Chem. Phys. 54, 1936–1947. doi:10.1063/1.1675122

CrossRef Full Text | Google Scholar

Cordiner, M. A., Teanby, N. A., Nixon, C. A., Vuitton, V., Thelen, A. E., and Charnley, S. B. (2019). ALMA Spectral Imaging of Titan Contemporaneous with Cassini's Grand Finale. A. j. 158 (2), 76. doi:10.3847/1538-3881/ab2d20

CrossRef Full Text | Google Scholar

Cunha de Miranda, B., Romanzin, C., Chefdeville, S., Vuitton, V., Žabka, J., Polášek, M., et al. (2015). Reactions of State-Selected Atomic Oxygen Ions O⁺(⁴S, ²D, ²P) with Methane. J. Phys. Chem. A. 119, 6082–6098. doi:10.1021/jp512846v

CrossRef Full Text | Google Scholar

Dickens, J. E., Irvine, W. M., DeVries, C. H., and Ohishi, M. (1997). Hydrogenation of Interstellar Molecules: A Survey for Methylenimine (CH₂NH). Ap. J. 479 (1), 307–312. doi:10.1086/303884

PubMed Abstract | CrossRef Full Text | Google Scholar

Ervin, K. M., and Armentrout, P. B. (1985). Translational Energy Dependence of Ar⁺+XY→ArX⁺+Y (XY=H₂,D₂,HD) from thermal to 30 eV c.M. J. Chem. Phys. 83, 166–189. doi:10.1063/1.449799

CrossRef Full Text | Google Scholar

Fathi, P., Geppert, W. D., Kaiser, A., and Ascenzi, D. (2016). Ion-neutral Reaction of the C₂H₂N⁺ Cation with C₂H₂: An Experimental and Theoretical Study. Mol. Astrophysics. 2, 1–11. doi:10.1016/j.molap.2015.09.002

CrossRef Full Text | Google Scholar

Frish, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2013). Gaussian 09, Revision D.01. Wallingford: Gaussian Inc. Available at: https://gaussian.com/

Gans, B., Garcia, G. A., Boyé-Péronne, S., Pratt, S. T., Guillemin, J.-C., Aguado, A., et al. (2019). Origin Band of the First Photoionizing Transition of Hydrogen Isocyanide. Phys. Chem. Chem. Phys. 21, 2337–2344. doi:10.1039/c8cp07737a

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldberg, N., Fiedler, A., and Schwarz, H. (1995). Cyanocarbene HCCN, Isocyanocarbene HCNC, Iminovinylidene HNCC, and Their Ionic Counterparts: A Combined Experimental and Theoretical Study. J. Phys. Chem. 99, 15327–15334. doi:10.1021/j100042a002

CrossRef Full Text | Google Scholar

Guillemin, J.-C., and Denis, J.-M. (1988). Synthese D'imines Lineaires Non-stabilisees Par Reactions Gaz-Solide Sous Vide(1). Tetrahedron 44 (14), 4431–4446. doi:10.1016/S0040-4020(01)86145-7

CrossRef Full Text | Google Scholar

Harland, P. W., and McIntosh, B. J. (1985). Enthalpies of Formation for the Isomeric Ions H_xCCN⁺ and H_xCNC⁺ (X = 0−3) by "monochromatic" Electron Impact on C₂N₂, CH₃CN and CH₃NC. Int. J. Mass Spectrom. Ion Process. 67 (1), 29–46. doi:10.1016/0168-1176(85)83035-4

CrossRef Full Text | Google Scholar

Heerma, W., Sarneel, M. M., and Dijkstra, G. (1986). Structure and Formation of Gaseous [C₃H₄N]⁺ Ions. Org. Mass. Spectrom. 21 (10), 681–687. doi:10.1002/oms.1210211013

CrossRef Full Text | Google Scholar

Holmes, J. L., Aubry, C., and Mayer, P. M. (2006). Assigning Structures to Ions in Mass Spectrometry. Boca Raton: Taylor & Francis.

Google Scholar

Holzmeier, F., Lang, M., Hader, K., Hemberger, P., and Fischer, I. (2013). H₂CN⁺ and H₂CNH⁺: New Insight into the Structure and Dynamics from Mass-Selected Threshold Photoelectron Spectra. J. Chem. Phys. 138 (21), 214310. doi:10.1063/1.4808050

CrossRef Full Text | Google Scholar

Hörst, S. M. (2017). Titan's Atmosphere and Climate. J. Geophys. Res. Planets. 122 (3), 432–482. doi:10.1002/2016je005240

CrossRef Full Text | Google Scholar

Imanaka, H., and Smith, M. A. (2007). Role of Photoionization in the Formation of Complex Organic Molecules in Titan’s Upper Atmosphere. Geophys. Res. Lett. 34 (2), L02204. doi:10.1029/2006gl028317

CrossRef Full Text | Google Scholar

Kim, J. K., Anicich, V. G., and Huntress, W. T. (1977). Product Distributions and Rate Constants for the Reactions of CH₃⁺, CH₄⁺, C₂H₂⁺, C₂H₃⁺, C₂H₄⁺, and C₂H₆⁺ Ions with CH₄, C₂H₂, C₂H₄, and C₂H₆. J. Phys. Chem. 81, 1798–1805. doi:10.1021/j100534a002

CrossRef Full Text | Google Scholar

Kramida, A., Ralchenko, Y., and Reader, J.NIST ASD Team (2020). NIST Atomic Spectra Database. version 5.8. Available at: http://physics.nist.gov/asd (Accessed July 25, 2021).

Google Scholar

Lago, A. F., and Baer, T. (2006). A Photoelectron Photoion Coincidence Study of the Vinyl Bromide and Tribromoethane Ion Dissociation Dynamics: Heats of Formation of C₂H₃⁺, C₂H₃Br, C₂H₃Br⁺, C₂H₃Br₂⁺, and C₂H₃Br₃. J. Phys. Chem. A. 110 (9), 3036–3041. doi:10.1021/jp053943x

PubMed Abstract | CrossRef Full Text | Google Scholar

Larsson, M., Geppert, W. D., and Nyman, G. (2012). Ion Chemistry in Space. Rep. Prog. Phys. 75 (6), 066901. doi:10.1088/0034-4885/75/6/066901

PubMed Abstract | CrossRef Full Text | Google Scholar

Ligterink, N. F. W., Calcutt, H., Coutens, A., Kristensen, L. E., Bourke, T. L., Drozdovskaya, M. N., et al. (2018). The ALMA-PILS Survey: Stringent Limits on Small Amines and Nitrogen-Oxides towards IRAS 16293-2422B. A. A 619, A28. doi:10.1051/0004-6361/201731980

CrossRef Full Text | Google Scholar

Martinez, O., Lattanzi, V., Thorwirth, S., and McCarthy, M. C. (2013). Detection of Protonated Vinyl Cyanide, CH₂CHCNH⁺, a Prototypical Branched Nitrile Cation. J. Chem. Phys. 138, 094316. doi:10.1063/1.4793316

CrossRef Full Text | Google Scholar

Mercier, B., Compin, M., Prevost, C., Bellec, G., Thissen, R., Dutuit, O., et al. (2000). Experimental and Theoretical Study of a Differentially Pumped Absorption Gas Cell Used as a Low Energy-Pass Filter in the Vacuum Ultraviolet Photon Energy Range. J. Vac. Sci. Technol. A. 18 (5), 2533–2541. doi:10.1116/1.1288196

CrossRef Full Text | Google Scholar

Milligan, D. B., Freeman, C. G., Maclagan, R. G. A. R., McEwan, M. J., Wilson, P. F., and Anicich, V. G. (2001). Termolecular Ion-Molecule Reactions in Titan's Atmosphere. II: The Structure of the Association Adducts of HCNH⁺ with C₂H₂ and C₂H₄. J. Am. Soc. Mass. Spectrom. 12, 557–564. doi:10.1016/S1044-0305(01)00215-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Minnhagen, L. (1973). Spectrum and the Energy Levels of Neutral Argon, Ar I. J. Opt. Soc. Am. 63, 1185–1198. doi:10.1364/JOSA.63.001185

CrossRef Full Text | Google Scholar

Nahon, L., de Oliveira, N., Garcia, G. A., Gil, J.-F., Pilette, B., Marcouillé, O., et al. (2012). DESIRS: a State-Of-The-Art VUV Beamline Featuring High Resolution and Variable Polarization for Spectroscopy and Dichroism at SOLEIL. J. Synchrotron Radiat. 19, 508–520. doi:10.1107/S0909049512010588

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, M. T., Rademakers, J., and Martin, J. M. L. (1994). Concerning the Heats of Formation of the [C,H₃,N]^+• Radical Cations. Chem. Phys. Lett. 221, 149–155. doi:10.1016/0009-2614(94)87031-4

CrossRef Full Text | Google Scholar

Nicolas, C., Alcaraz, C., Thissen, R., Žabka, J., and Dutuit, O. (2002). Effects of Ion Excitation on Charge Transfer Reactions of the Mars, Venus, and Earth Ionospheres. Planet. Space Sci. 50, 877–887. doi:10.1016/S0032-0633(02)00063-6

CrossRef Full Text | Google Scholar

Nixon, C. A., Lorenz, R. D., Achterberg, R. K., Buch, A., Coll, P., Clark, R. N., et al. (2018). Titan's Cold Case Files - Outstanding Questions after Cassini-Huygens. Planet. Space Sci. 155, 50–72. doi:10.1016/j.pss.2018.02.009

CrossRef Full Text | Google Scholar

Palmer, M. Y., Cordiner, M. A., Nixon, C. A., Charnley, S. B., Teanby, N. A., Kisiel, Z., et al. (2017). ALMA Detection and Astrobiological Potential of Vinyl Cyanide on Titan. Sci. Adv. 3, e1700022. doi:10.1126/sciadv.1700022

PubMed Abstract | CrossRef Full Text | Google Scholar

Polášek, M., Zins, E.-L., Alcaraz, C., Žabka, J., Křížová, V., Giacomozzi, L., et al. (2016). Selective Generation of the Radical Cation Isomers [CH₃CN]^•+ and [CH₂CNH]^•+ via VUV Photoionization of Different Neutral Precursors and Their Reactivity with C₂H₄. J. Phys. Chem. A. 120 (27), 5041–5052. doi:10.1021/acs.jpca.5b12757

CrossRef Full Text | Google Scholar

Richardson, V., Alcaraz, C., Geppert, W. D., Polášek, M., Romanzin, C., Sundelin, D., et al. (2021). The Reactivity of Methanimine Radical Cation (H₂CNH⁺) and its Isomer Aminomethylene (HCNH₂⁺) with Methane. Chem. Phys. Lett. 775, 138611. doi:10.1016/j.cplett.2021.138611

CrossRef Full Text | Google Scholar

Richardson, V., Tosi, P., Thissen, R., Alcaraz, C., Romanzin, C., Žabka, J., et al. unpublished data, manuscript in preparation.

Ruscic, B., and Bross, D. H. (2020). Active Thermochemical Tables (ATcT) Values Based on version 1.122p of the Thermochemical Network. Available at: https://atct.anl.gov/ (Accessed March 20, 2021).

Google Scholar

Salpin, J.-Y., Nguyen, M. T., Bouchoux, G., Gerbaux, P., and Flammang, R. (1999). Isomerization of Acetonitrile N-Methylide [CH₃CNCH₂]⁺ and N-Methylketenimine [CH₃NCCH₂]⁺ Radical Cations in the Gas Phase: Theoretical Study of the [C₃,H₅,N]⁺ Potential Energy Surface. J. Phys. Chem. A. 103, 938–946. doi:10.1021/jp983796l

CrossRef Full Text | Google Scholar

Sundelin, D., Ascenzi, D., Richardson, V., Alcaraz, C., Polášek, M., Romanzin, C., et al. (2021). The Reactivity of Methanimine Radical Cation (H₂CNH⁺) and its Isomer Aminomethylene (HCNH₂⁺) with C₂H₄. Chem. Phys. Lett. 777, 138677. doi:10.1016/j.cplett.2021.138677

CrossRef Full Text | Google Scholar

Suzuki, T., Ohishi, M., Hirota, T., Saito, M., Majumdar, L., and Wakelam, V. (2016). Survey Observations of a Possible Glycine Precursor, Methanimine (CH₂NH). Ap. J. 825 (1), 79. doi:10.3847/0004-637x/825/1/79

CrossRef Full Text | Google Scholar

Taquet, V., Wirström, E. S., and Charnley, S. B. (2016). Formation and Recondensation of Complex Organic Molecules during Protostellar Luminosity Outbursts. Ap. J. 821 (1), 46. doi:10.3847/0004-637x/821/1/46

CrossRef Full Text | Google Scholar

Teloy, E., and Gerlich, D. (1974). Integral Cross Sections for Ion-Molecule Reactions. I. The Guided Beam Technique. Chem. Phys. 4, 417–427. doi:10.1016/0301-0104(74)85008-1

CrossRef Full Text | Google Scholar

Theule, P., Borget, F., Mispelaer, F., Danger, G., Duvernay, F., Guillemin, J. C., et al. (2011). Hydrogenation of Solid Hydrogen Cyanide HCN and Methanimine CH₂NH at Low Temperature. Astron. Astrophysics 534, A64. ISSN 1432-0746. doi:10.1051/0004-6361/201117494

CrossRef Full Text | Google Scholar

Traeger, J. C., and McLoughlin, R. G. (1981). Absolute Heats of Formation for Gas-phase Cations. J. Am. Chem. Soc. 103, 3647–3652. doi:10.1021/ja00403a006

CrossRef Full Text | Google Scholar

Vigren, E., Hamberg, M., Zhaunerchyk, V., Kamińska, M., Thomas, R. D., Larsson, M., et al. (2009). The Dissociative Recombination of Protonated Acrylonitrile, CH₂CHCNH⁺, with Implications for the Nitrile Chemistry in Dark Molecular Clouds and the Upper Atmosphere of Titan. Ap. J. 695 (1), 317–324. doi:10.1088/0004-637x/695/1/317

CrossRef Full Text | Google Scholar

Vigren, E., Semaniak, J., Hamberg, M., Zhaunerchyk, V., Kaminska, M., Thomas, R. D., et al. (2012). Dissociative Recombination of Nitrile Ions with Implications for Titan's Upper Atmosphere. Planet. Space Sci. 60 (1), 102–106. doi:10.1016/j.pss.2011.03.001

CrossRef Full Text | Google Scholar

Vinatier, S., Bézard, B., Fouchet, T., Teanby, N., Dekok, R., Irwin, P., et al. (2007). Vertical Abundance Profiles of Hydrocarbons in Titan's Atmosphere at 15° S and 80° N Retrieved from Cassini/CIRS Spectra. Icarus 188 (1), 120–138. doi:10.1016/j.icarus.2006.10.031

CrossRef Full Text | Google Scholar

Vuitton, V., Yelle, R. V., Klippenstein, S. J., Hörst, S. M., and Lavvas, P. (2019). Simulating the Density of Organic Species in the Atmosphere of Titan with a Coupled Ion-Neutral Photochemical Model. Icarus 324, 120–197. doi:10.1016/j.icarus.2018.06.013

CrossRef Full Text | Google Scholar

Vuitton, V., Yelle, R. V., and McEwan, M. J. (2007). Ion Chemistry and N-Containing Molecules in Titan's Upper Atmosphere. Icarus 191 (2), 722–742. doi:10.1016/j.icarus.2007.06.023

CrossRef Full Text | Google Scholar

Wang, Z.-C., Cole, C. A., Snow, T. P., and Bierbaum, V. M. (2015). Experimental and Computational Studies of the Formation Mechanism of Protonated Interstellar Diazines. Ap. J. 798 (102), 102. doi:10.1088/0004-637X/798/2/102

CrossRef Full Text | Google Scholar

Westlake, J. H., Waite, J. H., Carrasco, N., Richard, M., and Cravens, T. (2014). The Role of Ion-Molecule Reactions in the Growth of Heavy Ions in Titan's Ionosphere. J. Geophys. Res. Space Phys. 119 (7), 5951–5963. doi:10.1002/2014ja020208

CrossRef Full Text | Google Scholar

Widicus Weaver, S. L., Laas, J. C., Zou, L., Kroll, J. A., Rad, M. L., Hays, B. M., et al. (2017). Deep, Broadband Spectral Line Surveys of Molecule-Rich Interstellar Clouds. Ap. JS. 232 (1), 3. doi:10.3847/1538-4365/aa8098

CrossRef Full Text | Google Scholar

Woods, R. C., Gudeman, C. S., Dickman, R. L., Goldsmith, P. F., Huguenin, G. R., Irvine, W. M., et al. (1983). The HCO(+)/HOC(+) Abundance Ratio in Molecular Clouds. Ap. J. 270, 583–588. doi:10.1086/161150

CrossRef Full Text | Google Scholar