Recognition: unknown
Eccentricity as a signature of hierarchical subsolar-mass mergers in collapsar disks
Pith reviewed 2026-05-07 11:25 UTC · model grok-4.3
The pith
Hierarchical mergers in collapsar disks build eccentricity that can survive to the final merger at levels around 0.1.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Numerical relativity simulations of hierarchical compact object mergers modeled as black holes in a disk-like geometry demonstrate the build-up of potentially large eccentricity for the final merger, of order e ≃ 0.6 initially, and show that, because of the short lifetime of the system, a substantial part of this eccentricity, up to e ≃ 0.1, can survive until merger in the general case. As a result, future detections of eccentricities in potential subsolar-mass gravitational-wave candidate events would be a strong indicator for a hierarchical formation scenario.
What carries the argument
Numerical relativity simulations of repeated black hole mergers in a simplified disk-like geometry that track the evolution of orbital eccentricity through successive captures.
If this is right
- Eccentricity of order 0.1 can persist until the final merger because the collapsar system lifetime is short.
- Detection of eccentricity in a subsolar-mass gravitational-wave event would indicate a hierarchical formation channel.
- Each step in the hierarchy can also produce its own electromagnetic counterpart.
- The mechanism supplies a concrete way to interpret sub-threshold gravitational-wave candidates linked to possible superkilonova transients.
Where Pith is reading between the lines
- Current and next-generation gravitational-wave detectors could search for this eccentricity signature in low-mass events.
- Absence of measurable eccentricity in multiple subsolar-mass candidates would constrain the fraction of mergers that proceed hierarchically.
- The same disk environment could be explored with different initial fragment masses or disk thicknesses to map how eccentricity retention depends on those parameters.
Load-bearing premise
The simplified disk-like geometry used in the simulations accurately represents the repeated capture and merger dynamics inside a real collapsar accretion disk.
What would settle it
A gravitational-wave detection of a subsolar-mass merger whose waveform shows measurable eccentricity near 0.1.
Figures
read the original abstract
In this work, we investigate gravitational-wave signatures of a proposed subsolar-mass merger scenario resulting from fragmentation inside a collapsar accretion disk. This scenario has gained recent interest with the electromagnetic transient AT2025ulz, a possible superkilonova counterpart candidate to the sub-threshold gravitational wave event S250818k. One prediction of fragmentation is the formation of multiple smaller neutron-star fragments, some of which might merge hierarchically. Such mergers are expected not only to produce individual electromagnetic counterparts, but also, because of their repeated capture and merger dynamics, to impart kicks to the system and thereby drive orbital eccentricity. By performing numerical relativity simulations of hierarchical compact object mergers modeled as black holes in a disk-like geometry consistent with this scenario, we demonstrate the build-up of potentially large eccentricity for the final merger, of order $e \simeq 0.6$ initially, and show that, because of the short lifetime of the system, a substantial part of this eccentricity , up to $e\simeq 0.1$, can survive until merger in the general case. As a result, future detections of eccentricities in potential subsolar-mass gravitational-wave candidate events would be a strong indicator for a hierarchical formation scenario.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates gravitational-wave signatures of hierarchical subsolar-mass mergers arising from fragmentation in collapsar accretion disks. It argues that repeated capture and merger dynamics impart kicks that build orbital eccentricity to e ≃ 0.6 initially, with a substantial fraction (up to e ≃ 0.1) surviving to merger owing to the short system lifetime. This is demonstrated via numerical relativity simulations that model the compact objects as black holes in a simplified disk-like geometry. The resulting eccentricity is proposed as a strong indicator for the hierarchical channel in future detections of subsolar-mass events, with reference to candidates such as S250818k and the associated transient AT2025ulz.
Significance. If the modeling assumptions hold, the work supplies a concrete, falsifiable prediction that could distinguish hierarchical formation in collapsar disks from other channels for subsolar-mass compact objects. The direct numerical-relativity demonstration of eccentricity build-up and partial survival constitutes a strength, providing quantitative guidance for gravitational-wave searches and linking to recent electromagnetic candidates.
major comments (2)
- [Numerical methods and setup] The central simulations model the fragments as black holes in a simplified disk-like geometry (described in the numerical methods and setup). Real fragments are neutron stars whose equation of state, tidal deformability, and interactions with a turbulent, magnetized, accreting disk (gas drag, dynamical friction) are omitted; these effects can alter kick velocities and damp eccentricity beyond pure gravitational-wave emission, directly affecting the claim that e ≃ 0.1 survives until merger in the general case.
- [Results] The reported eccentricity values (initial order 0.6, surviving up to 0.1) are obtained from the NR runs, yet no resolution studies, convergence tests, or error estimates on the eccentricity measurements are provided. Given that initial orbital parameters and disk geometry are free parameters, it is unclear how robust the survival result is across the explored space.
minor comments (2)
- [Abstract] The abstract states that eccentricity 'can survive until merger in the general case'; a brief summary of the range of initial conditions and disk parameters actually simulated would remove ambiguity.
- [Figures] Figure captions and axis labels for the eccentricity evolution plots should explicitly note the time coordinate (e.g., time to merger) and any averaging procedure used.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed comments on our manuscript. We address each major comment point by point below, providing clarifications on our modeling choices and numerical approach while indicating the revisions incorporated to strengthen the work.
read point-by-point responses
-
Referee: [Numerical methods and setup] The central simulations model the fragments as black holes in a simplified disk-like geometry (described in the numerical methods and setup). Real fragments are neutron stars whose equation of state, tidal deformability, and interactions with a turbulent, magnetized, accreting disk (gas drag, dynamical friction) are omitted; these effects can alter kick velocities and damp eccentricity beyond pure gravitational-wave emission, directly affecting the claim that e ≃ 0.1 survives until merger in the general case.
Authors: We agree that the modeling of fragments as black holes in a simplified disk-like geometry represents a deliberate approximation, as stated in the manuscript. This choice isolates the gravitational dynamics of hierarchical capture and kicks to demonstrate the potential for eccentricity build-up to ~0.6 and partial survival to ~0.1 due to short merger timescales. Neutron-star-specific effects and disk interactions could indeed introduce additional damping or alter kicks, and we do not claim quantitative universality. In the revised manuscript, we have expanded the discussion to explicitly acknowledge these limitations, revised the phrasing from 'in the general case' to 'within this simplified model', and emphasized the need for future work with full NS physics and magnetized disks to refine predictions. revision: partial
-
Referee: [Results] The reported eccentricity values (initial order 0.6, surviving up to 0.1) are obtained from the NR runs, yet no resolution studies, convergence tests, or error estimates on the eccentricity measurements are provided. Given that initial orbital parameters and disk geometry are free parameters, it is unclear how robust the survival result is across the explored space.
Authors: We appreciate this observation on numerical validation. The original submission presented results from our primary high-resolution runs without explicit convergence details. In the revised manuscript, we have added a new subsection detailing resolution studies across multiple grid resolutions, with convergence demonstrated for the extracted eccentricity values and error estimates provided based on inter-resolution differences. We have also included a brief analysis of variations in initial orbital parameters and disk geometry within the representative range explored, confirming that the qualitative survival of eccentricity persists, although exact values show some dependence. A comprehensive parameter survey remains beyond the current scope but is noted as future work. revision: yes
Circularity Check
No significant circularity; results from direct NR simulations
full rationale
The paper's central results on eccentricity build-up (e ≃ 0.6) and survival (up to e ≃ 0.1) are obtained as direct outputs from numerical relativity simulations of black holes in a simplified disk-like geometry. Initial conditions are chosen to represent the hierarchical capture scenario, but the eccentricity values are computed quantities, not fitted parameters renamed as predictions or self-defined by construction. No load-bearing self-citations, uniqueness theorems imported from prior author work, or ansatz smuggling via citation are present. The derivation chain consists of numerical experiments whose outputs are independent of the target claim by the paper's own description.
Axiom & Free-Parameter Ledger
free parameters (1)
- initial orbital parameters and disk geometry
axioms (1)
- domain assumption Fragmentation in a collapsar accretion disk leads to multiple smaller neutron-star fragments that merge hierarchically and impart kicks driving orbital eccentricity.
Reference graph
Works this paper leans on
-
[1]
Living Reviews in Relativity , keywords =
Abbott, B. P., et al. 2016, Living Rev. Rel., 19, 1, doi: 10.1007/s41114-020-00026-9
-
[2]
2017, PhRvL, 119, 161101, doi: 10.1103/PhysRevLett.119.161101
Abbott, B. P., et al. 2017a, Phys. Rev. Lett., 119, 161101, doi: 10.1103/PhysRevLett.119.161101
-
[3]
ApJ848(2), 13 (2017) https://doi.org/10.3847/2041-8213/aa920c arXiv:1710.05834 [astro-ph.HE]
Abbott, B. P., et al. 2017b, Astrophys. J. Lett., 848, L13, doi: 10.3847/2041-8213/aa920c
-
[4]
Exploring the Sensitivity of Next Generation Gravitational Wave Detectors
Abbott, B. P., et al. 2017c, Class. Quant. Grav., 34, 044001, doi: 10.1088/1361-6382/aa51f4
-
[5]
Abbott, B. P., et al. 2020, Astrophys. J. Lett., 892, L3, doi: 10.3847/2041-8213/ab75f5
-
[6]
Abbott, R., et al. 2020, Astrophys. J. Lett., 896, L44, doi: 10.3847/2041-8213/ab960f
-
[7]
Abbott, R., et al. 2021, Astrophys. J. Lett., 915, L5, doi: 10.3847/2041-8213/ac082e
-
[8]
Agarwal, A., Siegel, D. M., Metzger, B. D., & Nagele, C. 2026, Astrophys. J., 998, 57, doi: 10.3847/1538-4357/ae2856
-
[9]
and Jegou du Laz, Theophile and Coughlin, Michael W
Ahumada, T., et al. 2026, Publ. Astron. Soc. Pac., 138, 034101, doi: 10.1088/1538-3873/ae4539
-
[10]
Gauge conditions for long-term numerical black hole evolutions without excision
Alcubierre, M., Bruegmann, B., Diener, P., et al. 2003, Phys. Rev. D, 67, 084023, doi: 10.1103/PhysRevD.67.084023
-
[11]
D., Berger, E., Fong, W., et al
Alexander, K. D., et al. 2017, Astrophys. J. Lett., 848, L21, doi: 10.3847/2041-8213/aa905d Figure 5.Same as Fig. 2, but using Newtonian point-par- ticle dynamics
-
[12]
Anand, S., et al. 2021, Nature Astron., 5, 46, doi: 10.1038/s41550-020-1183-3 10 T able 1.Initial puncture coordinates (x, y, z) and coordinate velocities (v x, vy, vz) for all three simulated configurations. The central puncture has physical mass 3M ⊙ and all other initial punctures have physical mass 0.2M ⊙, and dimensionless black hole spins are denote...
-
[13]
Andreoni, I., et al. 2020, Astrophys. J., 890, 131, doi: 10.3847/1538-4357/ab6a1b
-
[14]
Bamber, J., Shapiro, S. L., Ruiz, M., & Tsokaros, A. 2025, Phys. Rev. D, 112, 024046, doi: 10.1103/8z1j-rzmx
-
[15]
Beloborodov, A. M. 2003, Astrophys. J., 588, 931, doi: 10.1086/374217
-
[16]
, year = 1952, month = jan, volume =
Bondi, H. 1952, Mon. Not. Roy. Astron. Soc., 112, 195, doi: 10.1093/mnras/112.2.195
-
[17]
, year = 1944, month = jan, volume =
Bondi, H., & Hoyle, F. 1944, Mon. Not. Roy. Astron. Soc., 104, 273, doi: 10.1093/mnras/104.5.273
-
[18]
Boschini, M., Loutrel, N., Gerosa, D., & Fumagalli, G. 2025, Phys. Rev. D, 111, 024008, doi: 10.1103/PhysRevD.111.024008
-
[19]
Brandt, S., & Bruegmann, B. 1997, Phys. Rev. Lett., 78, 3606, doi: 10.1103/PhysRevLett.78.3606
-
[20]
2023, , 957, 68, 10.3847/1538-4357/acfc1c
Burrows, A., Vartanyan, D., & Wang, T. 2023, Astrophys. J., 957, 68, doi: 10.3847/1538-4357/acfc1c
-
[21]
doi:10.1103/PhysRevLett.96.111101 , eid =
Zlochower, Y. 2006, Phys. Rev. Lett., 96, 111101, doi: 10.1103/PhysRevLett.96.111101
-
[22]
Cantiello, M., Yoon, S. C., Langer, N., & Livio, M. 2007, Astron. Astrophys., 465, L29, doi: 10.1051/0004-6361:20077115
-
[23]
Chambers, J. E. 2001, Icarus, 152, 205, doi: 10.1006/icar.2001.6639
-
[24]
V., Dietrich, T., Johnson-McDaniel, N
Chaurasia, S. V., Dietrich, T., Johnson-McDaniel, N. K., et al. 2018, Phys. Rev. D, 98, 104005, doi: 10.1103/PhysRevD.98.104005
-
[25]
Chen, W.-X., & Beloborodov, A. M. 2007, Astrophys. J., 657, 383, doi: 10.1086/508923
-
[26]
Chen, Y.-X., & Metzger, B. D. 2025, Astrophys. J. Lett., 991, L22, doi: 10.3847/2041-8213/ae045d
-
[27]
Chiang, E., & Youdin, A. N. 2010, Ann. Rev. Earth Planet. Sci., 38, 493, doi: 10.1146/annurev-earth-040809-152513
-
[28]
Chirenti, C., Gold, R., & Miller, M. C. 2017, Astrophys. J., 837, 67, doi: 10.3847/1538-4357/aa5ebb
-
[29]
Chornock, R., et al. 2017, Astrophys. J. Lett., 848, L19, doi: 10.3847/2041-8213/aa905c
-
[30]
Colpi, M., Shapiro, S. L., & Teukolsky, S. A. 1989, ApJ, 339, 318, doi: 10.1086/167299
-
[31]
Colpi, M., Shapiro, S. L., & Teukolsky, S. A. 1991, ApJ, 369, 422, doi: 10.1086/169771 11
-
[32]
Corman, M., East, W. E., & Read, J. S. 2026, https://arxiv.org/abs/2603.25102
-
[33]
Cowperthwaite, P. S., et al. 2017, Astrophys. J. Lett., 848, L17, doi: 10.3847/2041-8213/aa8fc7
-
[34]
Cromartie, H. T., et al. 2019, Nature Astron., 4, 72, doi: 10.1038/s41550-019-0880-2
-
[35]
De, S., & Siegel, D. M. 2021, Astrophys. J., 921, 94, doi: 10.3847/1538-4357/ac110b
-
[36]
Hessels, J. 2010, Nature, 467, 1081, doi: 10.1038/nature09466
-
[37]
Deppe, N., Throwe, W., Kidder, L. E., et al. 2026,SpECTRE v2026.04.01, 2026.04.01, 10.5281/zenodo.19373346 Zenodo, doi: 10.5281/zenodo.19373346
-
[38]
Drout, M. R., et al. 2017, Science, 358, 1570, doi: 10.1126/science.aaq0049
-
[39]
East, W. E., & Pretorius, F. 2012, Astrophys. J. Lett., 760, L4, doi: 10.1088/2041-8205/760/1/L4
-
[40]
E., Pretorius, F., & Stephens, B
East, W. E., Pretorius, F., & Stephens, B. C. 2012, Phys. Rev. D, 85, 124009, doi: 10.1103/PhysRevD.85.124009
-
[41]
Edgar, R. G. 2004, New Astron. Rev., 48, 843, doi: 10.1016/j.newar.2004.06.001
-
[42]
Eichler, D., Livio, M., Piran, T., & Schramm, D. N. 1989, Nature, 340, 126, doi: 10.1038/340126a0
-
[43]
A Horizon Study for Cosmic Explorer: Science, Observatories, and Community
Evans, M., et al. 2021, https://arxiv.org/abs/2109.09882
work page internal anchor Pith review arXiv 2021
-
[44]
M., Sravan, N., Cantrell, A., et al
Farr, W. M., Sravan, N., Cantrell, A., et al. 2011, Astrophys. J., 741, 103, doi: 10.1088/0004-637X/741/2/103
-
[45]
B., Havlickova, M., & Rasio, F
Ford, E. B., Havlickova, M., & Rasio, F. A. 2001, Icarus, 150, 303, doi: 10.1006/icar.2001.6588
-
[46]
Fryer, C. L. 1999, Astrophys. J., 522, 413, doi: 10.1086/307647
-
[47]
Fryer, C. L., & Kalogera, V. 2001, Astrophys. J., 554, 548, doi: 10.1086/321359
-
[48]
Fuller, J., Piro, A. L., & Jermyn, A. S. 2019, Mon. Not. Roy. Astron. Soc., 485, 3661, doi: 10.1093/mnras/stz514
-
[49]
Gillanders, J. H., Huber, M. E., Nicholl, M., et al. 2025, ApJL, 995, L27, doi: 10.3847/2041-8213/ae2125
-
[50]
1993, Icarus, 106, 247, doi: 10.1006/icar.1993.1169
Gladman, B. 1993, Icarus, 106, 247, doi: 10.1006/icar.1993.1169
-
[51]
Pretorius, F. 2012, Phys. Rev. D, 86, 121501, doi: 10.1103/PhysRevD.86.121501
-
[52]
Goldreich, P., Lithwick, Y., & Sari, R. 2004, Ann. Rev. Astron. Astrophys., 42, 549, doi: 10.1146/annurev.astro.42.053102.134004
-
[53]
Goldreich, P., & Tremaine, S. 1980, Astrophys. J., 241, 425, doi: 10.1086/158356
-
[54]
2023, ApJL, 952, L32, doi: 10.3847/2041-8213/ace779
Gottlieb, O., Jacquemin-Ide, J., Lowell, B., Tchekhovskoy, A., & Ramirez-Ruiz, E. 2023, ApJL, 952, L32, doi: 10.3847/2041-8213/ace779
-
[55]
2022a, MNRAS, 510, 4962, doi: 10.1093/mnras/stab3784
Tchekhovskoy, A. 2022a, MNRAS, 510, 4962, doi: 10.1093/mnras/stab3784
-
[56]
Gottlieb, O., Liska, M., Tchekhovskoy, A., et al. 2022b, ApJL, 933, L9, doi: 10.3847/2041-8213/ac7530
-
[57]
J., Busmann, M., Koehn, H., et al
Hall, X. J., Busmann, M., Koehn, H., et al. 2025, arXiv e-prints, arXiv:2510.24620, doi: 10.48550/arXiv.2510.24620
-
[58]
Hallinan, G., et al. 2017, Science, 358, 1579, doi: 10.1126/science.aap9855
-
[59]
Compact binary evolutions with the Z4c formulation
Hilditch, D., Bernuzzi, S., Thierfelder, M., et al. 2013, Phys. Rev. D, 88, 084057, doi: 10.1103/PhysRevD.88.084057
-
[60]
Hinder, I., Kidder, L. E., & Pfeiffer, H. P. 2018, PhRvD, 98, 044015, doi: 10.1103/PhysRevD.98.044015
-
[61]
A., Kumar, P., Agarwal, B., et al
Huerta, E. A., Kumar, P., Agarwal, B., et al. 2017, PhRvD, 95, 024038, doi: 10.1103/PhysRevD.95.024038
-
[62]
Kaaz, N., Murguia-Berthier, A., Chatterjee, K., Liska, M. T. P., & Tchekhovskoy, A. 2023, Astrophys. J., 950, 31, doi: 10.3847/1538-4357/acc7a1
-
[63]
Kasliwal, M. M., et al. 2017, Science, 358, 1559, doi: 10.1126/science.aap9455
-
[64]
Kasliwal, M. M., et al. 2025, Astrophys. J. Lett., 995, L59, doi: 10.3847/2041-8213/ae2000
-
[65]
Kidder, L. E., et al. 2017, J. Comput. Phys., 335, 84, doi: 10.1016/j.jcp.2016.12.059
-
[66]
Kim, Y., & Most, E. R. 2025, Phys. Rev. D, 111, 083025, doi: 10.1103/PhysRevD.111.083025
-
[67]
Kley, W., & Nelson, R. P. 2012, Ann. Rev. Astron. Astrophys., 50, 211, doi: 10.1146/annurev-astro-081811-125523
-
[68]
Kokkotas, K. D., & Schaefer, G. 1995, Mon. Not. Roy. Astron. Soc., 275, 301, doi: 10.1093/mnras/275.2.301
-
[69]
2006, ApJ, 642, 1131, doi: 10.1086/501448
Kokubo, E., Kominami, J., & Ida, S. 2006, ApJ, 642, 1131, doi: 10.1086/501448
-
[70]
Kreidberg, L., Bailyn, C. D., Farr, W. M., & Kalogera, V. 2012, Astrophys. J., 757, 36, doi: 10.1088/0004-637X/757/1/36
-
[71]
Theoretical and experimental constraints for the equation of state of dense and hot matter
Kumar, R., et al. 2024, Living Rev. Rel., 27, 3, doi: 10.1007/s41114-024-00049-6
-
[72]
Mon Not Roy Astron Soc 270:611
Lai, D. 1994, Mon. Not. Roy. Astron. Soc., 270, 611, doi: 10.1093/mnras/270.3.611
-
[73]
Lattimer, J. M., & Prakash, M. 2016, Phys. Rept., 621, 127, doi: 10.1016/j.physrep.2015.12.005
-
[74]
Lattimer, J. M., & Schramm, D. N. 1974, Astrophys. J. Lett., 192, L145, doi: 10.1086/181612 12
-
[75]
Lerner, Y., Stone, N. C., & Ofengeim, D. D. 2025, Mon. Not. Roy. Astron. Soc., 545, staf1835, doi: 10.1093/mnras/staf1835
-
[76]
Lin, D. N. C., & Papaloizou, J. 1986, Astrophys. J., 309, 846, doi: 10.1086/164653
-
[77]
Lissauer, J. J. 1993, Ann. Rev. Astron. Astrophys., 31, 129, doi: 10.1146/annurev.aa.31.090193.001021
-
[78]
Ma, L., & Fuller, J. 2019, Mon. Not. Roy. Astron. Soc., 488, 4338, doi: 10.1093/mnras/stz2009
-
[79]
MacFadyen, A. I., & Woosley, S. E. 1999, ApJ, 524, 262, doi: 10.1086/307790
-
[80]
Margalit, B., & Metzger, B. D. 2017, Astrophys. J. Lett., 850, L19, doi: 10.3847/2041-8213/aa991c
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.