Formation of the methyl cation by photochemistry in a protoplanetary disk.

Forty years ago, it was proposed that gas-phase organic chemistry in the interstellar medium can be initiated by the methyl cation CH3+ (refs. 1-3), but so far it has not been observed outside the Solar System4,5. Alternative routes involving processes on grain surfaces have been invoked6,7. Here we report James Webb Space Telescope observations of CH3+ in a protoplanetary disk in the Orion star-forming region. We find that gas-phase organic chemistry is activated by ultraviolet irradiation.

Diversity of protonated mixed pyrene-water clusters investigated by collision induced dissociation.

Protonated mixed pyrene-water clusters, (Py)m(H2O)nH+, where m = [1-3] and n = [1-10], are generated using a cryogenic molecular cluster source. Subsequently, the mass-selected mixed clusters undergo controlled collisions with rare gases, and the resulting fragmentation mass spectra are meticulously analyzed to discern distinct fragmentation channels. Notably, protonated water cluster fragments emerge for n ≥ 3, whereas they are absent for n = 1 and 2. The experimental results are complemented by theoretical calculations of structures and energetics for (Py)(H2O)nH+ with n = [1-4]. These calculations reveal a shift in proton localization, transitioning from the pyrene molecule for n = 1 and 2 to water molecules for n ≥ 3. The results support a formation scenario wherein water molecules attach to protonated pyrene PyH+ seeds, and, by extension, to (Py)2H+ and (Py)3H+ seeds. Various isomers are identified, corresponding to potential protonation sites on the pyrene molecule. Protonated polycyclic aromatic hydrocarbons are likely to be formed in cold, dense interstellar clouds and protoplanetary disks due to the high proton affinity of these species. Our findings show that the presence of protonated PAHs in these environments could lead to the formation of water clusters and mixed carbon-water nanograins, having a potential impact on the water cycle in regions of planet formation.

Near-infrared absorption and radiative cooling of naphthalene dimers (C10H8)2.

The radiative cooling of naphthalene dimer cations, (C10H8)2+ was studied experimentally through action spectroscopy using two different electrostatic ion-beam storage rings, DESIREE in Stockholm and Mini-Ring in Lyon. The spectral characteristics of the charge resonance (CR) band were observed to vary significantly with a storage time of up to 30 seconds in DESIREE. In particular, the position of the CR band shifts to the blue, with specific times (inverse of rates) of 0.64 s and 8.0 s in the 0-5 s and 5-30 s storage time ranges, respectively. These long-time scales indicate that the internal energy distribution of the stored ions evolves by vibrational radiative cooling, which is consistent with the absence of fast radiative cooling via recurrent fluorescence for (C10H8)2+. Density functional based tight binding calculations with local excitations and configuration interactions (DFTB-EXCI) were used to simulate the absorption spectrum for ion temperatures between 10 and 500 K. The evolution of the bandwidth and position with temperature is in qualitative agreement with the experimental findings. Furthermore, these calculations yielded linear temperature dependencies for both the shift and the broadening. Combining the relationship between the CR band position and the ion temperature with the results of the statistical model, we demonstrate that the observed blue shift can be used to determine the radiative cooling rate of (C10H8)2+.

Electronic effects in the dissociative ionisation of pyrene clusters.

We present a combined experimental and theoretical study on the dissociative ionisation of clusters of pyrene. We measured the experimental appearance energies in the photon energy range 7.2-12.0 eV of the fragments formed from neutral monomer loss for clusters up to the hexamer. The results obtained show a deviation from statistical dissociation. From electronic structure calculations, we suggest that the role of excited states must be considered in the interpretation of experimental results, even in these relatively large systems. Non-statistical effects in the dissociative ionization process of polycyclic aromatic hydrocarbon (PAH) clusters may have an impact on the assessment of mechanisms determining the stability of these clusters in astrophysical environments.

Threshold collision induced dissociation of protonated water clusters.

We report threshold collision induced dissociation experiments on protonated water clusters thermalized at low temperature for sizes n = 19-23. Fragmentation cross sections are recorded as a function of the collision energy and analyzed with a statistical model. This model allows us to account for dissociation cascades and provides values for the dissociation energies of each cluster. These values, averaging around 0.47 eV, are in good agreement with theoretical predictions at various levels of theory. Furthermore, the dissociation energies show a trend for the n = 21 magic and n = 22 anti-magic numbers relative to their neighbours, which is also in agreement with theory. These results provide further evidence to resolve the disagreement between previously published experimental values. A careful quantitative treatment of cascade dissociation in this model introduces interdependence between the dissociation energies of neighboring sizes, which reduces the number of free fitting parameters and improves both reliability and uncertainties on absolute dissociation energies deduced from experiments.

Formation of the acenaphthylene cation as a common C2H2-loss fragment in dissociative ionization of the PAH isomers anthracene and phenanthrene.

Polycyclic aromatic hydrocarbons (PAHs) are thought to be a major constituent of astrophysical environments, being the carriers of the ubiquitous aromatic infrared bands (AIBs) observed in the spectra of galactic and extra-galactic sources that are irradiated by ultraviolet (UV) photons. Small (2-cycles) PAHs were unambiguously detected in the TMC-1 dark cloud, showing that PAH growth pathways exist even at low temperatures. The processing of PAHs by UV photons also leads to their fragmentation, which has been recognized in recent years as an alternative route to the generally accepted bottom-up chemical pathways for the formation of complex hydrocarbons in UV-rich interstellar regions. Here we consider the C12H8+ ion that is formed in our experiments from the dissociative ionization of the anthracene and phenanthrene (C14H10) molecules. By employing the sensitive action spectroscopic scheme of infrared pre-dissociation (IRPD) in a cryogenic ion trap instrument coupled to the free-electron lasers at the FELIX Laboratory, we have recorded the broadband and narrow line-width gas-phase IR spectra of the fragment ions (C12H8+) and also the reference spectra of three low energy isomers of C12H8+. By comparing the experimental spectra to those obtained from quantum chemical calculations we have identified the dominant structure of the fragment ion formed in the dissociation process to be the acenaphthylene cation for both isomeric precursors. Ab initio molecular dynamics simulations are presented to elucidate the fragmentation process. This result reinforces the dominant role of species containing a pentagonal ring in the photochemistry of small PAHs.

Water Attachment onto Size-Selected Cationic Pyrene Clusters.

We report measurements of the attachment rates of water molecules onto mass-selected cationic pyrene clusters for size from n = 4 to 13 pyrene units and for different collision energies. Comparison of the attachment rates with the collision rates measured in collision-induced dissociation experiments provides access to the values of the sticking coefficient. The strong dependence of the attachment rates on size and collision energy is rationalized through a model in which we use a Langevin-type collision rate and adjust on experimental data the statistical dissociation rate of the water molecule from the cluster after attachment. This allows us to extrapolate our results to the conditions of isolation and long time scales encountered in astrophysical environments.

Isomer Differentiation of Trapped C16H10+ Using Low-Energy Collisions and Visible/VUV Photons.

Polycyclic aromatic hydrocarbons are major species in astrophysical environments, and this motivates their study in samples of astrochemical interest such as meteorites and laboratory analogues of stardust. Molecular analyses of carbonaceous matter in these samples show a dominant peak at m/z = 202.078 corresponding to C16H10. Obtaining information on the associated isomeric structures is a challenge for the molecular analysis of samples available in very small quantities (mg or less). Here we show that coupling laser desorption ionization mass spectrometry with ion trapping opens up the possibility of unraveling isomers by activating ion fragmentation via collisions or photon absorption. We report the best criteria for differentiating isomers with comparable dissociation energies, namely pyrene, fluoranthene, and 9-ethynylphenanthrene, on the basis of the parent dissociation curve and the ratio of dehydrogenation channels. Photoabsorption schemes (multiple photon absorption in the visible range and single photon absorption at 10.5 eV) are more effective in differentiating these isomers than activation by low energy collisions. The impact of the activation scheme on the fragmentation kinetics and dehydrogenation pathways is discussed. By analyzing the 10.5 eV photodissociation measurements with a simple kinetic model, we were able to derive a branching ratio for the H and 2H/H2 loss channels of the parent ions. The results suggest a role in the formation of H2 for bay hydrogens that are present in both fluoranthene and 9-ethynylphenanthrene. In addition, we suggest for the latter the presence of a highly competitive 2H loss channel, possibly associated with the formation of a pentagonal ring.

Compression and ablation of the photo-irradiated molecular cloud the Orion Bar.

The Orion Bar is the archetypal edge-on molecular cloud surface illuminated by strong ultraviolet radiation from nearby massive stars. Our relative closeness to the Orion nebula (about 1,350 light years away from Earth) means that we can study the effects of stellar feedback on the parental cloud in detail. Visible-light observations of the Orion Bar show that the transition between the hot ionized gas and the warm neutral atomic gas (the ionization front) is spatially well separated from the transition between atomic and molecular gas (the dissociation front), by about 15 arcseconds or 6,200 astronomical units (one astronomical unit is the Earth-Sun distance). Static equilibrium models used to interpret previous far-infrared and radio observations of the neutral gas in the Orion Bar (typically at 10-20 arcsecond resolution) predict an inhomogeneous cloud structure comprised of dense clumps embedded in a lower-density extended gas component. Here we report one-arcsecond-resolution millimetre-wave images that allow us to resolve the molecular cloud surface. In contrast to stationary model predictions, there is no appreciable offset between the peak of the H2 vibrational emission (delineating the H/H2 transition) and the edge of the observed CO and HCO+ emission. This implies that the H/H2 and C+/C/CO transition zones are very close. We find a fragmented ridge of high-density substructures, photoablative gas flows and instabilities at the molecular cloud surface. The results suggest that the cloud edge has been compressed by a high-pressure wave that is moving into the molecular cloud, demonstrating that dynamical and non-equilibrium effects are important for the cloud evolution.

Threshold collision induced dissociation of pyrene cluster cations.

We report threshold collision induced dissociation experiments on cationic pyrene clusters, for sizes n = 2-6. Fragmentation cross sections are recorded as a function of the collision energy and analyzed with a statistical model. This model can account for the dissociation cascades and provides values for the dissociation energies. These values, of the order of 0.7 eV-1 eV, are in excellent agreement with those previously derived from thermal evaporation. They confirm the charge resonance stability enhancement predicted by theoretical calculations. In addition, remarkable agreement is obtained with theoretical predictions for the two smaller sizes n = 2 and 3. For the larger sizes, the agreement remains good, although the theoretical values obtained for the most stable structures are systematically higher by 0.2 eV. This offset could be attributed to approximations in the calculations. Still, there is an indication in the results of an incomplete description of the role of isomerization and/or direct dissociation upon collisions. Finally, by-product clusters containing dehydrogenated species are found to dissociate at energies comparable to the non-dehydrogenated ones, which shows no evidence for covalent bonds within the clusters.

Experimental Approach to the Study of Anharmonicity in the Infrared Spectrum of Pyrene from 14 to 723 K.

Quantifying the effect of anharmonicity on the infrared spectrum of large molecules such as polycyclic aromatic hydrocarbons (PAHs) at high temperatures is the focus of a number of theoretical and experimental studies, many of them motivated by astrophysical applications. We recorded the IR spectrum of pyrene C16H10 microcrystals embedded in KBr pellets over a wide range of temperatures (14-723 K) and studied the evolution of band positions, widths, and integrated intensities with temperature. We identified jumps for some of the spectral characteristics of some bands in the 423-473 K range. These were attributed to a change of phase from crystal to molten in condensed pyrene, which appears to affect more strongly bands involving large CH motions. Empirical anharmonicity factors that quantify the linear evolution of band positions and widths with temperature for values larger than ∼150-250 K, depending on the band, were retrieved from both phases and averaged to provide recommended values for these anharmonicity factors. The derived values were found to be consistent with available gas phase data. We conclude about the relevance of the methodology to produce data that can be compared with calculated anharmonic IR spectra and provide input for models that simulate the IR emission of astro-PAHs.

Thermal evaporation of pyrene clusters.

This work presents a study of the thermal evaporation and stability of pyrene (C16H10)n clusters. Thermal evaporation rates of positively charged mass-selected clusters are measured for sizes in the range n = 3-40 pyrene units. The experimental setup consists of a gas aggregation source, a thermalization chamber, and a time of flight mass spectrometer. A microcanonical Phase Space Theory (PST) simulation is used to determine the dissociation energies of pyrene clusters by fitting the experimental breakdown curves. Calculations using the Density Functional based Tight Binding combined with a Configuration Interaction (CI-DFTB) model and a hierarchical optimization scheme are also performed in the range n = 2-7 to determine the harmonic frequencies and a theoretical estimation of the dissociation energies. The frequencies are used in the calculations of the density of states needed in the PST simulations, assuming an extrapolation scheme for clusters larger than 7 units. Using the PST model with a minimal set of adjustable parameters, we obtain good fits of the experimental breakdown curves over the full studied size range. The approximations inherent to the PST simulation and the influence of the used parameters are carefully estimated. The derived dissociation energies show significant variations over the studied size range. Compared with neutral clusters, significantly higher values of the dissociation energies are obtained for the smaller sizes and attributed to charge resonance in line with CI-DFTB calculations.

Direct Evidence of the Benzylium and Tropylium Cations as the Two Long-Lived Isomers of C7 H7.

Disentangling the isomeric structure of C7 H7 + is a longstanding experimental issue. We report here the full mid-infrared vibrational spectrum of C7 H7 + tagged with Ne obtained with infrared-predissociation spectroscopy at 10 K. Saturation depletion measurements were used to assign the contribution of benzylium and tropylium isomers and demonstrate that no other isomer is involved. Recorded spectral features compare well with density functional theory calculations. This opens perspectives for a better understanding and control of the formation paths leading to either tropylium or benzylium ions.

Unimolecular reaction energies for polycyclic aromatic hydrocarbon ions.

Imaging photoelectron photoion coincidence spectroscopy was employed to explore the unimolecular dissociation of the ionized polycyclic aromatic hydrocarbons (PAHs) acenaphthylene, fluorene, cyclopenta[d,e,f]phenanthrene, pyrene, perylene, fluoranthene, dibenzo[a,e]pyrene, dibenzo[a,l]pyrene, coronene and corannulene. The primary reaction is always hydrogen atom loss, with the smaller species also exhibiting loss of C2H2 to varying extents. Combined with previous work on smaller PAH ions, trends in the reaction energies (E0) for loss of H from sp2-C and sp3-C centres, along with hydrocarbon molecule loss were found as a function of the number of carbon atoms in the ionized PAHs ranging in size from naphthalene to coronene. In the case of molecules which possessed at least one sp3-C centre, the activation energy for the loss of an H atom from this site was 2.34 eV, with the exception of cyclopenta[d,e,f]phenanthrene (CPP) ions, for which the E0 was 3.44 ± 0.86 eV due to steric constraints. The hydrogen loss from PAH cations and from their H-loss fragments exhibits two trends, depending on the number of unpaired electrons. For the loss of the first hydrogen atom, the energy is consistently ca. 4.40 eV, while the threshold to lose the second hydrogen atom is much lower at ca. 3.16 eV. The only exception was for the dibenzo[a,l]pyrene cation, which has a unique structure due to steric constraints, resulting in a low H loss reaction energy of 2.85 eV. If C2H2 is lost directly from the precursor cation, the energy required for this dissociation is 4.16 eV. No other fragmentation channels were observed over a large enough sample set for trends to be extrapolated, though data on CH3 and C4H2 loss obtained in previous studies is included for completeness. The dissociation reactions were also studied by collision induced dissociation after ionization by atmospheric pressure chemical ionization. When modeled with a simple temperature-based theory for the post-collision internal energy distribution, there was reasonable agreement between the two sets of data.

Investigating the importance of edge-structure in the loss of H/H2 of PAH cations: the case of dibenzopyrene isomers.

We present a detailed study of the main dehydrogenation processes of two dibenzopyrene cation (C24H14+) isomers, namely dibenzo(a,e)pyrene (AE+) and dibenzo(a,l)pyrene (AL+). First, action spectroscopy under VUV photons was performed using synchrotron radiation in the 8-20 eV range. We observed lower dissociation thresholds for the non-planar molecule (AL+) than for the planar one (AE+) for the main dissociation pathways: H and 2H/H2 loss. In order to rationalize the experimental results, dissociation paths were investigated by means of density functional theory calculations. In the case of H loss, which is the dominant channel at the lowest energies, the observed difference between the two isomers can be explained by the presence in AL+ of two C-H bonds with considerably lower adiabatic dissociation energies. In both isomers the 2H/H2 loss channels are observed only at about 1 eV higher than H loss. We suggest that this is due to the propensity of bay H atoms to easily form H2. In addition, in the case of AL+, we cannot exclude a competition between 2H and H2 channels. In particular, the formation of a stable dissociation product with a five-membered ring could account for the low energy sequential loss of 2 hydrogens. This work shows the potential role of non-compact PAHs containing bay regions in the production of H2 in space.

Identification of the fragment of the 1-methylpyrene cation by mid-IR spectroscopy.

The fragment of the 1-methylpyrene cation, C17H11+, is expected to exist in two isomeric forms, 1-pyrenemethylium PyrCH2+ and the tropylium containing species PyrC7+. We measured the infrared (IR) action spectrum of cold C17H11+ tagged with Ne using a cryogenic ion trap instrument coupled to the FELIX laser. Comparison of the experimental data with density functional theory calculations allows us to identify the PyrCH2+ isomer in our experiments. The IR Multi-Photon Dissociation spectrum was also recorded following the C2H2 loss channel. Its analysis suggests combined effects of anharmonicity and isomerisation while heating the trapped ions, as shown by molecular dynamics simulations.

Anharmonic vibrational spectroscopy of polycyclic aromatic hydrocarbons (PAHs).

While powerful techniques exist to accurately account for anharmonicity in vibrational molecular spectroscopy, they are computationally very expensive and cannot be routinely employed for large species and/or at non-zero vibrational temperatures. Motivated by the study of Polycyclic Aromatic Hydrocarbon (PAH) emission in space, we developed a new code, which takes into account all modes and can describe all infrared transitions including bands becoming active due to resonances as well as overtone, combination, and difference bands. In this article, we describe the methodology that was implemented and discuss how the main difficulties were overcome, so as to keep the problem tractable. Benchmarking with high-level calculations was performed on a small molecule. We carried out specific convergence tests on two prototypical PAHs, pyrene (C16H10) and coronene (C24H12), aiming at optimising tunable parameters to achieve both acceptable accuracy and computational costs for this class of molecules. We then report the results obtained at 0 K for pyrene and coronene, comparing the calculated spectra with available experimental data. The theoretical band positions were found to be significantly improved compared to harmonic density functional theory calculations. The band intensities are in reasonable agreement with experiments, the main limitation being the accuracy of the underlying calculations of the quartic force field. This is a first step toward calculating moderately high-temperature spectra of PAHs and other similarly rigid molecules using Monte Carlo sampling.

Cationic Methylene-Pyrene Isomers and Isomerization Pathways: Finite Temperature Theoretical Studies.

This paper provides spectral characterizations of the two isomers of the 1-methylenepyrene cation, namely, the 1-pyrenemethylium and a pyrene-like isomer owing a tropylium cycle. Both are possible photodissociation products of the 1-methylpyrene cation and were proposed as potential contributors to the diffuse interstellar bands. In that respect, vibrational and electronic spectra are computed for the optimized structures at the density functional theory (DFT) and time-dependent (TD-)DFT levels. Finite temperature effects on these spectra are estimated from molecular dynamics simulations within the density functional-based tight-binding (DFTB) and TD-DFTB frameworks, these methods being first benchmarked against DFT and TD-DFT calculations. The computed spectra allow discrimination of the two isomers. When the temperature increases, bands are observed to redshift and merge. The isomerization mechanism is investigated with the metadynamics technique, a biased dynamics scheme allowing to probe reaction mechanisms with high energy barriers by investigating the free energy surface at various temperatures. Four pathways with similar barrier heights (3.5-4 eV) are found, showing that the interconversion process would only occur in interstellar clouds under photoactivation. The present study opens the way to simulations on larger methyl- and methylenePAHs of astrophysical interest and their experimental investigation.

Dissociation of the anthracene radical cation: a comparative look at iPEPICO and collision-induced dissociation mass spectrometry results.

The dissociation of the anthracene radical cation has been studied using two different methods: imaging photoelectron photoion coincidence spectrometry (iPEPCO) and atmospheric pressure chemical ionization-collision induced dissociation mass spectrometry (APCI-CID). Four reactions were investigated: (R1) C14H10(+•) → C14H9(+) + H, (R2) C14H9(+) → C14H8(+•) + H, (R3) C14H10(+•) → C12H8(+•) + C2H2 and (R4) C14H10(+•) → C10H8(+•) + C4H2. An attempt was made to assign structures to each fragment ion, and although there is still room for debate whether for the C12H8(+•) fragment ion is a cyclobuta[b]naphthalene or a biphenylene cation, our modeling results and calculations appear to suggest the more likely structure is cyclobuta[b]naphthalene. The results from the iPEPICO fitting of the dissociation of ionized anthracene are E0 = 4.28 ± 0.30 eV (R1), 2.71 ± 0.20 eV (R2), and 4.20 ± 0.30 eV (average of reaction R3) whereas the Δ(‡)S values (in J K(-1) mol(-1)) are 12 ± 15 (R1), 0 ± 15 (R2), and either 7 ± 10 (using cyclobuta[b]naphthalene ion fragment in reaction R3) or 22 ± 10 (using the biphenylene ion fragment in reaction R3). Modeling of the APCI-CID breakdown diagrams required an estimate of the postcollision internal energy distribution, which was arbitrarily assumed to correspond to a Boltzmann distribution in this study. One goal of this work was to determine if this assumption yields satisfactory energetics in agreement with the more constrained and theoretically vetted iPEPICO results. In the end, it did, with the APCI-CID results being similar.

Dynamics of hydrogen and methyl radical loss from ionized dihydro-polycyclic aromatic hydrocarbons: a tandem mass spectrometry and imaging photoelectron-photoion coincidence (iPEPICO) study of dihydronaphthalene and dihydrophenanthrene.

Ionized 1,2-dihydronaphthalene (C10H10(+)) and 9,10-dihydrophenanthrene (C14H12(+)) are homologous dihydrogenated polycyclic aromatic hydrocarbons containing adjacent sp(3) carbon sites. Tandem mass spectrometry involving kiloelectronvolt collision induced dissociation was employed to aid in the structural characterization of the products of the main dissociation channels, loss of H (and subsequent H and H2 losses in dihydronaphthalene) and CH3. Evident from both the CID mass spectra and the imaging photoelectron-photoion coincidence (iPEPICO) breakdown curves is the fact that there are two competitive routes to the loss of H. For 1,2-dihydronaphthalene these give activation energies of 2.22 ± 0.10 and 2.44 ± 0.05 eV, whereas only 2.37 ± 0.12 eV was obtained for 9,10-dihydrophenanthrene. The two parallel H-loss chaneels are believed to be the result of isomerization taking place to the methylindene ion and the 9-methylfluorene ion for 1,2-dihydronaphthalene and 9,10-dihydrophenanthren, respectively. Each newly formed isomer dissociates by H loss (one of the two competing H-loss reactions) and, of course, methyl loss. Methyl radical loss has similar kinetics for the two systems, E0 = 2.57 ± 0.12 eV, Δ(‡)S = 18 ± 11 J K(-1) mol(-1) for ionized dihydronaphthalene and E0 = 2.38 ± 0.15 eV, Δ(‡)S = -3 ± 15 J K(-1) mol(-1) for ionized dihydrophenanthrene, but as can be seen, the E0 and Δ(‡)S are slightly lower for the latter. The final bond rupture step in both H and CH3 loss is expected to be accompanied by a positive Δ(‡)S, thus the low energy H loss and CH3 loss originate from the isomer ion in both cases, with the entropy of activation being dominated by the isomerization step. In contrast, sp(3)-H loss from the dihydro-PAHs differ by little in both systems (E0 = 2.44 eV in ionized dihydronaphthalene and 2.37 eV in ionized dihydrophenanthrene and the Δ(‡)S values are 27 and 18 J K(-1) mol(-1), respectively). The presence of a second sp(3) carbon site has decreased the C-H bond dissociation energy relative to protonated naphthalene and protonated phenanthrene, possibly to facilitate the restoration of the unaltered PAH ion. The calculated dihedral angle is -34.3° in C10H10(+•) whereas C14H12(+•) has an angle of -49.6°, indicating that to restore the planar nature of the molecules, which is required for all reaction channels investigated, there is more rearrangement needed for 9,10-dihydrophenanthrene. Energetics and entropic values associated with H and H2 loss from [M - H](+) ions from ionized dihydronaphthalene were determined to be 2.72 eV, 9 ± 17 J K(-1) mol(-1), and 2.85 eV, 9 ± 7 J K(-1) mol(-1), respectively.