Deexcitation dynamics of superhydrogenated polycyclic aromatic hydrocarbon cations after soft-x-ray absorption

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.

Unimolecular fragmentation and radiative cooling of isolated PAH ions: A quantitative study

Time-resolved spontaneous and laser-induced unimolecular fragmentation of perylene cations (C₂₀H₁₂ ⁺) has been measured on timescales up to 2 s in a cryogenic electrostatic ion beam storage ring. We elaborate a quantitative model, which includes fragmentation in competition with radiative cooling via both vibrational and electronic (recurrent fluorescence) de-excitation. Excellent agreement with experimental results is found when sequential fragmentation of daughter ions co-stored with the parent perylene ions is included in the model. Based on the comparison of the model to experiment, we constrain the oscillator strength of the D₁ → D₀ emissive electronic transition in perylene (f^RF = 0.055 ± 0.011), as well as the absolute absorption cross section of the D₅ ← D₀ excitation transition (σ_abs > 670 Mb). The former transition is responsible for the laser-induced and recurrent fluorescence of perylene, and the latter is the most prominent in the absorption spectrum. The vibrational cooling rate is found to be consistent with the simple harmonic cascade approximation. Quantitative experimental benchmarks of unimolecular processes in polycyclic aromatic hydrocarbon ions like perylene are important for refining astrochemical models.

Ion Dissociation Dynamics of 1,2,3,4-Tetrahydronaphthalene: Tetralin as a Test Case For Hydrogenated Polycyclic Aromatic Hydrocarbons

The unimolecular dissociation of ionized tetralin was probed by tandem mass spectrometry, imaging photoelectron photoion coincidence (iPEPICO) spectroscopy, and theory. The major reactions observed were the loss of the hydrocarbons CH₃^•, C₂H₄, and C₃H₅^• together with H^•-atom loss. RRKM modeling of the iPEPICO data suggested a two-well potential energy surface. Ionized tetralin can lose all four neutrals via H-shift and ring-opening reactions or CH₃^• and C₂H₄ after interconversion to the 1-methylindane ion, a process similar to that found for ionized 1,2-dihydronaphthalene (isomerizing to form the 1-methylindene ion structure). This was confirmed at the B3LYP/6-31+G(d,p) level of theory, and potential mechanisms for all reactions are described. The ionization energy of tetralin was established from the threshold photoelectron spectrum to be 8.46 ± 0.01 eV.

Storage time dependent photodissociation action spectroscopy of polycyclic aromatic hydrocarbon cations in the cryogenic electrostatic storage ring DESIREE

The intrinsic absorption profile and radiative cooling rate of coronene cations are reported.

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.

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.

Structure and Dissociation Pathways of Protonated Tetralin (1,2,3,4-Tetrahydronaphthalene)

The infrared multiple-photon dissociation (IRMPD) spectrum of protonated tetralin (1,2,3,4-tetrahydronaphthalene, THN) has been recorded using an infrared free electron laser coupled to a Fourier transform ion cyclotron mass spectrometer. IR-induced fragmentation of the protonated parent [THN + H] ⁺, m/z 133, yielded a single fragment ion at m/z 91. No evidence for fragment ions at m/z 131 or 132 was observed, indicating that protonated THN ejects neither atomic H nor molecular H₂. Comparison of the experimental spectrum with density functional calculations (B3LYP/6-311++G(d,p)) of the two possible protonated isomers identifies a preference for the position of protonation. Possible decomposition pathways starting from both [THN + H(5)]⁺ and [THN + H(6)]⁺ are investigated. The potential energy profiles computed for these decomposition routes reveal that (1) the m/z 91 ionic product resembles the benzylium ion, but with the extra hydrogen and the methylene substituents in various ortho, meta, and para conformations around the aromatic ring and that (2) the decomposition process involving the [THN + H(6)]⁺ isomer is predominant, while the one involving the [THN + H(5)]⁺ may play a smaller role. Potential energy pathways from the initial decomposition product(s) to the benzylium and tropylium ions have also been computed. Given the relatively low barriers to these ions, it is concluded that the benzylium ion and, with sufficient activation, the tropylium ion plus neutral propene are the final products.

The sequence to hydrogenate coronene cations: A journey guided by magic numbers

The understanding of hydrogen attachment to carbonaceous surfaces is essential to a wide variety of research fields and technologies such as hydrogen storage for transportation, precise localization of hydrogen in electronic devices and the formation of cosmic H₂. For coronene cations as prototypical Polycyclic Aromatic Hydrocarbon (PAH) molecules, the existence of magic numbers upon hydrogenation was uncovered experimentally. Quantum chemistry calculations show that hydrogenation follows a site-specific sequence leading to the appearance of cations having 5, 11, or 17 hydrogen atoms attached, exactly the magic numbers found in the experiments. For these closed-shell cations, further hydrogenation requires appreciable structural changes associated with a high transition barrier. Controlling specific hydrogenation pathways would provide the possibility to tune the location of hydrogen attachment and the stability of the system. The sequence to hydrogenate PAHs, leading to PAHs with magic numbers of H atoms attached, provides clues to understand that carbon in space is mostly aromatic and partially aliphatic in PAHs. PAH hydrogenation is fundamental to assess the contribution of PAHs to the formation of cosmic H₂.

Threshold energies for single-carbon knockout from polycyclic aromatic hydrocarbons

We have measured absolute cross sections for ultrafast (femtosecond) single-carbon knockout from polycyclic aromatic hydrocarbon (PAH) cations as functions of He–PAH center-of-mass collision energy in the 10–200 eV range. Classical molecular dynamics (MD) simulations cover this range and extend up to 105 eV. The shapes of the knockout cross sections are well-described by a simple analytical expression yielding experimental and MD threshold energies of EthExp = 32.5 ± 0.4 eV and EthMD = 41.0 ± 0.3 eV, respectively. These are the first measurements of knockout threshold energies for molecules isolated in vacuo. We further deduce semiempirical (SE) and MD displacement energies, i.e., the energy transfers to the PAH molecules at the threshold energies for knockout, of TdispSE = 23.3 ± 0.3 eV and TdispMD = 27.0 ± 0.3 eV. The semiempirical results compare favorably with measured displacement energies for graphene (Tdisp = 23.6 eV).