---
title: 'MESA-QUEST: Quasi-star Simulation Toolkit'
url: https://www.emergentmind.com/topics/mesa-quasi-star-evolutionary-simulation-toolkit-mesa-quest
type: topic
---

# MESA-QUEST: Quasi-star Simulation Toolkit

Searching arXiv for the MESA-QUEST papers and related MESA context.
MESA Quasi-star Evolutionary Simulation Toolkit (MESA-QUEST) is an open, MESA-based implementation of quasi-star models designed to study the formation and growth of heavy black hole seeds via direct collapse. It extends the 1D stellar evolution “star” module of the Modules for Experiments in Stellar Astrophysics (MESA) to represent a quasi-star as a hydrostatic, radiation-supported envelope hosting a central, accreting black hole, with the core goal of reproducing the structural and evolutionary behavior of quasi-stars found with the Cambridge STARS code and of providing a flexible platform for exploring additional physics and parameter sensitivity relevant to the viability of quasi-stars as heavy seeds of supermassive black holes [2510.11772; 2507.12651]. The framework was introduced to address the rapid appearance of supermassive black holes at redshifts \(z \geq 9\), for which direct-collapse black hole (DCBH) scenarios and quasi-star intermediaries remain a leading heavy-seed pathway [2510.11772].

## 1. Astrophysical context and scientific motivation

The origin of the first supermassive black holes (SMBHs) observed at \(z \geq 9\) remains one of the most challenging open questions in astrophysics. Their rapid emergence suggests that massive “heavy seeds” must have formed early, possibly through the direct collapse of pristine gas clouds in the first galaxies [2510.11772]. In the DCBH scenario, pristine gas in early protogalaxies collapses monolithically without fragmenting, aided by strong Lyman–Werner radiation fields that suppress \( \mathrm{H}_2 \) cooling, low metallicity, and angular-momentum transport by non-axisymmetric instabilities [2510.11772].

Within this framework, the end-state of such a collapse can be a quasi-star: a massive, radiation-supported envelope that harbors a central, accreting black hole. Quasi-stars are compelling because the embedded black hole can accrete at a rate set by the envelope’s Eddington limit rather than the black hole’s own Eddington limit, enabling rapid growth to intermediate-mass black holes of \(\gtrsim 10^3\!-\!10^5\,M_\odot\) within the tight cosmic time budget implied by early SMBH detections [2510.11772]. The earlier MESA-QUEST implementation described the same scientific objective in terms of porting and extending the methodology previously realized in the Cambridge STARS code into MESA’s modular framework, thereby enabling tighter numerical control, general-relativistic stability corrections, and straightforward inclusion of additional physics [2507.12651].

The toolkit is therefore positioned at the intersection of DCBH formation theory, quasi-star structure theory, and numerical stellar evolution. Its scientific role is not merely to reproduce canonical quasi-star tracks, but to delineate the physical limits under which quasi-stars can remain stable, survive mass loss, and produce heavy seeds capable of evolving into the earliest SMBHs detected by JWST and Chandra [2510.11772].

## 2. Software architecture and integration with MESA

MESA-QUEST is built upon the Modules for Experiments in Stellar Astrophysics, specifically the 1D stellar evolution “star” module [2510.11772]. In the initial implementation, all new quasi-star physics—inner boundary treatment, accretion source term, Tolman–Oppenheimer–Volkoff (TOV) correction, sound-speed regularization, and Bondi-radius limiter—are implemented in `run_star_extras` within the MESA-QUEST repository [2507.12651]. The broader MESA substrate contributes the EOS, opacity, convection, atmosphere, timestep, mesh-refinement, and solver infrastructure on which quasi-star extensions are layered [2208.03651; 1009.1622].

The defining architectural modification is the replacement of the standard stellar center boundary conditions. Instead of the usual \(m(r=0)=0\) and \(l(r=0)=0\), MESA-QUEST imposes
\[
m(r=R_0)=M_{\rm BH}+M_{\rm cav}, \qquad l(r=R_0)=L_{\rm BH},
\]
where \(R_0\) is the inner boundary radius, \(M_{\rm BH}\) is the black hole mass, \(M_{\rm cav}\) is the cavity mass interior to \(R_0\), and \(L_{\rm BH}\) is the black-hole accretion luminosity [2510.11772]. This converts a stellar-center problem into an envelope-on-core problem in which the unresolved central region is represented by a point mass plus a prescribed luminosity source.

Two inner boundary definitions are available. One is the Bondi radius \(R_B\), which matches the quasi-hydrostatic envelope to a Bondi-like inflow. The other is the saturated-convection radius \(R_i\), following Coughlin (2024), which places the hydrostatic boundary deeper, where outward convective energy transport saturates [2510.11772]. The existence of these two boundary options is central to the framework because the boundary placement strongly affects the accessible black-hole growth regime.

The repository is public at `http://www.github.com/andysantarelli/MESA-QUEST`, and includes source modifications and example inlists illustrating how to enable quasi-star physics, choose inner boundary types, and select wind schemes [2510.11772]. Standard MESA build procedures are followed; the customized `star` module is compiled and executed via the included inlists. Output includes black-hole mass, envelope mass, luminosities, radii, inner-boundary radius, and mass-loss rates [2510.11772]. A typical 1D quasi-star track to \(\sim 10^3\,M_\odot\) black holes completes in hours on a modern workstation, with runtime depending on wind strength and whether the inner boundary approaches the photosphere [2510.11772].

## 3. Physical model, governing equations, and inner-boundary formalism

For \(r \geq R_0\), MESA-QUEST solves the standard stellar structure equations with quasi-star-specific modifications [2510.11772]. Hydrostatic equilibrium is written in MESA form with a GR correction,
\[
\frac{dP}{dr} = - \frac{G_{\rm eff} M(r)\rho}{r^2},
\]
where \(G_{\rm eff}\) includes a Tolman–Oppenheimer–Volkov correction appropriate to the deep potential well and is recalculated each timestep [2510.11772]. Mass conservation remains
\[
\frac{dM}{dr}=4\pi r^2 \rho,
\]
and radiative transport is
\[
\frac{dT}{dr} = - \frac{3\kappa \rho L}{16\pi a c T^3 r^2}.
\]
Convection is treated with the Ledoux criterion and mixing-length theory, with \(\alpha_{\rm MLT}=2\) in the fiducial runs [2510.11772].

The luminosity equation is modified so that the dominant energy source is black-hole accretion rather than nuclear burning:
\[
\frac{dL}{dr} = 4\pi r^2 \rho \left(\epsilon_{\rm nuc} - \epsilon_\nu - T\frac{ds}{dt}\right) + (\text{radiative + convective transport}),
\]
with \(\epsilon_{\rm nuc} \approx 0\) in fiducial quasi-star models and the primary energy source supplied through the inner boundary condition \(l(R_0)=L_{\rm BH}\) [2510.11772]. This formulation encodes the defining quasi-star assumption that photon transport and convection carry accretion power outward through the envelope.

The cavity mass \(M_{\rm cav}\) accounts for gas within the inner boundary. For Bondi-type models, the unresolved density field is taken as
\[
\rho(r)=\rho_0\left(\frac{r}{R_0}\right)^{-n/2},
\]
with \(n=1\) or \(3\) depending on whether angular momentum is transported inward or outward, and
\[
M_{\rm cav} = \int_{R_S}^{R_0} 4\pi r^2 \rho(r)\,dr,
\]
where the Schwarzschild radius \(R_S \ll R_0\) is typically negligible [2510.11772]. In the earlier Bondi-only implementation, \(M_{\rm cav}\) was estimated by extending an inner \(\rho \propto r^{-3/2}\) profile interior to the Bondi radius [2507.12651].

The Bondi inner boundary is defined by
\[
\frac{1}{2} m c_s^2 = \frac{G m M_{\rm BH}}{R_B}, \qquad
R_B = \frac{2 G M_{\rm BH}}{c_s^2},
\]
with \(c_s\) evaluated at the innermost resolved cell [2510.11772]. In the 2025 implementation, \(c_s\) is averaged over the innermost \(\lesssim 50\) zones and a safety limiter enforces \(|\Delta r_B|/r_B \leq 0.05\) per timestep to stabilize the boundary [2507.12651].

The saturated-convection boundary follows a deeper hydrostatic truncation based on dimensionless variables \(\xi=r/r_i\) and \(m_i = m(r)/M_i\), where \(M_i \equiv m(r_i)\), with governing ODE
\[
K_i \frac{d}{d\xi}\left[\xi^{-2}\left(\frac{dm_i}{d\xi}\right)^{1/3}\right]
= -\left(\frac{m_i}{\xi^4}\right)\left(\frac{dm_i}{d\xi}\right),
\]
and
\[
K_i = \left[\frac{L\sqrt{3}}{2\eta}\right]^{2/3}\left(\frac{r_i}{M_i}\right)^{5/3}\left(\frac{1}{G}\right).
\]
The associated convection-dominated accretion flow region satisfies \(\rho \propto r^{-1/2}\), and \(M_{\rm cav}\) is adjusted accordingly [2510.11772].

## 4. Accretion prescriptions, luminosity coupling, and microphysics

MESA-QUEST implements two accretion prescriptions, each coupled directly to the inner luminosity boundary [2510.11772]. The first is convection-limited Bondi accretion, following Ball et al. (2012), in which the black-hole growth rate is
\[
\dot{M}_{\rm BH} =
16\pi \eta \frac{(1-\epsilon)}{\epsilon}
\left[\frac{\rho}{c_s c^2 \Gamma}\right](G M_{\rm BH})^2,
\]
where \(\eta\) is the convective efficiency, \(\epsilon\) the radiative efficiency, \(\Gamma\) the adiabatic index, and the local \(\rho\) and \(c_s\) are evaluated at the inner boundary [2510.11772]. The associated luminosity is
\[
L = \frac{\epsilon}{1-\epsilon}\dot{M}_{\rm BH} c^2.
\]

The second is envelope-Eddington-scaled accretion, motivated by quasi-star theory. In this formulation,
\[
L_0 = \alpha L_{\rm E} = \alpha\,4\pi \left(\frac{c}{\kappa}\right) G M_*,
\]
where \(\alpha\) is a control parameter, \(\kappa\) is the opacity evaluated at the deep core, and \(M_*\) is the total quasi-star mass [2510.11772]. The implied accretion rate is
\[
\dot{M}_{\rm BH} =
\frac{1-\epsilon}{\epsilon}
\left(\frac{4\pi}{\kappa c}\right)\alpha G M_*.
\]
The fiducial Eddington-scaled calculations span \(\alpha \in \{0.5, 0.8, 1.0, 1.4\}\) with \(\epsilon = 0.1\) [2510.11772].

The relevant radiative limits are the envelope Eddington luminosity,
\[
L_{\rm E} = \frac{4\pi G M_* c}{\kappa},
\]
and the black-hole Eddington limit,
\[
L_{\rm Edd,BH} = \frac{4\pi G M_{\rm BH} c}{\kappa}.
\]
Core opacities show that \(\alpha \approx 1\) is a critical boundary between compact, high-opacity envelopes and expanded, lower-opacity envelopes [2510.11772]. This opacity sensitivity is significant because it controls whether the envelope enters runaway expansion or remains compact and stable.

The microphysics assumes primordial gas with fiducial composition \(X=0.75\) and \(Y=0.25\) [2510.11772]. Standard MESA EOS and opacity modules are used, with convection treated via the Ledoux criterion and mixing-length theory at \(\alpha_{\rm MLT}=2\) [2510.11772]. The earlier implementation specifies MESA’s “ML1” option with \(\alpha_{\rm MLT}=2\) to match Ball et al. (2011), and emphasizes that radiation pressure is substantial to dominant in the interior [2507.12651].

MESA itself contributes a more general infrastructure for time-dependent convection, improved energy conservation, automatic differentiation, and updated opacity and EOS handling [2208.03651]. The MESA VI description notes that MESA does not include GR radiation hydrodynamics, multi-D radiation transport, or a native black-hole accretion flow model; in quasi-star applications these are represented as 1D source terms injected near the inner boundary [2208.03651]. This contextualizes MESA-QUEST as a hydrostatic 1D quasi-star framework rather than a full GRRMHD treatment.

## 5. Winds, mass loss, parameter space, and practical workflow

A major extension of MESA-QUEST is the incorporation of multiple wind and mass-loss prescriptions [2510.11772]. Three regimes are explored.

| Regime | Prescription | Role in MESA-QUEST |
|---|---|---|
| Reimers | \(\dot{M}_{\rm R} = 4 \times 10^{-13}\eta (L/L_\odot)(R/R_\odot)(M_\odot/M)\,M_\odot\,{\rm yr}^{-1}\) with \(\eta=0.5\) | Lower bound |
| Dutch | van Loon et al. (2005) for \(T_{\rm eff} > 10^4\,{\rm K}\); de Jager et al. (1988) for \(T_{\rm eff} \le 10^4\,{\rm K}\) | Intermediate strength |
| Super-Eddington radiation-driven | Dotan & Shaviv (2011) form; Fiacconi et al. (2015) fit \(\dot{M}_{\rm wind}=1.4\times10^{-4} M_*^{0.96} M_{\rm BH}^{0.17}\,M_\odot\,{\rm yr}^{-1}\) | Upper bound |

The Reimers wind has minimal impact on quasi-star lifetimes in the reported models [2510.11772]. The Dutch scheme is stronger than Reimers and curtails growth when envelopes expand and cool; the Nugis & Lamers WR component is not applied because it requires \(X_H < 0.4\) [2510.11772]. The quasi-star-specific radiation-driven winds provide the strongest mass loss and serve as an upper bound on envelope erosion [2510.11772].

The fiducial initialization procedure is explicit. An envelope is built by accreting at \(0.1\,M_\odot\,{\rm yr}^{-1}\) onto a \(20\,M_\odot\) protostar until \(M_* = 10^4\,M_\odot\), following supermassive-star construction approaches [2510.11772]. A seed black hole of \(M_{{\rm BH},0}=10\,M_\odot\) is then embedded through the modified inner boundary conditions, using either a Bondi or saturated-convection inner boundary [2510.11772]. The practical, reproducible setup specified for the repository includes primordial composition, Ledoux convection, \(\alpha_{\rm MLT}=2.0\), TOV correction, a selectable inner boundary type, one of the two accretion schemes, and a chosen wind prescription [2510.11772].

Representative controls listed in the documentation include:
- `initial_mass = 20`
- `accumulate_mass_to = 1d4`
- `accretion_rate = 0.1`
- `use_ledoux_criterion = .true.`
- `mixing_length_alpha = 2.0`
- `use_TOV_correction = .true.`
- `bh_initial_mass = 10`
- `bh_inner_boundary_type = 'bondi' or 'saturated_convection'`
- `bh_accretion_scheme = 'bondi_convective' or 'alpha_edd'`
- `alpha_edd = 1.0`
- `radiative_efficiency = 0.1`
- `convective_efficiency = 0.1`
- `wind_scheme = 'none'/'reimers'/'dutch'/'radiation_driven'`
- `reimers_eta = 0.5` [2510.11772]

The recommended diagnostics are \(M_{\rm BH}(t)\), \(M_*(t)\), \(L(t)\), \(L_{\rm BH}(t)\), \(R_*(t)\), \(R_0(t)\), \(T_{\rm eff}(t)\), \(\kappa_{\rm core}(t)\), and \(\dot{M}_{\rm wind}(t)\) [2510.11772]. Best-practice guidance in the repository notes that \(\alpha \lesssim 1\) is preferred for stability-focused studies, the saturated-convection boundary is preferred for maximizing black-hole growth, and the three wind prescriptions may be used to bracket uncertainties from lower to upper mass-loss limits [2510.11772].

## 6. Evolutionary outcomes, validation, and interpretive issues

The principal result is that quasi-stars can grow central black holes to \(\gtrsim 10^3\,M_\odot\) under favorable conditions [2510.11772]. Without winds, saturated-convection inner boundaries yield substantially larger black-hole growth than Bondi boundaries. In saturated-convection models with \(\alpha \gtrsim 1\), MESA-QUEST reaches \(M_{\rm BH}/M_* \approx 0.55\), approximately five times higher than Bondi-limited cases, which yield \(\approx 0.11\) [2510.11772]. This agrees with the trend found by Coughlin (2024), which suggested ratios up to \(\sim 0.6\), and exceeds the \(\sim 0.167\) ratio reported by Ball et al. (2012) under Bondi-like conditions [2510.11772].

A central interpretive claim of the framework concerns termination. In both boundary schemes, growth halts when the inner boundary radius \(R_0\) approaches the photosphere; the integration fails not because of a physical instability but because the hydrostatic region vanishes [2510.11772]. On that basis, the reported black-hole-to-total-mass “limits” in earlier quasi-star studies are argued to be numerical or integration artifacts rather than physical caps [2510.11772]. This is one of the major controversies addressed by the toolkit, and it is presented as a boundary-placement and stopping-criterion issue rather than as a new physical instability.

Envelope behavior depends strongly on \(\alpha\). For \(\alpha < 1\), envelopes remain compact, hot, and dense; opacities remain high, preventing runaway expansion and maintaining stability over longer times [2510.11772]. For \(\alpha \geq 1\), lower core opacities and higher luminosities drive envelope expansion; in Bondi models, the Bondi radius can show step-like increases when separated convective zones merge and the star becomes fully convective, lowering \(c_s\) and increasing \(R_B\) [2510.11772].

With winds included, the qualitative hierarchy is clear. Reimers winds reduce final \(M_{\rm BH}\) by \(\lesssim 5\!-\!10\%\) for \(\alpha \geq 1\) in Bondi models and less in saturated-convection models [2510.11772]. Dutch winds significantly curtail growth in models with envelope expansion; for \(\alpha=1.4\), the envelope is ejected promptly after expansion disrupts equilibrium, while \(\alpha=1.0\) survives marginally longer but still loses mass rapidly [2510.11772]. Super-Eddington radiation-driven winds are most severe: quasi-stars survive only \(\sim 10^3\) yr, and black holes grow by at most a factor of \(\sim 2\) in mass across \(\alpha\) [2510.11772]. This suggests that, unless envelope accretion replenishes the lost mass at comparable rates, heavy-seed formation can be strongly quenched [2510.11772].

Validation against earlier quasi-star work is explicit. The earlier Bondi-based implementation reproduces the fiducial quasi-star of Ball et al. (2011) to within \(\approx 10\%\) in the final black-hole mass and closely tracks the same structural evolution [2507.12651]. Radial density profiles, inner and outer radii, luminosity evolution, and effective-temperature decline match the STARS results qualitatively, supporting the use of MESA as a validated quasi-star platform [2507.12651]. The later framework extends that validation by showing that Bondi inner-boundary models reproduce the qualitative behavior of Ball et al. (2012), while saturated-convection models recover the larger \(M_{\rm BH}/M_*\) trend associated with deeper hydrostatic boundaries [2510.11772].

## 7. Relation to direct-collapse and supermassive-star studies, limitations, and future directions

MESA-QUEST is closely related to MESA studies of supermassive primordial stars, which provide envelope structures analogous to quasi-star progenitors [2208.00008]. In those models, high accretion rates produce radiation-pressure-dominated, near-Eddington, cool red hypergiant envelopes with \(T_{\rm eff} \approx 10^4\,{\rm K}\), \(L \approx 10^9\!-\!10^{10}\,L_\odot\), and collapse driven by general-relativistic instability during central hydrogen burning [2208.00008]. MESA-QUEST adopts a comparable primordial-composition, Ledoux-convection, \(\alpha_{\rm MLT}=2\), TOV-corrected numerical environment, but replaces the nuclear-burning core with a black-hole inner boundary and accretion-powered luminosity [2510.11772; 2208.00008].

The astrophysical implications are correspondingly specific. Quasi-stars can plausibly produce heavy seeds of \(\gtrsim 10^3\,M_\odot\) under favorable conditions, especially with saturated convection near the black hole and modest winds [2510.11772]. Environmental requirements are consistent with DCBH formation: pristine gas, strong Lyman–Werner backgrounds maintaining \(T \gtrsim 10^4\,{\rm K}\), low angular momentum, and sustained inflows to maintain quasi-star envelopes against winds [2510.11772]. A plausible implication is that the viability of the quasi-star pathway depends less on a unique hydrostatic ceiling to \(M_{\rm BH}/M_*\) than on the mass-budget competition between black-hole feeding, envelope support, and wind-driven envelope loss.

Several limitations are explicit. Physics not yet included comprises rotation, magnetic fields, jets, photon trapping, GR-radiation hydrodynamics, and explicit envelope accretion during the quasi-star phase [2510.11772]. The choice of inner boundary remains uncertain: neither Bondi nor saturated convection is definitively “correct,” and the true hydrostatic breakdown radius depends on complex multi-D physics [2510.11772]. Likewise, the super-Eddington winds are implemented through semi-empirical fits rather than full radiation hydrodynamics, so mass-loss rates and critical thresholds may change in more complete treatments [2510.11772].

Future extensions identified by the framework are rotation, magnetic fields, GR-radiation hydrodynamics, and explicit envelope accretion [2510.11772]. The earlier implementation had already identified wind mass loss, envelope accretion, and photon trapping as natural next additions made possible by MESA’s modularity [2507.12651]. In that sense, MESA-QUEST is both a research code for current quasi-star calculations and an extensible numerical framework for testing whether quasi-stars can remain a viable heavy-seed channel for the earliest SMBHs.

Source: https://www.emergentmind.com/topics/mesa-quasi-star-evolutionary-simulation-toolkit-mesa-quest