pith. machine review for the scientific record. sign in

arxiv: 2512.09036 · v2 · submitted 2025-12-09 · 🌌 astro-ph.HE · astro-ph.GA

Probing the origin of the kilonova candidate GRB 230307A: analysis of host galaxy and offset

Pith reviewed 2026-05-16 23:20 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GA
keywords GRB 230307Akilonovabinary neutron starhost galaxy offsetnatal kicksglobular clusterNFW haloJWST
0
0 comments X

The pith

GRB 230307A's 40 kpc offset makes a globular cluster origin unlikely and requires fine-tuned natal kicks for a disk-born binary neutron star.

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

The paper tests why the kilonova candidate tied to GRB 230307A lies so far from its host galaxy at a projected distance of about 40 kpc. Two scenarios are considered: a merger inside a remote globular cluster or a binary neutron star born in the galactic disk and flung outward by natal kicks from the neutron stars' formation. JWST imaging shows no globular cluster matches the expected properties at the site, making that channel improbable. For the disk-kick channel the authors construct a mass model of the host from JWST and MUSE data assuming an NFW halo and use Bayesian inference to calculate the travel time needed for different kick speeds and starting positions. When these travel times are compared against population-synthesis simulations of binary neutron star mergers, only 0.1 percent of the simulated systems fall in the compatible region of kick-velocity and coalescence-time space.

Core claim

A globular-cluster origin is unlikely because JWST data reveal no suitable clusters at the transient location. The host galaxy's mass distribution is modeled via Bayesian inference under an NFW halo assumption using rotation-curve and stellar-mass constraints from MUSE and JWST. This model is used to compute the time for a disk-born BNS to reach the observed offset while marginalizing over kick velocity, direction, and birth radius. Population-synthesis runs supply the distribution of natal kicks and coalescence times; the two approaches overlap in parameter space, yet only 0.1 percent of simulated mergers satisfy the galaxy-mass model, implying that a disk origin is possible only under fine

What carries the argument

Bayesian mass modeling of the host galaxy under an NFW halo profile, combined with population-synthesis simulations that supply natal kick velocities and coalescence times for binary neutron star systems.

If this is right

  • A globular-cluster origin is disfavored by the absence of matching clusters in the JWST imaging.
  • Only a small fraction of disk-born BNS systems possess the kick velocity and merger time needed to reach the observed offset.
  • Most binary neutron star mergers are therefore expected to occur at smaller projected distances from their host galaxies.
  • Distant kilonovae like this one would be rare events requiring specific fine-tuned conditions on natal kicks and coalescence times.

Where Pith is reading between the lines

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

  • If the low compatibility fraction holds, large-offset kilonovae would be exceptional rather than representative of the typical BNS population, affecting rate estimates for such events.
  • Future proper-motion measurements of the merger site could directly test whether the required kick velocities are realized in nature.
  • Similar offset analyses on additional events could constrain the high-velocity tail of neutron-star natal kicks across different galaxy types.

Load-bearing premise

The NFW halo plus Bayesian mass model correctly describes the gravitational potential at 40 kpc, and the population-synthesis code plus stellar tracks fully sample the relevant natal-kick and coalescence-time distributions.

What would settle it

A direct kinematic measurement of the host's mass or rotation curve at large radii that deviates from the NFW Bayesian model, or an updated population-synthesis run that yields a substantially larger fraction of systems reaching 40 kpc.

Figures

Figures reproduced from arXiv: 2512.09036 by Amanda E. Araujo-Carvalho, Arianna Cortesi, Charles D. Kilpatrick, Christopher J. Conselice, Clecio R. Bom, Daniel Ruschel-Dutra, Davi C. Rodrigues, Duncan Austin, Fabricio Ferrari, Giuliano Iorio, Luidhy Santana-Silva, Luis Lomeli-Nu\~nez, Nathan Adams, Roberto Cid Fernandes, Thomas Harvey.

Figure 1
Figure 1. Figure 1: A color-magnitude diagram of a sample of globular clus￾ters from Blakeslee et al. (2012) in I- and H-band compared with limits on a counterpart to GRB 230307A in F070W and F150W. The red line corresponds to our joint detection threshold to point sources in F070W and F150W, with the grey region containing de￾tected and the white region undetected point sources. We rule out nearly all of the globular cluster… view at source ↗
Figure 2
Figure 2. Figure 2: Upper: Surface brightness profile obtained extracting the galaxy light in elliptical bins, following the galaxy axis ratio, of increasing radius, shown as blue dots with error bars. The black dashed line represents the 1D fit of these data, assuming an exponential profile. The red line shows the output of the 2D GALFIT for an exponential profile, which leads to a larger disk scale length (h) than the previ… view at source ↗
Figure 3
Figure 3. Figure 3: Left: Velocity map of the Hα emission (see Sec. 3). The solid black line represents the path of virtual extractions. A clear rotation pattern can be discerned in this map. Right: The rotation curve extracted from the map on the left with error bars indicating the 1σ interval of uncertainty [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Left. Marginalized 2D and 1D sections of the posterior distribution for the fitted parameters. In the 2D regions, the contours indicate the 1σ and 2σ credible regions, and the black dot marks the mode of each marginalized case. Right. The best-fit (MAP) rotation curve with a NFW dark matter halo (δ⋆ = 1.10, ∆i = 0.01, logc = 0.75, log(r200/kpc) = 2.64). In the z = 0 plane, the gravitational potential for t… view at source ↗
Figure 5
Figure 5. Figure 5: The time interval ∆t, as a function of the radial kick velocity component KR, for the binary system to travel from an initial radius R0 to the observed kilonova location. The solid curve labeled “mode” assumes the galaxy parameters (δ⋆,∆i,logc,logr200) are given by their best-fit values and that the initial radius is the most probable one, R0 = h. The shaded regions are the 1σ and 2σ credible intervals obt… view at source ↗
read the original abstract

We investigate the host galaxy of the long gamma-ray burst GRB 230307A, which is associated with a kilonova candidate likely produced by a binary neutron-star (BNS) merger. The transient occurred at a projected offset of ~40 kpc from its host. We consider two explanations for this large distance: (i) NSs that merge inside a remote globular cluster, or (ii) a BNS that formed in the disk whose orbit was strongly modified by the NS natal kicks. Using JWST data and comparisons with known globular clusters, we show that a globular-cluster origin is unlikely. Using JWST and MUSE data, we derive the host galaxy morphology, stellar mass, estimate the atomic gas (HI+He) contribution, and the host rotation curve. Assuming an NFW halo and applying Bayesian inference, we obtain a mass model for the host. From this model, we compute the time required for a disk-formed BNS, with a given natal kick, to reach the observed offset while marginalizing over uncertainties and over the initial position in the disk. We compare these results with BNS-merger simulations from a population-synthesis code combined with stellar evolutionary tracks, which provide the coalescence time and kick velocity for each realization. The two approaches have an overlap in the kick-time diagram, but only 0.1% of the simulated systems is compatible with the galaxy-mass model. This indicates that a disk origin is possible, but requires fine-tuned conditions for the kilonova to occur at such a large distance from the host galaxy.

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 / 3 minor

Summary. The manuscript analyzes the host galaxy of GRB 230307A and its associated kilonova candidate at a projected offset of ~40 kpc. JWST imaging is used to exclude a globular-cluster origin by direct comparison with known clusters. MUSE rotation-curve, JWST stellar-mass, and gas-mass data are combined to construct a Bayesian mass model assuming an NFW halo; from this model the authors compute the natal-kick and coalescence-time combinations required for a disk-born BNS to reach the observed offset, marginalizing over initial disk position and model uncertainties. These requirements are then compared with the joint distribution of kick velocities and merger times drawn from population-synthesis simulations, yielding an overlap of only 0.1 % and the conclusion that a disk origin is possible but demands fine-tuned conditions.

Significance. If the central result holds, the work quantifies the rarity of large-offset kilonovae from disk-born BNS systems and supplies a concrete observational test for natal-kick and merger-time distributions. Strengths include the direct JWST exclusion of a globular-cluster channel and the integration of MUSE/JWST constraints into a Bayesian mass model. The explicit comparison with population-synthesis realizations provides a falsifiable metric (the 0.1 % fraction) that can be revisited with improved halo constraints or larger simulation grids.

major comments (2)
  1. [Mass modeling / Bayesian inference] Mass-model section (Bayesian NFW fit): the MUSE rotation curve primarily constrains the inner disk and bulge; the halo concentration and total mass at 40 kpc are therefore set by the NFW functional form and chosen priors rather than direct dynamical tracers. No sensitivity tests are shown for plausible variations in concentration or halo mass, yet these parameters directly shift the escape speed and orbital timescale that define the allowed kick-time region and therefore the reported 0.1 % overlap fraction.
  2. [Comparison with population-synthesis simulations] Kick-time comparison section: the 0.1 % compatibility is obtained by solving for the minimum kick velocity and coalescence time needed to reach the observed offset under the derived potential and then counting the fraction of population-synthesis realizations that fall inside that region. The manuscript does not report the grid resolution, the sampling density of the natal-kick distribution, or the completeness of the stellar-evolution tracks at the relevant masses and metallicities; these details are required to assess whether the quoted fraction is robust or an artifact of under-sampling.
minor comments (3)
  1. [Methods] The text refers to 'marginalization over uncertainties' without specifying which parameters are marginalized (e.g., halo concentration, initial disk radius, inclination) or how the marginalization is performed numerically.
  2. [Figures] Figure captions and axis labels for the kick-time diagram should explicitly state the units of kick velocity (km s^{-1}) and coalescence time (Myr or Gyr) and indicate the contour levels used to define the 0.1 % region.
  3. [Simulations] The population-synthesis code and its version, together with the adopted initial-mass-function and common-envelope parameters, should be cited with a specific reference or repository link.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments highlight important aspects of the mass modeling and simulation comparison that require clarification and additional material. We address each point below and will revise the manuscript accordingly.

read point-by-point responses
  1. Referee: [Mass modeling / Bayesian inference] Mass-model section (Bayesian NFW fit): the MUSE rotation curve primarily constrains the inner disk and bulge; the halo concentration and total mass at 40 kpc are therefore set by the NFW functional form and chosen priors rather than direct dynamical tracers. No sensitivity tests are shown for plausible variations in concentration or halo mass, yet these parameters directly shift the escape speed and orbital timescale that define the allowed kick-time region and therefore the reported 0.1 % overlap fraction.

    Authors: We agree that the outer-halo parameters are largely set by the NFW functional form and priors once the inner rotation curve is fixed. In the revised manuscript we will add a dedicated sensitivity analysis in which we vary the concentration parameter over c = 5–15 and the halo mass within ±20 % of the fiducial value, recompute the escape-speed and orbital-time contours, and report the resulting range in the overlap fraction. This will quantify how robust the 0.1 % figure is to plausible changes in the outer potential. revision: yes

  2. Referee: [Comparison with population-synthesis simulations] Kick-time comparison section: the 0.1 % compatibility is obtained by solving for the minimum kick velocity and coalescence time needed to reach the observed offset under the derived potential and then counting the fraction of population-synthesis realizations that fall inside that region. The manuscript does not report the grid resolution, the sampling density of the natal-kick distribution, or the completeness of the stellar-evolution tracks at the relevant masses and metallicities; these details are required to assess whether the quoted fraction is robust or an artifact of under-sampling.

    Authors: We will include the requested technical details in the revised methods section: the grid resolution and sampling density used for the natal-kick distribution, together with the completeness of the stellar-evolution tracks at the relevant progenitor masses and metallicities. These additions will allow readers to evaluate whether the 0.1 % overlap is limited by sampling or by the underlying physics. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper derives a host mass model via Bayesian inference applied to MUSE rotation-curve data, JWST stellar-mass estimates, and gas contributions, under an explicit NFW halo assumption. It then computes the natal-kick and coalescence-time region required for a disk-born BNS to reach the observed 40 kpc offset by integrating orbits in that potential while marginalizing over initial disk position and model uncertainties. This allowed region is compared to the joint kick-time distribution obtained from an independent population-synthesis code plus stellar-evolution tracks. The reported 0.1 % overlap is therefore a direct numerical comparison between two externally sourced distributions rather than a quantity forced by construction from the same fitted parameters. No self-definitional steps, fitted-input predictions, or load-bearing self-citations appear in the central claim.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The analysis rests on the NFW dark-matter profile, standard stellar-evolution tracks inside the population-synthesis code, and the assumption that the observed rotation curve can be extrapolated to 40 kpc without additional components.

free parameters (2)
  • NFW concentration and scale radius
    Fitted via Bayesian inference to the rotation curve and stellar mass; central to the travel-time calculation.
  • natal kick velocity distribution parameters
    Taken from the population-synthesis code; the 0.1% fraction is sensitive to the high-velocity tail.
axioms (2)
  • domain assumption The host galaxy's gravitational potential is dominated by an NFW halo at 40 kpc
    Invoked when computing travel time from the fitted mass model.
  • domain assumption The population-synthesis code plus stellar tracks produce a representative sample of BNS systems with realistic kicks and coalescence times
    Used to generate the comparison distribution.

pith-pipeline@v0.9.0 · 5667 in / 1682 out tokens · 38756 ms · 2026-05-16T23:20:27.309352+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

72 extracted references · 72 canonical work pages · 21 internal anchors

  1. [1]

    [+s e۷/ /5EJJ 7e/ ˫& v 9eXk? _ Kl, (@ P 0P (@ P h^ta f f3OD G2РGp AX ? ˼ @

    thebibliography [1] 20pt to REFERENCES 6pt =0pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key on reference command Each re...

  2. [2]

    J., Conselice, C

    Adams, N. J., Conselice, C. J., Austin, D., et al. 2024, The Astrophysical Journal , 965, 169

  3. [3]

    P., Tollerud , E

    Astropy Collaboration , Robitaille , T. P., Tollerud , E. J., et al. 2013, http://dx.doi.org/10.1051/0004-6361/201322068 , 558, A33

  4. [4]

    M., Sip o cz , B

    Astropy Collaboration , Price-Whelan , A. M., Sip o cz , B. M., et al. 2018, http://dx.doi.org/10.3847/1538-3881/aabc4f , 156, 123

  5. [5]

    2013, PSFEx: Point Spread Function Extractor , Astrophysics Source Code Library, record ascl:1301.001

    Bertin , E. 2013, PSFEx: Point Spread Function Extractor , Astrophysics Source Code Library, record ascl:1301.001

  6. [6]

    1996, , 117, 393, 10.1051/aas:1996164

    Bertin , E. & Arnouts , S. 1996, http://dx.doi.org/10.1051/aas:1996164 Astronomy and Astrophysics Supplement , 117, 393

  7. [7]

    1991, http://dx.doi.org/10.1086/170058 Astrophys

    Bertola, F., Bettoni, D., Danziger, J., et al. 1991, http://dx.doi.org/10.1086/170058 Astrophys. J. , 373, 369

  8. [8]

    & Tremaine, S

    Binney, J. & Tremaine, S. 1988, Galactic Dynamics (Princeton Series in Astrophysics) (Princeton University Press)

  9. [9]

    Optical and IR Photometry of Globular Clusters in NGC1399: Evidence for Color-Metallicity Nonlinearity

    Blakeslee , J. P., Cho , H., Peng , E. W., et al. 2012, http://dx.doi.org/10.1088/0004-637X/746/1/88 , 746, 88 , arXiv: 1201.1031

  10. [10]

    B., van Dokkum P

    Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, http://dx.doi.org/10.1086/591786 The Astrophysical Journal , 686, 1503

  11. [11]

    S., Berger , E., Stevenson , S., et al

    Broekgaarden , F. S., Berger , E., Stevenson , S., et al. 2022, http://dx.doi.org/10.1093/mnras/stac1677 , 516, 5737 , arXiv: 2112.05763

  12. [12]

    2023, http://dx.doi.org/10.1093/mnras/stad2597 MNRAS , 526, 3273 , 2208.14974

    Cappellari , M. 2023, http://dx.doi.org/10.1093/mnras/stad2597 MNRAS , 526, 3273 , 2208.14974

  13. [13]

    R., Bailes , M., & Broekgaarden , F

    Chattopadhyay , D., Stevenson , S., Hurley , J. R., Bailes , M., & Broekgaarden , F. 2021, http://dx.doi.org/10.1093/mnras/stab973 , 504, 3682 , arXiv: 2011.13503

  14. [14]

    & Gunn, J

    Conroy, C. & Gunn, J. E. 2010, Astrophysics Source Code Library , ascl

  15. [15]

    Gravitational Waves from Black Holes in Merging Ultra-Dwarf Galaxies

    Conselice , C. J., Bhatawdekar , R., Palmese , A., & Hartley , W. G. 2020, , 890, arXiv:1907.05361, arXiv: 1907.05361

  16. [16]

    J., Adams, N., Harvey, T., et al

    Conselice, C. J., Adams, N., Harvey, T., et al. 2025, The Astrophysical Journal , 983, 30

  17. [17]

    G., Bressan , A., et al

    Costa , G., Shepherd , K. G., Bressan , A., et al. 2025, http://dx.doi.org/10.1051/0004-6361/202452573 , 694, A193 , arXiv: 2501.12917

  18. [18]

    & Fermi GBM Team

    Dalessi , S. & Fermi GBM Team . 2023, GRB Coordinates Network , 33407, 1

  19. [19]

    & Flachaire, E

    Davidson, R. & Flachaire, E. 2008, Journal of Econometrics , 146, 162

  20. [20]

    & Mandel , I

    Disberg , P. & Mandel , I. 2025, http://dx.doi.org/10.48550/arXiv.2505.22102 arXiv e-prints , arXiv:2505.22102 , arXiv: 2505.22102

  21. [21]

    Morfometryka -- A New Way of Establishing Morphological Classification of Galaxies

    Ferrari , F., de Carvalho , R. R., & Trevisan , M. 2015, http://dx.doi.org/10.1088/0004-637X/814/1/55 , 814, 55 , arXiv: 1509.05430

  22. [22]

    L., Madore, B

    Freedman, W. L., Madore, B. F., Hatt, D., et al. 2019, http://dx.doi.org/10.3847/1538-4357/ab2f73 , 882, 34

  23. [23]

    Compact Remnant Mass Function: Dependence on the Explosion Mechanism and Metallicity

    Fryer , C. L., Belczynski , K., Wiktorowicz , G., et al. 2012, http://dx.doi.org/10.1088/0004-637X/749/1/91 , 749, 91 , arXiv: 1110.1726

  24. [24]

    Gaia Collaboration , Vallenari , A., Brown , A. G. A., et al. 2023, http://dx.doi.org/10.1051/0004-6361/202243940 , 674, A1 , arXiv: 2208.00211

  25. [25]

    J., Chrimes , A

    Gaspari , N., Levan , A. J., Chrimes , A. A., & Nelemans , G. 2024, http://dx.doi.org/10.1093/mnras/stad3259 , 527, 1101 , arXiv: 2310.14773

  26. [26]

    The progenitors of compact-object binaries: impact of metallicity, common envelope and natal kicks

    Giacobbo , N. & Mapelli , M. 2018, http://dx.doi.org/10.1093/mnras/sty1999 , 480, 2011 , arXiv: 1806.00001

  27. [27]

    & Mapelli , M

    Giacobbo , N. & Mapelli , M. 2020, http://dx.doi.org/10.3847/1538-4357/ab7335 , 891, 141 , arXiv: 1909.06385

  28. [28]

    J., Adams , N

    Harvey , T., Conselice , C. J., Adams , N. J., et al. 2025, http://dx.doi.org/10.3847/1538-4357/ad8c29 , 978, 89 , arXiv: 2403.03908

  29. [29]

    2005, MNRAS, 364, 1105, doi: 10.1111/j.1365-2966.2005.09655.x

    Hobbs , G., Lorimer , D. R., Lyne , A. G., & Kramer , M. 2005, http://dx.doi.org/10.1111/j.1365-2966.2005.09087.x , 360, 974 , arXiv: astro-ph/0504584

  30. [30]

    Igoshev , A. P. 2020, http://dx.doi.org/10.1093/mnras/staa958 , 494, 3663 , arXiv: 2002.01367

  31. [31]

    P., Chruslinska, M., Dorozsmai, A., & Toonen, S

    Igoshev, A. P., Chruslinska, M., Dorozsmai, A., & Toonen, S. 2021, http://dx.doi.org/10.1093/mnras/stab2734 Mon. Not. Roy. Astron. Soc. , 508, 3345 , arXiv: 2109.10362

  32. [32]

    2023, http://dx.doi.org/10.1093/mnras/stad1630 , 524, 426 , arXiv: 2211.11774

    Iorio , G., Mapelli , M., Costa , G., et al. 2023, http://dx.doi.org/10.1093/mnras/stad1630 , 524, 426 , arXiv: 2211.11774

  33. [33]

    H., Guo Y., Hertzberg J., Katz N., Mo H

    Ivanova , N., Heinke , C. O., Rasio , F. A., Belczynski , K., & Fregeau , J. M. 2008, http://dx.doi.org/10.1111/j.1365-2966.2008.13064.x , 386, 553 , arXiv: 0706.4096

  34. [34]

    Explosion Mechanisms of Core-Collapse Supernovae

    Janka , H.-T. 2012, http://dx.doi.org/10.1146/annurev-nucl-102711-094901 Annual Review of Nuclear and Particle Science , 62, 407 , arXiv: 1206.2503

  35. [35]

    The ACS Virgo Cluster Survey. XII. The Luminosity Function of Globular Clusters in Early Type Galaxies

    Jord \'a n , A., McLaughlin , D. E., C \^o t \'e , P., et al. 2007, http://dx.doi.org/10.1086/516840 , 171, 101 , arXiv: astro-ph/0702496

  36. [36]

    Kalogera, V., Kolb, U., & King, A. R. 1998, http://dx.doi.org/10.1086/306119 Astrophys. J. , 504, 967 , arXiv: astro-ph/9803288

  37. [37]

    2023, http://dx.doi.org/10.1093/mnras/stad019 , 519, 5893 , arXiv: 2209.09252

    Kapil , V., Mandel , I., Berti , E., & M \"u ller , B. 2023, http://dx.doi.org/10.1093/mnras/stad019 , 519, 5893 , arXiv: 2209.09252

  38. [38]

    D., Fong , W.-f., Blanchard , P

    Kilpatrick , C. D., Fong , W.-f., Blanchard , P. K., et al. 2022, http://dx.doi.org/10.3847/1538-4357/ac3e59 , 926, 49 , arXiv: 2109.06211

  39. [39]

    Progenitors of gravitational wave mergers: Binary evolution with the stellar grid-based code ComBinE

    Kruckow , M. U., Tauris , T. M., Langer , N., Kramer , M., & Izzard , R. G. 2018, http://dx.doi.org/10.1093/mnras/sty2190 , 481, 1908 , arXiv: 1801.05433

  40. [40]

    A Globular Cluster Luminosity Function Distance to NGC 4993 Hosting a Binary Neutron Star Merger GW170817/GRB 170817A

    Lee , M. G., Kang , J., & Im , M. 2018, http://dx.doi.org/10.3847/2041-8213/aac2e9 , 859, L6 , arXiv: 1805.01127

  41. [41]

    SPARC: Mass Models for 175 Disk Galaxies with Spitzer Photometry and Accurate Rotation Curves

    Lelli , F., McGaugh , S. S., & Schombert , J. M. 2016, http://dx.doi.org/10.3847/0004-6256/152/6/157 , 152, 157 , arXiv: 1606.09251

  42. [42]

    J., Gompertz, B

    Levan , A. J., Gompertz , B. P., Salafia , O. S., et al. 2024, http://dx.doi.org/10.1038/s41586-023-06759-1 , 626, 737 , arXiv: 2307.02098

  43. [43]

    W., C \^o t \'e , P., et al

    Lim , S., Peng , E. W., C \^o t \'e , P., et al. 2024, http://dx.doi.org/10.3847/1538-4357/ad3444 , 966, 168 , arXiv: 2403.09926

  44. [44]

    M., Krajnovi \'c , D., Epinat, B., et al

    L \'o pez, C. M., Krajnovi \'c , D., Epinat, B., et al. 2025, Astronomy & Astrophysics , 702, A146

  45. [45]

    & Ferrari , F

    Lucatelli , G. & Ferrari , F. 2019, http://dx.doi.org/10.1093/mnras/stz2154 , 489, 1161 , arXiv: 1907.10188

  46. [46]

    H., Guo Y., Hertzberg J., Katz N., Mo H

    Macci\'o, A. V., Dutton, A. A., & Bosch, F. C. v. d. 2008, http://dx.doi.org/10.1111/j.1365-2966.2008.14029.x , 391, 1940 , arXiv: 0805.1926

  47. [47]

    2020, http://dx.doi.org/10.3847/1538-4357/ab584d , 888, 76 , arXiv: 1909.01371

    Mapelli , M., Spera , M., Montanari , E., et al. 2020, http://dx.doi.org/10.3847/1538-4357/ab584d , 888, 76 , arXiv: 1909.01371

  48. [48]

    2025, http://dx.doi.org/10.48550/arXiv.2510.06311 arXiv e-prints , arXiv:2510.06311 , arXiv: 2510.06311

    Marinacci , F., Baldi , M., Iorio , G., et al. 2025, http://dx.doi.org/10.48550/arXiv.2510.06311 arXiv e-prints , arXiv:2510.06311 , arXiv: 2510.06311

  49. [49]

    C., & de Almeida, \'A

    Marra, V., Rodrigues, D. C., & de Almeida, \'A . O. 2020, , 494, 2875

  50. [50]

    Meidt, S. E. et al. 2014, http://dx.doi.org/10.1088/0004-637X/788/2/144 , 788, 144 , arXiv: 1402.5210

  51. [51]

    2010, Galaxy Formation and Evolution (Cambridge University Press)

    Mo, H., van den Bosch, F., & White, S. 2010, Galaxy Formation and Evolution (Cambridge University Press)

  52. [52]

    A Universal Density Profile from Hierarchical Clustering

    Navarro, J. F., Frenk, C. S., & White, S. D. 1997, http://dx.doi.org/10.1086/304888 , 490, 493 , arXiv: astro-ph/9611107

  53. [53]

    Y., et al., 2002, @doi [ ] 10.1086/340952 , 124, 266

    Peng , C. Y., Ho , L. C., Impey , C. D., & Rix , H.-W. 2002, http://dx.doi.org/10.1086/340952 , 124, 266 , arXiv: astro-ph/0204182

  54. [54]

    D., Sivaramakrishnan, A., Lajoie, C.-P., et al

    Perrin, M. D., Sivaramakrishnan, A., Lajoie, C.-P., et al. 2014, in Space telescopes and instrumentation 2014: optical, infrared, and millimeter wave, Vol. 9143, SPIE, 1174

  55. [55]

    & Dall'Agnol De Oliveira , B

    Ruschel-Dutra , D. & Dall'Agnol De Oliveira , B. 2020, danielrd6/ifscube v1.0

  56. [56]

    2021, http://dx.doi.org/10.1093/mnras/stab2058 , 507, 74 , arXiv: 2107.07635

    Ruschel-Dutra , D., Storchi-Bergmann , T., Schnorr-M \"u ller , A., et al. 2021, http://dx.doi.org/10.1093/mnras/stab2058 , 507, 74 , arXiv: 2107.07635

  57. [57]

    M., Rom \'a n , J., Knapen , J

    S \'a nchez-Alarc \'o n , P. M., Rom \'a n , J., Knapen , J. H., et al. 2023, http://dx.doi.org/10.1051/0004-6361/202346719 , 677, A117 , arXiv: 2307.02527

  58. [58]

    2023, http://dx.doi.org/10.1093/mnras/stad2768 , 526, 2210 , arXiv: 2305.04955

    Sgalletta , C., Iorio , G., Mapelli , M., et al. 2023, http://dx.doi.org/10.1093/mnras/stad2768 , 526, 2210 , arXiv: 2305.04955

  59. [59]

    Singer, L. P. & Price, L. R. 2016, http://dx.doi.org/10.1103/PhysRevD.93.024013 Phys. Rev. D , 93, 024013

  60. [60]

    P., Chen, H.-Y., Holz, D

    Singer, L. P., Chen, H.-Y., Holz, D. E., et al. 2016, http://dx.doi.org/10.3847/2041-8205/829/1/l15 ApJ , 829, L15

  61. [61]

    Supplement: Going the Distance: Mapping Host Galaxies of LIGO and Virgo Sources in Three Dimensions Using Local Cosmography and Targeted Follow-up

    Singer , L. P., Chen , H.-Y., Holz , D. E., et al. 2016, http://dx.doi.org/10.3847/0067-0049/226/1/10 ApJS , 226, 10 , arXiv: 1605.04242

  62. [62]

    Merging black hole binaries with the SEVN code

    Spera , M., Mapelli , M., Giacobbo , N., et al. 2019, http://dx.doi.org/10.1093/mnras/stz359 , 485, 889 , arXiv: 1809.04605

  63. [63]

    Swaters , R. A. & Balcells , M. 2002, http://dx.doi.org/10.1051/0004-6361:20020449 , 390, 863 , arXiv: astro-ph/0204526

  64. [64]

    The rotation curves shapes of late-type dwarf galaxies

    Swaters , R. A., Sancisi , R., van Albada , T. S., & van der Hulst , J. M. 2009, http://dx.doi.org/10.1051/0004-6361:200810516 , 493, 871 , arXiv: 0901.4222

  65. [65]

    Taylor , M. B. 2005, in Astronomical Society of the Pacific Conference Series, Vol. 347, Astronomical Data Analysis Software and Systems XIV, ed. P. Shopbell , M. Britton , & R. Ebert , 29

  66. [66]

    2023, tommyod/KDEpy: v1.1.8

    Tommy, White, S., Alexandra, et al. 2023, tommyod/KDEpy: v1.1.8

  67. [67]

    L., O’Connor, B., et al

    Troja , E., Fryer , C. L., O'Connor , B., et al. 2022, http://dx.doi.org/10.1038/s41586-022-05327-3 , 612, 228 , arXiv: 2209.03363

  68. [68]

    van der Kruit, P. C. & Freeman, K. C. 2011, http://dx.doi.org/10.1146/annurev-astro-083109-153241 Ann. Rev. Astron. Astrophys. , 49, 301 , arXiv: 1101.1771

  69. [69]

    2016, Monthly Notices of the Royal Astronomical Society , 463, 3409

    Vazdekis, A., Koleva, M., Ricciardelli, E., R \"o ck, B., & Falc \'o n-Barroso, J. 2016, Monthly Notices of the Royal Astronomical Society , 463, 3409

  70. [70]

    On the formation history of Galactic double neutron stars

    Vigna-G \'o mez , A., Neijssel , C. J., Stevenson , S., et al. 2018, http://dx.doi.org/10.1093/mnras/sty2463 , 481, 4009 , arXiv: 1805.07974

  71. [71]

    2024, Nature, 626, 742, doi: 10.1038/s41586-023-06979-5

    Yang , Y.-H., Troja , E., O'Connor , B., et al. 2024, http://dx.doi.org/10.1038/s41586-023-06979-5 , 626, 742 , arXiv: 2308.00638

  72. [72]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.archive archivePrefix empty "" archivePrefix ":" * if FUNCTION format.primaryClass primaryClass empty "" " [" primaryClass * "]" * if FUNCTION format.eprint eprint duplicate empty 'skip " " archiveprefix empty 'skip " " * archiveprefix * ":" * if " " * swap * " " * if FUNCTION n.dashify 't := ...