Ashes of FIRE: Modeling Dust Grain Size Evolution in the Local Group with FIRE

This paper introduces a new discretized grain size evolution model within the FIRE-3 framework that reveals a bimodal dust size distribution and demonstrates how the interplay between dust growth, destruction, and coagulation processes shapes Local Group dust abundances and extinction curve properties, while highlighting the need for alternative mechanisms like top-down PAH formation to explain observed infrared emission features.

The Gist

Imagine the universe as a giant, cosmic kitchen. In this kitchen, stars are the chefs, and the "dust" they create is like the flour, sugar, and spices that make up the atmosphere of galaxies. This dust isn't just dirt; it's the essential ingredient that helps form new stars, regulates the temperature of space, and blocks light, creating the beautiful colors we see in the night sky.

For a long time, scientists trying to simulate how galaxies evolve treated this cosmic dust like a static ingredient—just a fixed amount of flour sitting in a bowl. They didn't really track how the flour changed, got bigger, or got smashed up.

This paper, "Ashes of FIRE," introduces a new, high-tech recipe book (a computer model) that finally tracks the life cycle of every single grain of dust, from its birth to its death. The authors used a powerful simulation code called GIZMO coupled with a physics engine called FIRE (Feedback in Realistic Environments) to see how dust behaves in galaxies like our own Milky Way and its smaller neighbors, the Large and Small Magellanic Clouds.

Here is the story of their findings, explained simply:

1. The Life Cycle of a Dust Grain

Think of a dust grain as a snowball in a blizzard. The paper tracks four main stages of its life:

• Birth: Stars explode (supernovae) or die gently (AGB stars), spewing out fresh "snow" (dust) into space.
• Growth: As this snowball drifts through cold, dense clouds, gas particles stick to it, making it grow bigger (like a snowball rolling down a hill).
• Shattering: If two snowballs collide at high speed in a turbulent wind, they smash into tiny fragments.
• Coagulation: If they collide gently in a calm, dense room, they stick together to form one giant clump.

2. The "Two-Peak" Discovery

The most exciting finding is that the authors discovered a bimodal distribution. Imagine a histogram (a bar chart) of dust sizes.

• Old models predicted a smooth slide: lots of tiny grains, fewer medium ones, and very few big ones (like a smooth hill).
• This new model shows two distinct peaks: a pile of tiny grains and a separate pile of medium-sized grains, with a "valley" in between where very few grains exist.

Why? Because the simulation is so detailed (high resolution) that it can see where things happen.

• The Small Peak: Created when big grains smash apart in violent supernova explosions.
• The Valley: Tiny grains are so eager to grow that as soon as they are born, they quickly grab onto gas and get bigger, skipping the "medium" size.
• The Large Peak: Once grains get big enough, they start sticking together gently in dense clouds, forming the second pile.

3. The "Local Group" Mystery

The authors tested their model on three different "galactic kitchens":

• The Milky Way (MW): A big, rich kitchen with lots of metals.
• The LMC & SMC: Smaller, poorer kitchens with fewer metals.

They found that the amount of dust in these kitchens is mostly determined by how fast new dust grows (accretion) versus how fast it gets destroyed by supernova shocks. This explains why simpler models (that didn't track grain sizes) still got the amount of dust right—they just got lucky by assuming a standard distribution.

However, the shape of the dust (how big the grains are) is different in each galaxy.

• In the smaller galaxies (LMC/SMC), the "gentle sticking" (coagulation) is less efficient because the gas is thinner.
• This means the grains stay small. Small grains block blue light very well, making the light from those galaxies look redder and steeper. The model successfully predicted this!

4. The Missing "PAH" Problem

Here is where the model hit a snag.

• The Expectation: We know from observations that there should be a huge population of ultra-tiny carbon grains (called PAHs, or Polycyclic Aromatic Hydrocarbons). These are like the "cherry on top" that glow in the infrared and create specific chemical signatures.
• The Reality: The model didn't produce enough of these tiny grains. Why? Because in the simulation, as soon as a tiny carbon grain is born, it grows up too fast, swallowing up gas and becoming a medium-sized grain before it can stay tiny.

The Solution? The authors suggest a "Top-Down" theory. Maybe these tiny grains aren't just growing up; maybe they are being made from larger grains that get broken down by UV radiation (like a giant cookie crumbling into crumbs). This "top-down" process might be the missing ingredient needed to keep the population of tiny grains alive.

5. The "Big Grain" Cutoff

The model also explains why we don't see dust grains larger than a certain size (about 0.3 microns).

• Theory A: Maybe big grains get trapped inside dense clouds and never escape to be seen.
• Theory B (The authors' favorite): When supernovae explode, they act like a giant blender, shattering any grain that gets too big. This "blender" effect creates a hard limit on how big dust can get.

The Big Picture

This paper is a major step forward because it treats dust not as a static background, but as a dynamic, living population that changes size and shape depending on its environment.

• Analogy: Previous models were like watching a time-lapse of a city and just counting the total number of cars. This new model is like having a GPS tracker on every single car, seeing which ones crash, which ones get bigger, and how traffic jams (dense clouds) change the flow.

The authors conclude that while their model explains how much dust exists and how big it generally is, there are still mysteries to solve—specifically regarding the tiniest grains (PAHs) and exactly how supernovae destroy the largest ones. But by resolving the "multiphase" nature of the galaxy (seeing the hot, cold, dense, and thin parts separately), they have unlocked a new level of understanding about the dusty universe.

Technical Summary

Here is a detailed technical summary of the paper "Ashes of FIRE: Modeling Dust Grain Size Evolution in the Local Group with FIRE".

1. Problem Statement

Interstellar dust plays a critical role in galaxy evolution (e.g., cooling, star formation, radiation pressure) and significantly affects observational data (extinction, emission). While observations of the Local Group (Milky Way, LMC, SMC) reveal complex variations in dust abundance, chemical composition, and grain size distributions, standard galaxy simulations have historically treated dust as a static quantity or used simplified "two-size" approximations.

• The Gap: Existing models often fail to self-consistently resolve the multiphase Interstellar Medium (ISM), leading to an inability to predict the specific mechanisms driving grain size evolution (shattering, coagulation, accretion) in different ISM phases. Consequently, many models cannot reproduce the observed bimodal grain size distributions or the specific variations in extinction curves (slopes and 2175 Å bump strengths) across different metallicity environments.
• The Challenge: There is a lack of rigorous, high-resolution models that track the evolution of specific dust species (silicates, carbonaceous, metallic iron) while accounting for all major lifecycle processes: creation, growth, destruction, shattering, and coagulation.

2. Methodology

The authors developed a new, discretized grain size evolution model integrated into the GIZMO meshless finite-mass code, coupled with the FIRE-3 stellar feedback and ISM physics suite.

• Discretized Size Distribution: Unlike "two-size" models, this model tracks the grain size distribution ($\partial n / \partial a$) across 16 logarithmically spaced bins ($a_{min} = 5$ Å to $a_{max} = 1$ $\mu$m) for each gas cell.
• Species Tracking: The model separately evolves silicates, carbonaceous dust, and metallic iron, accounting for their distinct physical properties (density, surface energy, charging, sputtering yields).
• Physical Processes Modeled:
  • Creation: Dust injection from Supernovae (SNe) and Asymptotic Giant Branch (AGB) stars with specific size distributions (log-normal for large grains, power-law for small fragments).
  • Growth: Gas-phase metal accretion onto dust grains, limited by temperature cutoffs and Coulomb enhancement factors.
  • Destruction: Thermal sputtering in hot gas, astration (incorporation into stars), and a novel prescription for SNR (Supernova Remnant) destruction that includes grain shattering.
  • Collisional Evolution:
    • Shattering: Occurs in turbulent gas (WIM) and SNRs when relative velocities exceed a threshold, transferring mass from large to small grains.
    • Coagulation: Occurs in dense molecular clouds when relative velocities are low, transferring mass from small to large grains.
  • Sub-resolution Physics: The model incorporates sub-resolution gas clumping (using Mach numbers) to enhance interaction rates and a specific enhancement factor for coagulation in dense gas.
• Simulations: A suite of idealized, non-cosmological disc galaxies representing the Milky Way (MW), Large Magellanic Cloud (LMC), and Small Magellanic Cloud (SMC) masses were run. The authors performed sensitivity tests by varying parameters for accretion, destruction, shattering, and coagulation.

3. Key Contributions

1. First High-Resolution Discretized Model in FIRE: This is the first implementation of a full, discretized grain size evolution model coupled with the high-fidelity FIRE-3 feedback physics, allowing for the resolution of the multiphase ISM.
2. Novel SNR Destruction Prescription: The authors introduce a new analytical prescription for dust destruction in SNRs that explicitly accounts for the shattering of large grains, a process often neglected or simplified in previous works.
3. Bimodal Distribution Prediction: The model uniquely predicts a bimodal grain size distribution (peaks at $\sim$5 nm and $\sim$0.1 $\mu$m) arising naturally from the spatial separation of ISM phases (shattering in hot/turbulent gas vs. coagulation in cold/dense gas), contrasting with the power-law distributions predicted by lower-resolution models.
4. Mechanistic Explanation of Local Group Variations: The study provides a physical explanation for the differences in dust properties between the MW, LMC, and SMC based on metallicity-dependent accretion and density-dependent coagulation efficiencies.

4. Key Results

• Dust Abundance & Composition:

  • Local Group dust abundances (Dust-to-Metal ratio, D/Z) are primarily determined by the balance between accretion (growth) and SNR destruction.
  • Variations in D/Z and elemental depletion (Si, Mg, C, Fe) across the MW, LMC, and SMC are driven by metallicity-dependent accretion rates.
  • Coagulation and shattering have minimal impact on total dust abundance, explaining why simpler models (assuming constant MRN distributions) can still match abundance observations.

• Grain Size Distribution (Bimodality):

  • The model produces a robust bimodal distribution: a small-grain peak ($\sim$5 nm) formed by shattering and subsequent accretion, and a large-grain peak ($\sim$0.1 $\mu$m) formed by coagulation.
  • This bimodality is a direct result of resolving the multiphase ISM, where shattering and coagulation occur in distinct environments.

• Extinction Curves:

  • Slope: The steepness of the UV-optical extinction curve is determined by the efficiency of coagulation. Inefficient coagulation (as in low-metallicity LMC/SMC) leaves more small silicate grains, resulting in steeper slopes.
  • 2175 Å Bump: The strength of the UV bump is determined by the mass fraction of small carbonaceous grains. The model predicts that inefficient coagulation leads to stronger bumps (more small carbonaceous grains).
  • Discrepancy: While the model reproduces the steeper slopes of the LMC/SMC, it fails to reproduce the observed weak/absent bumps in these galaxies. The model predicts too many small carbonaceous grains because they grow too efficiently via accretion.

• PAH Evolution & The "Missing" Small Grains:

  • The model fails to predict a population of very small ($<1$ nm) carbonaceous grains required for Mid-IR (MIR) emission features. In the simulation, these grains rapidly grow via accretion in the Cold Neutral Medium (CNM).
  • Implication: The authors suggest a "top-down" formation mechanism is necessary, where pre-existing grains are converted into PAHs by UV radiation in the diffuse ISM, inhibiting their growth via accretion.

• Large Grain Cutoff:

  • The sharp cutoff of grains larger than $0.3$$\mu$m is driven by efficient shattering in SNRs. Without this mechanism, large grains would survive and accumulate.

5. Significance and Implications

• Resolution Dependence: The paper highlights that the shape of the grain size distribution (bimodal vs. power-law) and the abundance of very small grains are highly dependent on simulation resolution. Low-resolution models that sub-grid the ISM cannot capture the spatial separation of dust processes, leading to incorrect predictions for extinction curves and PAH populations.
• PAH Formation: The inability of current "bottom-up" growth models to explain the lack of small carbonaceous grains in low-metallicity environments suggests that top-down PAH formation (grain fragmentation/aromatization) is a critical, missing ingredient in galaxy evolution models.
• Future Observations: The authors argue that extinction curves alone are degenerate regarding grain size distribution shapes. Future constraints must rely on dust emission spectra and polarization data to distinguish between bimodal and power-law distributions.
• High-Redshift Predictions: The findings suggest that dust properties in high-redshift, high-density galaxies may differ significantly from the Local Group due to enhanced coagulation, potentially altering the attenuation curves used to interpret early galaxy observations (e.g., with JWST).

In summary, this work represents a significant step forward in linking hydrodynamical galaxy simulations with detailed dust physics, revealing that the multiphase nature of the ISM is the key driver of the observed bimodal grain size distribution and highlighting the need for new physics (top-down PAH formation) to explain the smallest dust grains.