PAH Emission Power Spectra

• PAH emission power spectra are quantitative descriptions of PAH infrared features that trace molecular size, charge, and local interstellar medium conditions.
• Methodologies include spectral decomposition using quantum-statistical models or fitted Drude profiles and spatial analysis via Fourier transforms to determine power-law slopes.
• These spectra offer diagnostic insights into star formation, PAH survival in extreme environments, and feedback-driven structural scales in galaxies.

Polycyclic Aromatic Hydrocarbon (PAH) emission power spectra describe the distribution of infrared (IR) power emitted by PAH molecules following the absorption of UV/optical photons, generally expressed as a function of wavelength or spatial frequency in astrophysical environments. PAHs are the primary carriers of a family of broad and intense emission features at key mid-IR wavelengths (notably 3.3, 6.2, 7.7, 8.6, and 11.2 μm), which dominate the IR spectra of a wide range of galactic and extragalactic sources and encode information about ISM conditions, molecular size/charge, local radiation fields, and physical structure. Recent advances, notably with JWST and the PHANGS survey, have enabled direct measurement of both spectral-energy and spatial power spectra of PAH emission at tens-of-pc resolution, providing quantitative constraints on the structure, excitation, and evolution of PAH carriers in diverse environments (Lind-Thomsen et al., 1 Apr 2025, Draine et al., 2020, Richie et al., 19 Oct 2025).

1. Fundamental Properties of PAH Emission Spectra

PAH emission spectra are characterized by discrete vibrational bands at infrared wavelengths, originating from stretching and bending modes of C–H and C–C bonds. The primary bands are located at 3.3 μm (CH stretch), 6.2 μm (CC stretch), 7.7 and 8.6 μm (CC stretch and CH in-plane bend), and 11.2 μm (CH out-of-plane bend) (Peeters, 2011). Each band's relative intensity, central wavelength, and profile trace the PAH’s charge state, size, edge structure, and local excitation. Typical normalized ratios for “normal” ISM conditions are: $I_{3.3}/I_{6.2}$ ≈ 0.05–0.10, $I_{11.2}/I_{6.2}$ ≈ 0.4–0.6, $I_{7.7}/I_{6.2}$ ≈ 2–5. The aggregate of these bands contributes ~10–15% of a star-forming galaxy’s mid-IR power (Peeters, 2011, Lai et al., 2020).

PAH emission power spectra can be analyzed in two principal domains: (1) wavelength (spectral) space, describing the vibrational feature energy distribution per molecule or per unit volume, and (2) spatial-frequency space, describing how PAH emission is structured across different physical scales in galaxies (Lind-Thomsen et al., 1 Apr 2025).

2. Modeling and Measurement Methodologies

Spectral Decomposition and Single-Photon Approximation

PAH emission modeling employs either (a) quantum-statistical calculations for single-photon heating and cooling, or (b) empirical or fitted Drude/Lorentzian profiles for observed spectra (Richie et al., 19 Oct 2025, Draine et al., 2020, Lai et al., 2020). In the single-photon limit (SPA), appropriate for small PAH grains and interstellar radiation fields with intensity $U<100$, each absorption/emission event is treated independently. The time-integrated emission spectrum $P(\lambda)$ is constructed as a weighted sum over “basis spectra” $S(\lambda, E)$ for absorbed photon energies $E$, further integrated over the ambient photon flux and PAH size/charge distributions:

$$P(\lambda) = \int F(E)\, \sigma_{\rm abs}(E)\, E\, S(\lambda, E)\, dE$$

This formalism matches multi-photon treatments to better than 10% for 3–20 μm (Richie et al., 19 Oct 2025).

Spatial Power Spectra (Fourier Analysis)

PHANGS–JWST MIRI imaging enables spatial power-spectrum analysis at tens-of-pc scales. The 2D spatial power spectrum is computed via the Fourier transform of the intensity map $I(x, y)$, yielding:

$P(k_x, k_y) = |\mathcal{F}(k_x, k_y)|^2$

Radial averaging over annuli in ($k_x$, $k_y$) space produces a 1D $P(k_r)$. The spectrum is described by a power law $P(k_r) = b\,k_r^{-\alpha}$ across accessible scales, with the power-law slope $\alpha$ quantifying the relative dominance of large versus small-scale structure (Lind-Thomsen et al., 1 Apr 2025).

3. Physical Drivers and Band Ratio Diagnostics

Dependence on PAH Charge, Size, and Radiation Field

The distribution of emitted IR power among vibrational bands depends sensitively on PAH charge state, molecular size, and the hardness/intensity of the ambient radiation field. Cationic PAHs radiate strongly at 6.2, 7.7, and 8.6 μm, with the 6.2/11.2 or 7.7/11.2 band ratios serving as robust diagnostics of the “ionization parameter” ($\gamma = G_0 \sqrt{T_{\rm gas}}/n_e$) (Sidhu et al., 2023, Draine et al., 2020).

Smaller PAHs ($N_C \lesssim 10^2$) undergo larger temperature excursions and emit preferentially in short-λ bands (3.3 μm), whereas larger PAHs ($N_C \gtrsim 10^3$) favor long-λ bands (17 μm plateau). The ratio $P_{7.7}/P_{\rm IR}$ for a standard PAH population varies by a factor of ≈2.5 between hard (young starburst) and soft (M31 bulge) radiation fields; $P_{3.3}/P_{\rm IR}$ is even more sensitive and can fluctuate by a factor of 6 over the same range (Draine et al., 2020).

Observational Band-Power Fractions

In local star-forming galaxies, the 3.3 μm feature contributes ≈1.5–3% of the total PAH power, while the 7.7 μm and 11.3 μm bands together dominate the emission [$f_{3.3} \approx 0.015-0.03$, total $L_{\Sigma \rm PAH} / L_{\rm IR} \approx 5-10\%$; (Lai et al., 2020)]. Ratio trends are modified at high obscuration or low metallicity; in the dwarf galaxy II Zw 40, $f_{3.3}$ reaches ≈2.5% even as global PAH power is suppressed, indicating survival or effective reformation of small PAHs in hard radiation fields.

4. Spatial Power Spectrum Results from the JWST Era

MIRI observations of 13 PHANGS galaxies reveal that PAH emission (7.7 and 11.3 μm) displays a spatial power spectrum with a steeper slope ($\alpha_{7.7}=2.19^{+0.16}_{-0.15}$, $\alpha_{11.3}=1.88^{+0.25}_{-0.37}$) relative to thermal dust continuum bands ($\alpha_{10.0\mu m}=1.48^{+0.33}_{-0.47}$, $\alpha_{21.0\mu m}=0.94^{+0.23}_{-0.28}$) (Lind-Thomsen et al., 1 Apr 2025). Notably, only PAH-dominated bands display a statistically significant break in the power spectrum at a characteristic scale $\ell_0 = 160^{+110}_{-50}$ pc. Below $\ell_0$, power is suppressed, indicating a departure from simple scale-free/turbulent structure, attributed to the confinement of PAH emission to photo-dissociation regions (PDRs) where stellar feedback and UV fields truncate small-scale structure.

Furthermore, the distributions of $\alpha$ are narrower in the PAH bands, with $\sigma_\alpha \sim 0.2$, implying a uniform regulatory mechanism (e.g., stellar feedback) across local star-forming disks, in contrast to the diverse dust continuum morphologies ($\sigma_\alpha \sim 0.4$) (Lind-Thomsen et al., 1 Apr 2025).

5. Practical Implementation: Band Decomposition and Attenuation

PAH emission spectra are typically decomposed into a sum of Drude or Lorentzian components for each fundamental and overtone, as well as broad plateau features (Lai et al., 2020, Peeters, 2011). Accurate band fluxes require detailed continuum subtraction, careful accounting for absorption features (notably silicate at 9.7 μm), and consideration of attenuation geometry (mixed vs. obscured-continuum, which can alter recovered $L_{3.3}$ by up to a factor of ≈5 under high extinction).

In practice, the total PAH power is expressed as

$L_{\Sigma\,{\rm PAH}} = \sum_i L_{\rm PAH},i$

where $L_{\rm PAH},i$ is the luminosity integrated over band $i$. Band fractions $f_i = L_{\rm PAH},i / L_{\Sigma\,{\rm PAH}}$ provide an empirical “power-spectrum” view of the relative energetic weights among features (Lai et al., 2020).

Spatial power spectra require imaging data, Fourier transformation, PSF correction, and radial binning. The resultant $P(k_r)$ is fit to identify power-law regimes, locate breaks ($\ell_0$), and compare structure between emission components (Lind-Thomsen et al., 1 Apr 2025).

6. Applications and Astrophysical Implications

PAH emission power spectra function as sensitive diagnostics for ISM conditions, star-formation rates, and feedback processes across cosmic time. The 7.7 μm band serves as a robust estimate of the global PAH-to-dust mass ratio ($q_{\rm PAH}$), while the 3.3/11.2 and 6.2/11.2 μm ratios probe PAH size and ionization (Draine et al., 2020, Sidhu et al., 2023). In high-z galaxies and extreme starbursts, strong PAH emission indicates prodigious star formation even when AGN signatures dominate the continuum (Rawlings et al., 2012).

The spatially resolved suppression of PAH structure below ≈160 pc in PHANGS galaxies ties this scale to molecular cloud and PDR morphologies, establishing a feedback-regulated floor to structural self-similarity (Lind-Thomsen et al., 1 Apr 2025). Conversely, the absence of clear scale breaks in thermal continuum suggests a more diverse, integrated dust-emission morphology across galactic histories.

7. Extended Spectral Features and Future Directions

Spectroscopy in the 1–5 μm range predicts a rich, yet faint, suite of overtone, combination, and functional group signature bands. These “pseudo-continuum” structures encode elemental fractionation (e.g., D/H, N/C) and can now be accessed with JWST/NIRSpec (Allamandola et al., 2021). Band ratios of weak CD-stretch (4.3–4.8 μm) to CH-stretch (3.3 μm), or cyano-PAH C≡N features (4.46–4.50 μm), provide direct measures of the chemical evolution of the PAH population. The persistence of small PAH emission at high luminosity and low metallicity, as observed empirically (Lai et al., 2020), suggests robust mechanisms for PAH survival or reformation under intense ISM conditions.

Recent developments in computational tools (e.g., “pah_spec” package) support rapid generation and fitting of PAH spectra for arbitrary radiation fields and grain populations, which, combined with high-resolution imaging and spectroscopy, will enable detailed mapping of PAH evolution and ISM feedback across the local and distant universe (Richie et al., 19 Oct 2025).