{"id":"45f6034b-3aff-48a2-a521-c4b5a8898b1d","arxiv_id":"2506.18151","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"By matching LIGO-Virgo-KAGRA merger rates to gamma-ray burst rates, the authors infer that the dividing mass between short and long bursts is about 1.36 times the maximum neutron star mass.","lead":"This paper compares gravitational-wave merger rates with gamma-ray burst rates to infer the mass threshold separating short and long bursts from neutron star mergers. The inferred threshold, about 1.36 times the maximum neutron star mass, suggests heavy neutron stars can survive after merging and offers a new way to probe their interiors.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inferred k2=1.36 rests on a sharp total-mass cutoff between sGRB and lGRB progenitors; if the real classification is smooth or depends on mass ratio and spin, the quoted posterior is systematically biased.","rationale":"The reader's weakest assumption identifies the same core issue: the comparison of GW and GRB rates assumes a direct, sharp mapping from total binary mass to GRB class. My stress-test sharpens this into a concrete model misspecification: the step-function classification is not varied or marginalized over, so the quoted credible interval on k2 is conditional on an untested functional form. The paper's robustness checks are genuine and demonstrate insensitivity to the rate priors, q boundaries, and removal of M_VL, but they do not probe the sharpness of the mass-duration transition or allow for overlapping populations. This matters because the central scientific claim is not just about the data but about the physical interpretation: a high k2 implies heavy neutron stars survive long after merger. If the sGRB/lGRB division is actually governed by a smooth remnant-lifetime function (or by disk mass that depends on mass ratio and spin), then the inferred k2 is a convolution of that function with the mass distribution, not a direct measure of the survival threshold. The paper is transparent about this dependence in the Discussion, so it is not internally inconsistent; it is a limitation that the current evidence cannot resolve. Because the reader's verdict is already CONDITIONAL, my analysis does not move it. I agree with the conditional recommendation: the method is promising, but the headline constraint on k2 should be reframed as a sharp-threshold-model result until the mass-duration mapping is tested against smooth or multi-dimensional alternatives.","tokens_in":12409,"tokens_out":5667,"duration_ms":60787,"concrete_test":"Replace the sharp mass cuts in Eq. (2)-(3) with a smooth classification, e.g., P(sGRB|Mtot) = 1 / (1 + exp(-(M_LS - Mtot)/delta)) and P(lGRB|Mtot) = 1 - P(sGRB|Mtot) for Mtot < M_PC, with delta an unknown nuisance parameter (e.g., log-uniform from 0.02 to 0.5 Msun), and re-run the Bayesian inference. If the posterior median or 68% interval of k2 shifts by more than ~0.05 (compared to the quoted uncertainty of 0.08-0.09), the sharp-threshold assumption is load-bearing. Alternatively, calibrate P(sGRB|Mtot) from numerical relativity remnant-lifetime simulations over a grid of total masses, mass ratios, and equations of state, then repeat the inference and compare the resulting k2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The likelihood in Eq. (2)-(3) maps the observed sGRB and lGRB rates onto disjoint intervals in total binary mass: sGRBs from M_VL <= Mtot < M_LS, lGRBs from M_LS <= Mtot < M_PC, plus prompt-collapse mergers with mass ratios in [1.2,3]. This is a sharp step-function classification. The paper varies the boundaries q_low and q_high (Table I), but never varies the sharpness of the mass-duration transition, nor does it include a parameter for the fraction of mergers that produce the opposite burst class. The cited simulations (Izquierdo et al. 2025) indicate that BH-disk jets are always long-duration, with luminosity scaling with disk mass; the observed lGRB population is therefore a threshold in disk mass, which depends on total mass, mass ratio, and equation of state in a non-stepwise way. If the true classification probability P(sGRB|Mtot) is a smooth function with width comparable to the mass distribution scale, or if a non-negligible fraction of mergers below M_LS produce lGRBs (or vice versa), the likelihood in Eq. (3) is misspecified. The quoted 68% interval on k2 does not include this systematic uncertainty. The Discussion explicitly acknowledges that the results rely on the assumption that merger-induced lGRBs and sGRBs are mainly divided by total mass. Since the headline claim (k2 > 1.24 with 90% confidence) is essentially a statement about the location of this assumed boundary, this is the load-bearing assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a Bayesian framework that combines LVK merger-rate posterior samples with observed rates of short gamma-ray bursts and kilonova-associated long gamma-ray bursts to infer the ratios (k1, k2, k3) of characteristic binary total masses (very long-lived, long-lived/short-lived, and prompt-collapse) to the TOV mass M_TOV. The fiducial analysis yields k2 = 1.36^{+0.08}_{-0.09} (median and 68% confidence interval), with a 90% lower bound of 1.24, which the authors interpret as evidence that neutron stars substantially above the non-rotating maximum mass can survive for an extended period after merger. The inference is reported to be robust to variations in the assumed GRB rates and mass-ratio cuts.","tokens_in":12752,"tokens_out":9014,"duration_ms":86270,"significance":"If the assumed sharp mass-duration mapping is correct, the method provides a novel and independent constraint on neutron-star nuclear properties and identifies correlations among M_LS, M_PC, and M_TOV. The analysis is commendably anchored to external LVK rate posterior samples and public equation-of-state constraints, and it includes explicit robustness tests over rate priors and mass-ratio thresholds. The central quantitative claim, however, is entirely conditional on the assumed step-function classification of GRB duration by total mass; the paper acknowledges this reliance in the Discussion. The significance of the claimed constraint is therefore moderate until that assumption is systematically tested.","major_comments":[{"comment":"The likelihood models sGRB and lGRB production as a sharp partition of total mass into intervals (M_VL, M_LS] and (M_LS, M_PC]. There is no parameter for a smooth transition or for a fraction of mergers producing the 'wrong' burst class. Because k2 is defined as the location of the boundary, the quoted posterior on k2 is conditional on this step-function form. The Discussion states that the results rely on this assumption. I request either a smooth classification model P(sGRB|Mtot) with a width parameter, or a demonstration that the inference is invariant to such a width; without this, the 68% interval understates the systematic uncertainty.","section":"Method, Eq. (3)"},{"comment":"The normal likelihood for eta with mean 1 and sigma 1 is an ad hoc choice; eta is a positive quantity bounded by observations in [0.005,1], so the normal assigns roughly 16% probability to unphysical negative values and centers the likelihood at the upper edge of the observed range. The authors should justify the fiducial (eta, sigma_eta) = (1,1) or adopt a distribution with the correct support (e.g., lognormal) and confirm that the posterior is robust to this choice.","section":"Method, Eq. (3)"},{"comment":"The comparison equates the observed GRB rate ratio with the ratio of GW merger rates in the two mass ranges, with no explicit treatment of beaming or detection selection. If the literature rates are not already beaming-corrected, or if the beaming corrections differ between sGRBs and lGRBs, eta is not equal to the merger-rate ratio. Please state clearly whether the adopted rates are intrinsic (beaming-corrected) and quantify the impact of beaming uncertainties on k2.","section":"Method and Discussion"}],"minor_comments":[{"comment":"The fiducial row reads \"R_sGRB = 100100-100 Gpc^-1 yr^-1\", which appears to be a formatting error for 100^{+100}_{-100} Gpc^-3 yr^-1, and the units should be Gpc^-3 yr^-1, not Gpc^-1 yr^-1.","section":"Table I"},{"comment":"In the bullet list, \"t >> 1 s\" appears twice; please use the much-greater symbol (≫) for consistency with the rest of the text.","section":"Method, bullet list"},{"comment":"The sentence \"The variations in ηGW significantly surpass the observational uncertainty of GRBs with ≲10% fluctuations\" is ambiguous about whether the 10% changes apply to the parameters or to the uncertainty; please rephrase.","section":"Results, Figure 1"},{"comment":"The section title \"Agnostic to model\" is unidiomatic; consider \"Model-agnostic analysis\" or \"Robustness to model assumptions\".","section":"Section heading"},{"comment":"The y-axis label \"Long GRB Rate/Short GRB rate\" should be \"R_lGRB/R_sGRB\" for consistency with the text and equations.","section":"Figure 1"},{"comment":"The integral signs in Eq. (2) lack explicit limits; since the variables are rates, the integration domain should be specified.","section":"Method, Eq. (2)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the method is promising, but the central claim rests on the sharp mass-duration cut. The main risk is over-claiming robustness given the step-function classification; the authors should be encouraged to run a smooth-transition variant and to justify the normal likelihood for the rate ratio. No concerns about novelty disclosure or citation practice."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper compares LVK merger rates with sGRB and kilonova-associated lGRB rates to infer k2, the total binary mass dividing short from long GRBs, relative to M_TOV. The result k2 = 1.36^{+0.08}_{-0.09} is new, and the approach is a real step up from Perna et al. 2025: it uses the LVK observed mass distribution rather than the Galactic one, includes NSBH mergers, and fits several characteristic masses in one Bayesian framework. The robustness runs across rate priors and mass-ratio cuts are reassuring; k2 barely moves. The anti-correlation between k2 and M_TOV is a useful mapping for future joint constraints.\n\nThe soft spots are real but not fatal. The likelihood in Eq. (3) is a normal ansatz with wide, arbitrary sigmas; it stands in for beaming, selection, and rate systematics. That is acceptable for a first look, but the bigger issue is the step-function classification. The model puts sGRBs in one total-mass interval and lGRBs in the adjacent one; the posterior on k2 is essentially a statement about where that assumed boundary sits. The stress-test concern is valid: Izquierdo et al. (cited in the paper) find BH-disk jets are long-duration with luminosity scaling with disk mass, so the real sGRB/lGRB division is likely a smooth threshold in disk mass, not a sharp switch at one total mass. The paper acknowledges this in the Discussion, but the quoted 68% interval does not include it. Also, the M_TOV posterior basically mirrors the Legred et al. prior, so the new content is the correlation mapping, not a fresh TOV measurement. The abstract's \"robust against substantial ... model assumptions\" is overstated; it is robust within a specific class of models.\n\nWho is this for? Anyone working on NS equation of state, compact-object merger rates, or GRB central engines. It will get cited as a method paper. I would send it to peer review: it is clear, honest, and the core idea is worth having on the record, even though the headline number should be read as conditional on the assumed mass-duration mapping. A good referee will ask for a parametrized smooth transition or a reframing as a method demonstration.\n\nRecommendation: engage with it. Worth a careful read and a serious referee.","headline":"A genuinely new rate-based probe of the NS maximum mass, but the headline k2 rests on an untested step-function mapping between total mass and GRB duration.","tokens_in":13255,"tokens_out":2946,"would_cite":true,"duration_ms":29189,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"By comparing gravitational-wave and gamma-ray burst rates, this paper identifies a dividing total binary mass of 1.36 times the neutron-star TOV mass, implying heavy remnants survive long after merger.","keywords":["neutron star mergers","gamma-ray bursts","gravitational waves","TOV mass","equation of state","kilonova","merger remnant","GRB duration"],"falsifier":"A concrete falsifier is a gravitational-wave binary neutron star merger with total mass above the inferred $k_2 M_{\\mathrm{TOV}}$ threshold (about $3.0\\,M_\\odot$ if $M_{\\mathrm{TOV}}\\simeq 2.2\\,M_\\odot$) that still produces a short gamma-ray burst, or a kilonova-associated long burst from a merger well below that threshold. A statistical version would be a future joint sample in which burst duration varies continuously with mass ratio or spin rather than switching at a single total mass.","tokens_in":12193,"feed_emoji":"💥","tokens_out":15391,"duration_ms":124658,"temperature":0.7,"pith_summary":"The paper sets out to explain what determines whether a neutron-star merger produces a short or a long gamma-ray burst, proposing that total binary mass is the controlling factor. By comparing merger rates measured by gravitational-wave detectors with the observed rates of kilonova-associated long and short bursts, it finds the dividing total mass is $1.36^{+0.08}_{-0.09}$ times the maximum mass of a nonrotating neutron star (the TOV mass), with $k_2 > 1.24$ at 90% confidence. If correct, this gives a new, observation-based route to the neutron-star equation of state, independent of pulsar timing and radius measurements, and implies that neutron stars noticeably heavier than the TOV limit can survive for hundreds of milliseconds after merger. The inference is presented as insensitive to the large uncertainties in the gamma-ray burst rates and to the choice of mass-ratio cuts at the boundaries.","feed_headline":"Neutron-star merger mass of 1.36 TOV divides long and short GRBs","feed_subtitle":"The boundary is 1.36 times the maximum nonrotating neutron-star mass, so heavy remnants can survive after merger.","key_machinery":"The machinery is a ladder of three characteristic total masses, each expressed as a multiple of the TOV mass: $M_{\\mathrm{VL}}=k_1 M_{\\mathrm{TOV}}$ for the very-long-lived/long-lived boundary, $M_{\\mathrm{LS}}=k_2 M_{\\mathrm{TOV}}$ for the long-lived/short-lived boundary, and $M_{\\mathrm{PC}}=k_3 M_{\\mathrm{TOV}}$ for prompt collapse, with $k_3>k_2>k_1>1$. The paper assigns gravitational-wave-detected mergers to short-burst or long-burst categories by placing their total masses on this ladder and compares the resulting rate ratio to observed gamma-ray burst rates using a Bayesian likelihood whose two factors are normal distributions centered on the observed short-burst rate and the long/short rate ratio. The priors for $M_{\\mathrm{TOV}}$ and the neutron-star radius come from existing pulsar, gravitational-wave, and X-ray constraints, and $M_{\\mathrm{PC}}$ is tied to $M_{\\mathrm{TOV}}$ and radius through a semi-analytic fitting formula from numerical relativity. The key work of the ladder is to make the mass-to-duration mapping explicit, so that a small shift in any characteristic mass changes the predicted long/short burst ratio by much more than current observational uncertainty.","core_discovery":"The paper's central discovery is that the total binary mass separating short-burst and long-burst merger remnants is $k_2 M_{\\mathrm{TOV}}$ with $k_2 = 1.36^{+0.08}_{-0.09}$ (median and 68% confidence interval), and $k_2 > 1.24$ at 90% confidence, where $M_{\\mathrm{TOV}}$ is the maximum mass of a cold, nonrotating neutron star. Because $k_2$ sets the boundary between a remnant that lives on the hundreds-of-milliseconds timescale, powering a short burst, and one that collapses earlier and leaves a black-hole disk, powering a long burst, a large value means mergers with total mass well above the TOV mass do not collapse promptly. The same inference returns $k_1 = 1.11^{+0.11}_{-0.08}$ for the very-long-lived/long-lived boundary, $k_3 = 1.43^{+0.07}_{-0.08}$ for the prompt-collapse boundary, and $M_{\\mathrm{TOV}} = 2.21^{+0.20}_{-0.13} M_\\odot$, together with a strong correlation between $k_2$ and $k_3$ and an anti-correlation between $k_2$ and $M_{\\mathrm{TOV}}$.","pith_inferences":["The paper does not draw this conclusion, but if $k_2$ is robust the duration distribution of merger-induced bursts should show a sharp transition at a fixed total mass; future gravitational-wave detectors could test this by stacking events with measured total mass and electromagnetic counterparts.","A testable extension suggested by the anti-correlation between $k_2$ and $M_{\\mathrm{TOV}}$ is that an independent TOV-mass measurement from radio pulsars or nuclear theory could be inserted into this relation to predict where the short/long burst divide sits before more multi-messenger events arrive.","If the total-mass assumption is ever relaxed, the same rate-comparison method could be reversed to constrain the mass-ratio or spin distributions of merging neutron-star binaries instead of nuclear parameters."],"forward_implications":["If $k_2$ is really 1.36, neutron-star remnants with total masses up to about $1.36\\,M_{\\mathrm{TOV}}$ survive for hundreds of milliseconds before collapse, setting the timescale for short-burst emission and for ejecta that enriches the disk before the black hole forms.","The inferred anti-correlation between $k_2$ and $M_{\\mathrm{TOV}}$ means a sharper determination of either quantity sharpens the other, so an independent measurement of the TOV mass would translate directly into a constraint on the long-lived/short-lived mass boundary and the prompt-collapse mass.","A single joint gravitational-wave and gamma-ray burst detection of a neutron-star merger with measured total mass and burst duration would anchor $M_{\\mathrm{LS}}$ and, through the correlations, constrain $M_{\\mathrm{TOV}}$ and $M_{\\mathrm{PC}}$.","The result supports the picture in which kilonova-associated long bursts come from heavier mergers whose remnant collapses to a black hole surrounded by a massive disk, while short bursts come from lighter mergers with longer-lived neutron-star remnants."],"supporting_citations":[{"why":"Supplies the gravitational-wave population catalog from which the paper takes merger-rate posteriors as a function of component masses.","marker":"[29]"},{"why":"Adds the local neutron star-black hole merger rate from the GW230529 event, extending the comparison to neutron star-black hole mergers.","marker":"[30]"},{"why":"Numerical simulations showing that black-hole-powered jets with massive disks produce long-duration gamma-ray bursts, fixing one side of the duration mapping.","marker":"[12]"},{"why":"Argues from observed kilonova colors that neutron stars are the central engines of standard short gamma-ray bursts.","marker":"[11]"},{"why":"Parameter study of black-hole disks showing that even light disks produce long-duration jets, so short bursts require a neutron-star engine.","marker":"[24]"},{"why":"Numerical relativity result placing the long-lived/short-lived remnant transition near 1.3-1.4 TOV masses, used as the comparison point for the inferred $k_2$.","marker":"[36]"},{"why":"Provides the semi-analytic fitting formula connecting the prompt-collapse mass $M_{\\mathrm{PC}}$ to the TOV mass and radius.","marker":"[45]"},{"why":"Supplies the public posterior samples for neutron-star radius and TOV mass used as the prior for $M_{\\mathrm{TOV}}$ and $R_{\\mathrm{TOV}}$.","marker":"[47]"},{"why":"One of the kilonova-associated long gamma-ray burst observations constraining the long-burst rate used in the likelihood.","marker":"[5]"},{"why":"Kilonova-associated long burst GRB 230307A, contributing the long-burst rate and the evidence that mergers can produce such bursts.","marker":"[8]"}],"fun_headline_variants":["Neutron-star merger mass of 1.36 TOV splits GRB durations","GW and GRB data set neutron-star mass boundary at 1.36 TOV","Mergers exceeding 1.36 TOV produce longer gamma-ray bursts","Characteristic mass 1.36 TOV reveals neutron-star survival","New method constrains neutron-star maximum mass via bursts"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The inference assumes that the observed ratio of long to short gamma-ray burst rates directly equals the ratio of gravitational-wave merger rates in two total-mass ranges, with no correction for beaming, detector selection, or other burst-duration drivers such as spin, mass ratio, or viewing angle.","fun_headline_variants_meta":{"raw":{"variants":["Neutron-star merger mass of 1.36 TOV splits GRB durations","GW and GRB data set neutron-star mass boundary at 1.36 TOV","Mergers exceeding 1.36 TOV produce longer gamma-ray bursts","Characteristic mass 1.36 TOV reveals neutron-star survival","New method constrains neutron-star maximum mass via bursts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000308,"raw_usage":{"total_tokens":1847,"prompt_tokens":1121,"completion_tokens":726,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":737,"completion_tokens_details":{"reasoning_tokens":628}},"tokens_in":737,"tokens_out":726,"duration_ms":7476,"temperature":1.0,"reasoning_tokens":628,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:54:52.919678+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete falsifier is a gravitational-wave binary neutron star merger with total mass above the inferred $k_2 M_{\\mathrm{TOV}}$ threshold (about $3.0\\,M_\\odot$ if $M_{\\mathrm{TOV}}\\simeq 2.2\\,M_\\odot$) that still produces a short gamma-ray burst, or a kilonova-associated long burst from a merger well below that threshold. A statistical version would be a future joint sample in which burst duration varies continuously with mass ratio or spin rather than switching at a single total mass.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Numerical relativity result placing the long-lived/short-lived remnant transition near 1.3-1.4 TOV masses, used as the comparison point for the inferred $k_2$."}],"review_version":2}