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

Experimental radiative cooling rates of a polycyclic aromatic hydrocarbon cation

Several small Polycyclic Aromatic Hydrocarbons (PAHs) have been identified recently in the Taurus Molecular Cloud (TMC-1) using radio telescope observations. Reproducing the observed abundances of these molecules has been a challenge for astrochemical models. Rapid radiative cooling of PAHs by Recurrent Fluorescence (RF), the emission of optical photons from thermally populated electronically excited states, has been shown to efficiently stabilize small PAHs following ionization, augmenting their resilience in astronomical environments and helping to rationalize their observed high abundances. Here, we use a novel method to experimentally determine the radiative cooling rate of the cation of 1-cyanonaphthalene (C10H7CN, 1-CNN), the neutral species of which has been identified in TMC-1. Laser-induced dissociation rates and kinetic energy release distributions of 1-CNN cations isolated in a cryogenic electrostatic ion-beam storage ring are analysed to track the time evolution of the vibrational energy distribution of the initially hot ion ensemble as it cools. The measured cooling rate is in good agreement with the previously calculated RF rate coefficient. Improved measurements and models of the RF mechanism are needed to interpret astronomical observations and refine predictions of the stabilities of interstellar PAHs.

Cooling dynamics of energized naphthalene and azulene radical cations

Naphthalene and azulene are isomeric polycyclic aromatic hydrocarbons (PAHs) and are topical in the context of astrochemistry due to the recent discovery of substituted naphthalenes in the Taurus Molecular Cloud-1 (TMC-1). Here, the thermal- and photo-induced isomerization, dissociation, and radiative cooling dynamics of energized (vibrationally hot) naphthalene (Np+) and azulene (Az+) radical cations, occurring over the microsecond to seconds timescale, are investigated using a cryogenic electrostatic ion storage ring, affording "molecular cloud in a box" conditions. Measurement of the cooling dynamics and kinetic energy release distributions for neutrals formed through dissociation, until several seconds after hot ion formation, are consistent with the establishment of a rapid (sub-microsecond) Np+ ⇌ Az+ quasi-equilibrium. Consequently, dissociation by C2H2-elimination proceeds predominantly through common Az+ decomposition pathways. Simulation of the isomerization, dissociation, recurrent fluorescence, and infrared cooling dynamics using a coupled master equation combined with high-level potential energy surface calculations [CCSD(T)/cc-pVTZ], reproduce the trends in the measurements. The data show that radiative cooling via recurrent fluorescence, predominately through the Np+ D0 ← D2 transition, efficiently quenches dissociation for vibrational energies up to ≈1 eV above dissociation thresholds. Our measurements support the suggestion that small cations, such as naphthalene, may be more abundant in space than previously thought. The strategy presented in this work could be extended to fingerprint the cooling dynamics of other PAH ions for which isomerization is predicted to precede dissociation.

Cascade Infrared Thermal Photon Emission

The later stages of cooling of molecules and clusters in the interstellar medium are dominated by emission of vibrational infrared radiation. With the development of cryogenic storage it has become possible to experimentally study these processes. Recent storage ring results demonstrate that intramolecular vibrational redistribution takes place within the cooling process, and an harmonic cascade model has been used to interpret the data. Here we analyze this model and show that the energy distributions and the photon emission rates develop into near-universal functions that can be characterized with only a few parameters, irrespective of the precise vibrational spectra and oscillator strengths of the systems. We show that the photon emission rate and emitted power vary linearly with total excitation energy with a small offset. The time developments of ensemble internal energy distributions are calculated with respect to their first two moments. The excitation energy decreases exponentially with a rate constant which is the average of all k_1→0 Einstein coefficients, and the time development of the variance is also calculated.

Efficient stabilization of cyanonaphthalene by fast radiative cooling and implications for the resilience of small PAHs in interstellar clouds

After decades of searching, astronomers have recently identified specific Polycyclic Aromatic Hydrocarbons (PAHs) in space. Remarkably, the observed abundance of cyanonaphthalene (CNN, C₁₀H₇CN) in the Taurus Molecular Cloud (TMC-1) is six orders of magnitude higher than expected from astrophysical modeling. Here, we report unimolecular dissociation and radiative cooling rate coefficients of the 1-CNN isomer in its cationic form. These results are based on measurements of the time-dependent neutral product emission rate and kinetic energy release distributions produced from an ensemble of internally excited 1-CNN⁺ studied in an environment similar to that in interstellar clouds. We find that Recurrent Fluorescence - radiative relaxation via thermally populated electronic excited states - efficiently stabilizes 1-CNN⁺, owing to a large enhancement of the electronic transition probability by vibronic coupling. Our results help explain the anomalous abundance of CNN in TMC-1 and challenge the widely accepted picture of rapid destruction of small PAHs in space.

Radiative cooling rates of substituted PAH ions

The unimolecular dissociation and infrared radiative cooling rates of cationic 1-hydroxypyrene (OHPyr$^+$, \ce{C16H10O+}) and 1-bromopyrene (BrPyr$^+$, \ce{C16H9Br+}) are measured using a cryogenic electrostatic \rev{ion beam} storage ring. A novel numerical approach is developed to analyze the time dependence of the dissociation rate and to determine the absolute scaling of the radiative cooling rate coefficient. The model results show that radiative cooling competes with dissociation below the critical total vibrational energies \revv{$E_c=5.39(1)$}~eV for OHPyr$^+$ and \revv{5.90(1)}~eV for BrPyr$^+$. These critical energies and implications for radiative cooling dynamics are important for astrochemical models concerned with energy dissipation and molecular lifecycles. The methods presented extend the utility of storage ring experiments on astrophysically relevant ions.

Isochronous mass spectrometry in an electrostatic storage ring

For sensitive studies of molecular ions in electrostatic storage rings, the exact knowledge of the isobaric composition of stored beams from a variety of ion sources is essential. Conventional mass-filtering techniques are often inefficient to resolve the beam components. Here, we report the first isochronous mass spectrometry in an electrostatic storage ring, which offers a high mass resolution of Δm/m < 1 × 10^-5 even for heavy molecular species with m > 100 u and uncooled ion beams. Mass contaminations can be resolved and identified at relative fractions down to 0.02%.

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.

Statistical vibrational autodetachment and radiative cooling rates of para-benzoquinone

We report measurements of the statistical vibrational autodetachment (VAD, also called thermionic emission) and radiative cooling rates of isolated para-benzoquinone (pBQ, C₆H₄O₂) radical anions using the cryogenic electrostatic ion storage ring facility DESIREE. The results are interpreted using master equation simulations with rate coefficients calculated using statistical detailed balance theory. The VAD rate is determined by measuring the time-dependent yield of neutral pBQ due to spontaneous electron emission from a highly-excited ensemble of anions formed in an electron-attachment ion source. Competition with radiative cooling quenches the VAD rate after a critical time of τ_c = 11.00(5) ms. Master equation simulations which reproduce the VAD yield provide an estimate of the initial effective vibrational temperature of the ions of 1100(20) K, and provide insight into the anion formation scenario. A second measurement of the radiative cooling rate of pBQ^- stored for up to 0.5 s was achieved using time-dependent photodetachment action spectroscopy across the ²A_u ← ²B_2g and ²B_2u ← ²B_2g transitions. The rate at which hot-band contributions fade from the action spectrum is quantified by non-negative matrix factorisation. This is found to be commensurate with the average vibrational energy extracted from the simulations, with 1/e lifetimes of 0.16(3) s and 0.1602(7) s, respectively. Implications for astrochemistry are discussed.

Radiative relaxation in isolated large carbon clusters: Vibrational emission versus recurrent fluorescence

Recurrent fluorescence (RF) from isolated carbon clusters containing between 24 and 60 atoms is theoretically investigated as a function of internal energy, cluster size, and structural features. The vibrational relaxation kinetics and the associated IR emission spectra are determined by means of a Monte Carlo approach with vibrational density of states computed in the harmonic approximation. RF is generally found to be highly competitive with vibrational emission. The behaviors predicted for clusters of various sizes and archetypal structures indicate that the IR emission spectra are strongly influenced by RF, an energy gap law being obtained for the evolution of the RF rate constant depending on the electronic excitation state. The present results are relevant to the photophysics of the interstellar medium and could contribute to elucidating the carriers of the extended red emission bands and the continuum emission lying below the aromatic infrared bands believed to originate from mixed aromatic–aliphatic compounds.