Multiple ionization and hydrogen loss from neutral and positively-charged coronene

Addressing electronic and dynamical evolution of molecules and molecular clusters: DFTB simulations of energy relaxation in polycyclic aromatic hydrocarbons

We present a review of the capabilities of the density functional based Tight Binding (DFTB) scheme to address the electronic relaxation and dynamical evolution of molecules and molecular clusters following energy deposition via either collision or photoabsorption. The basics and extensions of DFTB for addressing these systems and in particular their electronic states and their dynamical evolution are reviewed. Applications to PAH molecules and clusters, carbonaceous systems of major interest in astrochemical/astrophysical context, are reported. A variety of processes are examined and discussed such as collisional hydrogenation, fast collisional processes and induced electronic and charge dynamics, collision-induced fragmentation, photo-induced fragmentation, relaxation in high electronic states, electronic-to-vibrational energy conversion and statistical versus non-statistical fragmentation. This review illustrates how simulations may help to unravel different relaxation mechanisms depending on various factors such as the system size, specific electronic structure or excitation conditions, in close connection with experiments.

Competition between Loss of H2 versus H+H in the Fragmentation of the Fluorene Cation

The fluorene cation is a frequently studied molecule in the context of fragmentation experiments. This is because of its potential role in interstellar chemistry, notably as a precursor of PAH cages. In this paper, we analyze H, H+ , H2 andH 2 + ${{\rm{H}}_2^ + }$ losses from the fluorene cation using the SMF (Statistical Molecular Fragmentation) model. We calculate the probabilities of all the 534 possible fragmentation channels as a function of the excitation energy, up to the loss of three hydrogens. Four different types of hydrogen atom pairings (from the same carbon, from the same ring, from different rings and across-the-bay) have been tested in order to determine which types contribute to the actual production of hydrogen molecules. The simulated breakdown curves are in very good agreement with different experimental results when same ring pairing is taken into account. It was possible to deduce from the model the locations of the emitted H, H+ , H2 andH 2 + ${{\rm{H}}_2^ + }$ species.

Multiple dehydrogenation of fluorene cation and neutral fluorene using the statistical molecular fragmentation model

The statistical molecular fragmentation (SMF) model was used to analyze the 306 fragmentation channels (containing 611 different species) that result from the fluorene (C13H10+) cation losing up to three hydrogen atoms (neutral radicals and/or a proton). Breakdown curves from such analysis permit one to extract experimentally inaccessible information about the fragmentation of the fluorene cation, such as the locations of the lost hydrogen atoms (or proton), yields of the neutral fragments, electronic states of the residues, and quantification of very low probability channels that would be difficult to detect. Charge localization during the fragmentation pathways was studied to provide a qualitative understanding of the fragmentation process. Breakdown curves for both the fluorene cation and neutral fluorene were compared. The SMF results match the rise and fall of the one hydrogen loss yield experimentally measured by imaging photoelectron-photoion coincidence spectroscopy using a VUV synchrotron.

Doubly charged coronene clusters-Much smaller than previously observed

The smallest doubly charged coronene cluster ions reported so far, Cor₁₅²⁺, were produced by charge exchange between bare coronene clusters and He²⁺ [H. A. B. Johansson et al., Phys. Rev. A 84, 043201 (2011)]. These dications are at least five times larger than the estimated Rayleigh limit, i.e., the size at which the activation barrier for charge separation vanishes. Such a large discrepancy is unheard of for doubly charged atomic or molecular clusters. Here we report the mass spectrometric observation of doubly charged coronene trimers, produced by electron ionization of helium nanodroplets doped with coronene. The observation implies that Cor₃²⁺ features a non-zero fission barrier too large to overcome under the present experimental conditions. The height of the barriers for the dimer and trimer has been estimated by means of density functional theory calculations. A sizeable barrier for the trimer has been revealed in agreement with the experimental findings.

Ion collision-induced chemistry in pure and mixed loosely bound clusters of coronene and C60 molecules

Ionization, fragmentation and molecular growth have been studied in collisions of 22.5 keV He2+- or 3 keV Ar+-projectiles with pure loosely bound clusters of coronene (C24H12) molecules or with loosely bound mixed C60-C24H12 clusters by using mass spectrometry. The heavier and slower Ar+ projectiles induce prompt knockout-fragmentation - C- and/or H-losses - from individual molecules and highly efficient secondary molecular growth reactions before the clusters disintegrate on picosecond timescales. The lighter and faster He2+ projectiles have a higher charge and the main reactions are then ionization by ions that are not penetrating the clusters. This leads mostly to cluster fragmentation without molecular growth. However, here penetrating collisions may also lead to molecular growth but to a much smaller extent than with 3 keV Ar+. Here we present fragmentation and molecular growth mass distributions with 1 mass unit resolution, which reveals that the same numbers of C- and H-atoms often participate in the formation and breaking of covalent bonds inside the clusters. We find that masses close to those with integer numbers of intact coronene molecules, or with integer numbers of both intact coronene and C60 molecules, are formed where often one or several H-atoms are missing or have been added on. We also find that super-hydrogenated coronene is formed inside the clusters.

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 (C₂₄H₁₄⁺) 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/H₂ 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/H₂ 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 H₂. In addition, in the case of AL⁺, we cannot exclude a competition between 2H and H₂ 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 H₂ in space.

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.

Dissociation of polycyclic aromatic hydrocarbons: molecular dynamics studies

We present dynamical studies of the dissociation of polycyclic aromatic hydrocarbon (PAH) radical cations in their ground electronic states with significant internal energy. Molecular dynamics simulations are performed, the electronic structure being described on-the-fly at the self-consistent-charge density functional-based tight binding (SCC-DFTB) level of theory. The SCC-DFTB approach is first benchmarked against DFT results. Extensive simulations are achieved for naphthalene [Formula: see text], pyrene [Formula: see text] and coronene [Formula: see text] at several energies. Such studies enable one to derive significant trends on branching ratios, kinetics, structures and hints on the formation mechanism of the ejected neutral fragments. In particular, dependence of branching ratios on PAH size and energy were retrieved. The losses of H and C₂H₂ (recognized as the ethyne molecule) were identified as major dissociation channels. The H/C₂H₂ ratio was found to increase with PAH size and to decrease with energy. For [Formula: see text], which is the most interesting PAH from the astrophysical point of view, the loss of H was found as the quasi-only channel for an internal energy of 30 eV. Overall, in line with experimental trends, decreasing the internal energy or increasing the PAH size will favour the hydrogen loss channels with respect to carbonaceous fragments.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

VUV photo-processing of PAH cations: quantitative study on the ionization versus fragmentation processes

Interstellar polycyclic aromatic hydrocarbons (PAHs) are strongly affected by the absorption of vacuum ultraviolet (VUV) photons in the interstellar medium (ISM), yet the branching ratio between ionization and fragmentation is poorly studied. This is crucial for the stability and charge state of PAHs in the ISM in different environments, affecting in turn the chemistry, the energy balance, and the contribution of PAHs to the extinction and emission curves. We studied the interaction of PAH cations with VUV photons in the 7 - 20 eV range from the synchrotron SOLEIL beamline, DESIRS. We recorded by action spectroscopy the relative intensities of photo-fragmentation and photo-ionization for a set of eight PAH cations ranging in size from 14 to 24 carbon atoms, with different structures. At photon energies below ~13.6 eV fragmentation dominates for the smaller species, while for larger species ionization is immediately competitive after the second ionization potential (IP). At higher photon energies, all species behave similarly, the ionization yield gradually increases, leveling off between 0.8 and 0.9 at ~18 eV. Among isomers, PAH structure appears to mainly affect the fragmentation cross section, but not the ionization cross section. We also measured the second IP for all species and the third IP for two of them, all are in good agreement with theoretical ones confirming that PAH cations can be further ionized in the diffuse ISM. Determining actual PAH dication abundances in the ISM will require detailed modeling. Our measured photo-ionization yields for several PAH cations provide a necessary ingredient for such models.

Formation of H2 from internally heated polycyclic aromatic hydrocarbons: excitation energy dependence

We have investigated the effectiveness of molecular hydrogen (H2) formation from Polycyclic Aromatic Hydrocarbons (PAHs) which are internally heated by collisions with keV ions. The present and earlier experimental results are analyzed in view of molecular structure calculations and a simple collision model. We estimate that H2 formation becomes important for internal PAH temperatures exceeding about 2200 K, regardless of the PAH size and the excitation agent. This suggests that keV ions may effectively induce such reactions, while they are unlikely due to, e.g., absorption of single photons with energies below the Lyman limit. The present analysis also suggests that H2 emission is correlated with multi-fragmentation processes, which means that the [PAH-2H](+) peak intensities in the mass spectra may not be used for estimating H2-formation rates.

Photodissociation of pyrene cations: structure and energetics from C16H10(+) to C14(+) and almost everything in between

The unimolecular dissociation of the pyrene radical cation, C16H10(+•), has been explored using a combination of computational techniques and experimental approaches, such as multiple photon absorption in the cold ion trap Piège à Ions pour la Recherche et l'Etude de Nouvelles Espèces Astrochimiques (PIRENEA) and imaging photoelectron photoion coincidence spectrometry (iPEPICO). In total, 22 reactions, involving the fragmentation cascade (H, C2H2, and C4H2 loss) from the pyrene radical cation down to the C14(+•) fragment ion, have been studied using PIRENEA. Branching ratios have been measured for reactions from C16H10(+•), C16H8(+•), and C16H5(+). Density functional theory calculations of the fragmentation pathways observed experimentally and postulated theoretically lead to 17 unique structures. One important prediction is the opening of the pyrene ring system starting from the C16H4(+•) radical. In the iPEPICO experiments, only two reactions could be studied, namely, R1 C16H10(+•) → C16H9(+) + H (m/z = 201) and R2 C16H9(+) → C16H8(+•) + H (m/z = 200). The activation energies for these reactions were determined to be 5.4 ± 1.2 and 3.3 ± 1.1 eV, respectively.