arxiv: 2604.24874 · v1 · submitted 2026-04-27 · 🌌 astro-ph.GA · astro-ph.SR

Recognition: unknown

N-body modelling of the ED-2 stream progenitor shows Gaia BH3's formation involved dynamical interactions

Cl\`audia Garcia-Diago, Daniel Mar\'in Pina, Giuliano Iorio, Marc Ard\`evol, Mark Gieles, Sara Rastello

Authors on Pith no claims yet

Pith reviewed 2026-05-08 02:08 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.SR

keywords Gaia BH3ED-2 streamN-body simulationsblack hole binariesdynamical interactionscluster evolutionstellar streams

0 comments

The pith

N-body simulations of the ED-2 stream progenitor show that Gaia BH3 formed through multiple dynamical interactions in the cluster.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tests whether the massive black hole binary Gaia BH3, located in the ED-2 stellar stream, could have formed through isolated binary evolution or required the dense environment of its parent cluster. The authors run N-body models of the progenitor cluster that include stellar evolution and compare the resulting stream structure and binary properties directly to observations. The models indicate that Gaia BH3 most probably arose as an exchange binary that experienced repeated strong encounters, making isolated or low-density formation unlikely. If correct, this means that cluster dynamics must be included when interpreting the growing sample of star-black hole binaries expected from future Gaia releases.

Core claim

The N-body simulations demonstrate that Gaia BH3 most likely formed as an exchange binary which underwent multiple strong dynamical interactions inside the progenitor cluster of the ED-2 stream, rather than evolving in isolation or a low-density cluster.

What carries the argument

Detailed N-body simulations of the progenitor cluster that incorporate single and binary stellar evolution, matched against the observed properties of the ED-2 stream and Gaia BH3.

If this is right

Cluster dynamics are required to explain the assembly of Gaia BH3 and should be considered for other star-black hole binaries in Gaia data.
Isolated binary evolution is disfavored for this system.
The progenitor cluster must have been dense enough to enable repeated strong encounters.
Future Gaia releases will need dynamical models to interpret newly discovered black hole binaries.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Other stellar streams may host similar dynamically assembled black hole binaries whose properties reflect their birth cluster rather than isolated evolution.
Population synthesis models that omit dense-cluster interactions will underpredict the number of massive black hole binaries with low-mass companions in the Milky Way halo.
Varying the initial density and binary fraction in the simulations could quantify how often such exchange binaries survive to the present day.

Load-bearing premise

The chosen initial conditions for the progenitor cluster produce a simulated stream and binary population that accurately match observations and can reliably separate isolated from dynamical formation channels.

What would settle it

Detection of a Gaia BH3 analogue in a low-density field population or failure of dense-cluster models to reproduce the ED-2 stream's observed kinematics and member properties would undermine the dynamical-formation conclusion.

Figures

Figures reproduced from arXiv: 2604.24874 by Cl\`audia Garcia-Diago, Daniel Mar\'in Pina, Giuliano Iorio, Marc Ard\`evol, Mark Gieles, Sara Rastello.

**Figure 1.** Figure 1: Positions of the stars of the ED-2 stream in Galactocentric cartesian coordinates. In black, view at source ↗

**Figure 2.** Figure 2: Probability distribution function (PDF) of the number view at source ↗

**Figure 3.** Figure 3: Positions of the stars of the ED-2 stream in Galactocentric cartesian coordinates. In black, view at source ↗

**Figure 4.** Figure 4: Probability distribution function of the number of S-BH binaries in our simulations as a function of the logarithm of the view at source ↗

**Figure 5.** Figure 5: Cumulative distribution of the properties of S-BH binaries in our simulations as a function of the logarithm of the period view at source ↗

**Figure 6.** Figure 6: Efficiency η, defined as the number of S-BH binaries NS−BH per unit of initial cluster mass, as a function of the initial cluster mass. In blue, models with a primordial binary population; in green, models without primordial binaries. The error bars represent the standard error of the mean (among multiple Nbody models), assuming NS−BH behaves like a Poissonian variable. 100 101 102 103 Number of strong… view at source ↗

**Figure 9.** Figure 9: Probability distribution function of the number of S-BH view at source ↗

**Figure 10.** Figure 10: Probability distribution function of the number of S view at source ↗

**Figure 11.** Figure 11: Probability distribution function of the number of S-BH view at source ↗

read the original abstract

Context. The Gaia collaboration announced the discovery of a binary of a massive black hole (33 M$_\odot$) with a low-mass giant star (Gaia BH3) in the ED-2 stellar stream. The properties of this binary, as well as its position in the stream, challenge a formation scenario invoking only isolated binary evolution. Aims. We aim to quantify the importance of dynamics in the formation of Gaia BH3 in the progenitor cluster of the ED-2 stream. Methods. We perform detailed N-body simulations of the progenitor cluster of the ED-2 stream, including the effects of single and binary stellar evolution. We compare these simulations to observations of the ED-2 stream and the properties of Gaia BH3. Results. We determine that Gaia BH3 most likely formed as an exchange binary which underwent multiple strong dynamical interactions. We highlight the importance of cluster dynamics in assembling Gaia BH3, and disfavour a formation scenario where it evolved in isolation and/or in a low-density cluster. Conclusions. The role of dynamics should be considered when interpreting properties of star-black hole binaries found in the next Gaia Data Release.

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.

Desk Editor's Note private letter to a colleague

N-body modeling makes a case for dynamical formation of Gaia BH3 but the initial conditions leave room for doubt on how strongly isolated evolution is ruled out.

read the letter

The main thing to know is that this paper uses N-body simulations to show Gaia BH3 likely formed through dynamical interactions in the ED-2 cluster progenitor, rather than evolving in isolation. They take standard N-body code with stellar evolution and run models of the progenitor to see what produces a stream like ED-2 with a binary like Gaia BH3 inside it. The result is a formation channel involving exchange binaries and multiple encounters. This is a solid step in applying these tools to a real observed system and highlighting dynamics over pure isolated evolution. One soft spot is the initial conditions. The stress-test note is right to flag that the chosen cluster masses, densities, and structures need to allow fair tests of both channels. If the setups were picked mainly to enable dense interactions, the isolated runs might not represent viable alternatives that could still match the observations. The abstract gives no quantitative metrics or details on avoiding tuning, so the disfavoring of isolated scenarios feels preliminary. This is for people in the subfield of stellar dynamics and compact object formation. A reader working on Gaia streams or black hole binaries will find the specific conclusion useful as a data point. The thinking is clear and it engages the literature on formation channels, so it should go to peer review.

Referee Report

2 major / 2 minor

Summary. The paper performs N-body simulations of the progenitor cluster of the ED-2 stellar stream, incorporating single and binary stellar evolution, and compares the results to observations of the stream and the Gaia BH3 binary (33 M⊙ black hole + low-mass giant). It concludes that Gaia BH3 most likely formed as an exchange binary undergoing multiple strong dynamical interactions, highlighting the role of cluster dynamics and disfavoring isolated or low-density evolution scenarios.

Significance. If the central claim holds after addressing the noted issues, the work would demonstrate the necessity of including dynamical interactions when modeling black hole binaries in stellar streams, with direct implications for interpreting future Gaia discoveries of star-BH systems. The forward-modeling approach using detailed N-body simulations with stellar evolution is a methodological strength that allows testing of formation channels.

major comments (2)

[Abstract and Methods] Abstract and Methods: The central claim that isolated/low-density formation is disfavoured requires demonstrating that the chosen initial conditions (mass, density, binary fraction) both reproduce the observed ED-2 stream properties and permit a genuine comparison of channels. No evidence is provided that a broader grid of lower-density initial conditions was tested to confirm that isolated formation cannot match the data under any plausible progenitor setup.
[Results] Results: The comparison of simulations to ED-2 observations and Gaia BH3 properties is described qualitatively but lacks quantitative metrics (e.g., kinematic match statistics, metallicity agreement with error bars, or goodness-of-fit values), making it difficult to assess how uniquely the dynamical channel is required or how post-hoc adjustments to initial conditions were avoided.

minor comments (2)

[Abstract] The abstract would be strengthened by briefly stating the N-body code and stellar evolution package employed.
[Introduction] Notation for binary parameters (e.g., masses, separations) should be defined consistently when first introduced.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments identify areas where the manuscript can be clarified and strengthened. We respond point-by-point to the major comments, indicating where revisions will be made.

read point-by-point responses

Referee: [Abstract and Methods] Abstract and Methods: The central claim that isolated/low-density formation is disfavoured requires demonstrating that the chosen initial conditions (mass, density, binary fraction) both reproduce the observed ED-2 stream properties and permit a genuine comparison of channels. No evidence is provided that a broader grid of lower-density initial conditions was tested to confirm that isolated formation cannot match the data under any plausible progenitor setup.

Authors: The initial conditions were selected to reproduce the observed mass, spatial extent, velocity dispersion, and metallicity of the ED-2 stream, consistent with independent constraints from Gaia data and prior stream modeling. Lower-density configurations disperse too rapidly to match these stream properties or to retain a massive black hole binary in the observed location and configuration. While a full grid of lower-density models was not performed owing to the substantial computational cost of N-body simulations with detailed stellar evolution, we will revise the Methods section to include an explicit justification of the chosen parameters, a discussion of why lower-density setups are inconsistent with the stream observations, and a statement on the limited but representative parameter space explored. This will support the comparison of formation channels without overstating the scope of the grid. revision: partial
Referee: [Results] Results: The comparison of simulations to ED-2 observations and Gaia BH3 properties is described qualitatively but lacks quantitative metrics (e.g., kinematic match statistics, metallicity agreement with error bars, or goodness-of-fit values), making it difficult to assess how uniquely the dynamical channel is required or how post-hoc adjustments to initial conditions were avoided.

Authors: We agree that quantitative metrics would improve the assessment of the results. In the revised manuscript we will add explicit quantitative comparisons, including root-mean-square residuals for positions and velocities, a Kolmogorov-Smirnov statistic for the kinematic distributions, and a table of simulated versus observed Gaia BH3 properties (component masses, orbital period, eccentricity) with observational uncertainties. We will also clarify in the Methods that the initial conditions were fixed in advance using independent literature values for the ED-2 progenitor and were not adjusted after inspecting the Gaia BH3 outcomes. These additions will better demonstrate the preference for the dynamical channel. revision: yes

Circularity Check

0 steps flagged

N-body forward modeling of ED-2 progenitor is independent of the Gaia BH3 formation conclusion

full rationale

The paper runs N-body simulations of a progenitor cluster whose initial conditions are selected to reproduce the observed ED-2 stream kinematics and metallicity; the binary formation channel (exchange vs isolated) then emerges as an output of the stellar evolution plus dynamical encounters. No equation or result is shown to be identical to an input parameter by construction, no fitted quantity is relabeled as a prediction, and no self-citation chain is invoked to justify uniqueness of the dense-cluster channel. The distinction between formation scenarios is therefore a genuine simulation outcome compared against external Gaia data rather than a tautology.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The modeling rests on standard assumptions about cluster origins and simulation physics plus free parameters for initial cluster conditions; no new entities are postulated.

free parameters (1)

Initial cluster mass, density, and structural parameters
Required to set up the N-body runs so that the simulated stream matches the observed ED-2 stream and Gaia BH3 properties.

axioms (2)

domain assumption The ED-2 stream originates from a single dissolved progenitor cluster
Invoked by the choice to model the stream as the remnant of one N-body cluster.
domain assumption N-body integration plus single/binary stellar evolution codes accurately capture the relevant dynamical and evolutionary processes
Standard assumption underlying all such simulations; stated in the methods description.

pith-pipeline@v0.9.0 · 5535 in / 1437 out tokens · 59525 ms · 2026-05-08T02:08:44.764379+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

40 extracted references · 3 canonical work pages · 2 internal anchors

[1]

Arca Sedda, M., Kamlah, A. W. H., Spurzem, R., et al. 2024, MNRAS, 528, 5119 Article number, page 9 of 11 A&A proofs:manuscript no. main

2024
[2]

2024, A&A, 687, L3

Balbinot, E., Dodd, E., Matsuno, T., et al. 2024, A&A, 687, L3

2024
[3]

2023, A&A, 678, A115

Balbinot, E., Helmi, A., Callingham, T., et al. 2023, A&A, 678, A115

2023
[4]

& Hut, P

Barnes, J. & Hut, P. 1986, Nature, 324, 446

1986
[5]

E., Bulik, T., & O’Shaughnessy, R

Belczynski, K., Holz, D. E., Bulik, T., & O’Shaughnessy, R. 2016, Nature, 534, 512

2016
[6]

& Hut, P

Casertano, S. & Hut, P. 1985, ApJ, 298, 80

1985
[7]

M., Casares, J., Muñoz-Darias, T., et al

Corral-Santana, J. M., Casares, J., Muñoz-Darias, T., et al. 2016, A&A, 587, A61 Di Carlo, U. N., Agrawal, P., Rodriguez, C. L., & Breivik, K. 2024, ApJ, 965, 22

2016
[8]

& Mayor, M

Duquennoy, A. & Mayor, M. 1991, A&A, 248, 485

1991
[9]

2024, The Open Journal of Astrophysics, 7, 38

El-Badry, K. 2024, The Open Journal of Astrophysics, 7, 38

2024
[10]

& Burdge, K

El-Badry, K. & Burdge, K. B. 2022, MNRAS, 511, 24

2022
[11]

Flower, P. J. 1996, ApJ, 469, 355

1996
[12]

L., Belczynski, K., Wiktorowicz, G., et al

Fryer, C. L., Belczynski, K., Wiktorowicz, G., et al. 2012, ApJ, 749, 91 Gaia Collaboration, Panuzzo, P., Mazeh, T., et al. 2024, A&A, 686, L2

2012
[13]

& Gnedin, O

Gieles, M. & Gnedin, O. Y . 2023, MNRAS, 522, 5340

2023
[14]

2018, MNRAS, 475, L15

Giesers, B., Dreizler, S., Husser, T.-O., et al. 2018, MNRAS, 475, L15

2018
[15]

2019, A&A, 632, A3

Giesers, B., Kamann, S., Dreizler, S., et al. 2019, A&A, 632, A3

2019
[16]

2024, A&A, 690, A144

Iorio, G., Torniamenti, S., Mapelli, M., et al. 2024, A&A, 690, A144

2024
[17]

M., et al

Jordi, C., Gebran, M., Carrasco, J. M., et al. 2010, A&A, 523, A48

2010
[18]

2024 [arXiv:2403.13579]

Kotko, I., Banerjee, S., & Belczynski, K. 2024 [arXiv:2403.13579]

work page arXiv 2024
[19]

2001, MNRAS, 322, 231 Küpper, A

Kroupa, P. 2001, MNRAS, 322, 231 Küpper, A. H. W., Maschberger, T., Kroupa, P., & Baumgardt, H. 2011, MNRAS, 417, 2300

2001
[20]

Y ., Lu, J

Lam, C. Y ., Lu, J. R., Udalski, A., et al. 2022, ApJ, 933, L23

2022
[21]

2026, A&A, 706, A105

Li, Z., Lu, X., Lü, G., et al. 2026, A&A, 706, A105

2026
[22]

2022, A&A, 664, A159

Mahy, L., Sana, H., Shenar, T., et al. 2022, A&A, 664, A159

2022
[23]

The role of accretion efficiency, natal kicks, and angular momentum transport in the formation of the Gaia black holes

Mapelli, M., Sgalletta, C., Müller-Horn, J., et al. 2026, arXiv e-prints, arXiv:2604.12839 Marín Pina, D., Gieles, M., Rastello, S., & Iorio, G. 2025, arXiv e-prints, arXiv:2510.07021 Marín Pina, D., Rastello, S., Gieles, M., et al. 2024, A&A, 688, L2

work page internal anchor Pith review Pith/arXiv arXiv 2026
[24]

J., & Hobbs, D

Mikkola, D., McMillan, P. J., & Hobbs, D. 2023, MNRAS, 519, 1989 Peñarrubia, J. 2021, MNRAS, 501, 3670

2023
[25]

B., Li, Z., Lombardi, Jr., J

Perets, H. B., Li, Z., Lombardi, Jr., J. C., & Milcarek, Jr., S. R. 2016, ApJ, 823, 113

2016
[26]

Peters, P. C. 1964, Physical Review, 136, 1224

1964
[27]

Plummer, H. C. 1911, MNRAS, 71, 460

1911
[28]

2026, A&A, 707, A217

Rastello, S., Iorio, G., Gieles, M., & Wang, L. 2026, A&A, 707, A217

2026
[29]

2023, MNRAS, 526, 740

Rastello, S., Iorio, G., Mapelli, M., et al. 2023, MNRAS, 526, 740

2023
[30]

E., de Koter, A., et al

Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444

2012
[31]

& Evans, C

Sana, H. & Evans, C. J. 2011, in IAU Symposium, V ol. 272, Active OB Stars:

2011
[32]

G., et al

Saracino, S., Kamann, S., Guarcello, M. G., et al. 2022, MNRAS, 511, 2914

2022
[33]

2023, MNRAS, 521, 3162

Saracino, S., Shenar, T., Kamann, S., et al. 2023, MNRAS, 521, 3162

2023
[34]

2022, Nature Astronomy, 6, 1085

Shenar, T., Sana, H., Mahy, L., et al. 2022, Nature Astronomy, 6, 1085

2022
[35]

Tanikawa, A., Cary, S., Shikauchi, M., Wang, L., & Fujii, M. S. 2024, MNRAS, 527, 4031

2024
[36]

S., et al

Tanikawa, A., Wang, L., Fujii, M. S., et al. 2025, The Open Journal of Astro- physics, 8, 79

2025
[37]

2020, MNRAS, 495, 4170

Tanikawa, A., Yoshida, T., Kinugawa, T., Takahashi, K., & Umeda, H. 2020, MNRAS, 495, 4170

2020
[38]

GWTC-4.0: Updating the Gravitational-Wave Transient Catalog with Observations from the First Part of the Fourth LIGO-Virgo-KAGRA Observing Run

Tanikawa, A., Yoshida, T., Kinugawa, T., et al. 2022, ApJ, 926, 83 The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Col- laboration, et al. 2025, arXiv e-prints, arXiv:2508.18082

work page internal anchor Pith review arXiv 2022
[39]

2010, AJ, 140, 1158

Torres, G. 2010, AJ, 140, 1158

2010
[40]

Wang, S., Zhao, X., Feng, F., et al. 2024, Nature Astronomy, 8, 1583 Article number, page 10 of 11 Daniel Marín Pina et al.:N-body modelling of the ED-2 progenitor showsGaiaBH3’s formation involved dynamics Appendix A: Dissolution time To find a combination of initial mass and half-mass density, (M0, ρh,0), that results in a given dissolution time,t dis, ...

2024