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

Coincidence mass spectrometry study of double ionization of pyrene by 70 eV electron impact

We have performed coincidence mass spectrometry of fragmentation of pyrene molecules by 70 eV electron impact. Ionized fragments have been mass selected using a reflectron time-of-flight mass spectrometer, and a field programmable gate array has been used for the timing of the electron and ion extraction pulses and for the event-by-event detection of the ions. Double ionization results in a number of prominent fragmentations producing two singly-ionized fragments with kinetic energies of up to a few eV. A number of fragmentations produce ions with four or more carbon atoms, which can only be formed by the breaking of at least three C-C bonds.

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.

Low-Energy Transformation Pathways between Naphthalene Isomers

The transformation pathways between low-energy naphthalene isomers are studied by investigating the topology of the energy landscape of this astrophysically relevant molecule. The threshold algorithm is used to identify the minima basins of the isomers on the potential energy surface of the system and to evaluate the probability flows between them. The transition pathways between the different basins and the associated probabilities were investigated for several lid energies up to 11 eV, this value being close to the highest photon energy in the interstellar medium (13.6 eV). More than a hundred isomers were identified and a set of 23 minima was selected among them, on the basis of their energy and probability of occurrence. The return probabilities of these 23 minima and the transition probabilities between them were computed for several lid energies up to 11 eV. The first connection appeared at 3.5 eV while all minima were found to be connected at 9.5 eV. The local density of state was also sampled inside their respective basins. This work gives insight into both energy and entropic barriers separating the different basins, which also provides information about the transition regions of the energy landscape.

VUV photoprocessing of oxygen-containing polycyclic aromatic hydrocarbons: iPEPICO study of the unimolecular dissociation of ionized benzofuran

Oxygen-containing polycyclic aromatic hydrocarbons (OPAHs) are potential contributors to the 11.3 m band in interstellar observations. To further explore their role in the interstellar medium, we have investigated their fate after photoprocessing by VUV radiation; in particular, we studied the dissociative photoionization of the simplest OPAH, benzofuran, with imaging photoelectron photoion coincidence spectroscopy, iPEPICO. Ionized benzofuran dissociates by loss of CO, followed by a sequential H atom loss. The parallel HCO-loss channel, leading to the same bicyclic C7H5+ fragment ion, is not competitive at low excess energies above the ionization threshold. However, the collision-induced dissociation tandem mass spectrometry results suggest that CO and HCO may be formed in parallel at higher energies. An RRKM fragmentation model reproduced the iPEPICO data well assuming the initial 1,2-H shift transition state to be rate determining to CO loss. The breakdown diagram and the measured dissociation rates agreed well at the CBS-QB3-calculated activation energy of 2.99 eV, which could be relaxed to 3.25 eV, and only a slight adjustment of the ab initio activation entropy. The model barrier to sequential H-loss is larger than the computed H-loss threshold and the breakdown diagram rises less steeply than predicted, which indicates suprastatistical kinetic energy release after the tight H-transfer transition state of the first step. HCO cleavage is possible after a ring-opening transition state, which is looser than and isoenergetic with the CO-loss transition state. However, a subsequent ring formation transition state at 3.85 eV is moderately tight, which suppresses HCO loss at low excess energies.

Survival of polycyclic aromatic hydrocarbon knockout fragments in the interstellar medium

Laboratory studies play a crucial role in understanding the chemical nature of the interstellar medium (ISM), but the disconnect between experimental timescales and the timescales of reactions in space can make a direct comparison between observations, laboratory, and model results difficult. Here we study the survival of reactive fragments of the polycyclic aromatic hydrocarbon (PAH) coronene, where individual C atoms have been knocked out of the molecules in hard collisions with He atoms at stellar wind and supernova shockwave velocities. Ionic fragments are stored in the DESIREE cryogenic ion-beam storage ring where we investigate their decay for up to one second. After 10 ms the initially hot stored ions have cooled enough so that spontaneous dissociation no longer takes place at a measurable rate; a majority of the fragments remain intact and will continue to do so indefinitely in isolation. Our findings show that defective PAHs formed in energetic collisions with heavy particles may survive at thermal equilibrium in the interstellar medium indefinitely, and could play an important role in the chemistry in there, due to their increased reactivity compared to intact or photo-fragmented PAHs.

Structural investigation of doubly-dehydrogenated pyrene cations

The vibrationally resolved spectra of the pyrene cation and doubly-dehydrogenated pyrene cation (C₁₆H₁₀˙⁺; Py⁺ and C₁₆H₈˙⁺; ddPy⁺) are presented. Infrared predissociation spectroscopy is employed to measure the vibrational spectrum of both species using a cryogenically cooled 22-pole ion trap. The spectrum of Py⁺ allows a detailed comparison with harmonic and anharmonic density functional theory (DFT) calculated normal mode frequencies. The spectrum of ddPy⁺ is dominated by absorption features from two isomers (4,5-ddPy⁺ and 1,2-ddPy⁺) with, at most, minor contributions from other isomers. These findings can be extended to explore the release of hydrogen from interstellar PAH species. Our results suggest that this process favours the loss of adjacent hydrogen atoms.

Dimerization of Polycyclic Aromatic Hydrocarbon Molecules and Radicals under Flame Conditions

This work presents a dynamic and kinetic study on the dimerization of polycyclic aromatic hydrocarbon (PAH) molecules and radicals under flame conditions using reactive force field (ReaxFF) molecular dynamics (MD) simulations. The accuracy of the ReaxFF force field is evaluated through comparing with quantum chemistry (QC) calculations of the barrier heights and species concentrations of PAHs reacting with H and OH radicals. A series of homobinary collisions between PAH molecules/radicals are performed to reveal the influence of temperature, molecular size, PAH composition, and the number of radical sites on the dynamics and kinetics of PAH dimerization. Instead of directly forming the strong covalent bonds, the majority of the binary collisions between PAH radicals are bound with weak intermolecular interactions. Effects of oxygen on PAH radical dimerization are also investigated, which indicates that the oxygenated PAH radicals are less likely to contribute to soot nucleation. In addition, the temperature, PAH characteristic, and radical site dependent collision efficiency for PAH radical-radical combinations is extracted from this study.

Atomic hydrogen interactions with gas-phase coronene cations: hydrogenation versus fragmentation

Sequential hydrogenation of polycyclic aromatic hydrocarbon (PAH) cations drives a gradual transition from a planar to a puckered geometry and from an aromatic to an aliphatic electronic structure. The resulting H-induced weakening of the molecular structure together with the exothermic nature of the consecutive H-attachment processes can lead to substantial molecular fragmentation. We have studied H attachment to gas-phase coronene cations in a radiofrequency ion trap using tandem mass spectrometry. With increasing hydrogenation, C2Hi loss and multifragmentation are identified as main de-excitation channels. To understand the dependence of both channels on H-exposure time, we have simulated the molecular stability and fragmentation channels of hydrogenated PAHs using a molecular dynamics approach employing potential energies determined by a density functional based tight binding method. As the coronene fragmentation patterns depend on the balance between energy deposition by H-attachment and the extent of cooling in between subsequent attachment processes, we investigate several scenarios for the energy distribution of hydrogenated PAHs. Good agreement between experiment and simulation is reached, when realistic energy distributions are considered.

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.

Dynamics and kinetics of reversible homo-molecular dimerization of polycyclic aromatic hydrocarbons

Physical dimerization of polycyclic aromatic hydrocarbons (PAHs) has been investigated via molecular dynamics (MD) simulation with the ReaxFF reactive force field that is developed to bridge the gap between the quantum mechanism and classical MD. Dynamics and kinetics of homo-molecular PAH collision under different temperatures, impact parameters, and orientations are studied at an atomic level, which is of great value to understand and model the PAH dimerization. In the collision process, the enhancement factors of homo-molecular dimerizations are quantified and found to be larger at lower temperatures or with smaller PAH instead of size independent. Within the capture radius, the lifetime of the formed PAH dimer decreases as the impact parameter increases. Temperature and PAH characteristic dependent forward and reverse rate constants of homo-molecular PAH dimerization are derived from MD simulations, on the basis of which a reversible model is developed. This model can predict the tendency of PAH dimerization as validated by pyrene dimerization experiments [H. Sabbah et al., J. Phys. Chem. Lett. 1(19), 2962 (2010)]. Results from this study indicate that the physical dimerization cannot be an important source under the typical flame temperatures and PAH concentrations, which implies a more significant role played by the chemical route.

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'.