The bound origin of low-mass stellar binaries
Pith reviewed 2026-05-19 20:46 UTC · model grok-4.3
The pith
Simulations show that 70-80% of low-mass binaries form already gravitationally bound when the second star appears.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the simulations, approximately 70-80% of binaries are gravitationally bound from the instant the second star forms. Binaries evolve and accrete together, affecting their planetary systems and chemical evolution. Half of the binaries are disrupted by the end of the star-formation epoch, so that about 40% of the final single stars belonged to a multiple at some point. Formation in multiples is the dominant mode of star formation, accounting for at least 57% of stars.
What carries the argument
High-resolution simulations of star cluster formation that include all key physics and stellar feedback mechanisms and that record the gravitational binding status of stellar pairs at the instant the second star appears.
Load-bearing premise
The simulations include all relevant physics and stellar feedback at sufficient resolution to correctly determine whether a newly formed pair is gravitationally bound at the instant the second star appears.
What would settle it
High-resolution observations of the youngest embedded clusters that measure the fraction of newly formed stellar pairs with negative total energy would directly test whether the bound fraction reaches 70-80%.
Figures
read the original abstract
Most main sequence stars, unlike our Sun, belong to multiple systems with two or more stars. How and when these multiples come together and become bound is uncertain, since the earliest stages of star formation are difficult to resolve. We analyze simulations of star cluster formation in Milky Way-like conditions, including all key physics and stellar feedback mechanisms, to understand how multiple systems form. We show that $\approx 70-80\%$ of binaries are gravitationally bound from the moment the second star forms. Binaries evolve and accrete together, which will affect their planetary systems and chemical evolution. Half of the binaries are disrupted by the end of the star-formation epoch, such that $\approx40\%$ of the final single stars belonged to a multiple at some point, with implications for the stellar initial mass function. Formation in multiples is the dominant mode of star formation, accounting for at least 57% of stars.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses hydrodynamical simulations of star cluster formation under Milky Way-like conditions that incorporate all key physics and stellar feedback. It reports that ≈70-80% of binaries are gravitationally bound from the instant the second star forms, that formation in multiples is the dominant channel (accounting for at least 57% of stars), and that roughly half of binaries are later disrupted, implying that ≈40% of final single stars were once members of multiples, with consequences for the IMF and planetary-system evolution.
Significance. If the central numerical result is robust, the work would establish that the majority of low-mass stars form already bound in multiples and that dynamical disruption during the embedded phase is common. This would directly inform the initial conditions for binary populations, the origin of the field IMF, and the chemical and dynamical history of planetary systems around stars that experienced prior multiplicity.
major comments (2)
- [§3 and §4] §3 (Numerical methods) and §4 (Results on binary binding): the headline 70-80% bound fraction is obtained by classifying pairs at the moment the second sink particle is inserted. No dedicated resolution series, variation of sink insertion criteria, or softening-length tests focused on this early-time diagnostic are presented. Because the classification occurs within a few dynamical times of sink creation, under-resolution of the local Jeans length or free-fall time could systematically bias the bound/unbound assignment; this directly affects the central claim.
- [§4.1] §4.1 (Definition of bound status): the precise numerical criterion used to decide whether a newly formed pair is gravitationally bound (e.g., total energy < 0, specific binding-energy threshold, or inclusion of gas potential) is not stated explicitly. Without this definition and a demonstration that it is insensitive to accretion radius or time-step artifacts, the reported percentage cannot be reproduced or assessed for numerical robustness.
minor comments (2)
- [Abstract and §2] The abstract states that “all key physics and stellar feedback mechanisms” are included, yet the main text does not tabulate the specific feedback channels (photoionization, winds, supernovae) or their numerical implementation parameters; a concise summary table would improve clarity.
- [Figures 3–5] Figure captions and axis labels for the time-evolution plots of binary fractions should explicitly note the time at which the second star forms and the exact epoch used for the “end of the star-formation epoch” to avoid ambiguity when comparing to observations.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major point below and indicate the revisions we will make.
read point-by-point responses
-
Referee: [§3 and §4] §3 (Numerical methods) and §4 (Results on binary binding): the headline 70-80% bound fraction is obtained by classifying pairs at the moment the second sink particle is inserted. No dedicated resolution series, variation of sink insertion criteria, or softening-length tests focused on this early-time diagnostic are presented. Because the classification occurs within a few dynamical times of sink creation, under-resolution of the local Jeans length or free-fall time could systematically bias the bound/unbound assignment; this directly affects the central claim.
Authors: We agree that a dedicated resolution study focused on the early-time binding diagnostic would strengthen the result. The simulations resolve the local Jeans length by a minimum of four cells at sink insertion (see §3), and the binding classification is performed once the second sink has formed and the local gravitational potential is dominated by the stellar components. Nevertheless, we did not present a targeted convergence series varying refinement criteria or softening length specifically for this metric. In the revised manuscript we will add a paragraph in §4 justifying the robustness on the basis of the existing resolution and the short dynamical time between sink creation and classification; we will also note that a full higher-resolution suite would be a valuable extension for future work. revision: partial
-
Referee: [§4.1] §4.1 (Definition of bound status): the precise numerical criterion used to decide whether a newly formed pair is gravitationally bound (e.g., total energy < 0, specific binding-energy threshold, or inclusion of gas potential) is not stated explicitly. Without this definition and a demonstration that it is insensitive to accretion radius or time-step artifacts, the reported percentage cannot be reproduced or assessed for numerical robustness.
Authors: We thank the referee for identifying this omission in clarity. Pairs are classified as bound when the total mechanical energy (kinetic energy computed in the center-of-mass frame plus the mutual gravitational potential energy of the two sinks) is negative; the gas potential is omitted at the instant of classification because the stellar potential dominates locally. In the revised manuscript we will state this criterion explicitly in §4.1 and add a short sensitivity test, performed on the existing data, showing that the bound fraction remains stable under modest variations in measurement timing and accretion radius. revision: yes
Circularity Check
No circularity: results are direct simulation diagnostics
full rationale
The paper reports measured fractions (≈70-80% bound at second-star formation, 57% of stars in multiples) obtained by post-processing outputs from hydrodynamical cluster-formation simulations. These quantities are counted from the simulation state at sink insertion times using gravitational binding criteria applied to the hydrodynamical fields; they are not obtained by fitting parameters to the target statistic, redefining inputs in terms of outputs, or invoking self-citations whose validity depends on the present result. The derivation chain is therefore self-contained against external benchmarks and does not reduce to its own inputs by construction.
Axiom & Free-Parameter Ledger
free parameters (2)
- Initial cloud mass, density, and turbulence spectrum
- Stellar feedback efficiency parameters
axioms (1)
- domain assumption The included physics (self-gravity, hydrodynamics, magnetic fields, and stellar feedback) are sufficient to determine gravitational binding at the instant of second-star formation.
Reference graph
Works this paper leans on
-
[1]
The Evolution of Protostars: Insights from Ten Years of Infrared Surveys with Spitzer and Herschel
Dunham, M. M.et al.Beuther, H., Klessen, R. S., Dullemond, C. P. & Henning, T. (eds)The Evolution of Protostars: Insights from Ten Years of Infrared Surveys with Spitzer and Herschel. (eds Beuther, H., Klessen, R. S., Dullemond, C. P. & Henning, T.)Protostars and Planets VI, 195–218 (2014). 1401.1809
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[2]
Protostars and Planets VII , year = 2023, editor =
Pineda, J. E.et al.Inutsuka, S., Aikawa, Y., Muto, T., Tomida, K. & Tamura, M. (eds)From Bubbles and Filaments to Cores and Disks: Gas Gathering and Growth of Structure Leading to the Formation of Stellar Systems. (eds Inutsuka, S., Aikawa, Y., Muto, T., Tomida, K. & Tamura, M.)Protostars and Planets VII, Vol. 534 ofAstronomical Society of the Pacific Con...
-
[3]
Offner, S. S. R.et al.Inutsuka, S., Aikawa, Y., Muto, T., Tomida, K. & Tamura, M. (eds)The Origin and Evolution of Multiple Star Systems. (eds Inutsuka, S., Aikawa, Y., Muto, T., Tomida, K. & Tamura, M.)Astronomical Society of the Pacific Conference Series, Vol. 534 ofAstronomical Society of the Pacific Conference Series, 275 (2023)
work page 2023
-
[4]
Chevance, M.et al.The Molecular Cloud Lifecycle.Space Sci. Rev.216, 50 (2020)
work page 2020
-
[5]
Bastian, N.et al.Evidence for environmentally dependent cluster disruption in M83.Mon. Not. R. Astron. Soc.417, L6–L10 (2011)
work page 2011
-
[6]
Offner, S. S. R.et al.Beuther, H., Klessen, R. S., Dullemond, C. P. & Henning, T. (eds)The Origin and Universality of the Stellar Initial Mass Function. (eds Beuther, H., Klessen, R. S., Dullemond, C. P. & Henning, T.)Protostars and Planets VI, 53–75 (2014). 1312.5326
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[8]
Guszejnov, D.et al.Effects of the environment on the multiplicity properties of stars in the STARFORGE simulations.Mon. Not. R. Astron. Soc.518, 4693–4712 (2023)
work page 2023
-
[9]
Iben, J., I. & Tutukov, A. V. Supernovae of type I as end products of the evolution of binaries with components of moderate initial mass.Astrophys. J. Suppl. Ser. 54, 335–372 (1984)
work page 1984
-
[10]
Li, L.-X. & Paczy´ nski, B. Transient Events from Neutron Star Mergers. Astrophys. J. Lett.507, L59–L62 (1998)
work page 1998
-
[11]
Portegies Zwart, S. F. & McMillan, S. L. W. Black Hole Mergers in the Universe. Astrophys. J. Lett.528, L17–L20 (2000)
work page 2000
-
[12]
Sana, H.et al.Binary Interaction Dominates the Evolution of Massive Stars. Science337, 444 (2012)
work page 2012
-
[13]
Conroy, C. & Kratter, K. M. Runaway Stars and the Escape of Ionizing Radiation from High-redshift Galaxies.Astrophys. J.755, 123 (2012)
work page 2012
-
[14]
Duchˆ ene, G. & Kraus, A. Stellar Multiplicity.Annu. Rev. Astron. Astrophys. 51, 269–310 (2013)
work page 2013
-
[15]
Rafikov, R. R. Building Tatooine: Suppression of the Direct Secular Excitation in Kepler Circumbinary Planet Formation.Astrophys. J. Lett.764, L16 (2013)
work page 2013
-
[16]
Moe, M. & Kratter, K. M. Impact of binary stars on planet statistics - I. Planet 16 occurrence rates and trends with stellar mass.Mon. Not. R. Astron. Soc.507, 3593–3611 (2021)
work page 2021
-
[17]
Schnee, S.et al.An Observed Lack of Substructure in Starless Cores.Astrophys. J.718, 306–313 (2010)
work page 2010
-
[18]
Bate, M. R. Predicting the properties of binary stellar systems: the evolution of accreting protobinary systems.Mon. Not. R. Astron. Soc.314, 33–53 (2000)
work page 2000
-
[19]
Bate, M. R. Stellar, brown dwarf and multiple star properties from a radiation hydrodynamical simulation of star cluster formation.Mon. Not. R. Astron. Soc. 419, 3115–3146 (2012)
work page 2012
-
[20]
Krumholz, M. R., Klein, R. I. & McKee, C. F. Radiation-hydrodynamic Simula- tions of the Formation of Orion-like Star Clusters. II. The Initial Mass Function from Winds, Turbulence, and Radiation.Astrophys. J.754, 71 (2012)
work page 2012
-
[21]
Li, P. S., Klein, R. I. & McKee, C. F. Formation of stellar clusters in magnetized, filamentary infrared dark clouds.Monthly Notices of the Royal Astronomical Society473, 4220–4241 (2017). URL https://doi.org/10.1093/mnras/stx2611
-
[24]
Kuruwita, R. L. & Haugbølle, T. The contribution of binary star formation via core fragmentation on protostellar multiplicity.Astron. Astrophys.674, A196 (2023)
work page 2023
-
[25]
On the Fragmentation of Gas Clouds Into Galaxies and Stars
Hoyle, F. On the Fragmentation of Gas Clouds Into Galaxies and Stars. Astrophys. J.118, 513–+ (1953)
work page 1953
-
[26]
Tsuribe, T. & Inutsuka, S.-I. Criteria for Fragmentation of Rotating Isothermal Clouds. I. Semianalytic Approach.Astrophys. J.526, 307–313 (1999)
work page 1999
-
[27]
Boss, A. P., Fisher, R. T., Klein, R. I. & McKee, C. F. The Jeans Condition and Collapsing Molecular Cloud Cores: Filaments or Binaries?Astrophys. J.528, 325–335 (2000)
work page 2000
-
[28]
Padoan, P. & Nordlund, ˚A. The Stellar Initial Mass Function from Turbulent Fragmentation.Astrophys. J.576, 870–879 (2002). 17
work page 2002
-
[29]
T.Single and multiple star formation in turbulent molecular cloud cores
Fisher, R. T.Single and multiple star formation in turbulent molecular cloud cores. Ph.D. thesis, AA(UNIVERSITY OF CALIFORNIA, BERKELEY) (2002)
work page 2002
-
[30]
Offner, S. S. R., Kratter, K. M., Matzner, C. D., Krumholz, M. R. & Klein, R. I. The Formation of Low-mass Binary Star Systems Via Turbulent Fragmentation. Astrophys. J.725, 1485–1494 (2010)
work page 2010
-
[31]
Kuffmeier, M., Calcutt, H. & Kristensen, L. E. The bridge: a transient phe- nomenon of forming stellar multiples. Sequential formation of stellar companions in filaments around young protostars.Astron. Astrophys.628, A112 (2019)
work page 2019
-
[32]
Rozner, M., Generozov, A. & Perets, H. B. Binary formation through gas-assisted capture and the implications for stellar, planetary, and compact object evolution. Mon. Not. R. Astron. Soc.521, 866–880 (2023)
work page 2023
-
[33]
URL https: //doi.org/10.1093/mnras/staa3902
Cournoyer-Cloutier, C.et al.Implementing primordial binaries in simulations of star cluster formation with a hybrid mhd and direct n-body method.Monthly Notices of the Royal Astronomical Society501, 4464–4478 (2020). URL https: //doi.org/10.1093/mnras/staa3902
-
[34]
Grudi´ c, M. Y.et al.The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert.Mon. Not. R. Astron. Soc.512, 216–232 (2022)
work page 2022
-
[35]
Guszejnov, D.et al.Effects of the environment and feedback physics on the initial mass function of stars in the STARFORGE simulations.Mon. Not. R. Astron. Soc.515, 4929–4952 (2022)
work page 2022
-
[36]
Lada, C. J. & Dame, T. M. The Mass-Size Relation and the Constancy of GMC Surface Densities in the Milky Way.Astrophys. J.898, 3 (2020)
work page 2020
-
[37]
Larson, R. B. Turbulence and star formation in molecular clouds.Monthly Notices of the Royal Astronomical Society194, 809–826 (1981). URL https: //doi.org/10.1093/mnras/194.4.809
-
[38]
Protostars and Planets VII , year = 2023, editor =
Chevance, M.et al.Inutsuka, S., Aikawa, Y., Muto, T., Tomida, K. & Tamura, M. (eds)The Life and Times of Giant Molecular Clouds. (eds Inutsuka, S., Aikawa, Y., Muto, T., Tomida, K. & Tamura, M.)Protostars and Planets VII, Vol. 534 ofAstronomical Society of the Pacific Conference Series, 1 (2023). 2203.09570
-
[39]
M., Wandelt, B., Heiles, C., Falgarone, E
Crutcher, R. M., Wandelt, B., Heiles, C., Falgarone, E. & Troland, T. H. Magnetic Fields in Interstellar Clouds from Zeeman Observations: Inference of Total Field Strengths by Bayesian Analysis.Astrophys. J.725, 466–479 (2010)
work page 2010
-
[40]
Grudi´ c, M. Y.et al.A model for the formation of stellar associations and clusters from giant molecular clouds.Mon. Not. R. Astron. Soc.506, 3239–3258 (2021). 18
work page 2021
-
[41]
A.et al.The Spatial Structure of Young Stellar Clusters
Kuhn, M. A.et al.The Spatial Structure of Young Stellar Clusters. I. Subclusters. Astrophys. J.787, 107 (2014)
work page 2014
-
[42]
Gouliermis, D. A. Unbound Young Stellar Systems: Star Formation on the Loose. Publ. Astron. Soc. Pac.130, 072001 (2018)
work page 2018
-
[43]
Ward, J. L., Kruijssen, J. M. D. & Rix, H.-W. Not all stars form in clusters – gaia-dr2 uncovers the origin of ob associations.Monthly Notices of the Royal Astronomical Society495, 663–685 (2020). URL https://doi.org/10.1093/mnras/ staa1056
-
[44]
Guszejnov, D., Hopkins, P. F. & Krumholz, M. R. Protostellar feedback in tur- bulent fragmentation: consequences for stellar clustering and multiplicity.Mon. Not. R. Astron. Soc.468, 4093–4106 (2017)
work page 2017
-
[45]
Kristensen, L. E. & Dunham, M. M. Protostellar half-life: new methodology and estimates.Astron. Astrophys.618, A158 (2018)
work page 2018
-
[46]
Moe, M. & Kratter, K. M. Dynamical Formation of Close Binaries during the Pre-main-sequence Phase.Astrophys. J.854, 44 (2018)
work page 2018
-
[47]
Farias, J. P., Offner, S. S. R., Grudi´ c, M. Y., Guszejnov, D. & Rosen, A. L. Stellar populations in STARFORGE: the origin and evolution of star clusters and associations.Mon. Not. R. Astron. Soc.527, 6732–6751 (2024)
work page 2024
-
[48]
Farias, J. P., Offner, S. S. R., Kerr, R. & Grudi´ c, M. Y. Stellar populations in STARFORGE II: Comparison with observations.arXiv e-printsarXiv:2506.00240 (2025)
-
[49]
Cournoyer-Cloutier, C.et al.Massive Star Cluster Formation with Binaries. I. Evolution of Binary Populations.Astrophys. J.977, 203 (2024)
work page 2024
-
[51]
G.et al.The Solar Neighborhood
Winters, J. G.et al.The Solar Neighborhood. XLV. The Stellar Multiplicity Rate of M Dwarfs Within 25 pc.Astron. J.157, 216 (2019)
work page 2019
-
[52]
Offner, S. S. R.et al.Inutsuka, S., Aikawa, Y., Muto, T., Tomida, K. & Tamura, M. (eds)The Origin and Evolution of Multiple Star Systems. (eds Inutsuka, S., Aikawa, Y., Muto, T., Tomida, K. & Tamura, M.)Protostars and Planets VII, Vol. 534 ofAstronomical Society of the Pacific Conference Series, 275 (2023). 2203.10066
work page internal anchor Pith review Pith/arXiv arXiv 2023
- [53]
-
[54]
Heggie, D. C. Binary evolution in stellar dynamics.Mon. Not. R. Astron. Soc. 173, 729–787 (1975)
work page 1975
-
[55]
Padoan, P., Pan, L., Juvela, M., Haugbølle, T. & Nordlund, ˚A. The Origin of Massive Stars: The Inertial-inflow Model.Astrophys. J.900, 82 (2020)
work page 2020
-
[56]
Moe, M. & Di Stefano, R. Mind Your Ps and Qs: The Interrelation between Period (P) and Mass-ratio (Q) Distributions of Binary Stars.Astrophys. J. Suppl. Ser.230, 15 (2017)
work page 2017
-
[57]
On the variation of the initial mass function.Mon
Kroupa, P. On the variation of the initial mass function.Mon. Not. R. Astron. Soc.322, 231–246 (2001). URL http://adsabs.harvard.edu/cgi-bin/ nph-bib query?bibcode=2001MNRAS.322..231K&db key=AST
work page 2001
-
[58]
Hopkins, P. F. & Raives, M. J. Accurate, meshless methods for magnetohydro- dynamics.Mon. Not. R. Astron. Soc.455, 51–88 (2016)
work page 2016
-
[60]
Klessen, R. S. & Burkert, A. The formation of stellar clusters: Gaussian cloud conditions. i.The Astrophysical Journal Supplement Series128, 287 (2000). URL https://dx.doi.org/10.1086/313371
-
[61]
Bate, M. R., Bonnell, I. A. & Bromm, V. The formation of a star cluster: pre- dicting the properties of stars and brown dwarfs.Mon. Not. R. Astron. Soc.339, 577–599 (2003)
work page 2003
-
[62]
Bate, M. R., Bonnell, I. A. & Price, N. M. Modelling accretion in protobinary systems.Mon. Not. R. Astron. Soc.277, 362–376 (1995)
work page 1995
-
[63]
Grudi´ c, M. Y., Guszejnov, D., Hopkins, P. F., Offner, S. S. R. & Faucher-Gigu` ere, C.-A. STARFORGE: Towards a comprehensive numerical model of star cluster formation and feedback.Mon. Not. R. Astron. Soc.506, 2199–2231 (2021)
work page 2021
-
[64]
Offner, S. S. R., Klein, R. I., McKee, C. F. & Krumholz, M. R. The Effects of Radiative Transfer on Low-Mass Star Formation.Astrophys. J.703, 131–149 (2009)
work page 2009
-
[65]
Grudi´ c, M. Y. & Guszejnov, D. MakeCloud (2021). URL https://github.com/ mikegrudic/MakeCloud
work page 2021
-
[66]
Ostriker, E. C., Stone, J. M. & Gammie, C. F. Density, Velocity, and Magnetic Field Structure in Turbulent Molecular Cloud Models.Astrophys. J.546, 980– 1005 (2001). 20
work page 2001
-
[67]
Bate, M. R. Stellar, brown dwarf and multiple star properties from hydrody- namical simulations of star cluster formation.Mon. Not. R. Astron. Soc.392, 590–616 (2009)
work page 2009
-
[68]
Lee, A. T., Offner, S. S. R., Kratter, K. M., Smullen, R. A. & Li, P. S. The Formation and Evolution of Wide-orbit Stellar Multiples In Magnetized Clouds. Astrophys. J.887, 232 (2019)
work page 2019
-
[69]
Grishin, E., Perets, H. B., Zenati, Y. & Michaely, E. Generalized Hill-stability criteria for hierarchical three-body systems at arbitrary inclinations.Mon. Not. R. Astron. Soc.466, 276–285 (2017)
work page 2017
-
[70]
Offner, S. S. R., Taylor, J. & Grud´ ıc, M. Y. The Life and Times of Star-forming Cores: An Analysis of Dense Gas in the STARFORGE Simulations.Astrophys. J.982, 138 (2025). 21 Extended data alpha2 alpha1 alpha4 alpha2,z0.1 seeds=3 seeds=2 seeds=2 seeds=1 Count Frac Count Frac Count Frac Count Frac All Total 1792 545 600 647 1.00 1370 551 819 1.00 994 544 ...
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.