The Lifecycle and Emission Properties of PAHs in Cosmological Hydrodynamic Galaxy Formation Simulations
summary
The gist
The gist The first cosmological model for the lifecycle and luminous properties of polycyclic aromatic hydrocarbons (PAHs) in galaxies as they evolve from z = 6 → 0 is presented by coupling
In short
The model tracks PAH evolution from high redshift (z=6) to today by simulating galaxy growth and dust grain shattering in the interstellar medium. It shows that as galaxies become more diffuse, PAH production increases, linking metallicity to PAH abundance. This suggests grain-grain shattering is a primary driver of cosmic PAH evolution.
Key concepts
- PAHs
- Polycyclic Aromatic Hydrocarbons are ultrasmall (less than 13 Angstrom) carbonaceous dust grains. The model assumes these grains form in situ through the shattering of larger dust grains within the interstellar medium, which is driven by changes in gas density over cosmic time.
- Grain-Grain Shattering
- This process describes how larger dust particles break apart into smaller ones. In this model, the rate of shattering increases as galaxies evolve; as the interstellar medium becomes more diffuse at lower redshifts, collision velocities rise, leading to a greater fraction of ultrasmall PAHs.
- PAH Mass Fraction (qPAH)
- This is the relative mass of PAHs compared to total dust. The model finds that qPAH rises significantly from about 5 x 10^-4 at z~4 to around 10^-2 at z~0, indicating a substantial increase in PAH abundance as the universe ages.
- PAH-Metallicity Relationship (PZR)
- This relationship connects the metallicity of a galaxy to its PAH abundance. The model naturally reproduces this by showing that galaxies enrich their diffuse interstellar medium and grow their PAH fraction simultaneously, linking rising metallicity directly to higher qPAH.
Terminology used across episodes
This episode discusses
- The Lifecycle and Emission Properties of PAHs in Cosmological Hydrodynamic Galaxy Formation Simulations · Paper Radio
- Ashes of FIRE: Modeling Dust Grain Size Evolution in the Local Group with FIRE
- Galaxies with grains: unraveling dust evolution and extinction curves with hydrodynamical simulations
- Silicate cosmic dust grain collisions in the interstellar medium: A molecular dynamics study · Paper Radio
- The Astrodust+PAH Model: A Unified Description of the Extinction, Emission, and Polarization from Dust in the Diffuse Interstellar Medium
- SPHEREx mapping of diffuse PAH and H II emission in the Galactic plane
- Detection of the 2175 UV Bump at z>7: Evidence for Rapid Dust Evolution in a Merging Reionisation-Era Galaxy
- Dust enrichment from core-collapse supernovae and extinction curves in the high-redshift universe
- Dust and Grain Size Evolution in Galaxy Simulations: What Matters and What Does Not
- PAH Emission Spectra and Band Ratios for Arbitrary Radiation Fields with the Single Photon Approximation
- CALIMA: On-the-fly dust and PAH evolution for radiation-hydrodynamics galaxy formation simulations
- The Metallicity Dependence of PAH Emission in Galaxies II: Insights from JWST/NIRCam Imaging of the Smallest Dust Grains in M101
The paper
The Lifecycle and Emission Properties of PAHs in Cosmological Hydrodynamic Galaxy Formation Simulations · Read on arXiv
DESIKA NARAYANAN, PAUL TORREY, MASSIMILIANO PARENTE, GRANT P. DONNELLY, DHRUV T. ZIMMERMAN, ALEX M. GARCIA, HELENA M. RICHIE, J.-D. T. SMITH, BRANDON S. HENSLEY, FEDERICO MARINACCI, JED MCKINNEY, ALEXANDRA POPE, GERGO POPPING, LAURA V. SALES, KARIN SANDSTROM, ETHAN SAVITCH
Department of Astronomy, University of Florida · Cosmic Dawn Center at the Niels Bohr Institute, University of Copenhagen and DTU-Space Technical University of Denmark Department of Astronomy, University of Virginia Virginia Institute for Theoretical Astronomy, University of Virginia The NSF-Simons AI Institute for Cosmic Origins Ritter Astrophysical Research Center Department of Physics and Astronomy University of Toledo Physics and Astronomy Department University of Pittsburgh Jet Propulsion Laboratory California Institute of Technology Department of Physics and Astronomy Augusto Righi University Bologna INAF Astrophysics and Space Science Observatory Bologna Department of Astronomy University of Texas at Austin Department of Astronomy University of Massachusetts Amherst European Southern Observatory European Southern Observatory Karl-Schwarzschild-Str. Garching Germany
We present the first cosmological model for the lifecycle and luminous properties of PAHs in galaxies as they evolve from z=6-->0. We model 40 zoom-in galaxies, coupled with an on-the-fly model for the evolution of dust grains in the ISM. We approximate PAHs as ultrasmall (a < 13 Angstrom) carbonaceous dust grains, and couple this model with single-photon excitation calculations to compute the emergent mid-infrared spectra. (1) If we assume that dust is large upon formation, then PAHs are naturally able to form in situ in the ISM via grain-grain shattering. Interstellar collision velocities increase in low density, diffuse gas in our model; as galaxies evolve, the increase in fractional mass of diffuse gas drives an increase in grain-grain collision velocities and a corresponding rise in the PAH mass fraction (qPAH) from 5 x 10-4 at z 4 to 10-2 at z 0. (2) Increased PAH production in the diffuse ISM results in an inverse relationship between qPAH and the molecular gas fraction. (3) The PAH light-to-mass ratio scales linearly with the radiation field intensity (LPAH/MPAH G 0) but anti-correlates with qPAH, because high-Sigma SFR galaxies have a denser ISM that suppresses shattering. This means the physical qPAH and observed LPAH/LFIR do not evolve in lockstep. (4) The PAH-metallicity relationship (PZR) arises naturally in this framework: galaxies enrich and grow their diffuse ISM fraction simultaneously, linking rising metallicity to rising qPAH. Our models represent the first to reproduce the PZR observed across z=0-2. (5) The LPAH-SFR and LPAH-M mol relations emerge from two effects: more massive galaxies have larger PAH reservoirs, and higher-SFR galaxies excite their PAHs more efficiently per unit mass. Taken together, these results suggest that grain-grain shattering in the diffuse ISM is the main driver behind the evolution of cosmic PAH abundances.
Transcript
Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.
Vera: We're looking at the paper: "The Lifecycle and Emission Properties of PAHs in Cosmological Hydrodynamic Galaxy Formation Simulations".
Jocelyn: Here's the gist of the paper. The gist The first cosmological model for the lifecycle and luminous properties of polycyclic aromatic hydrocarbons (PAHs) in galaxies as they evolve from z = 6 → 0 is presented by coupling cosmological hydrodynamic zoom-in simulations with an on-the-fly model for dust grain evolution and single-photon excitation calculations.
Vera: How it works The research employs a framework where PAHs are assumed to be ultrasmall (a < 13 A) carbonaceous dust grains that form predominantly via grain-grain shattering.
Jocelyn: This process is driven by the evolution of ISM physical conditions, specifically the local gas density, which is the single most important factor controlling shattering rates.
Vera: As galaxies evolve across cosmic time, typical ISM conditions shift from dense and molecular at high redshift to increasingly diffuse at low redshift.
Paper summary: Vera: To recap, this paper models the lifecycle of PAHs from z=six to zero using cosmological hydrodynamic simulations and an on-the-fly model for dust grain evolution <ref:2606.20809#pg1>. Their central thesis is that grain-grain shattering in the diffuse ISM is the main driver behind the evolution of cosmic PAH abundances across cosmic time.
Jocelyn: They claim that by assuming PAHs are ultrasmall carbonaceous grains, smaller than thirteen Angstroms, they can naturally explain why we see a rise in PAH mass fraction from about five times ten to the minus four at redshift four to about ten to the minus two at redshift zero.
Subrahmanyan: The paper claims this shattering mechanism provides a unified physical framework that accounts for three different scaling relations: the PZR, the LPAH over SFR relation, and the LPAH over molecular gas relations simultaneously across all redshifts.
Vera: They also find that increased PAH production in diffuse gas results in an inverse relationship between qPAH and the molecular gas fraction of a galaxy. This is because more material is being pushed into that diffuse phase where shattering is more active.
Jocelyn: The paper also highlights a decoupling between the physical mass fraction of PAHs and their observed luminosity-to-mass ratio, showing that they don't evolve in lockstep due to the radiation field intensity.
Subrahmanyan: It matters because it connects how galaxies build up their chemical content, through metallicity, to how they emit light, through PAH properties. If this model holds up across redshifts like z=zero to two, it gives us a strong physical mechanism for why we see these specific scaling relations in the data.
Vera: The authors show that more massive galaxies have larger PAH reservoirs and higher star formation rate galaxies excite their PAHs more efficiently per unit mass. But they ultimately point to the shattering as the main driver of evolution.
Jocelyn: So, it's not just one thing—it’s the interplay between gas dynamics, radiation fields, and grain evolution that controls how these molecules appear in cosmic history. That's what this work is about.
Conclusion: Vera: Looking at the full title, "The Lifecycle and Emission Properties of PAHs in Cosmological Hydrodynamic Galaxy Formation Simulations," it really tells you that this paper is about tracking the whole story of these molecules from their formation through their emission properties across cosmic time.
Jocelyn: And it’s important to remember that the authors are Desika Narayanan and his team, working on this complex simulation framework. Their work provides a physical picture linking galaxy structure formation directly to PAH abundance evolution.
Subrahmanyan: The implication for us is that we have a mechanism—grain-grain shattering in diffuse gas—that could explain why we see the observed PAH-metallicity relationship across redshifts zero to two. It moves this relationship from being just an observation to being linked to a physical process in how galaxies build their interstellar medium.
Vera: So, for someone listening who only watches the show, what this means is that when you look at galaxy data from different epochs, you shouldn't just see correlations; you should be thinking about the underlying physical process of grain shattering happening in the diffuse gas.
Jocelyn: It means future observations need to be sensitive enough to probe those diffuse phases better, because if this model is right, we need to see how that change in gas density and collision velocity affects the PAH light we observe.
Subrahmanyan: And for confirming this whole picture, Subrahmanyan would say that a key observation needed is one that can directly measure the changes in grain size distribution or collision velocities within those diffuse ISM phases.
Vera: That seems to be the direction they are pushing—connecting the simulation output back to observable quantities like the observed PAH emission features. That's what makes this kind of modeling useful for us as observational astronomers.
More episodes
- 2508.15337-Hierarchical Black Hole Mergers in AGN Disks: Tracing Massive Black Hole Growth Across Cosmic Time
- 2512.04600-The dynamical memory of tidal stellar streams: Joint inference of the Galactic potential and the progenitor of GD-1 with flow matching
- 2602.22204-Recovering 3D Magnetic Turbulence from Single-frequency Faraday Screens
- 2603.01472-Diversity in Lyman Continuum Escape at z about0.3 Revealed by WISE Infrared Observations
- 2603.08513-ML in Astrophysical Turbulence I: Predicting Prestellar Cores in Magnetized Molecular Clouds using eXtreme Gradient Boosting
- 2604.09539-The missing ultra-faint satellites of the Milky Way
- 2604.11904-JOYS: Launch and destruction of dust in protostellar jets. The case of BHR71-IRS1 with JWST MIRI
- 2605.09829-Stardust Galaxies at z>9: A Dust-Origin Transition Behind the Excess of UV-Bright Galaxies
- 2605.22163-A Unified HI Rotation Curve Database for 122 Dwarf and Irregular Galaxies
- 2606.04081-The role of major mergers in triggering super-Eddington accretion