Modeling the emission spectra of polycyclic aromatic hydrocarbons by recurrent fluorescence

This paper presents a theoretical statistical model incorporating Herzberg-Teller and Duschinsky rotation effects to calculate recurrent fluorescence rates in PAH cations, revealing that symmetry-forbidden transitions may significantly enhance their cooling efficiency and stability in the interstellar medium.

The Gist

The Cosmic Glow: How Space Molecules Stay Cool

Imagine the vast, cold emptiness of space (the Interstellar Medium). Floating in this darkness are tiny, flat, ring-shaped molecules called Polycyclic Aromatic Hydrocarbons (PAHs). You can think of them as microscopic, multi-layered pancakes made of carbon and hydrogen.

For decades, astronomers have seen a mysterious "glow" in the infrared light coming from these regions. They call these Aromatic Infrared Bands (AIBs). The big question has always been: How do these tiny molecules survive the harsh environment of space, and how do they glow?

Space is full of high-energy ultraviolet (UV) light from stars. When a PAH molecule gets hit by a UV photon, it's like getting punched by a freight train. The molecule absorbs the energy and gets incredibly hot and excited. If it can't cool down fast enough, it might shatter (fragment) and disappear.

The "Recurrent Fluorescence" Escape Trick

This paper focuses on a specific survival trick these molecules use called Recurrent Fluorescence (RF).

Think of the molecule as a hot, bouncing ball.

1. The Hit: A UV photon hits the ball, making it vibrate wildly (internal energy).
2. The Bounce: Usually, the ball just vibrates and slowly releases heat as infrared light (like a cooling cup of coffee).
3. The Trick (RF): Sometimes, the ball is so hot that it doesn't just vibrate; it actually jumps up to a higher "floor" (an excited electronic state) and then immediately jumps back down. When it jumps back down, it releases a flash of visible or near-infrared light. This is Recurrent Fluorescence.

It's like a trampoline. You jump up (absorb energy), land, and then immediately bounce back up to a different height before landing again. This "bounce" helps the molecule dump its excess energy quickly, preventing it from breaking apart.

The New Discovery: The "Silent" Helpers

For a long time, scientists thought only the "loud" jumps (bright, easy-to-see transitions) mattered for this cooling process. They ignored the "silent" jumps because, according to the rules of symmetry, those transitions were supposed to be forbidden (impossible).

The authors of this paper built a new, super-detailed computer model to see what's really happening.

They looked at three specific "pancakes": Naphthalene, Anthracene, and Pyrene.

Here is what they found, using a simple analogy:

• The Old View: Imagine a crowded room where only people wearing bright neon jackets (symmetry-allowed transitions) are allowed to dance and cool the room down.
• The New View: The authors discovered that the people in the dark clothes (symmetry-forbidden transitions) are actually dancing just as hard, if not harder, when the room gets hot enough.

Why?
Because the "dark" dancers start on a lower floor. Even though they are harder to see (they have a weaker "dance signal"), they are much closer to the ground. When the molecule is hot, it's much easier to get to that lower floor. Once there, they can jump back down and release energy.

The paper shows that for these molecules, these "forbidden" jumps might actually be responsible for more than half of the cooling power in certain energy ranges. It's like realizing that the quiet, shy people in the corner are actually the ones doing the heavy lifting to keep the party from exploding.

The Tools They Used

To figure this out, the scientists didn't just guess; they built a statistical model.

• The Microcanonical Ensemble: Imagine trying to count the energy of a single, isolated molecule. It's like trying to predict the exact path of a single leaf in a storm. It's very hard math.
• The Canonical Approximation: To make the math workable, they treated the single molecule as if it were part of a huge crowd of molecules at a specific temperature. It's like predicting the average behavior of a crowd to understand one person. They proved that for these space molecules, this "crowd" math works perfectly, even for a single molecule.

They also included Herzberg-Teller effects and Duschinsky rotation.

• Analogy: Imagine the molecule isn't a rigid plastic toy, but a squishy jelly. When it vibrates, it wobbles and changes shape. This shape-shifting changes how it interacts with light. The authors' model accounts for this "jelly wobble," which previous models ignored. This wobble is what allows the "forbidden" dark transitions to happen at all.

Why This Matters for the Universe

1. Stability: If these molecules can cool down efficiently via these "forbidden" jumps, they are much less likely to shatter when hit by starlight. This explains why we see so many of them in space.
2. The Glow: This new understanding helps explain the specific colors and shapes of the infrared bands astronomers see. The "dark" transitions add a specific glow to the long-wavelength (red/infrared) side of the spectrum that previous models couldn't explain.

The Bottom Line

This paper is like upgrading the map of a city. Previously, scientists thought only the main highways (bright transitions) were used for traffic (energy release). This research shows that the quiet backroads (forbidden transitions) are actually major thoroughfares, especially during rush hour (high energy).

By understanding these backroads, we finally understand how these tiny cosmic pancakes survive the heat of the stars and continue to glow, painting the universe with its famous infrared colors.

Technical Summary

1. Problem Statement

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous in the interstellar medium (ISM) and are hypothesized to be the source of the observed Aromatic Infrared Bands (AIBs). A critical challenge in astrochemistry is understanding how small PAHs survive the harsh ISM environment, particularly against high-energy Vacuum Ultraviolet (VUV) photons (up to 13.6 eV) that can cause fragmentation.

Recurrent Fluorescence (RF), also known as Poincaré fluorescence, is a proposed relaxation mechanism where a PAH cation, after absorbing a photon and undergoing internal conversion (IC) to a hot vibrational ground state, undergoes Inverse Internal Conversion (IIC) to a low-lying electronic excited state, followed by radiative decay (fluorescence) back to the ground state. This process competes with fragmentation and IR emission.

The Gap: Existing theoretical models for RF have historically relied on simplified statistical approaches (e.g., Nitzan and Jortner) that:

• Use the Condon approximation (ignoring geometry-dependent transition dipoles).
• Neglect Herzberg-Teller (HT) effects.
• Ignore vibrational progressions and Duschinsky rotations (mode mixing between electronic states).
• Often assume that symmetry-forbidden transitions (dark states) do not contribute significantly, focusing only on "bright" symmetry-allowed transitions.

The authors aim to develop a rigorous statistical model that accounts for these missing physical effects to accurately predict RF rates and spectra, specifically investigating the role of symmetry-forbidden low-lying states in highly symmetric PAH cations (Naphthalene, Anthracene, Pyrene).

2. Methodology

The authors developed a new theoretical framework based on time-dependent correlation functions within the harmonic approximation, bridging microcanonical and canonical ensembles.

• Theoretical Framework:

  • The RF differential rate constant $k_{RF}(\omega; E)$ is derived from the transition dipole moment time autocorrelation function.
  • Herzberg-Teller (HT) Expansion: The transition dipole moment is expanded linearly around the equilibrium geometry to include both Franck-Condon (FC) and HT contributions. This allows symmetry-forbidden transitions to become active via vibronic coupling.
  • Duschinsky Rotation: The model explicitly accounts for the mixing of normal modes between the ground and excited electronic states via the Duschinsky rotation matrix ($J$) and displacement vector ($K$).
  • Ensemble Equivalence: While the physical system is microcanonical (fixed energy $E$), the authors justify using a canonical ensemble with an effective temperature ($T_{eff}$) to compute the correlation functions, as this is computationally tractable for large systems. They tested two definitions for $T_{eff}$: matching the peak of the distribution ($T_{max}$) or matching the average energy ($T_{av}$), finding $T_{av}$ to be more accurate for small molecules like Naphthalene.

• Computational Inputs:

  • Quantum Chemistry: Calculations were performed using Time-Dependent Density Functional Theory (TD-DFT) with the $\omega$B97X-D functional and cc-pVTZ basis set.
  • Molecules Studied: Cations of Naphthalene ($Np^+$), Anthracene ($An^+$), and Pyrene ($Py^+$), all belonging to the $D_{2h}$ point group.
  • Parameters: Harmonic frequencies, normal mode coordinates, Duschinsky matrices, and transition dipole moments (and their derivatives) were extracted from TD-DFT calculations.
  • Spectral Simulation: The correlation functions were computed via Fast Fourier Transform (FFT) to generate absorption, photoelectron (PE), and RF emission spectra.

3. Key Contributions

1. Novel Statistical Model: The first application of a full correlation-function-based approach to RF that simultaneously includes HT effects, Duschinsky rotations, and full vibrational progressions.
2. Quantification of Forbidden Transitions: The study rigorously demonstrates that low-lying, symmetry-forbidden electronic states (specifically the $D_1$ state) can dominate the RF cooling rate in specific energy regimes, contrary to previous assumptions that only "bright" states ($D_2$) matter.
3. Validation of Canonical Approximation: The paper validates the use of canonical ensemble statistics with an effective temperature ($T_{av}$) to approximate microcanonical RF rates even for relatively small PAHs (like Naphthalene), provided the temperature is calibrated to the internal energy.
4. Spectral Prediction: The model successfully reproduces experimental RF emission spectra (Stokes shifts and broadening) and predicts the existence of significant emission in the near-IR/NIR region (1000–2000 nm) due to the $D_1 \to D_0$ transition.

4. Key Results

• Dominance of the $D_1$ State:

  • For $Np^+$, $An^+$, and $Py^+$, the lowest excited state ($D_1$) is symmetry-forbidden (dark) in the Condon limit. However, HT coupling activates it.
  • Despite having a much lower oscillator strength than the allowed $D_2$ transition, the $D_1 \to D_0$ transition contributes more than 50% to the total RF rate at internal energies below ~6.5 eV.
  • This is because the $D_1$ state has a significantly lower adiabatic energy, leading to a much higher population probability via IIC, which compensates for the weak transition dipole.
  • At higher energies (>6.5 eV), the $D_2$ state eventually dominates due to the $\omega^3$ dependence of the Einstein coefficient favoring higher frequencies.

• Spectral Characteristics:

  • Broadening: The emission spectra are broad and lack distinct vibrational structure due to the high density of vibrational states at the relevant internal energies.
  • Stokes Shift: The model predicts a Stokes shift of ~1200–2000 cm$^{-1}$ between absorption and emission, consistent with experimental observations by Saito et al. and Kusuda et al.
  • Peak Positions: The $D_2 \to D_0$ emission peaks around 13,000–14,000 cm$^{-1}$ (visible/NIR), while the $D_1 \to D_0$ emission peaks at lower frequencies (5,000–7,000 cm$^{-1}$), extending into the near-IR.

• HT Contributions:

  • The HT contribution to the oscillator strength increases with temperature (internal energy).
  • For symmetry-forbidden transitions, the HT mechanism is the only source of emission. The analysis shows that specific low-frequency ungerade modes (breaking inversion symmetry) drive this coupling.

• Comparison with Experiment:

  • The simulated RF spectra for $Np^+$ and $An^+$ match experimental bandpass filter data regarding the shape and the red-shifted maximum relative to absorption.
  • The model explains the "long-wavelength tail" observed in experiments as the contribution of the $D_1$ state and the broadening of the $D_2$ state.

5. Significance and Implications

• PAH Stability in the ISM: The findings suggest that RF is a more efficient cooling mechanism than previously thought. The significant contribution of low-lying forbidden states means PAHs can dissipate energy more effectively via fluorescence rather than fragmentation, enhancing their survival in the harsh, ionized ISM environment.
• Interpretation of AIBs: The model provides a more accurate theoretical basis for interpreting the AIBs and the near-IR continuum, suggesting that forbidden transitions play a non-negligible role in the radiative cooling of interstellar dust.
• Methodological Advancement: The developed framework offers a robust tool for studying relaxation dynamics in large molecular systems where full quantum dynamics are intractable. It bridges the gap between simple statistical models and high-level quantum chemistry.
• Future Directions: The authors propose extending this model to include anharmonic effects (crucial for fragmentation thresholds) and competing processes like isomerization and IR emission to build a complete kinetic model of PAH relaxation cascades.

In conclusion, this work fundamentally revises the understanding of recurrent fluorescence in PAHs by demonstrating that symmetry-forbidden transitions, mediated by Herzberg-Teller coupling, are critical for the thermal stability of interstellar PAHs.