---
title: 'A-SLOTH: Tracing Early Stars & Local Observables'
url: https://www.emergentmind.com/topics/a-sloth
type: topic
---

# A-SLOTH: Tracing Early Stars & Local Observables

A-SLOTH, short for **Ancient Stars and Local Observables by Tracing Halos**, is a semi-analytical model of galaxy formation focused on the early Universe and on linking the first stars and galaxies to present-day and high-redshift observables. It is built on dark-matter merger trees generated either from Extended Press–Schechter theory or imported from numerical simulations, and it applies analytical baryonic prescriptions to model gas cooling, star formation, feedback, chemical enrichment, and stellar populations, including both metal-free Population III and metal-enriched Population II stars. Public releases describe it as the first public code that connects the formation of the first stars and galaxies to observables, while emphasizing moderate computational requirements, modularity, and applications ranging from Milky Way stellar archaeology to reionization, transients, compact remnants, and high-redshift galaxy populations [2206.00223][2209.07339].

## 1. Origins, scope, and scientific niche

A-SLOTH was developed to address the dynamic-range problem inherent in first-star studies: the relevant physics spans minihaloes of \(10^5\)–\(10^6\,M_\odot\) at \(z\gtrsim 20\), atomic-cooling haloes, Milky Way progenitors, and present-day dwarf satellites, across \(\sim 13.8\) Gyr of cosmic evolution. In this regime, semi-analytic models are computationally advantageous because they can sit on top of dark-matter-only merger trees and replace unresolved baryonic processes with parametrized prescriptions, making systematic parameter studies feasible where full hydrodynamical calculations remain prohibitive [2209.07339].

The code’s stated purpose is not generic galaxy formation alone, but the reconstruction of “ancient stars and local observables” from hierarchical halo assembly. That emphasis gives A-SLOTH a distinctive position among semi-analytic models: it is optimized for Population III star formation, the Pop III-to-Pop II transition, the build-up of extremely metal-poor stellar populations, and the use of the Milky Way and its satellites as a fossil record of early-Universe physics. The software description in the Journal of Open Source Software notes that A-SLOTH can “build up a Milky-Way-like galaxy star by star in only a few minutes,” and that its modular design is intended to let users add new physics without constructing a new framework from scratch [2209.07339].

## 2. Computational framework and merger-tree architecture

At its core, A-SLOTH evolves baryons along halo merger trees. The original implementations used Extended Press–Schechter trees generated with the Parkinson et al. algorithm, while later applications also ingested trees from numerical simulations such as the Caterpillar Milky Way zooms and a high-resolution \(8\,h^{-1}\,\mathrm{cMpc}\) cosmological box. Across these modes, the workflow is consistent: read or generate a merger tree, step through haloes in time order, update gas reservoirs, form stars, apply feedback, track metals, and assemble observables from the resulting stellar and halo histories [2209.07339][2508.09331][2506.06139].

Later baryon-cycle implementations make the state variables explicit by evolving cold gas, hot gas, cumulative stellar mass, and unbound outflows. A representative star-formation law used in current A-SLOTH applications is  
$$
\dot{M}_\star \simeq \eta\,\frac{M_{\rm cold}}{t_{\rm cold,ff}},
$$
where \(\eta\) is a star-formation efficiency per mean free-fall time, \(M_{\rm cold}\) is the cold gas mass, and \(t_{\rm cold,ff}\) is the cold-gas free-fall time. In this sense, the model is semi-analytic in the strict SAM sense: halo growth is inherited from the merger tree, while baryonic evolution is encoded in coupled, timestep-based reservoir equations [2506.06139].

A-SLOTH also samples stars individually in several of its major implementations. That feature is central to its treatment of stochastic feedback, because ionizing luminosities, stellar lifetimes, and supernova events can be attached to discrete stellar masses rather than only to IMF-averaged populations. In Milky Way satellite applications, this is extended to individually sampled Pop II stars; in first-star applications, the code likewise samples individual Pop III stars from an explicit IMF [2206.00223][2202.01220].

## 3. Physical prescriptions: cooling, stellar populations, and feedback

A-SLOTH distinguishes between molecular-cooling minihaloes and atomic-cooling haloes, and it explicitly separates metal-free Population III star formation from metal-enriched Population II star formation. In multiple applications, Pop III formation occurs in gas below a critical metallicity threshold of order \(Z \lesssim 10^{-5}\,Z_\odot\), while Pop II star formation takes over once enrichment crosses that threshold. The Pop III IMF is generally represented as a power law,
$$
\frac{dN}{dM_\star} \propto M_\star^{-\alpha_{\rm III}},
$$
with \(M_{\rm min}\), \(M_{\rm max}\), and \(\alpha_{\rm III}\) treated as tunable parameters in the calibrated framework [2208.01673][2410.05393].

The feedback channels are radiative, chemical, and mechanical. Radiative feedback includes ionizing photons and Lyman–Werner suppression of \(\mathrm{H}_2\) cooling; chemical feedback follows stellar yields and the resulting Pop III-to-Pop II transition; mechanical feedback is dominated by supernova-driven gas removal and halo outflows. In the high-redshift galaxy implementation, the outflow prescription is explicitly energy-driven, with mass loading proportional to supernova energy injection divided by halo binding energy. That scaling captures the expected strong suppression of star formation in shallow potential wells and is one of the key control knobs in reproducing both dwarf-galaxy properties and high-redshift metallicities [2206.00223][2506.06139].

Chemical evolution is a persistent theme across A-SLOTH papers. The code tracks metals through halo growth, accretion, enrichment, and outflows, and later versions incorporate a stochastic description of inhomogeneous mixing by assigning star-forming gas a metallicity offset relative to the halo-average metallicity. This allows the framework to model not only global metallicity growth but also the low-metallicity tails relevant for extremely metal-poor stars, second-generation enrichment signatures, and Population III relic diagnostics [2209.07339][2506.06139].

## 4. Calibration strategy and inferred parameter space

The public-release description states that A-SLOTH was originally calibrated to six observables, including the optical depth to Thomson scattering, the stellar mass of the Milky Way and its satellite galaxies, the number of extremely-metal poor stars, and the cosmic star-formation rate density at high redshift. Those observables were chosen to tie together local stellar archaeology and cosmological reionization-era constraints within a single forward model [2206.00223].

A later full statistical calibration expanded this to a likelihood based on **nine independent observables** and produced posterior constraints for **11 parameters** describing early star formation and feedback. The nine observables span Milky Way-specific and cosmologically representative quantities, and the inferred best-fit model gives a Pop III IMF slope of \(dN/dM \propto M^{-1.77}\), with \(M_{\rm min}=13.6\,M_\odot\) and \(M_{\rm max}=197\,M_\odot\). The same calibration also reports posterior medians \(f_{\rm esc,II}=0.175\), \(f_{\rm esc,III}=0.525\), \(v_{\rm sv}/\sigma_{\rm sv}=1.75\), and \(c_{\rm ZIGM}=3.32\), while emphasizing that the IMF-generating parameters remain poorly constrained and that the IMF slope could range from log-flat to Salpeter within the allowed posterior volume [2410.05393].

| Parameter | Posterior median | Source |
|---|---:|---|
| \(\alpha_{\rm III}\) | \(1.77\) | [2410.05393] |
| \(M_{\rm min}\) | \(13.6\,M_\odot\) | [2410.05393] |
| \(M_{\rm max}\) | \(197\,M_\odot\) | [2410.05393] |
| \(f_{\rm esc,II}\) | \(0.175\) | [2410.05393] |
| \(f_{\rm esc,III}\) | \(0.525\) | [2410.05393] |
| \(v_{\rm sv}/\sigma_{\rm sv}\) | \(1.75\) | [2410.05393] |
| \(c_{\rm ZIGM}\) | \(3.32\) | [2410.05393] |

This calibration architecture is one of the reasons A-SLOTH is used as a forward model rather than only as an exploratory toy model. A plausible implication is that the framework’s value lies as much in its parameter posteriors and degeneracy structure as in any single best-fit realization, because multiple later application papers reuse these calibrated parameter ranges to test new observables without abandoning consistency with earlier Milky Way and reionization constraints [2410.05393].

## 5. Major application domains

A-SLOTH has been extended into a broad suite of application-specific models, all retaining the merger-tree backbone while modifying the baryonic modules or the observables.

| Application | Representative result | Source |
|---|---|---|
| Milky Way satellites and ultra-faints | The \(z=0\) stellar mass–halo mass relation develops a plateau in the ultra-faint regime; the plateau stellar mass is set by how many stars form before supernovae regulate further star formation | [2202.01220] |
| Pop III contribution to high-\(z\) galaxies | Pop III fractions are highest in \(10^7\)–\(10^9\,M_\odot\) haloes; Pop III-dominated galaxies are predicted to be too faint for direct JWST detection | [2208.01673] |
| Multidimensional comparison to MW satellites | Unsupervised clustering applied to A-SLOTH satellites identified mean stellar metallicity as the dominant mismatch between the fiducial model and observations | [2209.11248] |
| White dwarf luminosity function | A-SLOTH star-formation histories imply that the faint end of the white dwarf luminosity function can carry a Pop III remnant signature under a bottom-heavy Pop III IMF | [2409.00868] |
| High-\(z\) chemical evolution and IMF constraints | Reproducing the observed mass–metallicity–SFR relation at \(z\sim4\)–10 favored a Kroupa-like Pop II IMF with \(m_{\max}\gtrsim 200\,M_\odot\) under the adopted yield and outflow model | [2506.06139] |
| Dust in early galaxies | \(z\sim7\) ALMA galaxies are consistent with normal dust growth, whereas dust-poor JWST galaxies at \(z\gtrsim10\) require radiation-pressure-driven dust clearing during recent highly efficient bursts | [2309.02415] |
| Balmer-series observability of Pop III | Default Pop III \(H\alpha\) fluxes peak near \(\sim10^{-21}\,\mathrm{erg\,s^{-1}\,cm^{-2}}\), remaining 2–3 orders of magnitude below the quoted JWST/NIRSpec threshold without strong lensing | [2508.09331] |
| Little Red Dots and direct-collapse black holes | DCBH-seeded SMBH populations agree better than stellar-remnant seeds with observed LRD demographics and with the spectral properties of the most extreme LRD candidates | [2508.14155] |

Taken together, these studies show that A-SLOTH has evolved from a Milky Way archaeology model into a general early-Universe inference platform. Its observables now include halo stars, satellite mass functions, white dwarfs, cosmic SFRD, nebular and dust emission, black-hole populations, and prospective gravitational-wave and transient signatures [2206.00223][2409.00868].

## 6. Limitations, uncertainties, and methodological status

A-SLOTH remains a semi-analytic model, and its public software description explicitly notes that it uses simplified, parametrized baryonic physics rather than resolving small-scale hydrodynamics. Predictions at very low halo masses or extreme metallicities therefore rely on approximate prescriptions and extrapolations, and inference about the first stars remains model-dependent. This limitation is structural, not incidental: the framework gains speed and parameter-space coverage by replacing direct fluid evolution with analytic or statistical subgrid rules [2209.07339].

Later application papers make the main uncertainties explicit. In the white-dwarf study, the faint-end Pop III signature is sensitive to the assumed Pop III initial–final mass relation, core composition, and the neglect of binary evolution; the paper states directly that adopting a different Pop III IFMR can erase the faint-end hump [2409.00868]. In the high-redshift chemical-evolution study, the preferred \(m_{\max}\gtrsim 200\,M_\odot\) result depends on stellar-evolution tracks up to \(600\,M_\odot\), metallicity calibrations, a small simulation volume, and phenomenological outflow parameters [2506.06139]. In the dust extension, the conclusions depend on uncertain supernova dust yields, dust geometry, grain properties, and the assumption of tight gas–dust coupling [2309.02415]. In the LRD study, DCBH seeding criteria, Bondi-like gas structure, instantaneous BH mergers, and non-self-consistent dust treatment remain important caveats [2508.14155].

Even with those limitations, A-SLOTH’s significance lies in the combination of public availability, modularity, and cross-domain calibration. It can function as a forward model for Bayesian inference, a generator of synthetic galaxy populations, a subgrid layer for larger cosmological analyses, and a bridge between first-star physics and local fossil records. In that sense, A-SLOTH is best understood not as a single static model, but as an extensible semi-analytic research program for tracing halo assembly into stellar, chemical, radiative, and compact-remnant observables across cosmic time [2206.00223][2209.07339][2410.05393].

Source: https://www.emergentmind.com/topics/a-sloth