# Reactivity of the Ethenium Cation (C2H5+) with Ethyne (C2H2): A Combined Experimental and Theoretical Study

Vincent Richardson, Miroslav Polášek, Claire Romanzin, Paolo Tosi, Roland

Thissen, Christian Alcaraz, Ján Žabka, Daniela Ascenzi

► To cite this version:

Vincent Richardson, Miroslav Polášek, Claire Romanzin, Paolo Tosi, Roland Thissen, et al.. Reactivity of the Ethenium Cation (C2H5+) with Ethyne (C2H2): A Combined Experimental and Theoretical Study. Molecules, 2024, 29 (4), pp.810. 10.3390/molecules29040810. [https://doi.org/10.3390/molecules29040810](https://doi.org/10.3390/molecules29040810). [hal-04455989](https://HAL.inistre.fr/hal-04455989)

HAL Id: hal-04455989

https://hal.science/hal-04455989v1

Submitted on 13 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Reactivity of the Ethenium Cation ($C_2H_5^+$) with Ethyne ($C_2H_2$): A Combined Experimental and Theoretical Study

Vincent Richardson <fnref n="1" />,<fnref n="2" /> ID, Miroslav Polášek <fnref n="3" /> ID, Claire Romanzin <fnref n="4" />,<fnref n="5" />, Paolo Tosi <fnref n="1" /> ID, Roland Thissen <fnref n="4" />,<fnref n="5" /> ID, Christian Alcaraz <fnref n="4" />,<fnref n="5" /> ID, Ján Žabka <fnref n="3" /> ID and Daniela Ascenzi <fnref n="1" />,* ID

**Abstract:** The gas-phase reaction between the ethyl cation ($C_2H_5^+$) and ethyne ($C_2H_2$) is re-investigated by measuring absolute reactive cross sections (CSs) and branching ratios (BRs) as a function of collision energy, in the thermal and hyperthermal energy range, via tandem-guided ion beam mass spectrometry under single collision conditions. Dissociative photoionization of $C_2H_5Br$ using tuneable VUV radiation in the range 10.5–14.0 eV is employed to generate $C_2H_5^+$, which has also allowed us to explore the impact of increasing (vibrational) excitation on the reactivity. Reactivity experiments are complemented by theoretical calculations, at the G4 level of theory, of the relative energies and structures of the most relevant stationary points on the reactive potential energy hypersurface (PES) and by mass-analyzed ion kinetic energy (MIKE) spectrometry experiments to probe the metastable decomposition from the $[C_4H_7]^+$ PES and elucidate the underlying reaction mechanisms. Two main product channels have been identified at a centre-of-mass collision energy of ~0.1 eV: (a) $C_3H_3^+$+CH$_4$, with BR = 0.76 ± 0.05 and (b) $C_4H_5^+$+H$_2$, with BR = 0.22 ± 0.02. A third channel giving $C_2H_3^+$ in association with $C_2H_4$ is shown to emerge at both high internal excitation of $C_2H_5^+$ and high collision energies. From CS measurements, energy-dependent total rate constants in the range $4.3 \times 10^{-11} - 5.2 \times 10^{-10} \text{ cm}^3\cdot\text{molecule}^{-1}\cdot\text{s}^{-1}$ have been obtained. Theoretical calculations indicate that both channels stem from a common covalently bound intermediate, $CH_3CH_2CHCH^+$, from which barrierless and exothermic pathways exist for the production of both cyclic c-$C_3H_3^+$ and linear H$_2CCCH^+$ isomers of the main product channel. For the minor $C_4H_5^+$ product, two isomers are energetically accessible: the three-member cyclic isomer c-$C_3H_2(CH_3)^+$ and the higher energy linear structure $CH_2CHCCH_2^+$, but their formation requires multiple isomerization steps and passages via transition states lying only 0.11 eV below the reagents' energy, thus explaining the smaller BR. Results have implications for the modeling of hydrocarbon chemistry in the interstellar medium and the atmospheres of planets and satellites as well as in laboratory plasmas (e.g., plasma-enhanced chemical vapor deposition of carbon nanotubes and diamond-like carbon films).

**Keywords:** ion-molecule reactions; Titan; planetary atmospheres; interstellar medium; plasma; acetylene; Cassini mission

## 1. Introduction

The gas-phase reactivity of ions is believed to play a significant role in the formation of complex organic molecules (COMs) in both the interstellar medium (ISM) and the upper atmospheres of planets and their satellites [1,2]. Though ions are typically less abundant,

### check for updates

**Citation:** Richardson, V.; Polášek, M.; Romanzin, C.; Tosi, P.; Thissen, R.; Alcaraz, C.; Žabka, J.; Ascenzi, D. Reactivity of the Ethenium Cation ($C_2H_5^+$) with Ethyne ($C_2H_2$): A Combined Experimental and Theoretical Study. *Molecules* **2024**, *29*, 810. https://doi.org/10.3390/molecules29040810

Academic Editor: Yongle Li

Received: 22 December 2023

Revised: 24 January 2024

Accepted: 29 January 2024

Published: 9 February 2024

### CC BY

**Copyright:** © 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Department of Physics, University of Trento, 38123 Trento, Italy; vincent.richardson@liverpool.ac.uk (V.R.); paolo.tosi@unitn.it (P.T.)</footnote>

<footnote n="2">² Department of Physics, University of Liverpool, Oxford Street, Liverpool L69 7ZE, UK</footnote>

<footnote n="3">³ J. Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences, Dolejšškova 2155/3, 182 23 Prague, Czech Republic; miroslav.polasek@jh-inst.cas.cz (M.P.); jan.zabka@jh-inst.cas.cz (J.Ž.)</footnote>

<footnote n="4">⁴ Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France; claire.romanzin@universite-paris-saclay.fr (C.R.); roland.thissen@universite-paris-saclay.fr (R.T.); christian.alcaraz@universite-paris-saclay.fr (C.A.)</footnote>

<footnote n="5">⁵ Synchrotron Soleil, L'Orme des Merisiers, 91190 Saint-Aubin, France</footnote>

<footnote>* Correspondence: daniela.ascenzi@unitn.it</footnote>
and therefore harder to detect, in the ISM than neutrals, their greater reactivity makes them more significant to the overall chemistry than their number densities indicate. Notably, the rapid reaction of charged species with neutral partners can lead to the formation of complex ions that act as important chemical intermediates since they can recombine with electrons and/or anions to produce increasingly complex neutral species under astrophysical conditions [3,4].

In our Solar System, Titan, the largest moon of Saturn, is an environment of particular interest to astrochemistry as, similar to that of the early Earth, the atmosphere of Titan is dense and dominated by $N_2$, with methane ($CH_4$) as the second most abundant species [5,6]. Under such conditions, a rich ion chemistry is driven by ionization and dissociative ionization of its primary components by a combination of extreme ultraviolet radiation from the Sun [7], energetic magnetospheric electrons and ions [8], and cosmic rays. The ions and radicals that are formed are then able to go on to react with surrounding neutrals, thereby leading to the formation of complex molecular species, including some with prebiotic potential [3,9–12]. Thus, “complex” ions with m/z ratios as large as 350 have been detected [9,13–17] by the ion neutral mass spectrometer (INMS) and Cassini plasma spectrometer-ion beam sensor (CAPS-IBS) aboard the Cassini spacecraft that made several targeted flybys of Titan [18,19].

With an estimated density of $2 \times 10^2 \text{ cm}^{-3}$, $C_2H_5^+$ is thought to be the second most abundant ion in Titan's ionosphere, where it is expected to be the major contributor to the intense m/z 29 peak observed in the mass spectra of Titan's ionosphere measured by the Cassini INMS [14,20]. Outside Titan's atmosphere, the presence of $C_2H_5^+$ has been postulated on the basis of results from the Cassini Grand Finale mission [21–24], where the INMS was used to measure the composition of Saturn's upper atmosphere and its chemical interactions with material originating in the rings [22]. Further away from our Sun, although the New Horizons mission to Pluto did not make a direct measurement of its ionosphere, $C_2H_5^+$ has been suggested as one of the most abundant ions on the basis of the presence of hydrocarbons [25,26]. In addition, recent ALMA (Atacama large millimeter/submillimeter array) interferometric measurements have detected the presence of N-containing organics (HCN and its isomer HNC, $CH_3CN$) and propyne ($CH_3CCH$) in Pluto's atmosphere [27], thus suggesting the presence of a Titan-like chemistry. $C_2H_5^+$ has also been unambiguously identified among the diverse range of ions detected by the ROSINA high-resolution mass spectrometer employed in the Rosetta mission to comet 67P/Churyumov-Gerasimenko [28]. Though $C_2H_5^+$ has yet to be identified in the ISM, it could feasibly be formed via the protonation of $C_2H_4$, which has been detected [29,30]. $C_2H_2$ (ethyne, also known as acetylene) is not only one of the most abundant neutral species in Titan's ionosphere [14], but has also been detected in a range of astrochemical environments, including, but not limited to, the ionosphere of Pluto [31,32] and the ISM [33–37].

The reaction of $C_2H_5^+$ with $C_2H_2$ is of particular interest as it has been shown to be one of the four most deterministic ion-neutral reactions for models of Titan's ionosphere [38]. From previous experimental studies [39–41], it is well established that the main reaction products are $C_3H_3^+$ plus CH$_4$ and $C_4H_5^+$ plus H$_2$. A signature of the presence of $C_3H_3^+$ and $C_4H_5^+$ on Titan's upper atmosphere has been demonstrated by the ion densities at m/z 39 and 53 obtained from the Cassini INMS data, and assigned, at least partly, to the $C_3H_3^+$ and $C_4H_5^+$ ions, respectively [42]. More recently, the small cyclic molecule cyclopropenylidene (c-C$_3H_2$) has been detected on Titan's atmosphere using spectroscopic observations made with the ALMA telescope [43], and dissociative electron recombination of the cyclopropenylidene cation (c-C$_3H_3^+$) is among the mechanisms that have been proposed for the synthesis of c-C$_3H_2$. The other ionic product of the reaction, $C_4H_5^+$, has previously been identified as a key intermediate in the formation of benzene (c-C$_6H_6$) in Titan's ionosphere [44] and, with both $C_2H_5^+$ and $C_2H_2$ having been identified in a range of other astrochemical environments, this reaction could well play a significant role in the formation of benzene in such environments.
In addition to the astrochemical significance of the title reaction, C₂H₅⁺ and C₂H₂ are important species in laboratory plasmas, such as those used in microwave or radiofrequency discharges in methane for plasma-enhanced chemical vapor deposition of carbon nanotubes and diamond-like carbon films [45,46], atmospheric pressure plasma reactors for non-oxidative conversion of CH₄ into higher hydrocarbons and H₂ [47], or for the generation of value-added fuels from CH₄ and CO₂ [48–50]. For instance, chemical kinetics models aimed at exploring the plasma chemistry operative in nanosecond pulsed discharges or dielectric barrier discharges for the partial oxidation and dry reforming of methane predict the presence of C₂H₅⁺ among the relevant hydrocarbon ions [51–54].

Though the reactivity of C₂H₅⁺ with C₂H₂ has been studied previously [39–41], our recent experimental observations [55,56] have found some inconsistencies in the relative BRs of the two main ionic products (C₃H₃⁺ and C₄H₅⁺), thereby prompting this re-evaluation of the existing literature values. This is required because accurate kinetics data, in terms of both rate constants and BRs, are essential to improve the reliability and predictive power of chemical models used both in natural and laboratory plasmas.

In this work, we present experimental BRs and reaction cross sections (CSs), as well as the results of computational calculations on the potential Energy Surface (PES) for the reaction of C₂H₅⁺ with C₂H₂. Additionally, we present data on the CSs as a function of both the collision energy and the photon energy used to generate the ions (a proxy of the internal energy) in order to probe their respective impacts, the latter being of particular relevance to the laboratory plasmas involved in the dry reforming of methane [57]. From this, we are able to identify an alternative major product at threshold from that reported in previous literature [39,40], along with BRs and reaction CSs for other minor products, including a previously unobserved channel that predominates at the higher energies relevant to conversion plasmas. These experimental observations have been complemented with computational work on the reaction PES and further experimental measurements via the technique of mass-analyzed ion kinetic energy (MIKE) spectrometry to probe the underlying dynamics behind the observed products and the implications for both the neutral fragments and the ionic isomers formed.

## 2. Experimental Results

### 2.1. Generation of C₂H₅⁺

The dissociative photoionization of C₂H₅Br has been studied extensively [58–63] with three potential pathways being relevant to the E<sub>phot</sub> range used as part of this work. These three pathways, along with their respective lower lying dissociation limits relative to the respective ground state C₂H₅Br in its zero-point vibrational level from [58], are given by Equations (1)–(3).

$$C_2H_5Br + h\nu \rightarrow C_2H_5^+ + Br^\bullet \qquad + 11.13 \text{ eV} \quad (1)$$

$$\rightarrow C_2H_4^{+\bullet} + HBr \qquad + 11.21 \text{ eV} \quad (2)$$

$$\rightarrow C_2H_3^+ + H_2 + Br^\bullet \qquad + 13.24 \text{ eV} \quad (3)$$

The fragmentation into C₂H₅⁺ and Br is the only channel for which experimental appearance energy (AE) values are available in the literature with the most recent values being 11.132 ± 0.008 eV [59] and 11.130 ± 0.005 eV [60]. Due to the presence of other fragments in addition to the desired C₂H₅⁺ ion, we conclude that we also need to consider the potential isobaric contaminant arising from ¹³CCH₄⁺• as this has a similar formation enthalpy to our target ion. In order to characterize the fragmentation process, we have recorded photoionization efficiency (PIE) curves (i.e., ion yields as a function of photon energy) for pathways (1)–(3), and the results are shown in Figure 1.

For the C₂H₅⁺ fragment ion (m/z 29) an AE of 11.02 ± 0.20 eV has been measured, which is consistent with the above-mentioned literature values [59,60] as well as older determinations [61,62]. As the intensity of the m/z 28 signal is substantially less than the predominant m/z 29 fragment, at least up to E<sub>phot</sub> ~13 eV, we conclude that the ¹³CCH₄⁺•
impurity should be limited to < 1% of the m/z 29 signal. However, as the m/z 28 signal increases above $E_{phot} = 13.0$ eV, reactivity data with the C₂H₅⁺ ion have not been measured at $E_{phot} > 13.0$ eV in order to limit the impact from the isobaric contaminant ¹³CCH₄⁺•. Further consideration on the appearance energy of the m/z 28 fragment is given in the Supplementary Material.

For completeness, we have also conducted a PIE scan for the m/z 27 fragment (C₂H₃⁺), also shown in the bottom panel of Figure 1. Though the small signal intensity does not allow for a precise determination of the AE, for which no data exists in the literature, this pathway shows a sharp AE at ~ 13.2 eV, in excellent agreement with the literature thermodynamic dissociation limit of 13.24 eV [58]. As a result, when $E_{phot} + E_{CM} > 13.2$ eV, we may expect to observe a m/z 27 product corresponding to the collisional fragmentation of C₂H₅⁺.

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Photoionization efficiency (PIE) curves for the dissociative photoionization of C₂H₅Br in the 10.5–14 eV photon energy range: filled red dots are for the C₂H₅⁺ ion (m/z 29), open blue squares are for the C₂H₄⁺ ion (m/z 28), and green line is for the C₂H₃⁺ ion (m/z 27). Insets (a) and (b) show enlarged plots of the threshold regions for the corresponding m/z values. The solid black lines indicate the linear extrapolations used to extract the experimental appearance energy (AE) (reported with a vertical arrow).</figcaption></figure>

The other observation of note from Figure 1 is the m/z 29 trend as a function of $E_{phot}$, which shows two approximately linear increases interspersed with a plateau between ~11.7 and ~12.0 eV. The large increase in the PIE signal at about 12 eV might be due to the opening of a dissociation channel in which the C₂H₅⁺ ion is produced in combination with a Br atom in the excited spin-orbit state Br*(²P₁/₂) that lies 0.457 eV above the ground state [64]. From calculations of the spin-orbit-coupled potential energy curves (PEC) for the fragmentation of the ethylbromide cation [58], it appears that the electronic ground state of C₂H₅Br⁺ is split into 1A' and 1A'' states, both dissociating to ground-state C₂H₅⁺ plus Br, while the PEC of the electronically excited 2A' state is dissociative along the C-Br stretch coordinate and correlates with C₂H₅⁺ + Br*. The calculated energetics for the two processes are 11.13 eV and 11.59 eV, respectively. This is in agreement with the He(I) photoelectron spectrum of ethyl bromide [61,65] from which the onset of the electronically excited state 2A' is around 12 eV. In the same work, direct fragmentation from the electronically excited state of C₂H₅Br⁺ was observed without a prior radiationless transition to the ground electronic state. From this, we can speculate that the increase in the C₂H₅⁺ PIE signal around $E_{phot} = 12$ eV might be related to the dissociative ionization
from the first electronically excited state of C₂H₅Br⁺. Here, data as a function of E<sub>CM</sub> have been collected at E<sub>phot</sub> = 11.3, 12, and 13 eV in order to explore the different features of the C₂H₅⁺ PIE curve.

## 2.2. Reaction with C₂H₂

At low photon energies, the reaction of C₂H₅⁺ with C₂H₂ yields products at *m/z* 39 and 53; meanwhile, at higher E<sub>phot</sub>, an additional product is observed at *m/z* 27. The most probable product channels are indicated as Equations (4)–(9), and the reaction enthalpies estimated from the literature values are summarized and compared with results from our computational calculations in Table 1.

$$C_2H_5^+ + C_2H_2 \rightarrow c-C_3H_3^+ + CH_4 \quad m/z 39 \quad (4)$$

$$\rightarrow H_2CCCH^+ + CH_4 \quad m/z 39 \quad (5)$$

$$\rightarrow c-C_3H_2(CH_3)^+ + H_2 \quad m/z 53 \quad (6)$$

$$\rightarrow CH_2CHCCH_2^+ + H_2 \quad m/z 53 \quad (7)$$

$$\rightarrow C_2H_3^+ + C_2H_4 \quad m/z 27 \quad (8)$$

$$\rightarrow C_2H_3^+ + H_2 + C_2H_2 \quad m/z 27 \quad (9)$$

Previous calculations [66] indicate that there are two isomers of C₃H₃⁺ that are thermodynamically feasible for this reaction, namely c-C₃H₃⁺ and H₂CCCH⁺, with c-C₃H₃⁺ being the lowest energy isomer by 1.2–1.3 eV. The energetics of C₄H₅⁺ isomers have also been studied previously [67,68], with the most stable structure being the cyclic one reported in reaction (6). For completeness, as part of the computational work presented in Section 3, we have also looked for pathways leading to the higher energy isomers of [C₄H₅]⁺, and the only other isomer for which a submerged formation pathway has been found (at the G4 level of theory) is CH₂CHCCH₂⁺, reported in reaction (7). We have, therefore, limited our discussions of the *m/z* 53 channel to consideration of these two ions.

<table><caption>**Table 1.** Reaction enthalpies (ΔH°₂₉₸) for the identified products of the reaction of C₂H₅⁺ with C₂H₂ both from literature<fnref n="1" /> and our calculations²,³. Values are given in eV.</caption><thead><tr><th>Reaction Products</th><th>Equation</th><th>ΔH°<sub>298</sub>, lit. Values<sup>1</sup></th><th>ωB97X-D/cc-pVTZ<sup>2</sup></th><th>G4<sup>3</sup></th></tr></thead><tbody><tr><td>c-C<sub>3</sub>H<sub>3</sub><sup>+</sup> + CH<sub>4</sub></td><td>(4)</td><td>-1.36 ± 0.10<sup>4</sup></td><td>-1.57</td><td>-1.36</td></tr><tr><td>H<sub>2</sub>CCCH<sup>+</sup> + CH<sub>4</sub></td><td>(5)</td><td>-0.16 ± 0.10<sup>4</sup></td><td>-0.25</td><td>-0.19</td></tr><tr><td>c-C<sub>3</sub>H<sub>2</sub>(CH<sub>3</sub>)<sup>+</sup> + H<sub>2</sub></td><td>(6)</td><td>-1.68 ± 0.10<sup>5</sup></td><td>-1.77</td><td>-1.57</td></tr><tr><td>CH<sub>2</sub>CHCCH<sub>2</sub><sup>+</sup> + H<sub>2</sub></td><td>(7)</td><td>-1.01 ± 0.10<sup>5</sup></td><td>-</td><td>-0.92</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub><sup>+</sup> + C<sub>2</sub>H<sub>4</sub></td><td>(8)</td><td>+0.39 ± 0.08<sup>6</sup></td><td>+0.37</td><td>+0.36</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub><sup>+</sup> + H<sub>2</sub> + C<sub>2</sub>H<sub>2</sub></td><td>(9)</td><td>+2.21 ± 0.08<sup>6</sup></td><td>+2.29</td><td>+2.09</td></tr></tbody></table>

The rate constant and product distributions (BRs) of the reaction between C₂H₅⁺ and C₂H₂ have been measured previously using an ion cyclotron resonance (ICR) mass spectrometer, where C₂H₅⁺ ions were generated using electron ionization of C₂H₅Cl and C₂H₆ with a nominal electron energy of 16 eV [39]. The measured total rate constant was 1.9 × 10⁻¹⁰ cm³·molecule⁻<fnref n="1" />·s⁻<fnref n="1" /> (somewhat larger than an older value of 1.0 × 10⁻¹⁰ cm³·molecule⁻<fnref n="1" />·s⁻<fnref n="1" /> reported in [73]), and the two products at *m/z* 39 and 53 were observed with BRs of 0.36 and 0.64, respectively. A comparison of the BRs collected as part of this work (at E<sub>phot</sub> = 11.3 eV, E<sub>CM</sub> ~ 0.1 eV) with those from Kim et al. [39] is

<footnote n="1">¹ All values in this column have been estimated using the following formation enthalpies: ΔfH°₂₉₸(C₂H₅⁺) = 9.36 ± 0.03 eV [69] and ΔfH°₂₉₸(C₂H₂) = 2.35 ± 0.01 eV [70]. ² Present work: calculations at the ωB97X-D/cc-pVTZ level of theory have been performed using Gaussian 16 [71]. ³ Present work: calculations at the G4 level have been performed using Gaussian 16 [71]. ⁴ Formation enthalpy for CH₄ is from NIST [70], while the values for c-C₃H₃⁺ and H₂CCCH⁺ are from Ref. [69]. ⁵ Formation enthalpies of c-C₃H₂(CH₃)⁺ and CH₂CHCCH₂⁺ are from calculations at the MP4SDTQ/6-311++G(2df,p) level of theory from Cunje et al. [67]. ⁶ Formation enthalpy for C₂H₄ is from NIST [70], while that of C₂H₃⁺ is from Lago et al. [72].</footnote>
made in Table 2. It should be noted that, while previous experiments found the reaction channel [C₄H₅]⁺ plus H₂ to be the dominant one with a BR = 0.64, present results show that the main channel is the production of [C₃H₃]⁺ in association with a CH₄ molecule, with a BR = 0.76 ± 0.05. A more detailed consideration on this discrepancy is given in Sections 4.2 and 4.3.

The data presented here therefore not only involve the identification of a new product channel at *m/z* 27 (C₂H₃⁺) that is particularly relevant for higher energy plasmas, but also a significant update to the BRs at low internal energies. This latter point, in particular, could potentially lead to a marked change in the predicted abundances of higher mass ions in Titan's ionosphere [38], with further modeling work required for a quantitative assessment to be made.

<table><caption>**Table 2.** Branching ratios (BRs) for the products for the reaction of C₂H₅⁺ with C₂H₂ measured at *E*<sub>phot</sub> = 11.3 eV and *E*<sub>CM</sub> ~ 0.1 eV.</caption><thead><tr><th>Ionic Reaction Products</th><th>Product Mass (m/z)</th><th>Equations</th><th>BR, This Work</th><th>Literature BR <sup>1</sup></th></tr></thead><tbody><tr><td>C<sub>2</sub>H<sub>3</sub><sup>+</sup></td><td>27</td><td>(8)/(9)</td><td>0.02 ± 0.02</td><td>0.00</td></tr><tr><td>C<sub>3</sub>H<sub>3</sub><sup>+</sup></td><td>39</td><td>(4)/(5)</td><td>0.76 ± 0.05</td><td>0.36</td></tr><tr><td>C<sub>4</sub>H<sub>5</sub><sup>+</sup></td><td>53</td><td>(6)/(7)</td><td>0.22 ± 0.02</td><td>0.64</td></tr></tbody></table>

Table 3 reports the energy-dependent total rate constants $k_{tot}(E_{ave})$ obtained from the CSs measured as a function of the collision energy at $E_{phot} = 11.3$ eV (to be presented in Section 2.2.2). The energy-dependent total rate constants have been estimated using the expression:

$$k_{tot}(E_{ave}) = \langle v \rangle \cdot \sigma_{tot} \qquad (10)$$

where $\sigma_{tot}$ is the total reaction CS (summed over all the product channels), and $\langle v \rangle$ is the average relative velocity that can be estimated from the collision energy $E_{CM}$ (see [74,75] for a more detailed treatment). It should be noted that the rates obtained at the hyperthermal velocities probed in our experiment are comparable with the value $1.9 \times 10^{-10} \text{ cm}^3 \cdot \text{molecule}^{-1} \cdot \text{s}^{-1}$ determined via ICR mass spectrometry [39] and are smaller than the rate estimated from a simple Langevin capture model ($1.18 \times 10^{-9} \text{ cm}^3 \cdot \text{molecule}^{-1} \cdot \text{s}^{-1}$, see last column of Table 3).

The significant under-observation with respect to the Langevin rate is unsurprising, as the model assumes that if the reactants have enough energy to surmount the centrifugal barrier, the reaction forms from the resultant ion-molecule complex with unit probability, with no accounting for any dissociation back to the reactants. Not only has this assumption been shown to be false in a number of previous studies [55,56,76], but the limitations of capture theories in providing quantitative assessments of reaction rates, especially for complex systems such as in the current study, have been explored extensively in literature [77].

<table><caption>**Table 3.** Energy-dependent total rate constants ($k_{tot}$) measured at $E_{phot} = 11.3$ eV and at various values of the average collision energy ($E_{ave}$).</caption><thead><tr><th><i>E</i><sub>ave</sub>, eV</th><th><i>k</i><sub>tot</sub><sup>1</sup></th><th><i>k</i><sub>tot</sub>/<i>k</i><sub>L</sub><sup>2</sup></th></tr></thead><tbody><tr><td>0.096</td><td>1.7 × 10<sup>-10</sup></td><td>0.14</td></tr><tr><td>0.12</td><td>1.5 × 10<sup>-10</sup></td><td>0.13</td></tr><tr><td>0.52</td><td>4.3 × 10<sup>-11</sup></td><td>0.04</td></tr><tr><td>2.02</td><td>1.4 × 10<sup>-10</sup></td><td>0.12</td></tr><tr><td>4.52</td><td>5.2 × 10<sup>-10</sup></td><td>0.44</td></tr></tbody></table>

<footnote n="1">¹ Kim et al. [39].</footnote>

<footnote n="1">¹ $k_{tot}$ is the total (i.e., summed over all the product channels) rate constant (in cm³·molecule⁻¹·s⁻¹) at the specified average collision energy ($E_{ave}$), estimated as detailed in the text. ² $k_L$ is the Langevin rate constant for C₂H₅⁺ plus C₂H₂ equal to $1.18 \times 10^{-9} \text{ cm}^3 \cdot \text{molecule}^{-1} \cdot \text{s}^{-1}$.</footnote>
### 2.2.1. Cross Sections as a Function of the Photon Energy

Absolute CSs as a function of $E_{phot}$ have been collected at low collision energy for the m/z 27 ($C_2H_3^+$), 39 ($C_3H_3^+$) and 53 ($C_4H_5^+$) products, with results presented in Figure 2. The two products that are most abundant at low $E_{phot}$, namely m/z 39 and 53, show a marked decrease in the CS with increasing $E_{phot}$, while the m/z 27 product shows an increase above its AE of ~11.4 eV. Importantly, we note that the m/z 39 and 53 signals decrease below the AE of the m/z 27 product and so, while we expected the competition between the channels to contribute to the decrease as $E_{phot}$ is increased, we must also consider what leads to the initial decrease in CS for the m/z 39 and 53 products.

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Cross sections (CSs) as a function of $E_{phot}$ for the m/z 27, 39, and 53 products of the reaction of $C_2H_5^+$ with $C_2H_2$ at $E_{CM}$ = 0.12 eV. Dashed lines represent the photon energies at which CSs have been measured as a function of the collision energy for the various product channels.</figcaption></figure>

As the m/z 27 product ($C_2H_3^+$) can correspond with either the endothermic proton transfer, represented by reaction (8), or the collision-induced dissociation (CID) of the parent ion, reaction (9), further data are required to be able to elucidate the relative contributions of these two channels. To do this, we have also performed a more limited study on the reactivity of $C_2D_5^+$ with $C_2H_2$, where the direct proton (deuteron in this case) transfer product $C_2H_2D^+$ will appear at m/z 28, while the CID of the parent ion to give $C_2D_3^+$ will be shifted to m/z 30, with complex- and adduct-mediated processes leading to a level of H/D scrambling to give products at m/z 28, 29, and 30. The trends of these three products as a function of the photon energy are shown in Figure 3, where a clear preference is shown for the m/z 28 product, indicating that the observed signal at low collision energies predominately corresponds to the proton transfer.

The observation of minor products at m/z 29 and 30 at low collision energies, as well as the relative intensities of the three products with $C_2H_5^+$ under different conditions are discussed in more detail in Section 2.2.2.

### 2.2.2. Cross Sections as a Function of the Collision Energy

For each of the three products, data as a function of $E_{CM}$ have been collected at $E_{phot}$ = 11.3, 12.0, and 13.0 eV, in order to explore the combined effects of $E_{phot}$ and $E_{CM}$ on the reaction. Here, we present both the CS of the three channels as a function of the collision energy at the lowest photon energy of 11.3 eV (Figure 4) and a comparison of the CS as a function of the collision energy for the m/z 27 product at the three different photon energies (Figure 5).
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Cross sections (CSs) as a function of $E_{phot}$ for the m/z 28, 29, and 30 products of the reaction of $C_2D_5^+$ with $C_2H_2$ at $E_{CM}$ = 0.08 eV. Dashed lines represent the photon energies at which CSs have been measured as a function of the collision energy for the various product channels.</figcaption></figure>

The formation of the m/z 39 product ($[C_3H_5]^+$) is the major pathway at low collision energies (as observed in Figure 2), but its CS decreases with increasing collision energy, consistent with a pathway proceeding via formation of an adduct, with a subsequent increase with increasing collision energy that emerges at $E_{cm} > 2$ eV. This is a general trend since the collision energy shortens the precursor lifetime and drastically reduces the reaction efficiency. The inhibition of the low collision energy feature by increasing photon energy has already been shown in Figure 2, while the high collision energy feature is approximately independent of the photon energy, and so, a comparison of the collision energy dependence of this channel at different photon energies is limited to Figure S1 in the Supplementary Material. However, we briefly note here that, as can be seen through comparison with Figure S2 in the Supplementary Material, the high-collision energy rise in the m/z 39 signal is matched by a corresponding decrease in m/z 53 signal, indicating that this is a real effect and not an instrumental artefact.

The m/z 53 product ($[C_4H_5]^+$) exhibits a similar decrease with increasing collision energy to the m/z 39 product, again consistent with a complex-mediated pathway due to the reduction of precursor lifetime with increasing collision energy, as outlined in the previous paragraph. However, unlike the m/z 39 product, this is not accompanied by an increase in CS at high collision energies. As with the m/z 39 product, the inhibition of the low collision energy feature by increasing photon energy has already been shown in Figure 2, and so, a comparison of the collision energy dependence of this channel at different photon energies is limited to Figure S2 in the Supplementary Material.

The m/z 27 product ($C_2H_3^+$) shows a significant increase in CS with increasing collision energy, starting from around 0.2–0.3 eV. When we compare this trend at different photon energies (Figure 5), we note that not only does the measured CS increase with increasing photon energy over the entire collision energy range, but that this increase is especially pronounced at low collision energies. It is important to note that the internal energy variations of the $C_2H_5^+$ with changes in the photon energy are not easily quantifiable since, in the dissociative ionization process, a fraction of the photon energy is taken by both the ejected electron and the Br atom fragment. However, the observed trends with photon energy point to an increase in the amount of internal energy available to the $C_2H_5^+$ reagent ion.
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Cross sections as a function of *E*<sub>CM</sub> (in log–log scale) for the m/z 27, 39, and 53 products of the reaction of C<sub>2</sub>H<sub>5</sub><sup>+</sup> with C<sub>2</sub>H<sub>2</sub> at *E*<sub>phot</sub> = 11.3 eV. The dashed line represents the collision energy at which the corresponding scan as a function of the photon energy (Figure 2) was conducted.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Cross sections as a function of *E*<sub>CM</sub> (in log–log scale) for the m/z 27 product (C<sub>2</sub>H<sub>3</sub><sup>+</sup>) of the reaction of C<sub>2</sub>H<sub>5</sub><sup>+</sup> with C<sub>2</sub>H<sub>2</sub>, at *E*<sub>phot</sub> = 11.3, 12.0, and 13.0 eV. The dashed line represents the collision energy at which the corresponding scan as a function of the photon energy (Figure 2) was conducted.</figcaption></figure>
As with the photon energy dependence, we have collected data for the reaction of C₂D₅⁺ with C₂H₂ as a function of the collision energy for the *m/z* 28, 29, and 30 products, at both *E*<sub>phot</sub> = 11.3 and 12.0 eV, with the results shown in Figure 6. As already mentioned, the deuteron transfer product (C₂H₂D⁺) will appear at *m/z* 28, while the CID of the parent ion to give C₂D₃⁺ will be shifted to *m/z* 30. From Figure 6 it can be noted that, at *E*<sub>phot</sub> = 11.3 eV, the *m/z* 30 product is only observed at collision energies greater than ~ 1 eV, while the deuteron transfer at *m/z* 28 shows an appearance energy of about 0.3 eV that is not so dissimilar from the expected reaction endothermicity of Equation (8) (See Table 1). It is therefore possible to assign the observed *m/z* 27 CS in Figure 5 below a collision energy of ~1 eV, at *E*<sub>phot</sub> = 11.3 eV, to the endothermic proton transfer.

However, while a direct deuteron transfer mechanism can explain the observed *m/z* 28 signal, the *m/z* 29 product in Figure 6 can be rationalized only by assuming a contribution from H/D scrambling from a complex mediated mechanism, that would give C₂HD₂⁺ in addition to the other two channels:

$$C_2D_5^+ + C_2H_2 \rightarrow [C_4D_5H_2]^+ \rightarrow C_2H_2D^+ + C_2D_4 \quad \text{m/z 28} \quad (11)$$

$$\rightarrow C_2HD_2^+ + C_2HD_3 \quad \text{m/z 29} \quad (12)$$

$$\rightarrow C_2D_3^+ + C_2H_2D_2 \quad \text{m/z 30} \quad (13)$$

By focusing only on the data at low collision energies (i.e., below ~ 1 eV), the larger CS values observed at the larger photon energy *E*<sub>phot</sub> = 12.0 eV (right panel of Figure 6) compared to *E*<sub>phot</sub> = 11.3 eV (left panel of Figure 6) for the *m/z* 30, 29, and 28 products can therefore be rationalized as due to the direct and complex mediated deuteron transfer reactions (11)–(13), the probability of which increases due to the larger amount of internal energy available to the reacting cation. The red dashed line is a smoothed sum of the CSs for the three channels, to be compared with the collision energy trend for the *m/z* 27 product of Figure 5 measured at the same *E*<sub>phot</sub> = 12.0 eV: the similarity in absolute value and trend of the CSs is an indication that the interpretation in terms of scrambling in complex mediated deuteron transfer is feasible.

## 2.3. MIKE Spectra

In order to probe the metastable decomposition from the [C₄H₇]⁺ PES, MIKE spectra have been recorded for CH₂C(CH₃)CH₂⁺ ions generated via dissociative ionization of 2-methylpropene (CH₂C(CH₃)₂), with the results shown in Figure 7. Full KER distributions of the *m/z* 39 and 53 products have been determined from the metastable peak shapes using the META algorithm described in [78].

The peak at *m/z* 53, corresponding to the loss of H₂, has a dish-topped shape characteristic of a process involving a large kinetic energy release (KER). Contrastingly, the *m/z* 39 peak, corresponding to the loss of CH₄, has a markedly sharper profile, characteristic of a process that not only involves a lower KER but also a narrower KER distribution. The obtained KER distributions of the *m/z* 39 and 53, products are given in Figures S3 and S4, respectively, in the Supplementary Material. The peak at *m/z* 29 corresponds to the ejection of C₂H₅⁺, which is the entrance channel in the title reaction.

An equivalent MIKE spectrum of [C₄D₇]⁺ has also been recorded, in this case starting from d₈-2-methylpropene (CD₂C(CD₃)₂), where, despite the shifting of peak masses due to deuteration, very similar peak shapes are observed. This is shown in Figure S5 in the Supplementary Material, with the corresponding KER distributions for the *m/z* 42 and 58 fragments given by Figures S6 and S7, respectively. A further discussion of these results is given in Section 4.
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Cross sections (CSs) as a function of $E_{CM}$ (in log-log scale) for the m/z 28, 29, and 30 products of the reaction of $C_2D_5^+$ with $C_2H_2$ at $E_{phot} = 11.3$ eV (left panel) and 12.0 eV (right panel). The red dashed line is a smoothed sum of the CSs for the three channels. The black dashed lines represent the collision energy at which the corresponding scans as a function of the photon energy (Figure 3) were conducted.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 7.** Mass-analyzed ion kinetic energy (MIKE) spectrum of the m/z 55 ion formed via dissociative ionization of 2-methylpropene. Insets show shapes of the m/z 39 and 53 peaks, while the peak at m/z 54 is an artefact arising from leaking of m/z 56 ions into the precursor ion signal.</figcaption></figure>
### 3. Computational Results

In order to rationalize the observed products and trends as a function of both $E_{phot}$ and $E_{CM}$, the PES of the reaction has been explored through computational calculations, as described in Section 5.2. The PES was initially mapped at the ωB97X-D/cc-pVTZ level of theory, but a comparison of the relative energies obtained with those calculated from literature formation enthalpies (as shown in Table 1) indicated that a higher level of theory is required. The PES shown here is, therefore, that calculated at the G4 level of theory.

Results of the calculations are in terms of the stationary points along the reaction pathways and are shown in Figures 8 and 9, where the zero energy is given by the sum of the energies of the reactants, and the relative energies for species (reported in parentheses in the following) are in eV. Structural diagrams and atomic coordinates of the intermediates and transition states are given in Figures S9–S45 and Tables S4–S40 of the Supplementary Material. An overview of the reaction pathways and rearrangements of the intermediates is summarized in Scheme 1.

<figure><img src="image_9.png" /><figcaption>**Scheme 1.** A sketch of the rearrangements and reaction pathways for the reaction of $C_2H_5^+$ with $C_2H_2$ leading to product channels (4), (5), (6), and (7). The lettered labels of the different structures (e.g., **W1**) and transition states (e.g., **M1-W2**) correspond to those described in the text and in Figures 8 and 9).</figcaption></figure>

For all products at low $E_{CM}$, the reaction proceeds via the entrance channel ion-induced dipole complex **W1** (−0.44), from which the covalently-bound adduct $CH_3CH_2CHCH^+$ (−1.16), hereafter referred to as **M1**, can form via the transition state **W1-M1** (−0.28). Though endothermic, the *m/z* 27 product $C_2H_3^+$ can form in combination with $C_2H_4$ (0.36) either directly from the reagents or via **W1**. In both cases, this proceeds without a barrier.

As already mentioned and shown in Table 1, for the *m/z* 39 product, both the cyclic (c-$C_3H_3^+$) and linear ($H_2CCCH^+$) isomers are energetically allowed, with the former being more stable by ~1.2 eV, see reactions (4) and (5). Similarly, for the *m/z* 53 product, both the cyclic (c-CHCHC($CH_3^+$)) and linear ($CH_2CHCCH_2^+$)) isomers are energetically allowed, with the former being more stable by ~0.6–0.7 eV, see reactions (6) and (7).
<figure><img src="image_10.png" /><figcaption>**Figure 8.** The potential energy surface for the reaction of C₂H₅⁺ with C₂H₂: the most relevant reaction pathways, energetically accessible in the experiments and including intermediates and transition states (TSs), leading to the formation of *m/z* 27, 39, and 53 products are reported. Energies are computed, as detailed in the text, at the G4 level of theory. The zero of the energy scale is given by the sum of the reagent energies. For clarity reasons, only some of the intermediates' structures are reported. All the structures for intermediates and TSs are available in the Supplementary Material.</figcaption></figure>

From **M1**, the lower energy c-C₃H₃⁺ isomer of [C₃H₃]⁺ can form following a simultaneous [1,3] cyclization and [1,2] H-shift to give **W2** (−1.47), an ion-induced dipole complex of c-C₃H₃⁺ and CH₄, via **M1-W2** (−1.06). **W2** is then able to fragment into the separated products (−1.36) without a barrier. IRC calculations on the transformation from **M1** into products have been performed and are presented in the Supplementary Material (as animated gif file *TS M1-W2.gif*).

The lower energy c-C₃H₂(CH₃)⁺ isomer of [C₄H₅]⁺ can form from **M1** via an initial [1,3] H-shift to give the *anti*-conformer of CH₃CHCHCH₂⁺ (−2.77), hereafter referred to as **M2**, via **M1-2** (−0.96). **M2** can then undergo a [1,2] H-shift to give CH₃CHCCH₃⁺ (−2.05), hereafter referred to as **M3**, via **M2-3** (−1.27). **M3** can then undergo a simultaneous [1,3] cyclization and H₂-ejection via **M3-W3** (−0.11) to give **W3** (−1.58), an ion-induced dipole complex of c-C₃H₂(CH₃)⁺ and H₂ that can fragment into the separated products (−1.57) without a barrier. IRC calculations on the transformation from **M3** into products have been performed and are presented in the Supplementary Material (as animated gif file *TS M3-W3.gif*): the methyl group migrates to the carbene C atom, thus closing the three-member ring, and, as a consequence of such closure, two H atoms are eliminated from the migrating methyl group in the form of an H₂ molecule.

The higher energy *m/z* 39 and 53 isomers, H₂CCCH⁺ and CH₂CHCCH₂⁺, can both form from **M2** via an initial interconversion into the higher energy *syn*-conformer of CH₃CHCHCH₂⁺ (−2.63), hereafter referred to as **M4**, via **M2-4** (−1.84). It should be noted that **M4** can also be accessed from **M3** via a [1,2] H-shift over **M3-4** (−1.29). **M4** can then undergo a [1,2] H-shift between the two central C atoms to give CH₃CH₂CCH₂⁺ (−1.84), hereafter referred to as **M5**, via **M4-5** (−1.63). **M5** can then undergo a simultaneous C-C
cleavage and [1,2] H-shift between the CH₂ and CH₃ moieties via **M5-W4** (−0.26) to give **W4** (−0.31), an ion-induced dipole complex of CH₂CCH⁺ and CH₄. **W4** can then fragment into the separated products (−0.19) without a barrier. IRC calculations on the transformation from **M5** into such products have been performed and are presented in the Supplementary Material (as animated gif file TS M5-W4.gif).

Alternatively, **M5** can undergo a [1,2] H₂-ejection via **M5-W5** (−0.11) to give **W5** (−0.93), a weakly bound adduct of CH₂CHCCH₂⁺ and H₂, that can again fragment into the separated products (−0.92) without a barrier. IRC calculations on the transformation from **M5** into such products have been performed and are presented in the Supplementary Material (as animated gif file TS M5-W5.gif).

Finally, when sufficient internal energy is available, the endothermic proton transfer channel to give C₂H₃⁺ (*m/z* 27) and C₂H₄ (+0.36) can proceed either directly from the reactants or from the initial **W1** complex. All of these reaction pathways are shown graphically in Figure 8. As discussed earlier, at higher collision energies, there is also the possibility for collision-induced dissociation (CID) of the parent ion into C₂H₃⁺ and H₂, with the calculated reaction enthalpy for this process (Table 1) showing a good level of agreement, at both levels of theory, with the literature value.

<figure><img src="image_11.png" /><figcaption>**Figure 9.** Structure and energies of the intermediates and transition states (TSs) relevant to the inter-conversion of [C₄H₇]⁺ adducts formed by the reaction of C₂H₅⁺ with C₂H₂. Energies are computed, as detailed in the text, at the G4 level of theory. The zero of the energy scale is given by the sum of the reagent energies. For clarity reasons, only some of the intermediate structures are reported. All the structures, for intermediates and TSs are available in the Supplementary Material.</figcaption></figure>

In addition to calculations on the pathways leading to products, we have also probed the [C₄H₇]⁺ PES more widely in order to explore minima and transition states relevant to H/D scrambling, with the results shown graphically in Figure 9, where the zero energy is again given by the sum of the energies of the reactants.

Starting from **M2**, the first rearrangement to note is a [1,4] H-shift along the C atom chain that produces a mirrored **M2** intermediate via **M2-2** (−1.24). A similar [1,4] H-shift is also possible from **M4** via **M4-4** (−1.11). From **M3**, there is the potential for the CH hydrogen to adopt a bridged position between the two central C atoms via **M3-6** (−2.11).
to give **M6** (−2.11). From **M4**, a partial delocalization of one of the methyl hydrogens between the two terminal C atoms to give the strained **M7** structure (−1.26) is possible via **M4-7** (−1.18).

Two further interconversion pathways proceed from **M5**, with the lowest energy being a [1,2] CH₃-shift to give CH₂C(CH₃)CH₂⁺ (−2.35), hereafter referred to as **M8**, via **M5-8** (−1.33). The resultant structure is then able to isomerize into itself via a [1,3] H-shift over **M8-8** (−1.43). **M5** is also able to undergo a simultaneous cyclization and [1,3] H-shift via **M5-9** (−1.20) to give c-CH₂CH₂CH(CH₂)⁺ (−2.37), hereafter referred to as **M9**.

From **M9** two four-membered ring structures, **M10** (−1.72) and **M11** (−1.71), can be formed. The C<sub>s</sub> symmetric **M11** is, in fact, a π-protonated cyclobutene, which easily isomerizes by H-migration via **M10-11** (−1.71) to the C<sub>s</sub>-symmetric cyclobutylium ion **M10**. The ring in **M10** is slightly puckered with a C(H)-C(H<sub>2</sub>)-C(H<sub>2</sub>)-C(H<sub>2</sub>) dihedral angle of 6.8°. A slight increase of this angle while maintaining C<sub>s</sub> symmetry leads to **M9-10** (−1.72) from which the system goes downhill, reaches the valley-ridge inflection point, and proceeds further to another TS. This type of PES describes a reaction mechanism that has been referred to as a two-step-no-intermediate [79]. The latter TS is also a barrier for isomerization of **M9** into itself, i.e., **M9-9** (−2.38). From **M10**, planarizing the ring leads to a C<sub>2v</sub> symmetric structure representing the ring inversion transition state **M10-10** (−1.72).

Though the scrambling is discussed in more detail in Section 4 and in Section 5 of the Supplementary Material, we draw attention to the fact that almost all the additional minima identified in Figure 9 are symmetric, allowing for efficient exchange of H atoms.

## 4. Discussion

From the calculations presented in Section 3, it is possible to rationalize the experimental CSs and their trends with collision energy for the various products.

### 4.1. Product at m/z 27: C₂H₃⁺

As already mentioned, two different pathways are possible for the formation of C₂H₃⁺ ions (m/z 27): the endothermic proton transfer channel (8) and the CID of the reagent ion (9). Through the CSs measured for the reactivity of C₂D₅⁺ with C₂H₂ (see Figure 3 and 6), we have been able to assess that both pathways are operative in different collision energy regimes. In particular, the production of C₂H₃⁺ at high collision energies (larger than ~1 eV) is assigned to the CID of the parent ion, while the product observed at lower energy is the result of the endothermic proton transfer mediated by the electrostatically bound ion-molecule complex W1.

### 4.2. Products at m/z 39 ([C₃H₃]⁺) and m/z 53 ([C₄H₅I]⁺)

For the [C₃H₃]⁺ (m/z 39) and [C₄H₅]⁺ (m/z 53) products, calculations indicate that they both proceed via the **M1** covalently-bound intermediate formed from **W1**. The decrease in CS with increasing *E*<sub>phot</sub> for both channels (see Figure 2) can be rationalized by a combination of two factors: (i) the increased competition due to the emergence of the endothermic channels (8) and (9) giving C₂H₃⁺ (m/z 27) and (ii) the inhibition of adduct formation when the internal energy content of the reagent ion increases.

#### 4.2.1. Rationalization of m/z 39 and 53 Branching Ratios

To rationalize the BRs for the m/z 39 and 53 products, several aspects should be considered. From a thermochemical point of view, the calculated energies at the G4 level (see Table 1) indicate that the production of the lowest energy isomer c-C₃H₂(CH₃)⁺ (m/z 53) plus H₂, reaction (6), with ΔH°₂₉₸ = -1.57 eV, is more exothermic than the CH₄ elimination to give either c-C₃H₃⁺ via reaction (4) (by 0.21 eV), or l-C₃H₃⁺ via reaction (5) (by 1.38 eV). On the other hand, the production of the higher energy isomer CH₂CHCCH₂⁺ (m/z 53) plus H₂, reaction (7) with an exothermicity ΔH°₂₉₸ = -0.92 eV, is thermodynamically favoured over reaction (5) but not over reaction (4).
From a kinetic perspective, we must consider both the relative heights of the transition states and the number of rearrangements required prior to fragmentation. As all four relevant product channels stem from the initial adduct **M1**, we only need to consider the kinetics starting from this point.

For the *m/z* 39 products, the formation of c-C₃H₃⁺ proceeds via a single interconversion over **M1-W2**, which has a relative energy of -1.06 eV. In contrast, the formation of CH₂CCH⁺ proceeds via four interconversions over **M1-2** (-0.96), **M2-4** (-1.84), **M4-5** (-1.63), and **M5-W4** (-0.26). The formation of the cyclic isomer is therefore not only more exothermic than the formation of the linear isomer, but is also significantly more favored kinetically.

For the *m/z* 53 products, both isomers proceed following the interconversion to **M2** via **M1-2**, and so the relative kinetic feasibility can be considered from this point. For the c-C₃H₂(CH₃)⁺ isomer, the formation can then proceed following an isomerization to **M3** via **M2-3** (-1.27) and rearrangement via **M3-W3** (-0.11). For the CH₂CHCCH₂⁺ isomer, however, formation requires three successive rearrangements via **M2-4**, **M4-5**, and **M5-W5** (-0.11). While the maximum barrier heights for these two processes are the same, the much greater length of the formation pathway for the CH₂CHCCH₂⁺ isomer means that this pathway is much less favorable from a kinetic perspective.

### 4.2.2. Comparison of Computed Energy Pathways with MIKE Spectra

These conclusions are supported by the MIKE spectra collected as part of this work. Starting from **M8** (see Figure 9), the [C₄H₇]⁺ ions formed via the dissociative ionization of the 2-methylpropane can then isomerize to give **M5**, from which all five fragmentation pathways (including that to give C₂H₅⁺ plus C₂H₂) are then accessible (see Figure 8). From the broad shape of the *m/z* 53 peak (Figure 7), we are able to conclude that the [C₄H₅]⁺ fragment ions are formed with a large KER, consistent with fragmentation products that are much lower in energy than the adjacent transition state. This is consistent with the formation of the c-C₃H₂(CH₃)⁺ isomer, where the energy of the products (-1.57) is significantly lower than the transition state, **M3-W3** (-0.11), leading to the fragmentation. Should the major formation pathway be that to give the linear CH₂CHCCH₂⁺, we would instead observe a narrow distribution consistent with the higher energy of the products (-0.92) relative to the barrier to fragmentation via **M5-W5** (-0.11).

By contrast, for the *m/z* 39 peak, the sharp distribution (Figure 7) is consistent with a fragmentation pathway where the energy of the products lies only slightly lower than the barrier to fragmentation. Again comparing this to the PES, we note that it is consistent with the lower-energy pathway to form c-C₃H₃⁺, where the products (-1.36) are only 0.3 eV lower than the barrier to fragmentation (**M1-W2**, -1.06), with the formation of the CH₂CCH⁺ isomer involving products which are higher in energy than the barrier to fragmentation.

The total KER distributions obtained from the *m/z* 39 and 53 peaks of the MIKE spectra are reported in Figures S3 and S4 of the Supplementary Material, respectively, and show a clear difference in the distributions for the two channels. For the H₂ elimination channel, there is a tail that extends up to ~ 1.1 eV, while that for CH₄ elimination is much smaller, extending only up to ~ 0.30 eV. The KER distributions for the equivalent *m/z* 42 (C₃D₃⁺) and 57 (C₄D₅⁺) peaks in the MIKE spectrum of the deuterated analog (see Figure S5 of the Supplementary Material) are reported in Figures S6 and S7 of the Supplementary Material, with the high level of agreement between the hydrogenated and deuterated samples indicating that there is little to no contribution to the fragmentation from tunneling. While a comparison of the KER distributions with the calculated energies indicates that a fraction of the available energy is transferred into internal excitation of the fragments, a quantification of the effect is beyond the scope of the work.

Given the observed peak profiles from the MIKE spectra and the calculated PES, we are therefore able to conclude that the *m/z* 39 and 53 products correspond to c-C₃H₃⁺ and c-C₃H₂(CH₃)⁺, respectively.
When focusing on the kinetics of the formation of c-C₃H₃⁺ and c-C₃H₂(CH₃)⁺, i.e., the predominant isomer for the two main mass channels, we note that the formation of c-C₃H₃⁺ (*m/z* 39) not only involves a much lower reaction barrier than the formation of c-C₃H₂(CH₃)⁺ (*m/z* 53), but also requires one fewer rearrangement prior to fragmentation. Therefore, while the formation of the *m/z* 53 (c-C₃H₂(CH₃)⁺) product is the most exothermic pathway, the calculated PES is consistent with the experimental observation of the *m/z* 39 product, c-C₃H₃⁺, as the dominant pathway.

A confirmation of the proposed reaction pathways can be obtained from the data collected on the reactivity of C₂D₅⁺ with C₂H₂, by analyzing the extent of the H/D scrambling patterns in the mass spectra of [C₃(H, D)₃]⁺ and [C₄(H, D)₅]⁺ products. To avoid lengthening the main paper, this discussion is reported in Section 5 of the Supplementary Material.

### 4.3. *Comparison with Previous Results*

A direct comparison between our results and those in [39] should be made with care since reaction rates and BRs obtained using different mass spectrometric techniques may show some discrepancies, as previously outlined among Fourier Transform-ICR, selected ion flow drift tube methods and ion trap mass spectrometry [80–82]. One potential source of discrepancy is the different ranges of ion kinetic energies explored using different techniques. While in drift tube set-ups the ions have thermal translational and rotational energy distributions due to collisions with the buffer gas, a critical issue with early ICR experiments was the determination of the kinetic energy of the ions entering the ICR cell [83]. Due to the higher potentials at the trapping plates than at the center of the ICR cell, ions formed at or near the trapping plates were translationally excited with velocities characteristic of temperatures above thermal, but a quantitative assessment of the reactant ion kinetic energy was not always possible. In our guided ion beam technique, the kinetic energy of the ions entering the scattering cell is well defined and measurements are performed as a function of collision energy.

Another source of discrepancy may be related to an unknown amount of internal energy in the reactant ions as, in ICR experiments, ions can have different internal energies according to the ion generation and isolation method. In the previous study from literature [39], the C₂H₅⁺ ions were generated by electron ionization of C₂H₅Cl and C₂H₆ with a nominal electron energy of 16 eV, but inefficient internal relaxation prior to the reaction may have led to some portion of the reactant ions being generated in vibrationally excited states. In our experiment, the use of dissociative photoionization at well-defined photon energies allows for the control of internal excitation. As shown in Figure 2, BRs of the reaction products are strongly affected by the photon energy, it is therefore possible that the different BRs measured in that previous work [39] compared to this one are due to the presence of different amounts of internally excited C₂H₅⁺ ions. However, regardless of the source of the deviation between this work and previous literature, the calculated PES detailed in Section 3 and the MIKE spectra presented in Section 2.3, allow us to fully rationalize our experimental BRs.

## 5. Materials and Methods

### 5.1. *Experimental Set-Up*

The experimental set-up used for reactivity measurements has been described in detail elsewhere [55,56,84], and so, only a brief description is given here. The data on the reactivity of C₂H₅⁺ has been collected using the CERISES apparatus [85,86] in combination with the DESIRS beamline [87] of the SOLEIL synchrotron. CERISES is a guided ion beam tandem mass spectrometer consisting of two octopoles located between two quadrupole mass filters. This allows for the mass selection of reagent and product ions of ion-neutral reactions, with the neutral reagent being introduced in a scattering cell surrounding the final part of the first octopole.

The C₂H₅⁺/C₂D₅⁺ reagent ions used in this work were generated via VUV photoionization and fragmentation of C₂H₅Br (Sigma-Aldrich, 98% purity)/C₂D₅Br (Sigma-Aldrich, 99%
purity with 99% deuteration), introduced into the source at approximately 1.8 × 10⁻⁶ mbar. The wide tuneability and intensity of the DESIRS VUV photon source allows for data to be collected as a function of the photon energy which, in turn, acts as a proxy for the internal energy of the reagent ions. The photon energies used were in the range of 10.4–14 eV, with a resolution of 20–40 meV defined by the monochromator slit setting. Photons of energies greater than 15.7 eV were removed by an Argon gas filter [88]. Photon energies in the absolute scale were obtained using the absorption lines of Argon around 11.823 and 14.304 eV [64,89], with systematic shifts of 10–20 meV above the tabulated values.

Neutral reagent (C₂H₂), Air Liquide, 99.7% purity) pressures used were of the order of 1 × 10⁻⁵ mbar, in order to ensure operation close to the single collision regime. In this way, we were able to reduce the contribution from secondary collisions and limit the attenuation of the parent beam to under 10%. As the C₂H₂ was dissolved in acetone (CH₃C(O)CH₃) as a stabilizer, the sample was first passed through a cryogenic trap containing dry ice in order to remove any acetone present before introduction into the experiment. Absolute pressures were measured using a MKS 398H differential manometer.

The collision energy available to the reactants depends on both the ionic charge (in this case +1) and the potential difference between the ion source and reaction cell. The retarding potential method [90] has been used to determine the maximum of the first derivative of the parent ion yield, the corresponding voltage of which defines the zero of the kinetic energy in the laboratory frame, which can then be converted into center-of-mass collision energies (*E*<sub>CM</sub>). By changing the potentials of the reaction cell and all subsequent elements on the path of ions, we can scan a collision energy range from ≈ 0.07 to ≈ 6 eV in the center of mass frame. The FWHM of the collision energy in the scanned range is 0.16 eV in the *E*<sub>CM</sub> frame.

Data as a function of both photon and collision energies have been collected in the "multi-scan" mode where the signals for all ionic species of interest are collected at a given set of conditions before moving to the next one. In this way, we are able to drastically reduce any potential effects from drifts in source or reaction cell pressure. Absolute CS data rely on a calibration of the effective length of the reaction cell performed using the ArD⁺ product channel of the well-known reaction between Ar⁺ and D₂ measured with a similar guided ion beam tandem mass spectrometer [74]. Cross sections are typically given with 10% errors resulting from the uncertainty in the absolute pressure measurement. However, for experimental conditions where ionic counting rates are small (close to appearance energy and/or channels with very small branching ratio), the uncertainty can exceed 10%. In such cases, a propagation of uncertainty is performed for all counting rates involved in the cross section measurement, applying a square root function on each counting measurement to derive its uncertainty.

Mass-analyzed ion kinetic energy (MIKE) spectra have been recorded using a ZAB2-SEQ tandem mass spectrometer of reverse double-focusing geometry, in which the magnetic sector precedes the electrostatic sector. Precursor ions (CH₂C(CH₃)CH₂⁺/CD₂C(CD₃)CH₂⁺) are prepared by dissociative electron ionization (70 eV, 150°C) of CH₂C(CH₃)₂ (Merck, 99.3% purity) and CD₂C(CD₃)₂ (Merck, 99% purity and 99% deuteration), respectively. Following acceleration of 8 keV, precursor ions are mass-selected by the magnetic sector and allowed to dissociate spontaneously in the second field-free region of the instrument. The resultant ions are then analyzed by scanning the electrostatic sector. For peak profile measurements, the main beam width is reduced to 4 eV, corresponding to an energy resolution of 2000. For the *m/z* 55 ions (CH₂C(CH₃)CH₂⁺), the metastable time window for ion fragmentation under these conditions is 9.4–18.3 µs, with the *m/z* 62 ions (CD₂C(CD₃)CD₂⁺) having a metastable time window of 10.0–19.4 µs.

## 5.2. Computational Methodology

Electronic structure calculations to identify the underlying reaction mechanisms (i.e., stable intermediates and transition states) and aid in the interpretation of experimental results were performed with the Gaussian 16 suite of programs [71]. The ωB97X-D/cc-pVTZ
level of theory, where ωB97X-D is a long-range corrected hybrid density functional that includes empirical atom-atom dispersion corrections [91] and cc-pVTZ refers to the correlation consistent polarized Valence Triple Zeta basis set (i.e., a triple-zeta basis set intended specifically for use in correlated calculations) was initially used for a preliminary mapping of the most relevant stationary points of the PES. Intrinsic reaction coordinate (IRC) calculations have then been performed to ensure that each TS connects to the relevant minima.

The accuracy of the relative energies obtained at the ωB97X-D/cc-pVTZ level of theory was tested against available literature values for the reaction enthalpies of the various product channels (as shown in Table 1 and in Table S1 of the Supplementary Material). Further accuracy tests were performed by calculating proton affinities of various organic molecules (ethyne, ethene, propene, butadiene, cyclobutene), as reported in Table S2 of the Supplementary Material. From such a comparison, it was concluded that the accuracy of the relative energies at the ωB97X-D/cc-pVTZ level of theory is not satisfactory, and thus, the relative energies were all recalculated using a composite method denoted as G4 [92] to obtain more accurate energies, and only the latter results are presented in Section 3.

Reaction enthalpies (ΔH°₂₉₸) have also been calculated for comparison with literature values using both ωB97X-D/cc-pVTZ and G4 levels of theory.

**6. Conclusions**

Absolute cross sections and branching ratios for the reaction of the ethyl cation with ethyne have been measured as a function of collision energies, in the hyperthermal energy range (from ≈ 0.07 to ≈ 5 eV in the CM frame), under a single collision regime using a guided ion beam tandem mass spectrometer. The C₂H₅⁺ reactant ion is generated via dissociative ionization of C₂H₅Br using tunable VUV radiation in the energy range 10.4–14.0 eV. The two main product channels are [C₃H₃]⁺ plus CH₄ (BR = 0.76 ± 0.05) and [C₄H₅]⁺ plus H₂ (BR = 0.22 ± 0.02), where the reported values have been determined at a photon energy Eₚhöt = 11.3 eV and a collision energy in the CM frame ECM = 0.08 eV. At a higher internal excitation of C₂H₅⁺ and/or high collision energy, we observed an additional endothermic channel giving C₂H₃⁺ in association with either C₂H₄ or C₂H₂ plus H₂.

Ab initio calculations of the most relevant stationary points on the potential energy surfaces allow for the rationalization of experimental results. Both channels stem from a common covalently bound intermediate, CH₃CH₂CHCH⁺, from which barrierless and exothermic pathways exist for the production of both cyclic and linear isomers of both main channels. From a comparison between the MIKE spectra of C(CH₃)(CH₂)₂⁺ and C(CD₃)(CD₂)₂⁺ and the calculated PES, we are able to conclude that, for both mass channels, the predominant isomer is the lower energy cyclic one, and that the contribution from the tunneling mechanisms is minimal.

However, the main result of this study is the re-evaluation of the product BRs for the two main reaction channels. Present experimental values, supported by theoretical calculations, indicate that c-C₃H₃⁺ plus CH₄ channel is the dominant one, at odds with previous literature data [39] where [C₄H₅]⁺ plus H₂ was indicated as the primary product channel. Though the formation of c-C₃H₂(CH₃)⁺ plus H₂ is the most exothermic pathway, the barrier to fragmentation (−0.11 eV) is significantly higher than that required for c-C₃H₃⁺ formation (−1.06 eV), with the additional rearrangement required prior to fragmentation further enhancing the relative kinetic inhibition for this channel.

To assess the implications of our study on the modeling of Titan's atmosphere, a sensitivity test was performed using the Titan photochemical model of V. Vuitton and collaborators [6] to estimate the changes in the density of ions and neutrals when the BRs for the title reaction are replaced with the present values. The photochemical model was run using the two different sets of BRs, as reported in Table 2, with the results summarized as follows. In each case, the densities that are reported are given as an integration over the atmospheric column, peaking around 1000 km and negligible below 700 km, as illustrated in Figures 39–42 of the cited work [6].
* The density of $[C_3H_3]^+$ ions is increased by ~30% when our new BRs are used, despite the fact that the main formation reaction for $[C_3H_3]^+$ in the model is $C_2H_4^+ + C_2H_2 \rightarrow [C_3H_3]^+ + CH_3$ and not the title reaction. This is a relevant difference that should be considered in light of the fact that a change in the abundance of $[C_3H_3]^+$ can induce changes in other species, such as the recently detected $c-C_3H_2$ [43].

* The density of the $[C_4H_5]^+$ ion decreases by about a factor of 2 when our new BRs are used. The repercussion on the production of $[C_6H_7]^+$ is lower (with a maximum change of ~20%), as the main formation reaction for this ion is $[C_3H_5]^+ + CH_3CCH \rightarrow [C_6H_7]^+ + H_2$ rather than $[C_4H_5]^+$ plus $C_2H_4$. However, it still marks a significant change to the predicted density of this key astrochemical species in Titan's environment.

Beyond the atmosphere of Titan, present results have implications for the chemical evolution of a range of different astrochemical environments due to the near-ubiquitous observation of both reactants. Notably, while the radioastronomical detection of the symmetric ion $c-C_3H_3^+$ is difficult, its linear isomer $H_2CCCH^+$ has been very recently detected in a cold dense molecular cloud core [93], while both $c-C_3H_2$ and $H_2CCC$ isomers have long been detected in the ISM [94,95]. Finally, the observation of a third product, $C_2H_3^+$, at high internal and collision energies, has implications for the modeling of high-energy plasmas.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29040810/s1, Supplementary Material pdf file containing: Figures S1–S45 and Tables S1–S40; animated gif files: TS M1-W2.gif, TS M3-W3.gif, TS M5-W4.gif, TS M5-W5.gif 18 [58,69,70,96–98].

**Author Contributions:** Conceptualization: D.A., V.R., M.P.; data curation: R.T., M.P., C.A.; formal analysis: V.R., R.T., D.A., C.A.; funding acquisition: D.A., P.T.; investigation: C.R., V.R., D.A., C.A., R.T., J.Ž.; methodology: D.A., R.T., C.R., C.A., M.P., P.T., J.Ž.; project administration: D.A., R.T., C.A.; resources: C.R., C.A., R.T., J.Ž.; software: R.T., M.P., P.T.; supervision: D.A., R.T., C.R., C.A.; visualization: V.R., D.A.; writing—original draft: V.R. and D.A.; all authors contributed to writing, review, and editing. All authors have read and agreed to the published version of the manuscript.

**Funding:** This project has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska Curie grant agreement No 811312 for the project "Astro-Chemical Origins" (ACO) and from MUR PRIN 2020 project n. 2020AFB3FX "Astrochemistry beyond the second period elements". MP and JZ acknowledge support from the Ministry of Education, Youth and Sports of the Czech Republic (grant No. LTC20062) and the Czech Science Foundation (grant No. 21-11931J). VR acknowledges funding for a PhD fellowship from the Department of Physics, University of Trento. This article is based upon work from COST Action CA18212—Molecular Dynamics in the GAS phase (MD-GAS), supported by COST (European Cooperation in Science and Technology).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data are contained within the article and Supplementary Materials.

**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 20180118 and 20190249. We thank Veronique Vuitton for running the Titan photochemical model.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Larsson, M.; Geppert, W.D.; Nyman, G. Ion chemistry in space. *Rep. Prog. Phys.* **2012**, *75*, 066901. [CrossRef]

2. Agúndez, M.; Wakelam, V. Chemistry of Dark Clouds: Databases, Networks, and Models. *Chem. Rev.* **2013**, *113*, 8710–8737. [CrossRef]

3. Vigren, E.; Semaniak, J.; Hamberg, M.; Zhaunerchyk, V.; Kaminska, M.; Thomas, R.D.; af Ugglas, M.; Larsson, M.; Geppert, W.D. Dissociative recombination of nitrile ions with implications for Titan's upper atmosphere. *Planet. Space Sci.* **2012**, *60*, 102–106. [CrossRef]
4. Yuen, C.H.; Ayouz, M.A.; Balucani, N.; Ceccarelli, C.; Schneider, I.F.; Kokoouline, V. Dissociative recombination of CH₂NH₂⁺: A crucial link with interstellar methanimine and Titan ammonia. *Mon. Notices Royal Astron. Soc.* **2019**, *484*, 659–664. [CrossRef]

5. Hörst, S.M. Titan's atmosphere and climate. *J. Geophys. Res. Planets* **2017**, *122*, 432–482. [CrossRef]

6. Vuitton, V.; Yelle, R.; Klippenstein, S.; Hörst, S.; Lavvas, P. Simulating the density of organic species in the atmosphere of Titan with a coupled ion-entral photochemical model. *ICARUS* **2019**, *324*, 120–197. [CrossRef]

7. Imanaka, H.; Smith, M.A. Role of photoionization in the formation of complex organic molecules in Titan's upper atmosphere. *Geophys. Res. Lett.* **2007**, *34*. [CrossRef]

8. Coates, A.J.; Wellbrock, A.; Waite, J.H.; Jones, G.H. A new upper limit to the field-aligned potential near Titan. *Geophys. Res. Lett.* **2015**, *42*, 4676–4684. [CrossRef] [PubMed]

9. Westlake, J.; Waiter, J.; Carrasco, N.; Richard, M.; Cravens, T. The role of ion-molecule reactions in the growth of heavy ions in Titan's ionosphere. *J. Geophys. Res. Space Phys.* **2014**, *119*, 5951–5963. [CrossRef]

10. Vinatier, S.; Bézard, B.; Fouchet, T.; Teanby, N.A.; de Kok, R.; Irwin, P.G.; Conrath, B.J.; Nixon, C.A.; Romani, P.N.; Flasar, F.M.; et al. Vertical abundance profiles of hydrocarbons in Titan's atmosphere at 15S and 80N retrieved from Cassini/CIRS spectra. *ICARUS* **2007**, *188*, 120–138. [CrossRef]

11. Nixon, C.; Lorenz, R.; Achterberg, R.; Buch, A.; Coll, P.; Clark, R.; Courtin, R.; Hayes, A.; Iess, L.; Johnson, R.; et al. Titan's cold case files-Oustanding questions after Cassini-Huygens. *Planet. Space Sci.* **2018**, *155*, 50–72. [CrossRef]

12. Sittler, Jr., E.C.; Cooper, J.F.; Sturner, S.J.; Ali, A. Titan's ionospheric chemistry, fullerenes, oxygen, galactic cosmic rays and the formation of exobiological molecules on and within its surfaces and lakes. *ICARUS* **2020**, *344*, 113246. [CrossRef]

13. Waite, J.; Niemann, H.; Yelle, R.; Kasprzak, W.; Cravens, T.; Luhmann, J.; McNutt, R.; Ip, W.; Gell, D.; De La Haye, V.; et al. Ion Neutral Mass Spectrometer results from the first flyby of Titan. *Science* **2005**, *308*, 982–986. [CrossRef] [PubMed]

14. Vuitton, V.; Yelle, R.; McEwan, M. Ion chemistry and N-containing molecules in Titan's upper atmosphere. *ICARUS* **2007**, *191*, 722–742. [CrossRef]

15. Vuitton, V.; Yelle, R.V.; Lavvas, P. Composition and chemistry of Titan's thermosphere and ionosphere. *Philos. Trans. Royal Soc. A-Math. Phys. & Engin. Sci.* **2009**, *367*, 729–741. [CrossRef]

16. Cui, J.; Yelle, R.; Vuitton, V.; Waite, J.H. nad Kasprzak, W.; Gell, D.; Niemann, H.; Müller-Wodarg, I.; Borgeeren, N.; Fletcher, G.; Patrick, E.; et al. Analysis of Titan's neutral upper atmosphere from Cassini Ion Neutral Mass Spectrometer measurements. *ICARUS* **2009**, *200*, 581–615. [CrossRef]

17. Ali, A.; Sittler Jr., E.C.; Chornay, D.; Rowe, B.R.; Puzzarini, C. Organic chemistry in Titan's upper atmosphere and its astrobiological consequences: I. Views towards Cassini plasma spectrometer (CAPS) and ion neutral mass spectrometer (INMS) experiments in space. *Planet. Space Sci.* **2015**, *109*, 46–63. [CrossRef]

18. Spilker, L. Cassini-Huygens' exploration of the Saturn system: 13 years of discovery. *Science* **2019**, *364*, 1046–1051. [CrossRef] [PubMed]

19. Haythornthwaite, R.P.; Coates, A.J.; Jones, G.H.; Wellbrock, A.; Waite, J.H.; Vuitton, V.; Lavvas, P. Heavy Positive Ion Groups in Titan's Ionosphere from Cassini Plasma Spectrometer IBS Observations. *Planet. Sci. J.* **2021**, *2*, 26. [CrossRef]

20. Cravens, T.E.; Robertson, I.P.; Waite Jr., J.H.; Yelle, R.V.; Kasprzak, W.T.; Keller, C.N.; Ledvina, S.A.; Niemann, H.B.; Luhmann, J.G.; McNutt, R.L.; et al. Composition of Titan's ionosphere. *Geophys. Res. Lett.* **2006**, *33*, L07105. [CrossRef]

21. Waite, Jr., J.H.; Perryman, R.S.; Perry, M.E.; Miller, K.E.; Bell, J.; Cravens, T.E.; Glein, C.R.; Grimes, J.; Hedman, M.; Cuzzi, J.; et al. Chemical interactions between Saturn's atmosphere and its rings. *Science* **2018**, *362*, eaat2382. [CrossRef]

22. Moore, L.; Cravens, T.E.; Mueller-Wodarg, I.; Perry, M.E.; Waite, Jr., J.H.; Perryman, R.; Nagy, A.; Mitchell, D.; Persoon, A.; Wahlund, J.E.; et al. Models of Saturn's Equatorial Ionosphere Based on In Situ Data From Cassini's Grand Finale. *Geophys. Res. Lett.* **2018**, *45*, 9398–9407. [CrossRef]

23. Cravens, T.E.; Moore, L.; Waite Jr., J.H.; Perryman, R.; Perry, M.; Wahlund, J.E.; Persoon, A.; Kurth, W.S. The Ion Composition of Saturn's Equatorial Ionosphere as Observed by Cassini. *Geophys. Res. Lett.* **2019**, *46*, 6315–6321. [CrossRef]

24. Serigano, J.; Hörst, S.M.; He, C.; Gautier, T.; Yelle, R.V.; Koskinen, T.T.; Koskinen, M.G.; Trainer, M.J.R. Compositional Measurements of Saturn's Upper Atmosphere and Rings from Cassini INMS: An Extended Analysis of Measurements from Cassini's Grand Finale Orbits. *J. Geophys. Res. Planets* **2022**, *127*, e2022JE007238. [CrossRef]

25. Bagenal, F.; Horanyi, M.; McComas, D.J.; McNutt, R.L., Jr.; Elliott, H.A.; Hill, M.E.; Brown, L.E.; Delamere, P.A.; Kollmann, P.; Krimigis, S.M.; et al. Pluto's interaction with its space environment: Solar wind, energetic particles, and dust. *Science* **2016**, *351*, aad9045. [CrossRef]

26. Gladstone, G.R.; Stern, S.A.; Ennico, K.; Olkin, C.B.; Weaver, H.A.; Young, L.A.; Summers, M.E.; Strobel, D.F.; Hinson, D.P.; Kammer, J.A.; et al. The atmosphere of Pluto as observed by New Horizons. *Science* **2016**, *351*, aad8866. [CrossRef]

27. Lellouch, E.; Butler, B.; Moreno, R.; Gurwell, M.; Lavvas, P.; Bertrand, T.; Fouchet, T.; Strobel, D.F.; Moullet, A. Pluto's atmosphere observations with ALMA: Spatially-resolved maps of CO and HCN emission and first detection of HNC. *ICARUS* **2022**, *372*, 114722. [CrossRef]

28. Beth, A.; Altwegg, K.; Balsiger, H.; Berthelier, J.J.; Combi, M.R.; De Keyser, J.; Fiethe, B.; Fuselier, S.A.; Galand, M.; Gombosi, T.I.; et al. ROSINA ion zoo at Comet 67P. *Astron. Astroph.* **2020**, *642*. [CrossRef]

29. Betz, A.L. Ethylene in IRC+10216. *Astrophys. J.* **1981**, *244*, L103–L105. [CrossRef]

30. Hinkle, K.H.; Wallance, L.; Richter, M.J.; Cernicharo, J. Ethylene in the circumstellar envelope of IRC+10216. *Proc. Int. Astron. Union* **2008**, *4*, 161–162. [CrossRef]
31. Steffl, A.J.; Young, L.A.; Strobel, D.F.; Kammer, J.A.; Evans, J.S.; Stevens, M.H.; Schindhelm, R.N.; Parker, J.W.; Stern, S.A.; Weaver, H.A.; et al. Pluto's Ultraviolet Spectrum, Surface Reflectance, and Airglow Emissions. *Astron. J.* **2020**, *159*, 274. [CrossRef]

32. Kammer, J.A.; Gladstone, G.R.; Young, L.A.; Steffl, A.J.; Parker, J.W.; Greathouse, T.K.; Retherford, K.D.; Versteeg, M.H.; Strobel, D.F.; Summers, M.E.; et al. New Horizons Observations of an Ultraviolet Stellar Occultation and Appulse by Pluto's Atmosphere. *Astron. J.* **2020**, *159*, 26. [CrossRef]

33. Lacy, J.H.; Evans, N.J.; Achtermann, J.M.; Bruce, D.E.; Arens, J.F.; Carr, J.. Discovery of Interstellar Acetylene. *Astrophys. J. Lett.* **1989**, *342*, L43. [CrossRef]

34. Ridgway, S.T.; Hall, D.N.; Kleinmann, S.G.; Weinberger, D.A.; Wojslaw, R.S. Circumstellar acetylene in the infrared spectrum of IRC +10216. *Nature* **1976**, *264*, 345–346. [CrossRef]

35. Kostiuk, T.; Espenak, F.; Deming, D.; Mumma, M.J.; Zipoy, D.; Keady, J. Study of Velocity Structure in IRC +10216 Using Acetylene Line Profiles. *Bull. Am. Astron. Soc.* **1985**, *17*, 570.

36. Keady, J.J.; Ridgway, S.T. The IRC +10216 Circumstellar Envelope. III—Infrared Molecular Line Profiles. *Astrophys. J.* **1993**, *406*, 199–214. [CrossRef]

37. Evans, N.J.; Lacy, J.H.; Carr, J.S. Infrared Molecular Spectroscopy toward the Orion IRc2 and IRc7 Sources—A new Probe of Physical Conditions and Abundances in Molecular Clouds. *Astrophys. J.* **1989**, *342*, L43–L46. [CrossRef]

38. Carrasco, N.; Alcaraz, C.; Dutuit, O.; Plessis, S.; Thissen, R.; Vuitton, V.; Yelle, R.; Pernot, P. Sensitivity of a Titan ionospheric model to the ion-molecule reaction parameters. *Planet. Space Sci.* **2008**, *56*, 1644–1657. [CrossRef]

39. Kim, J.; Anicich, V.; Huntress, W. Product distributions and rate constants for the reactions of CH<sub>3</sub><sup>+</sup>, CH<sub>4</sub><sup>+</sup>, C<sub>2</sub>H<sub>2</sub><sup>+</sup>, C<sub>2</sub>H<sub>3</sub><sup>+</sup>, C<sub>2</sub>H<sub>4</sub><sup>+</sup>, and C<sub>2</sub>H<sub>6</sub><sup>+</sup> ions with CH<sub>4</sub>, C<sub>2</sub>H<sub>2</sub>, C<sub>2</sub>H<sub>4</sub>, and C<sub>2</sub>H<sub>6</sub>. *J. Phys. Chem.* **1977**, *81*, 1798–1805. [CrossRef]

40. Anicich, V. Evaluated Bimolecular Ion-Molecule Gas Phase Kinetics of Positive Ions for Use in Modelling Planetary Atmospheres Cometary Comae/Interstellar Clouds. *J. Phys. Chem. Ref. Data* **1993**, *22*, 1469–1569. [CrossRef]

41. Anicich, V. *An Index of the Literature for Bimolecular Gas Phase Cation-Molecule Reaction Kinetics*; JPL Open Repository: Pasadena, CA, USA, 2003. Available online: https://hdl.handle.net/2014/7981 (accessed on 1 May 2023).

42. Mandt, K.E.; Gell, D.A.; Perry, M.; Waite, Jr., J.H.; Crary, F.A.; Young, D.; Magee, B.A.; Westlake, J.H.; Cravens, T.; Kasprzak, W.; et al. Ion densities and composition of Titan's upper atmosphere derived from the Cassini Ion Neutral Mass Spectrometer: Analysis methods and comparison of measured ion densities to photochemical model simulations. *J. Geophys. Res. Planets* **2012**, *117*, E10006. [CrossRef]

43. Nixon, C.A.; Thelen, A.E.; Cordiner, M.A.; Kisiel, Z.; Charnley, S.B.; Molter, E.M.; Serigano, J.; Irwin, P.G.J.; Teanby, N.A.; Kuan, Y.J. Detection of Cyclopropenylidene on Titan with ALMA. *Astron. J.* **2020**, *160*, 205. [CrossRef]

44. Vuitton, V.; Yelle, R.V.; Cui, J. Formation and distribution of benzene on Titan. *J. Geophys. Res. Space Phys.* **2008**, *113*, E05007. [CrossRef]

45. Misina, M.; Pokorny, P. Energy and composition of the ion flux in microwave electron cyclotron resonance/radio frequency methane plasma. *Surf. Coat. Technol.* **2003**, *173–174*, 914–917. [CrossRef]

46. Okita, A.; Suda, Y.; Ozeki, A.; Sugawara, H.; Sakai, Y.; Oda, A.; Nakamura, J. Predicting the amount of carbon in carbon nanotubes grown by CH<sub>4</sub> rf plasmas. *J. Appl. Phys.* **2006**, *99*, [CrossRef]

47. Heijkers, S.; Aghaei, M.; Bogaerts, A. Plasma-Based CH<sub>4</sub> Conversion into Higher Hydrocarbons and H<sub>2</sub>: Modelling to Reveal the Reaction Mechanisms of Different Plasma Sources. *J. Phys. Chem. C* **2020**, *124*, 7016–7030. [CrossRef] [PubMed]

48. Yao, S.; Okumoto, M.; Nakayama, A.; Suzuki, E. Plasma Reforming and Coupling of Methane with Carbon Dioxide. *Energy Fuels* **2001**, *15*, 1295–1299. [CrossRef]

49. Inayat, A.; Tariq, R.; Khan, Z.; Ghenai, C.; Kamil, M.; Jamil, F.; Shanableh, A. A comprehensive review on advanced ther-mochemical processes for bio-hydrogen production via mircowave and plasma technologies. *Biomass Convers. Biorefin.* **2020**, *13*, 8593–8602. [CrossRef]

50. Khoja, A.H.; Tahir, M.; Amin, N.A.S. Recent developments in non-thermal catalytic DBD plasma reactor for dry reforming of methane. *Energy Convers. Manag.* **2019**, *183*, 529–560. [CrossRef]

51. De Bie, C.; van Dijk, J.; Bogaerts, A. The Dominant Pathways for the Conversion of Methane into Oxygenates and Syngas in an Atmospheric Pressure Dielectric Barrier Discharge. *J. Phys. Chem. C* **2015**, *119*, 22331–22350. [CrossRef]

52. Wang, W.; Berthelot, A.; Zhang, Q.; Bogaerts, A. Modelling of plasma-based dry reforming: How do uncertainties in the input data affect the calculation results? *J. Phys. D-Appl. Phys.* **2018**, *51*, 204003. [CrossRef]

53. Chen, T.Y.; Taneja, T.S.; Rousso, A.C.; Yang, S.; Kolemen, E.; Ju, Y. Time-resolved in situ measurements and predictions of plasma-assisted methane reforming in a nanosecond-pulsed discharge. *Proc. Combust. Inst.* **2021**, *38*, 6533–6540. [CrossRef]

54. Zhang, L.; Heijkers, S.; Wang, W.; Martini, L.M.; Tosi, P.; Yang, D.; Fang, Z.; Bogaerts, A. Dry reforming of methane in a nanosecond repetitively pulsed discharge: Chemical kinetics modeling. *Plasma Sources Sci. Technol.* **2022**, *31*, 055014. [CrossRef]

55. Richardson, V.; Ascenzi, D.; Sundelin, D.; Alcaraz, C.; Romanzin, C.; Thissen, R.; Guillemin, J.C.; Polášek, M.; Tosi, P.; Žabka, J.; et al. Experimental and Computational Studies on the Reactivity of Methanimine Radical Cation (H<sub>2</sub>CNH<sup>+</sup>) and its Isomer Aminomethulene (HCNH<sub>2</sub><sup>+</sup>) with C<sub>2</sub>H<sub>2</sub>. *Front. Astron. Space Sci.* **2021**, *8*, 158. [CrossRef]

56. Sundelin, D.; Ascenzi, D.; Richardson, V.; Alcaraz, C.; Polášek, M.; Romanzin, C.; Thissen, R.; Tosi, P.; Žabka, J.; Geppert, W.D. The reactivity of methanimine radical cation (H<sub>2</sub>CNH<sup>+</sup>) and its isomer aminomethylene (HCNH<sub>2</sub><sup>+</sup>) with C<sub>2</sub>H<sub>4</sub>. *Chem. Phys. Lett.* **2021**, *777*, 138677. [CrossRef]

57. Maitre, P.A.; Bieniek, M.S.; Kechagiopoulos, P.N. Modelling excited species and their role on kinetic pathways in the non-oxidative coupling of methane by dielectric barrier discharge. *Chem. Eng. Sci.* **2021**, *234*, 116399. [CrossRef]
58. Gardiner, S.H.; Karsili, T.N.V.; Lipciuc, M.L.; Wilman, E.; Ashfold, M.N.R.; Vallance, C. Fragmentation dynamics of the ethyl bromide and ethyl iodide cations: a velocity-map imaging study. *Phys. Chem. Chem. Phys.* **2014**, *16*, 2167–2178. [CrossRef]

59. Borkar, S.; Sztáray, B. Self-Consistent Heats of Formation for the Ethyl Cation, Ethyle Bromide, and Ethyl Iodide from Threshold Photoelectron Photoion Coincidence Spectroscopy. *J. Phys. Chem. A.* **2010**, *114*, 6117–6123. [CrossRef]

60. Baer, T.; Song, Y.; Liu, J.; Chen, W.; Ng, C.Y. Pulsed field ionization-photoelectron photoion coincidence spectroscopy with synchrotron radiation: The heat of formation of the C₂H₅⁺ ion. *Faraday Discuss.* **2000**, *115*, 137–145. [CrossRef]

61. Miller, B.E.; Baer, T. Kinetic energy release distribution in the fragmentation of energy-selected vinyl and ethyl bromide ions. *Chem. Phys.* **1984**, *85*, 39–45. [CrossRef]

62. Traeger, J.C.; McLoughlin, R.G. Absolute heats of formation for gas-phase cations. *J. Am. Chem. Soc.* **1981**, *103*, 3647–3652. [CrossRef]

63. Xu, D.; Price, R.J.; Huang, J.; Jackson, W.M. Photodissociation of the Ethyl Bromide Cation at 355 nm by Means of TOF-MS and Ion Velocity Imaging Techniques. *Z. Für Phys. Chem.* **2001**, *215*, 253–271. [CrossRef]

64. Kramida, A.; Ralchenko, Y.; Reader, J.; NIST ASD Team. NIST Atomic Spectra Database (Version 5.11). 2023. Available online: https://physics.nist.gov/asd (accessed on 1 May 2023).

65. Kimura, K.; Katsumatu, S.; Y., A.; Yamaziki, T.; Iwata, S. *Handbook of He(I) Photoelectron Spectra of Fundamental Organic Molecules*; Japan Scientific Societies Press und Halstead Press: Tokyo, Japan; New York, NY, USA, 1981.

66. Holmes, J.; Aubry, C.; Mayer, P. *Assigning Structures to Ions in Mass Spectrometry*; CRC Press, Taylor and Francis Group: Boca Raton, FL, USA, 2006. [CrossRef]

67. Cunje, A.; Lien, M.H.; Hopkinson, A.C. The C₄H₅⁺ Potential Energy Surface. Structure, Relative Energies, and Enthalpies of Formation of Isomers of C₄H₅⁺. *J. Org. Chem.* **1996**, *61*, 5212–5220. [CrossRef]

68. Catani, K.J.; Muller, G.; da Silva, G.; Bieske, E.J. Electronic spectrum and photodissociation chemistry of the linear methyl propargyl cation H₂C₄H₃⁺. *J. Chem. Phys.* **2017**, *146*, 044307. [CrossRef]

69. Ruscic, B.; Bross, D. Active Thermochemical Tables (ATcT) Values Based on ver. 1.130 of the Thermochemical Network. 2022. Available online: https://doi.org/10.17038/CSE/1997229 (accessed on 1 May 2023).

70. Linstrom, P.J.; Mallard, W.G. NIST Chemistry WebBook-Standard Reference Database n. 69. Available online: https://doi.org/10.18434/T4D303 (accessed on 1 May 2023).

71. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision C.01, 2016, Gaussian Inc. Wallingford CT. Available online: https://gaussian.com/(accessed on 1 May 2023).

72. Lago, A.; Baer, T. 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₂⁺, and C₂H₃Br₃. *J. Phys. Chem. A.* **2006**, *110*, 3036–3041. [CrossRef]

73. Blair, A.S.; Heslin, E.J.; Harrison, A.G. Biomolecular reactions of trapped ions. IV. Reactions in gaseous ethane and mixtures with acetylene and methane-d4. *J. Am. Chem. Soc.* **1972**, *94*, 2935–2939. [CrossRef]

74. Ervin, K.M.; Armentrout, P. Translational energy dependence of Ar⁺ + XY - ArX⁺ + Y (XY = H₂, D₂, HD) from thermal to 30 eV c.m. *J. Chem. Phys.* **1985**, *83*, 166–189. [CrossRef]

75. Nicolas, C.; Alcaraz, C.; Thissen, R.; Žabka, J.; Dutuit, O. Effects of ion excitation on ion-molecule reactions of the Mars, Venus, and Earth ionospheres. *Plan. Space Sci.* **2002**, *50*, 877–887. [CrossRef]

76. Richardson, V.; Valença Ferreira de Aragão, E.; He, X.; Pirani, F.; Mancini, L.; Faginas-Lago, N.; Rosi, M.; Martini, L.M.; Ascenzi, D. Fragmentation of interstellar methanol by collisions with He⁺: an experimental and computational study. *Phys. Chem. Chem. Phys.* **2022**, *24*, 22437–22452. [CrossRef]

77. Tsikritea, A.; Diprose, J.A.; Softley, T.P.; Heazlewood, B.R. Capture theory models: An overview of their development, experimental verification, and applications to ion-molecule reactions. *J. Chem. Phys.* **2022**, *157*, 060901. [CrossRef]

78. Szilágyi, Z.; Vékey, K. A simple algorithm for the calculation of kinetic energy release distributions. *Eur. Mass Spectrom.* **1995**, *1*, 507–518. [CrossRef]

79. Ess, D.H.; Wheeler, S.E.; Iafe, R.G.; Xu, L.; Çelebi Ölçüm, N.; Houk, K.N. Bifurcations on Potential Energy Surfaces of Organic Reactions. *Angew. Chem. Int. Ed.* **2008**, *47*, 7592–7601. [CrossRef]

80. Bruce, J.E.; Eyler, J.R. Probing Trapped Ion Energies Via Ion-Molecule Reaction Kinetics: Fourier Transform Ion Cyclotron Resonance Mass Spectrometry. *J. Am. Soc. Mass Spectrom.* **1992**, *3*, 727–733. [CrossRef]

81. Grover, R.; Decouzon, M.; Maria, P.C.; Gal, J.F. Reliability of Fourier transform-ion cyclotron resonance determinations of rate constants for ion/molecule reactions. *Eur. J. Mass Spectrom.* **1996**, *2*, 213–223. [CrossRef]

82. Carpignano, R.; Operti, L.; Rabezzana, R.; Vaglio, G.A. Chemometric Approach to Evaluate the Parameters Affecting the Determination of Reaction Rate Constants by Ion Trap Mass Spectrometry. *J. Am. Soc. Mass Spectrom.* **1998**, *9*, 938–944. [CrossRef]

83. Kemper, P.R.; Bowers, M.T. An improved tandem mass spectrometer-ion-cyclotron-resonance spectrometer. *Int. J. Mass Spectrom.* *Ion Phys.* **1983**, *52*, 1–24. [CrossRef]

84. Richardson, V.; Alcaraz, C.; Geppert, W.D.; Polášek, M.; Romanzin, C.; Sundelin, D.; Thissen, R.; Tosi, P.; Žabka, J.; Ascenzi, D. The reactivity of methanimine radical cation (H₂CNH⁺) and its isomer aminomethylene (HCNH₂⁺) with methane. *Chem. Phys. Lett.* **2021**, *775*, 138611. [CrossRef]

85. Alcaraz, C.; Nicolas, C.; Thissen, R.; Žabka, J.; Dutuit, O. ¹⁵N⁺ + CD₄ and O⁺ + ¹³CO₂ State-Selected Ion-Molecule Reactions Relevant to the Chemistry of Planetary Ionospheres. *J. Phys. Chem. A.* **2004**, *108*, 9998–10009. [CrossRef]
86. Cunha de Miranda, B.; Romanzin, C.; Chefdeville, S.; Vuitton, V.; Žabka, J.; Polášek, M.; Alcaraz, C. Reactions of State-Selected Atomic Oxygen Ions O⁺ (⁴S, ²D, ²P) with Methane. *J. Phys. Chem. A* **2015**, *119*, 6082–6098. [CrossRef] [PubMed]

87. Nahon, L.; de Oliveira, N.; Garcia, G.; Gil, J.F.; Pilette, B.; Marcouilleé, O.; Lagarde, B.; Polack, F. DESIRS: A state-of-the-art VUV beamline featuring high resolution and variable polarization for spectroscopy and dichroism at SOLEIL. *J. Synchrotron Rad.* **2012**, *19*, 508–520. [CrossRef] [PubMed]

88. Mercier, B.; Complin, M.; Prevost, C.; Bellec, G.; Thissen, R.; Dutuit, O.; Nahon, L. 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* **2000**, *18*, 2533–2541. [CrossRef]

89. Minnhagen, L. Spectrum and the energy levels of neutral argon, Ar I. *J. Opt. Soc. Am.* **1973**, *63*, 1185–1198. [CrossRef]

90. Teloy, E.; Gerlich, D. Integral cross sections for ion-molecule rections. I. The guided beam technique. *Chem. Phys.* **1974**, *4*, 417–427. [CrossRef]

91. Mardirossian, N.; Head-Gordon, M. ωB97X-V: A 10-parameter, range-separated hybrid, generalized gradient approximation density functional with nonlocal correlation, designed by a survival-of-the-fittest strategy. *Phys. Chem. Chem. Phys.* **2014**, *16*, 9904–9924. [CrossRef]

92. Curtiss, L.A.; Redfern, P.C.; Raghavachari, K. Gaussian-4 theory. *J. Chem. Phys.* **2007**, *126*, 084108. [CrossRef]

93. Silva, W.G.D.P.; Cernicharo, J.; Schlemmer, S.; Marcelino, N.; Loison, J.C.; Agúndez, M.; Gupta, D.; Wakelam, V.; Thorwirth, S.; Cabezas, C.; et al. Discovery of H₂CCCH⁺ in TMC-1. *Astron. Astrophys.* **2023**, *676*, L1. Erratum in *Astron. Astrophys.* **2023**, *676*, C5. [CrossRef]

94. Cernicharo, J.; Gottlieb, C.A.; Guelin, M.; Killian, T.C.; Paubert, G.; Thaddeus, P.; Vrtilek, J.M. Astronomical Detection of H₂CCC. *Astrophys. J. Lett.* **1991**, *368*, L39. [CrossRef]

95. Thaddeus, P.; Vrtilek, J.M.; Gottlieb, C.A. Laboratory and astronomical identification of cyclopropenylidene, C₃H₂. *Astrophys. J.* **1985**, *299*, L63–L66. [CrossRef]

96. Tsai, B.P.; Werner, A.S.; Baer, T. A photoion-photoelectron coincidence (PIPECO) study of fragmentation rates and kinetic energy release in energy selected metastable ions. *J. Chem. Phys.* **1975**, *63*, 4384–4392. [CrossRef]

97. Booze, J.A.; Weitzel, K.M.; Baer, T. The rates of HCl loss from energy-selected ethylchloride ions: A case of tunneling through an H-atom transfer barrier. *J. Chem. Phys.* **1991**, *94*, 3649–3656. [CrossRef]

98. Kuck, D. Half a century of scrambling in organic ions: Complete, incomplete, progressive and composite atom interchange. *Int. J. Mass Spectrom.* **2002**, *213*, 101–144. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.