pith. sign in

arxiv: 2605.16657 · v1 · pith:EEWGAKIUnew · submitted 2026-05-15 · 🌌 astro-ph.SR · astro-ph.GA

The bound origin of low-mass stellar binaries

Pith reviewed 2026-05-19 20:46 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.GA
keywords stellar binariesstar formationmultiple systemsgravitational bindingstellar feedbackinitial mass functionnumerical simulationscluster formation
0
0 comments X

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.

The paper uses detailed numerical simulations of star cluster formation in Milky Way-like conditions, incorporating stellar feedback and other key physics, to track how multiple star systems arise. It establishes that most binaries become gravitationally bound at the exact moment the second star forms rather than through later capture. These pairs then accrete material and evolve together until dynamical interactions disrupt roughly half of them by the end of the star-formation period. As a result, formation inside multiples becomes the main channel for stars, involving at least 57% of all stars, and about 40% of today's single stars once belonged to such systems. This picture directly shapes expectations for the stellar initial mass function and for the early environments of planetary systems.

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

Figures reproduced from arXiv: 2605.16657 by A. Generozov, D. Guszejnov, H.B. Perets, K.M. Kratter, M.Y. Grudi\'c, S.S.R. Offner.

Figure 1
Figure 1. Figure 1: Starforge uses a numerical framework implemented in the Gizmo code. The public version of the Gizmo code, which includes self-gravity, MHD, radiation transfer and various other physics modules is available at https://bitbucket.org/phopkins/ gizmo-public/src/master/. Acknowledgements. We thank the anonymous referees for constructive feedback that dramatically improved the quality of the paper. We thank Juan… view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [§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.
  2. [§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)
  1. [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.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

2 free parameters · 1 axioms · 0 invented entities

The central percentages rest on the assumption that the chosen hydrodynamical code and sub-grid feedback prescriptions faithfully capture the binding dynamics at the moment of second-star formation.

free parameters (2)
  • Initial cloud mass, density, and turbulence spectrum
    Chosen to represent Milky Way-like conditions; values are not derived from first principles within the paper.
  • Stellar feedback efficiency parameters
    Sub-grid parameters controlling radiation, winds, and supernovae that affect gas dispersal and binary survival.
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.
    Invoked when the authors state that the simulations contain 'all key physics'.

pith-pipeline@v0.9.0 · 5711 in / 1291 out tokens · 35165 ms · 2026-05-19T20:46:16.607947+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

65 extracted references · 65 canonical work pages · 3 internal anchors

  1. [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

  2. [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. [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)

  4. [4]

    Rev.216, 50 (2020)

    Chevance, M.et al.The Molecular Cloud Lifecycle.Space Sci. Rev.216, 50 (2020)

  5. [5]

    Bastian, N.et al.Evidence for environmentally dependent cluster disruption in M83.Mon. Not. R. Astron. Soc.417, L6–L10 (2011)

  6. [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

  7. [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)

  8. [9]

    & Tutukov, A

    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)

  9. [10]

    & Paczy´ nski, B

    Li, L.-X. & Paczy´ nski, B. Transient Events from Neutron Star Mergers. Astrophys. J. Lett.507, L59–L62 (1998)

  10. [11]

    Portegies Zwart, S. F. & McMillan, S. L. W. Black Hole Mergers in the Universe. Astrophys. J. Lett.528, L17–L20 (2000)

  11. [12]

    Science337, 444 (2012)

    Sana, H.et al.Binary Interaction Dominates the Evolution of Massive Stars. Science337, 444 (2012)

  12. [13]

    & Kratter, K

    Conroy, C. & Kratter, K. M. Runaway Stars and the Escape of Ionizing Radiation from High-redshift Galaxies.Astrophys. J.755, 123 (2012)

  13. [14]

    & Kraus, A

    Duchˆ ene, G. & Kraus, A. Stellar Multiplicity.Annu. Rev. Astron. Astrophys. 51, 269–310 (2013)

  14. [15]

    Rafikov, R. R. Building Tatooine: Suppression of the Direct Secular Excitation in Kepler Circumbinary Planet Formation.Astrophys. J. Lett.764, L16 (2013)

  15. [16]

    & Kratter, K

    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)

  16. [17]

    J.718, 306–313 (2010)

    Schnee, S.et al.An Observed Lack of Substructure in Starless Cores.Astrophys. J.718, 306–313 (2010)

  17. [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)

  18. [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)

  19. [20]

    R., Klein, R

    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)

  20. [21]

    S., Klein, R

    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

  21. [24]

    Kuruwita, R. L. & Haugbølle, T. The contribution of binary star formation via core fragmentation on protostellar multiplicity.Astron. Astrophys.674, A196 (2023)

  22. [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)

  23. [26]

    & Inutsuka, S.-I

    Tsuribe, T. & Inutsuka, S.-I. Criteria for Fragmentation of Rotating Isothermal Clouds. I. Semianalytic Approach.Astrophys. J.526, 307–313 (1999)

  24. [27]

    P., Fisher, R

    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)

  25. [28]

    & Nordlund, ˚A

    Padoan, P. & Nordlund, ˚A. The Stellar Initial Mass Function from Turbulent Fragmentation.Astrophys. J.576, 870–879 (2002). 17

  26. [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)

  27. [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)

  28. [31]

    & Kristensen, L

    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)

  29. [32]

    & Perets, H

    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)

  30. [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

  31. [34]

    Y.et al.The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert.Mon

    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)

  32. [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)

  33. [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)

  34. [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

  35. [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

  36. [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)

  37. [40]

    Y.et al.A model for the formation of stellar associations and clusters from giant molecular clouds.Mon

    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

  38. [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)

  39. [42]

    Gouliermis, D. A. Unbound Young Stellar Systems: Star Formation on the Loose. Publ. Astron. Soc. Pac.130, 072001 (2018)

  40. [43]

    Rajabi and M

    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

  41. [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)

  42. [45]

    Kristensen, L. E. & Dunham, M. M. Protostellar half-life: new methodology and estimates.Astron. Astrophys.618, A158 (2018)

  43. [46]

    & Kratter, K

    Moe, M. & Kratter, K. M. Dynamical Formation of Close Binaries during the Pre-main-sequence Phase.Astrophys. J.854, 44 (2018)

  44. [47]

    P., Offner, S

    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)

  45. [48]

    P., Offner, S

    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)

  46. [49]

    Cournoyer-Cloutier, C.et al.Massive Star Cluster Formation with Binaries. I. Evolution of Binary Populations.Astrophys. J.977, 203 (2024)

  47. [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)

  48. [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

  49. [53]

    & Moe, M

    Tokovinin, A. & Moe, M. Formation of close binaries by disc fragmentation and migration, and its statistical modelling.Mon. Not. R. Astron. Soc.491, 19 5158–5171 (2020)

  50. [54]

    Heggie, D. C. Binary evolution in stellar dynamics.Mon. Not. R. Astron. Soc. 173, 729–787 (1975)

  51. [55]

    & Nordlund, ˚A

    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)

  52. [56]

    & Di Stefano, R

    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)

  53. [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

  54. [58]

    Hopkins, P. F. & Raives, M. J. Accurate, meshless methods for magnetohydro- dynamics.Mon. Not. R. Astron. Soc.455, 51–88 (2016)

  55. [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

  56. [61]

    R., Bonnell, I

    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)

  57. [62]

    R., Bonnell, I

    Bate, M. R., Bonnell, I. A. & Price, N. M. Modelling accretion in protobinary systems.Mon. Not. R. Astron. Soc.277, 362–376 (1995)

  58. [63]

    Y., Guszejnov, D., Hopkins, P

    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)

  59. [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)

  60. [65]

    Grudi´ c, M. Y. & Guszejnov, D. MakeCloud (2021). URL https://github.com/ mikegrudic/MakeCloud

  61. [66]

    C., Stone, J

    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

  62. [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)

  63. [68]

    T., Offner, S

    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)

  64. [69]

    B., Zenati, Y

    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)

  65. [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 ...