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

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.

Gas-phase C60Hn+q (n = 0-4, q = 0,1) fullerenes and fulleranes: spectroscopic simulations shed light on cosmic molecular structures

The discovery of C60, C60+, and C70 in the interstellar medium has ignited a profound interest in the astrochemistry of fullerene and related systems. In particular, the presence of diffuse interstellar bands and their association with C60+ has led to the hypothesis that hydrogenated derivatives, known as fulleranes, may also exist in the interstellar medium and contribute to these bands. In this study, we systematically investigated the structural and spectroscopic properties of C60Hn+q (n = 0-4, q = 0,1) using an automated global minimum search and density functional theory calculations. Our results revealed novel global minimum structures for C60H2 and C60H4, distinct from previous reports. Notably, all hydrogenated fullerenes exhibited lower ionization potentials and higher proton affinities compared to C60. From an astrochemical perspective, our results exposed the challenges in establishing definitive spectroscopic criteria for detecting fulleranes using mid-infrared and UV-Vis spectroscopies. However, we successfully identified distinct electronic transitions in the near-infrared range that serve as distinctive signatures of cationic fulleranes. We strongly advocate for further high-resolution experimental studies to fully explore the potential of these transitions for the interstellar detection of fulleranes.

A Spectroscopic View on Cosmic PAH Emission

ConspectusPolycyclic aromatic hydrocarbon molecules (PAHs) are ubiquitously present at high abundances in the Universe. They are detected through their infrared (IR) fluorescence UV pumped by nearby massive stars. Hence, their infrared emission is used to determine the star formation rate in galaxies, one of the key indicators for understanding the evolution of galaxies. Together with fullerenes, PAHs are the largest molecules found in space. They significantly partake in a variety of physical and chemical processes in space, influencing star and planet formation as well as galaxy evolution.Since the IR features from PAHs originate from chemical bonds involving only nearest neighbor atoms, they have only a weak dependence on the size and structure of the molecule, and it is therefore not possible to identify the individual PAH molecules that make up the cosmic PAH family. This strongly hampers the interpretation of their astronomical fingerprints. Despite the lack of identification, constraints can be set on the characteristics of the cosmic PAH family thanks to a joint effort of astronomers, physicists, and chemists.This Account presents the spectroscopic properties of the cosmic PAH emission as well as the intrinsic spectroscopic properties of PAHs and astronomical modeling of the PAH evolution required for the interpretation of the cosmic PAH characteristics. We discuss the observed spectral signatures tracing PAH properties such as charge, size, and structure and highlight the related challenges. We discuss the recent success of anharmonic calculations of PAH infrared absorption and emission spectra and outline the path forward. Finally, we illustrate the importance of models on PAH processing for the interpretation of the astronomical data in terms of the charge balance and PAH destruction.Throughout this Account, we emphasize that huge progress is on the horizon on the astronomical front. Indeed, the world is eagerly awaiting the launch of the James Webb Space Telescope (JWST). With its incredible improvement in spatial resolution, combined with its complete spectral coverage of the PAH infrared emission bands at medium spectral resolution and superb sensitivity, the JWST will revolutionize PAH research. Previous observations could only present spectra averaged over regions with vastly different properties, thus greatly confusing their interpretation. The amazing spatial resolution of JWST will disentangle these different regions. This will allow us to quantify precisely how PAHs are modified by the physical conditions of their host environment and thus trace how PAHs evolve across space. However, this will only be achieved when the necessary (and still missing) fundamental properties of PAHs, outlined in this Account, are known. We strongly encourage you to join this effort.

Interaction of Aromatic Molecules with Forsterite: Accuracy of the Periodic DFT-D4 Method

Density functional theory (DFT) has provided deep atomic-level insights into the adsorption behavior of aromatic molecules on solid surfaces. However, modeling the surface phenomena of large molecules on mineral surfaces with accurate plane wave methods (PW) can be orders of magnitude more computationally expensive than localized atomic orbitals (LCAO) methods. In the present work, we propose a less costly approach based on the DFT-D4 method (PBE-D4), using LCAO, to study the interactions of aromatic molecules with the {010} forsterite (Mg2SiO4) surface for their relevance in astrochemistry. We studied the interaction of benzene with the pristine {010} forsterite surface and with transition-metal cations (Fe2+ and Ni2+) using PBE-D4 and a vdW-inclusive density functional (Dion, Rydberg, Schröder, Langreth, and Lundqvist (DRSLL)) with LCAO methods. PBE-D4 shows good agreement with coupled-cluster methods (CCSD(T)) for the binding energy trend of cation complexes and with PW methods for the binding energy of benzene on the forsterite surface with a difference of about 0.03 eV. The basis set superposition error (BSSE) correction is shown to be essential to ensure a correct estimation of the binding energies even when large basis sets are employed for single-point calculations of the optimized structures with smaller basis sets. We also studied the interaction of naphthalene and benzocoronene on pristine and transition-metal-doped {010} forsterite surfaces as a test case for PBE-D4. Yielding results that are in good agreement with the plane wave methods with a difference of about 0.02-0.17 eV, the PBE-D4 method is demonstrated to be effective in unraveling the binding structures and the energetic trends of aromatic molecules on pristine and transition-metal-doped forsterite mineral surfaces. Furthermore, PBE-D4 results are in good agreement with its predecessor PBE-D3(BJM) and with the vdW-inclusive density functionals, as long as transition metals are not involved. Hence, PBE-D4/CP-DZP has been proven to be a robust theory level to study the interaction of aromatic molecules on mineral surfaces.

Superhydrogenation of pentacene: the reactivity of zigzag-edges

Investigating the hydrogenation of carbonaceous materials is of interest in a wide range of research areas including electronic device development, hydrogen storage, and, in particular, astrocatalytic formation of molecular hydrogen in the universe. Polycyclic Aromatic Hydrocarbons (PAHs) are ubiquitous in space, locking up close to 15% of the elementary carbon. We have used thermal desorption measurements to study the hydrogenation sequence of pentacene from adding one additional H to the fully hydrogenated pentacene species. The experiments reveal that hydrogenated species with an even number of excess H atoms are highly preferred over hydrogenated species with an odd number of H atoms. In addition, the experiments show that specific hydrogenation states of pentacene with 2, 4, 6, 10, 16 and 22 extra H atoms are preferred over other even numbers. We have investigated the structural stability and activation energy barriers for the superhydrogenation of pentacene using Density Functional Theory. The results reveal a preferential hydrogenation pattern set by the activation energy barriers of the hydrogenation steps. Based on these studies, we formulate simple concepts governing the hydrogenation that apply equally well for different PAHs.

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.

Laboratory evidence for the formation of hydrogenated fullerene molecules

Experimental evidence for the formation of hydrogenated fullerene molecules is presented. Films of C60 were grown on a highly oriented pyrolytic graphite (substrate) and exposed to a beam of deuterium atoms. Thermal desorption combined with mass spectrometry was used to determine the deuterated fullerene products formed, revealing a maximum degree of deuteration corresponding to C60D36. Release of D2 from the deuterated C60 film occurs at a much higher temperature than for D-saturated graphite.

Insight into The Nature of Rim Bonds in Coronene

Coronene is known in chemistry as an aromatic or even superaromatic molecule while it has 24 π-electrons which does not conform to the 4 n + 2 Huckel's rule. Chemical bonding description of it is not settled in chemistry and five models were reported in the literature. According to our model, coronene has two concentric π-systems responsible for aromaticity inside of the molecule. In addition to that there are six peripheral 2c-2e π-bonds, which makes coronene aliphatic/aromatic at the same time. However, recent experiments and calculations put in question the presence of peripheral π-bonds. In order to resolve this issue, we computationally studied reaction mechanism of the Cl₂ molecule with C₂H₄, C₆H₆, and C₂₄H₁₂. As it turned out, coronene behaves in a way similar to ethylene by adding Cl₂ molecule. Thus, it proves that coronene indeed has peripheral double bonds which is also consistent with its experimental geometrical features. Our chemical bonding model allows to identify the most reactive atoms in coronene and other PAHs; therefore, it is a matter of importance for combustion and soot formation studies.

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.

The structure of coronene cluster ions inferred from H2 uptake in the gas phase

Mass spectra of helium nanodroplets doped with H₂ and coronene feature anomalies in the ion abundance that reveal anomalies in the energetics of adsorption sites. The coronene monomer ion strongly adsorbs up to n = 38 H₂ molecules indicating a commensurate solvation shell that preserves the D_6h symmetry of the substrate. No such feature is seen in the abundance of the coronene dimer through tetramer complexed with H₂; this observation rules out a vertical columnar structure. Instead we see evidence for a columnar structure in which adjacent coronenes are displaced in parallel, forming terraces that offer additional strong adsorption sites. The experimental value for the number of adsorption sites per terrace, approximately six, barely depends on the number of coronene molecules. The displacement estimated from this number exceeds the value reported in several theoretical studies of the bare, neutral coronene dimer.