---
title: 'StarTrack: Modular Binary-Population Synthesis'
url: https://www.emergentmind.com/topics/startrack
type: topic
---

# StarTrack: Modular Binary-Population Synthesis

StarTrack is a modular, Monte Carlo binary-population synthesis code for evolving large ensembles of zero-age main-sequence binaries through stellar winds, Roche-lobe overflow, common-envelope evolution, supernovae, compact-object formation, and, where relevant, gravitational-wave inspiral. In the cited literature it is used to predict populations of double compact objects, X-ray binaries, ultraluminous X-ray sources, Galactic black holes, non-interacting black-hole binaries, and wide double white dwarfs, and to connect those populations to gravitational-wave, electromagnetic, and astrometric observables [1202.4901, 1409.8360, 1306.5265, 2008.04890, 1909.04435, 2006.08317, 1209.0175].

## 1. Code identity and scope

Published descriptions characterize StarTrack as a binary-evolution population-synthesis framework that follows either \(2\times 10^6\) zero-age binaries or, in other applications, ensembles of \(\sim 10^6\!-\!10^7\) binaries, depending on the study design and the quantity being modeled [1202.4901, 1409.8360]. The code is repeatedly described as modular, and its core remit is isolated stellar and binary evolution rather than \(N\)-body cluster dynamics or purely phenomenological rate fitting [1202.4901].

The underlying workflow is consistent across applications. Primordial binaries are drawn from prescribed initial distributions, evolved through mass loss and binary interaction, subjected to compact-remnant formation prescriptions and natal-kick models, and then filtered according to the scientific target: compact binaries that merge within a Hubble time, accreting systems that contribute to X-ray luminosity functions, black-hole binaries detectable by Gaia or LAMOST, or wide white-dwarf pairs that survive Galactic perturbations [1306.5265, 2006.08317, 1209.0175].

StarTrack is also routinely embedded in larger population models. In cosmological calculations it is combined with metallicity-dependent star-formation histories and detector selection functions to produce redshift-dependent merger-rate densities and detectable-event rates [1409.8360, 2008.04890]. In Milky Way applications it is convolved with disk, bulge, and halo star-formation histories and chemical-evolution models to produce synthetic present-day Galactic catalogs [1908.08775, 2006.08317]. This dual role—as an evolutionary engine and as the kernel of a synthetic-universe pipeline—is one of the defining features of the code.

## 2. Evolutionary prescriptions and governing equations

Across the cited studies, StarTrack uses analytic single-star evolution fits, metallicity-dependent wind prescriptions, stability criteria for Roche-lobe overflow, an \(\alpha\)-\(\lambda\) treatment of common-envelope evolution, supernova remnant prescriptions of the “rapid” and “delayed” type, and natal kicks drawn from Maxwellian distributions with fallback-dependent suppression for black holes in many model families [1306.5265, 1202.4901, 1902.07718]. Several papers also include pair-instability and pulsational pair-instability prescriptions in the remnant-mass calculation [1902.07718, 1908.08775].

Initial conditions are not fixed by a single universal StarTrack setup; instead they vary by application. Recurring choices include Kroupa-like or three-segment Kroupa-type initial mass functions, flat mass-ratio distributions, orbital separations distributed either flat in \(\log a\) or via \(p(\log P)\propto(\log P)^{-0.55}\), and eccentricity distributions that are either thermal, \(f(e)=2e\), or \(p(e)\propto e^{-0.42}\) [1202.4901, 1010.0511, 1908.08775]. The presence of multiple input families is not incidental: it reflects the code’s role as a parametric synthesis framework rather than a single fixed binary-evolution realization.

A central equation is the Webbink-style common-envelope energy balance. In one standard form used in StarTrack studies,
$$
\alpha_{\rm CE}
\Bigl[
\frac{G\,M_{\rm d,f}\,M_{\rm a}}{2\,a_f}
-
\frac{G\,M_{\rm d,i}\,M_{\rm a}}{2\,a_i}
\Bigr]
=
\frac{G\,M_{\rm d,i}\,M_{\rm d,env}}{\lambda\,R_{\rm L}} ,
$$
where \(M_{\rm d,i}\) and \(M_{\rm d,f}\) are the donor masses before and after envelope ejection, \(M_{\rm a}\) is the accretor mass, \(a_i\) and \(a_f\) are the pre- and post-common-envelope separations, \(M_{\rm d,env}\) is the envelope mass, and \(\lambda\) is the binding-energy parameter [1202.4901].

Natal kicks are likewise encoded in compact formulas. For fallback-modulated black-hole kicks,
$$
V_k = V_{\rm max}(1-f_{\rm fb}),
$$
with \(V_{\rm max}\) drawn from a Maxwellian and \(f_{\rm fb}\) the fallback fraction [1202.4901]. Once a double compact object forms, gravitational-wave merger times can be computed with the Peters prescription,
$$
t_{\rm GW}
=
\frac{5}{256}\,
\frac{c^5\,a_f^4}{G^3\,M_1\,M_2\,(M_1+M_2)}
(1-e^2)^{7/2},
$$
which links post-interaction orbital architecture directly to merger observability [1202.4901].

These prescriptions give StarTrack its characteristic structure: binary-population outputs emerge from the coupling of mass loss, remnant formation, binary survival, and orbital shrinkage. The code’s scientific utility derives from this coupling, but so does its sensitivity to uncertain physics.

## 3. Common-envelope treatment and donor-structure criteria

Common-envelope evolution is the dominant structural uncertainty in much of the StarTrack literature. Earlier StarTrack-based compact-binary studies often bracketed the problem with two contrasting assumptions about Hertzsprung-gap donors: either a Hertzsprung-gap common envelope always leads to merger, or such donors are allowed to survive via the usual \(\alpha\)-\(\lambda\) formalism [1409.8360, 1202.4901]. This single modeling choice already shifts predicted merger rates by factors of a few to orders of magnitude, especially for BH–BH and BH–NS systems [1202.4901].

A later revision made the criterion more restrictive for massive donors. In the implementation described in "The impact of common envelope development criteria on the formation of LIGO/Virgo sources" [2102.05649], common-envelope onset for evolved H-rich giants with \(M_{\rm ZAMS}>18\,M_\odot\) depends on radius thresholds \(R_S(M_{\rm don},Z)\) and \(R_U(M_{\rm don},Z)\), together with a mass-ratio cutoff \(q_{\rm CE}(M_{\rm don},Z)\). In the three illustrative models reported there, the local BH–BH merger rate changes from \(62\) to \(88\) to \(18\ {\rm Gpc}^{-3}\,{\rm yr}^{-1}\), while NS–NS rates change from \(148\) to \(148\) to \(322\ {\rm Gpc}^{-3}\,{\rm yr}^{-1}\) [2102.05649]. The point is not merely numerical spread; it is that modest modifications to the common-envelope development rule can move the dominant formation channel from common-envelope–assisted evolution to channels with no common-envelope phase at all.

An even more explicit donor-structure criterion appears in "Development of convective envelopes in massive stars: Implications for gravitational wave sources" [2410.17315]. There, a star is defined to have developed an outer convective envelope once the mass in outer convection zones exceeds \(10\%\) of the total H-rich envelope mass. StarTrack then replaces the earlier evolutionary-type precheck in `preCEOutcome()` with an envelope-type precheck based on whether
\[
T_{\rm eff,donor}\le T_{\rm eff,conv}(L_{\rm donor},Z_{\rm donor}) .
\]
If the condition is satisfied, the donor is treated as having a convective envelope and the system proceeds to the \(\alpha\lambda\)-formalism; otherwise the radiative-envelope case is assigned an immediate merger, producing a Thorne–Żytkow object or quasi-star when the accretor is a neutron star or black hole [2410.17315]. The same study adds new data tables in `/data/stars/`, modifies `src/common_envelope.F90` and `src/binary_evolve.F90`, and records the implementation under the commit tag `conv_env_2024` [2410.17315].

The astrophysical consequences are substantial. In the \(M_{\alpha{\rm ML}1.5\_RMAX}\) model, the local BH–BH merger rate for total masses \(20\!-\!50\,M_\odot\) is reduced by \(\sim 20\times\), from \(\sim 66\) to \(3\ {\rm Gpc}^{-3}\,{\rm yr}^{-1}\), and the total-mass distribution becomes bimodal, with peaks at \(\sim 10\!-\!15\,M_\odot\) and \(\sim 30\!-\!60\,M_\odot\) [2410.17315]. Under the strong red-supergiant pulsation scenario, both the TŻO and quasi-star populations in the Galaxy collapse to \(<1\) object in the present epoch, making detection essentially unfeasible [2410.17315]. These results strongly suggest that StarTrack’s common-envelope channel is no longer well represented by any single Hertzsprung-gap rule; donor envelope structure has become a code-level state variable of primary importance.

## 4. Synthetic-universe construction and statistical inference

StarTrack outputs are often post-processed into cosmological rate models. In one formulation, each simulation yields a discrete merger population
$$
\rho_{\Lambda}(\theta)=\sum_i s_i\,\delta(\theta-\theta_i),
$$
where \(\theta\) denotes binary parameters and \(s_i\) are weights accounting for metallicity bins and cosmic assembly. Detector selection is then applied through
$$
R_{\Lambda}(\theta)
=
\rho_{\Lambda}(\theta)\,
p_{\rm det}(\theta)\,
\frac{dV_c}{dz}\,
(1+z)^{-1},
$$
and expected counts follow by integration over masses and redshift [2303.05436]. This framework turns StarTrack from a forward binary-evolution code into a generative model for observed compact-binary catalogs.

Earlier work uses an analogous construction for strong-lensing forecasts with the Einstein Telescope. There the redshift-dependent intrinsic merger density \(\dot n_0(z_s)\) from StarTrack is combined with the detector selection function and a singular isothermal sphere lens population to obtain the yearly rate of strongly lensed, detectable mergers [1409.8360]. The same logic reappears in stochastic-background calculations, where metallicity-resolved StarTrack outputs are convolved with cosmic star-formation and metallicity histories and then integrated over source redshift to compute \(\Omega_{\rm GW}(f)\) [2008.04890].

StarTrack has also been coupled directly to Bayesian population inference. In the inhomogeneous-Poisson formulation used for comparison with gravitational-wave catalogs,
$$
p(d|\Lambda)
=
e^{-\mu_\Lambda}
\prod_{\alpha}
\frac{\mu_{\Lambda,\alpha}^{N_\alpha}}{N_\alpha!}
\prod_j
\int
\bar{\rho}_\Lambda(\theta)\,
\mathcal{L}_j(\theta)\,d\theta ,
$$
where \(\mu_{\Lambda,\alpha}\) are expected counts by merger type and \(\mathcal{L}_j(\theta)\) are event-level likelihoods [2303.05436]. Closely related work replaces direct posterior-sample handling with bounded multivariate normal “Normal Approximate Likelihood” fits, enabling fast evaluation of event likelihoods across a large StarTrack simulation bank [2209.03790].

These developments are methodologically significant because they transform StarTrack from a source of tabulated rates into a parameterized synthetic-universe engine. Once detector selection, cosmology, and event likelihoods are included, code parameters such as kick dispersion, mass-transfer efficiency, and wind scaling can be compared directly with gravitational-wave observations rather than only with population-level summary statistics.

## 5. Scientific applications

StarTrack’s most visible application is the prediction of compact-object merger populations. In strong-lensing forecasts for the Einstein Telescope, the expected rate is about \(50\!-\!100\) strongly lensed inspiral events per year, with the BH–BH channel dominating the lensed yield and the high-BH-kicks scenario reducing the prediction to only a few events per year [1409.8360]. In stochastic-background calculations, the population I/II contribution reaches \(\Omega_{\rm GW}\sim 1.0\times10^{-9}\) at \(25\) Hz, making the background detectable at \(3\sigma\) after about \(7\) years of observation with current-generation ground-based detectors at design sensitivity; population III contributes about one order of magnitude less to the total background but dominates the residual background in 3G detectors after subtraction of resolvable sources [2008.04890].

A second major application is X-ray binary and ultraluminous X-ray source modeling. In the SINGS comparison, StarTrack-based theoretical X-ray luminosity functions are convolved with galaxy star-formation histories, and the best global models are consistent with \(\alpha_{\rm CE}\lambda\approx 0.1\) and a \(50\%\) uniform–\(50\%\) “twins” initial mass-ratio distribution [1306.5265]. In the ULX–double-compact-object study, typically \(50\%\) of merging BH–BH progenitor binaries are found to have evolved through a ULX phase, while the fraction of observed ULXs that will form merging double compact objects in the future varies between \(5\%\) and \(40\%\), depending on common-envelope model and metallicity [1909.04435].

StarTrack is also used for Galactic census problems. A Milky Way synthetic catalog predicts that the present-day Galaxy contains about \(1.2\times 10^8\) single black holes with average mass about \(14\,M_\odot\), and quotes current Galactic merger rates of \(3\!-\!81\ {\rm Myr}^{-1}\) for BH–BH, \(1\!-\!9\ {\rm Myr}^{-1}\) for BH–NS, and \(14\!-\!59\ {\rm Myr}^{-1}\) for NS–NS systems across two common-envelope models [1908.08775]. In astrometric and spectroscopic detectability studies, Gaia is predicted to observe \(\sim 41\!-\!340\) non-interacting black-hole binaries, falling to \(\sim 10\!-\!70\) if the recent thin-disk star formation is low, while LAMOST is expected to detect \(\lesssim 14\) such systems [2006.08317]. For white-dwarf binaries, StarTrack-based calculations indicate a significant observable population of wide WDWDs with orbital separations \(10^2\!-\!10^5\) AU, and the corresponding SDSS search identified twelve high-confidence wide WDWD pairs [1209.0175].

A plausible implication is that StarTrack’s scientific range is unusually broad not because the code solves one narrowly defined rate problem, but because its outputs can be reweighted, projected, and selection-filtered into many different observational spaces. The same binary-evolution core can therefore feed gravitational-wave cosmology, X-ray population synthesis, Galactic stellar-remnant cartography, and survey-yield forecasting.

## 6. Uncertainties, controversies, and name collisions

The most persistent conclusion across the StarTrack literature is that predictions are dominated by a small number of uncertain physical ingredients. Varying the common-envelope binding parameter \(\lambda\) between \(0.01\) and \(10\) shifts merger rates by up to \(\sim 2\) orders of magnitude, and adopting full black-hole kicks suppresses BH–BH rates by \(2\!-\!3\) dex [1202.4901]. Later work on very massive binaries shows that different assumptions about common-envelope survival, mass and angular-momentum loss, stellar mixing, pair-instability mass loss, and supernova outbursts can send the same observed system toward a close BH–BH merger, a wide BH–BH binary, a Thorne–Żytkow object, or total disruption by pair-instability supernovae [2108.10885]. The convective-envelope revision goes further, arguing that the common-envelope channel for BH–BH mergers with \(M_{\rm tot}>20\,M_\odot\) has been considerably overestimated and that predictions for systems above \(50\,M_\odot\) hinge heavily on limited understanding of stellar structure and mass loss close to the Eddington limit [2410.17315].

Observation-driven inference does not eliminate this uncertainty. One StarTrack-based analysis that varies \(\sigma_{\rm eff}\), \(f_a\), \(\beta\), and \(f_{\rm wind1}\) finds a maximum-likelihood model \(K0559\) with \(f_a=0.922\pm 0.019\), \(\beta=0.768\pm 0.18\), \(\sigma_{\rm eff}=108.3\pm 10.8\ {\rm km\,s}^{-1}\), and \(f_{\rm wind1}=0.328\pm 0.022\) [2303.05436]. Another isolated-binary analysis, using a different simulation family and likelihood approximation, finds a 4D likelihood peak at model \(D411\) with \(f_a\approx 0.19\), \(\beta\approx 0.94\), \(\sigma_{\rm kick}\approx 206\ {\rm km\,s}^{-1}\), and \(f_{\rm wind}\approx 1.00\) [2209.03790]. These are not directly commensurate constraints, but together they show that StarTrack-based inference remains conditional on model family, parameterization, and post-processing assumptions.

A separate source of confusion is lexical rather than astrophysical. "S.T.A.R.-Track" [2306.17602] is an object-centric, transformer-based framework for end-to-end 3D object tracking; on the nuScenes validation split it reports AMOTA \(=37.9\%\) and IDS \(=372\), and on the test split with a stronger backbone it reports AMOTA \(=43.9\%\). "The LSPE-Strip Pointing Reconstruction and Star Tracker" [2501.05604] describes an optical pointing instrument for a CMB telescope, with a prototype RMS accuracy of approximately \(3\) arcseconds and systematic errors below \(10\) arcseconds. These systems are unrelated to the astrophysical StarTrack population-synthesis code despite the near-identical naming.

Taken together, the literature presents StarTrack less as a finished predictive theory than as a highly structured computational hypothesis space for isolated binary evolution. Its continuing importance lies precisely in that role: it makes uncertain stellar and binary physics explicit, parameterized, and testable against increasingly heterogeneous data.

Source: https://www.emergentmind.com/topics/startrack