---
title: PNe in Cosmological Simulations (PICS)
url: https://www.emergentmind.com/topics/pne-in-cosmological-simulations-pics
type: topic
---

# PNe in Cosmological Simulations (PICS)

PICS, abbreviated from **PNe In Cosmological Simulations**, is a modeling framework for predicting planetary nebula (PN) populations directly within hydrodynamical cosmological simulations, thereby linking PN formation to realistic stellar populations with simulation-derived star formation histories, ages, metallicities, spatial distributions, and kinematics [2412.08702]. Its central scientific target is the planetary nebula luminosity function (PNLF), especially the long-standing problem of the apparently universal bright-end cutoff in \(\mathrm{[O\,III]}\,\lambda 5007\), which had previously been modeled mainly with artificial stellar populations and often fixed metallicity assumptions [2412.08702]. Within this framework, PN populations occur naturally in galaxies with diverse evolutionary pathways, and the resulting models are used both to interpret extragalactic distance measurements and to connect PN observables to galaxy formation and stellar population physics [2501.17926].

## 1. Definition, scope, and nomenclature

PICS was introduced to move PN population synthesis from idealized stellar populations to the self-consistent stellar content of cosmological simulations and observations [2501.17926]. In this usage, each stellar particle is treated as a single stellar population (SSP), and PN properties are derived from the particle’s age, mass, metallicity, and initial mass function, then placed back into the galaxy with inherited positional and kinematic information [2412.08702]. This makes PICS simultaneously a population-synthesis framework and a bridge between post-AGB stellar evolution, nebular emission, and galaxy assembly histories.

The framework is explicitly motivated by two deficiencies in earlier PNLF models: the limited realism of artificially constructed or fixed-metallicity parent populations, and the lack of a solid theoretical explanation for the bright-end universality of the PNLF across galaxy types [2501.17926]. PICS is therefore designed to be modular, allowing substitution and testing of different lifetime prescriptions, initial-to-final mass relations (IFMRs), post-AGB tracks, nebular models, and extinction treatments [2501.17926].

The acronym is potentially ambiguous. An unrelated astrophysical pipeline, also called PICS, denotes the **Pipeline for Images of Cosmological Strong lensing**, which simulates strong gravitational lensing in galaxy clusters rather than planetary nebulae [1511.03673]. In PN studies, however, PICS refers specifically to **PNe In Cosmological Simulations** [2412.08702].

## 2. Modeling architecture and population-synthesis workflow

The core PICS workflow is a staged mapping from an SSP to a synthetic PN population. The sequence is stated as

\[
\text{SSP Properties} \rightarrow \text{Lifetime Function} \rightarrow \text{IFMR} \rightarrow \text{Post-AGB Track} \rightarrow \text{Nebular Model}.
\]

This chain converts a parent stellar population into PN brightnesses and, by aggregation over stellar particles, into a galaxy-wide PNLF [2412.08702].

For each stellar particle, PICS uses age, total mass, metallicity, and IMF as the SSP input. The stellar lifetime function maps the age of the population to the masses of stars reaching the post-AGB phase; the IFMR then determines the central-star mass from the initial mass and metallicity; post-AGB tracks provide the luminosity and effective temperature evolution of the central star; and the nebular model converts these stellar properties into observable \(\mathrm{[O\,III]}\,\lambda 5007\) emission [2412.08702]. In the fiducial implementations described in the project papers, the framework adopts metallicity-sensitive lifetimes, IFMRs, and post-AGB tracks from Miller Bertolami (2016), and an empirical PN emission model from Valenzuela et al. (2019) that includes both optically thick and optically thin nebulae [2412.08702].

The principal observable is the absolute \(\lambda 5007\) magnitude,
\[
M(5007) = -2.5 \log I(5007) - 13.74,
\]
where \(I(5007)\) is the intensity at 10 pc distance [2501.17926]. The PNLF is then constructed as the distribution
\[
\frac{dN}{dM(5007)},
\]
with \(N\) the number of PNe [2412.08702].

A compact summary of the main components is as follows.

| Component | Source or prescription | Role |
|---|---|---|
| SSP input | Age, mass, metallicity, IMF | Parent population |
| Lifetime function | Miller Bertolami (2016) fits/interpolation | Time to PN phase |
| IFMR | Miller Bertolami (2016) and alternatives | Central-star mass |
| Post-AGB tracks | Miller Bertolami (2016) | \(L\) and \(T_\mathrm{eff}\) evolution |
| PN model | Valenzuela et al. (2019) | \(\mathrm{[O\,III]}\,\lambda5007\) output |
| Extinction treatment | Later project extensions | Emergent flux modification |

Because PNe are assigned to stellar particles, they inherit the stellar particle’s location and velocity, so the synthetic populations are not only photometric but also spatially and kinematically embedded in the host galaxy [2412.08702].

## 3. Metallicity, helium abundance, and the theoretical origin of the bright end

The central physical result of PICS is that realistic stellar populations, especially their metallicity distributions and metallicity-dependent lifetimes, are required to reproduce the bright end of the PNLF across galaxy types [2412.08702]. The framework finds that the bright-end cutoff is reproduced only when the dependence of stellar lifetimes on metallicity is included; in old, metal-rich systems, neglecting metallicity yields too few high-mass central stars and hence too few very bright PNe [2412.08702]. The stated theoretical implication is that the universality of the PNLF bright end is **not fundamental**, but depends on whether the relevant stellar physics is modeled within each galaxy’s true evolutionary context [2412.08702].

The project’s first installment generalizes this result by isolating the roles of metallicity, helium abundance, and the IFMR [2501.17926]. Metallicity affects the time to reach the PN phase, the mapping from initial to final mass, and the post-AGB luminosity and temperature evolution. A key outcome is that **old metal-rich populations can harbor much brighter PNe than old metal-poor ones**, and that the observed bright end can be reached even for stellar populations of 10 Gyr at high metallicity [2501.17926]. This directly addresses the tension between the old stellar populations of massive early-type galaxies and their apparently invariant \(M^*\) values.

Helium abundance enters at high metallicity. PICS explores a linear increase of \(Y\) with \(Z\) and a “cutoff He” case in which helium saturates above solar metallicity. The study identifies helium abundance as a vital ingredient at high metallicities and examines how helium saturation changes lifetime behavior and thus the availability of bright PN progenitors [2501.17926]. This suggests that the bright end cannot be understood from metallicity alone when the stellar population reaches super-solar regimes.

The IFMR is likewise a controlling parameter. The project finds that the PNLFs of old stellar populations are sensitive to the IFMR, allowing for the production of bright PNe, and that the bright end for old populations is especially sensitive to the low-initial-mass portion of the relation [2501.17926]. In practical terms, small shifts in the IFMR can materially change whether a nominally old population still produces central stars massive enough to populate \(M^*\).

## 4. Comparison with observations in the Milky Way and early-type galaxies

PICS was directly compared with the local Milky Way PNLF by simulating a region of a Milky Way-like disk galaxy and contrasting it with the empirical PNLF for PNe within 2 kpc of the Sun, described as statistically complete down to six orders of magnitude below the bright end [2412.08702]. The modeled PNLF closely matches the observed shape, including normalization over six magnitudes, with only minor differences of approximately \(0.2\) dex at the dim end and slight discrepancies at the bright end that are attributed to expected galaxy-to-galaxy variation and observational uncertainties such as circumstellar dust not modeled in the Milky Way data comparison [2412.08702]. This establishes that the same framework can reproduce both the bright-end behavior and the overall PN population in a realistic galactic environment.

The framework was also tested against early-type galaxies (ETGs) observed with VLT-MUSE, in a study that compared \(M^*\) and the luminosity-specific PN number at the top 0.5 mag, \(\alpha_{0.5}\), to stellar population parameters and to matched Magneticum Pathfinder analogue galaxies processed with PICS [2509.10175]. The observed ETGs have average mass-weighted ages and metallicities that are typically old,
\[
9 < \mathrm{Age} < 13.5~\mathrm{Gyr},
\]
and rather metal-rich,
\[
-0.4 < \mathrm{[M/H]} < +0.2,
\]
and within these intervals the study finds \(M^*\) to be independent of age and metallicity [2509.10175].

The same study reports a positive correlation between \(\alpha_{0.5}\) and the mass fraction of stellar populations with ages of \(2\)–\(10\) Gyr, and a negative correlation with the fraction older than \(10\) Gyr [2509.10175]. Similar trends are found in the PICS analogue galaxies, and the authors state that the presence of at least \(\sim 2\%\) of stellar mass younger than 10 Gyr is, in principle, sufficient to form the \(M^*\) PNe in ETGs [2509.10175]. They further argue that sufficiently massive PN central stars can be produced without invoking exotic binary evolution channels, using the updated post-AGB models incorporated into PICS [2509.10175]. A common misconception in the ETG literature is therefore that bright PNe in old, metal-rich systems necessarily require non-standard progenitors; the PICS-based result is that a small intermediate-age component can be sufficient.

For observational fitting, the ETG study uses the analytic PNLF form
\[
\phi_1(M) \propto e^{0.307 M} \left\{ 1 - e^{3 (M^* - M)} \right\},
\]
and defines
\[
\alpha_{0.5} = \text{Number of PNe in top 0.5 mag} \,/\, L_{\rm bol}
\]
[2509.10175]. Within the PICS interpretation, \(\alpha_{0.5}\) acts as a stellar population tracer, while \(M^*\) remains relatively invariant over the old, metal-rich parameter space probed in those data [2509.10175].

## 5. Circumnebular extinction and the bright-end cutoff

A subsequent extension of the framework addresses circumnebular extinction as a mechanism shaping the observed PNLF bright end [2509.25306]. The theoretical problem is that young stellar populations, which host more massive PN central stars, are expected to produce intrinsically brighter \(\mathrm{[O\,III]}\) emission than is observed, whereas the empirical bright-end cutoff remains near a nearly constant \(M^*(5007)\) [2509.25306]. The proposed solution is self-extinction by dust produced in AGB mass loss, especially for higher-mass central stars, which reduces the emergent \(\mathrm{[O\,III]}\) flux and compresses intrinsic differences toward the observed cutoff [2509.25306].

In this implementation, PICS adds a circumnebular extinction module to the existing birth distribution, IFMR, post-AGB, and nebular-emission machinery. The treatment starts from empirical fits for \(c(\mathrm{H}\beta)\) as a function of final mass, expressed as
\[
c(\mathrm{H}\beta) = S \times M_{\rm final,bright} + I,
\]
and then converts \(c(\mathrm{H}\beta)\) to \(A_{5007}\) [2509.25306]. Six extinction recipes were tested, differing in sample choice and fitting method, including Orthogonal Regression, Ordinary Least Squares, and an optimized fit for PNLF flatness [2509.25306].

The principal result is differential sensitivity with stellar age. Steeper extinction-mass slopes lead to large changes in the bright-end cutoff for **young** stellar populations, while older PNe are much less affected because the extinction recipes converge to little or no extinction at long post-AGB ages [2509.25306]. Thus, circumnebular extinction is described as crucial primarily for the youngest PN populations [2509.25306]. At the same time, direct comparison with observed PNe in nearby galaxies shows that the modeled extinctions do not fully reproduce the large observed spread or the highest extinction values, even after adding scatter and applying a global shift of \(\Delta c(\mathrm{H}\beta)=+0.08\) to some relations [2509.25306].

This leaves an unresolved tension. The project states that improvements on both the modeling and observational sides are needed, and lists possible causes including unmodeled PN formation channels, binary evolution, observational bias, or incomplete physics regarding dust formation and ejection [2509.25306]. Relative to the earlier phase of the project, where internal dust extinction had not yet been explicitly incorporated [2501.17926], the extinction study represents a move from omission of the effect to a parameterized exploration of how strongly it sculpts the PNLF.

## 6. Dynamical applications, simulation context, and open limitations

Because PICS assigns PNe to star particles with inherited positions and velocities, the framework can be used not only for luminosity functions but also for dynamical studies in low-surface-brightness environments [2604.14272]. In an application to the Magneticum Box4 (uhr) simulation, the method was used to model PNe in a massive elliptical galaxy of stellar mass \(4\times10^{11}\,M_\odot\) hosting a prominent tidal stream from a progenitor of \(2.4\times10^{10}\,M_\odot\) [2604.14272]. The study finds that PNe in the brightest 1.5 mag of the luminosity function are sufficient to recover the velocity gradient and dispersion of the massive stream, whereas restricting the tracer set to only the brightest 1.0 mag yields sparse sampling and stronger contamination by halo PNe [2604.14272]. In that context, PNe are presented as an attractive alternative to expensive deep IFU observations for constraining the gravitational potential of galaxy outskirts [2604.14272].

This dynamical use follows directly from the cosmological embedding of the model. The same properties that make PICS useful for explaining the PNLF bright end—self-consistent star formation histories, chemical enrichment histories, and realistic mixtures of SSPs—also make it suitable for mock tracer catalogues in streams, halos, and other faint structures [2412.08702]. A plausible implication is that PNLF modeling and PN kinematics are no longer separable topics once the tracer population is generated within a resolved cosmological galaxy.

Several limitations remain explicit in the project papers. For very high metallicities, \(Z>0.04\), some stellar models are extrapolated rather than directly simulated [2412.08702]. In the original Milky Way comparison, circumstellar extinction was not yet included and was identified as a source of discrepancy at the dim end [2412.08702]. The simulated galaxies used for local comparisons are similar to, but not exact analogues of, the Milky Way, with differences in bulge and disk populations [2412.08702]. In the dynamical application, stochastic sampling, contamination by halo PNe in projection, and the need to verify assumptions in unusual stellar populations are all noted as practical constraints [2604.14272].

Taken together, the PICS program establishes a cosmologically grounded PN population-synthesis framework in which the PNLF is interpreted as an emergent consequence of realistic stellar ages, metallicities, helium abundances, IFMR choices, nebular physics, and extinction, rather than as a phenomenon to be reproduced with artificial populations alone [2412.08702]. Its later extensions show that this same framework can be used for distance-scale work, stellar-population inference, and the kinematics of faint galactic substructure [2509.10175].

Source: https://www.emergentmind.com/topics/pne-in-cosmological-simulations-pics