{"id":"ebf5ea44-8983-4a09-a2d3-98cd320f16f2","arxiv_id":"2507.09175","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For a GW200115-compatible BHNS with a 3.6 solar mass black hole and 1.5 solar mass neutron star, the neutron star may fill its Roche lobe before the ISCO, hinting at a catastrophic mass-transfer channel that could hide electromagnetic counterparts.","lead":"This paper builds quasi-equilibrium models of black hole neutron star binaries matching the masses and spins of the GW200105 and GW200115 gravitational wave events. For one GW200115-like configuration, the neutron star appears to overflow its Roche lobe before reaching the last stable orbit, suggesting a third possible outcome beyond plunge and tidal disruption.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Roche-lobe overflow claim for q=0.4167 conflicts with the paper's own κ-based mass-shedding estimate; the third-fate conclusion rests on the less GR-informed diagnostic.","rationale":"The reader correctly identifies the Newtonian/2PN Roche-lobe comparison as the weakest assumption. I agree with that diagnosis and with the conditional verdict: the paper is worth publishing as a parameter study, but the third-fate claim is not established. My stress-test sharpens the concern into an internal inconsistency: the same paper's κ-based mass-shedding estimate, extracted from the full GR QE solution, places r_MS below r_ISCO for the very sequence used to claim pre-ISCO Roche-lobe overflow. This makes the central claim vulnerable to a direct computational check, not merely to a general worry about Newtonian approximations. If a future time evolution shows no mass transfer before ISCO, the conclusion collapses; if it does show mass transfer, the paper needs a more careful justification of why the Newtonian lobe is the right criterion. The conditional verdict remains appropriate because the paper is transparent and the question is well posed, but the current evidence is insufficient to promote the third fate from hypothesis to result.","tokens_in":16830,"tokens_out":5365,"duration_ms":71562,"concrete_test":"For the q=0.4167, d=24.3 M_sun FUKA solution, compute κ and locate the saddle/cusp point of the actual XCTS enthalpy potential along the line joining the two centers, rather than the Newtonian potential of Eq. (40). If the stellar surface does not reach the saddle and κ stays around 0.68, the Roche-lobe overflow is not realized in the full GR solution. An even more decisive check is to evolve this initial data with the Einstein Toolkit, as the paper promises, and test whether any mass transfer begins before the binary reaches r_ISCO=23.34 M_sun.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that for q=0.4167 the NS fills its Roche lobe before reaching ISCO, based on R_NS=7.76 M_sun exceeding the Eggleton radius 7.44 M_sun and the 2PN radius 6.90 M_sun at d=24.3 M_sun. The problem is not merely that the Roche-lobe formulas are Newtonian; it is that the paper's own GR-informed mass-shedding diagnostic contradicts the inferred overflow. Table III for the same q=0.4167 sequence gives 1/κ=1.46 at d=24.3, i.e. κ≈0.68, far from the cusp value κ→0 that is the accepted signature of mass shedding in QE sequences. Inserting that κ into the paper's Eq. (38) gives r_MS=19.11 M_sun, which is smaller than r_ISCO=23.34 M_sun. Thus, by the κ-based estimate, the star does not shed mass before ISCO and would plunge first. The paper reports both numbers without reconciling them. The Roche-lobe comparison, by contrast, treats the deformed relativistic star as a sphere of radius R_NS in a Newtonian point-mass potential (Eq. 40), ignoring the actual tidal flattening measured by κ. The conclusion therefore depends on preferring the less GR-informed estimate over the sequence's own deformation indicator. This is the load-bearing weak spot: if κ is the better diagnostic, the third-fate claim fails; if the Roche-lobe comparison is better, the paper must explain why κ remains far from zero at an overflowing configuration. The paper's own admission that the analysis is 'only indicative' does not resolve this inconsistency.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper constructs quasi-equilibrium (QE) initial-data sequences for black hole-neutron star (BHNS) binaries with component masses and spins compatible with the gravitational-wave events GW200105 and GW200115, using the FUKA framework and the Togashi equation of state. For each sequence the authors locate the ISCO by minimizing the binding energy and estimate the mass-shedding radius using the Kyutoku et al. formula with a kappa-based elongation factor. They then add a Newtonian and post-Newtonian Roche-lobe analysis, comparing the neutron star proper radius to the Eggleton and 2PN Roche-lobe radii. For the GW200115-compatible configuration with q=0.4167, they find that at the last converged separation d=24.3 solar masses the proper radius R_NS=7.76 solar masses exceeds both Roche-lobe estimates, while the extrapolated ISCO radius is r_ISCO=23.34 solar masses. The paper interprets this as evidence that the neutron star may fill its Roche lobe before reaching ISCO, proposing a third fate beyond plunge or tidal disruption: catastrophic mass transfer that suppresses electromagnetic emission.","tokens_in":17198,"tokens_out":7435,"duration_ms":81156,"significance":"If established, the proposed third fate would be a new physical channel for BHNS mergers and would affect how GW200115-like events with no electromagnetic counterparts are interpreted. The paper also provides a systematic set of QE-sequence data (ISCO fits, mass-shedding radii, kappa values, and Roche-lobe radii) for many mass-ratio and spin combinations, which is a useful resource, and it includes a resolution study supporting the numerical accuracy of the sequences. However, the central third-fate conclusion rests on a Newtonian/2PN Roche-lobe comparison that is not reconciled with the paper's own kappa-based mass-shedding diagnostic, and on an extrapolated ISCO value with no uncertainty estimate. As presented, the claim is therefore provisional rather than established.","major_comments":[{"comment":"The central claim for q=0.4167 is not reconciled with the paper's own GR-informed mass-shedding diagnostic. At the last converged separation d=24.3 solar masses, Table III lists 1/kappa=1.46 (kappa about 0.68), far from the cusp value kappa to 0 that is the standard signature of mass shedding in QE sequences. Inserting this c_R into Eq. (38) with R_NS=7.76 solar masses and M_BH/M_NS=2.4 gives r_MS about 19.11 solar masses, which is smaller than the extrapolated r_ISCO=23.34 solar masses. The kappa-based estimate therefore predicts plunge before mass shedding, in direct conflict with the overflow conclusion drawn from Table IV. The paper reports both numbers without explaining why the Roche-lobe comparison should be preferred over the sequence's own deformation indicator; if the Roche-lobe comparison is correct, the paper must explain why kappa remains at about 0.68 at a configuration that is already overflowing. This tension is load-bearing for the third-fate claim.","section":"IV.C-IV.D, Tables III-IV, Eqs. (38)-(40)"},{"comment":"The overflow comparison compares a GR proper radius with Newtonian and 2PN Roche-lobe radii for spherical stars. R_NS is the proper circumferential radius of a deformed relativistic star, whereas R_LE (Eggleton) and R_LPN (Ratkovic et al.) are spherical approximations for point-mass binaries; the paper's own statements in Section V acknowledge this but do not quantify its effect. In particular, the Roche-lobe calculation uses the Newtonian potential Eq. (40) and ignores the tidal flattening that kappa measures. Without a GR-informed filling criterion (for example, comparing the effective-potential cusp from the FUKA solution with the stellar surface), the conclusion that the NS fills its Roche lobe before ISCO remains an unquantified approximation rather than a supported result.","section":"IV.D and Section V, Eq. (40), Table IV"},{"comment":"The ordering r_ISCO < d_min that underlies the 'before reaching ISCO' claim depends on an extrapolated value. For q=0.4167, d_ISCO=23.34 solar masses is obtained by fitting E_b(Omega) with a parabola and then inverting the relation M_infinity Omega = a/d + b, yet d_min=24.3 solar masses is the last converged point; the extrapolation distance is less than one solar mass. No error estimates, residual plots, or robustness checks are given for the fit parameters or for d_ISCO. Since the margin is small, modest extrapolation uncertainty could reverse the ordering and eliminate the third-fate scenario. Please quantify the uncertainty or extend the sequences, or soften the claim accordingly.","section":"IV.B, Eqs. (35)-(37), Table II"}],"minor_comments":[{"comment":"The sentence 'RN S > RLLE t' contains a typo; it should read R_NS > R_LE, and the table header should distinguish R_LP from R_LPN more clearly.","section":"IV.D and Table IV"},{"comment":"The abstract and introduction contain grammatical errors, including 'we are suggestion a third fate' and 'will be smaller of ISCO'; a careful language edit is needed.","section":"Abstract and Section I"},{"comment":"Reference [5] is malformed; the author list should be 'M. A. Abramowicz, M. Calvani, and L. Nobili'.","section":"References"},{"comment":"Table VII labels R_NS as the 'Areal radius' while Table III calls it the 'proper circumferential radius'; the same definition should be used consistently throughout.","section":"Tables III and VII"},{"comment":"Eq. (38) is typeset as 'rM S= 21/3cR'; the intended expression 2^{1/3} c_R (M_BH/M_NS)^{1/3} R_NS should be written unambiguously.","section":"Eq. (38)"}],"recommendation":"major_revision","confidential_remarks":"The main obstacle is the internal inconsistency between the kappa-based mass-shedding radius and the Roche-lobe overflow claim. I would not recommend acceptance until the authors either reconcile the two diagnostics with a GR-informed criterion or substantially soften the third-fate conclusion. The paper otherwise fits the journal's scope, and the QE-sequence data are useful."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper constructs new FUKA quasi-equilibrium sequences for BHNS configurations matched to GW200105 and GW200115, then applies Newtonian and post-Newtonian Roche-lobe diagnostics to ask whether the neutron star fills its Roche lobe before reaching ISCO. That's a genuinely new way to frame the problem, and the sequence data itself looks like solid work: convergence is checked, the Togashi EOS is specified, and the paper is candid that time evolution is needed to confirm any instability. Credit where it's due: the QE sequences and the comparison of GR circumferential radii with Roche-lobe estimates are not in the earlier literature, and the paper does not oversell its evidence.\n\nThe soft spot is exactly where the stress-test note lands. For q=0.4167, the paper's own κ-based mass-shedding radius from Eq. (38) is r_MS=19.11 M⊙, while r_ISCO=23.34 M⊙. That says the star reaches ISCO first and plunges. But the Roche-lobe comparison says R_NS=7.76 M⊙ exceeds the Eggleton radius 7.44 M⊙ and the 2PN radius 6.90 M⊙ at d=24.3 M⊙, implying overflow before ISCO. Both numbers appear in the same tables, and the paper never reconciles them. The κ value at that point is 0.68, far from the cusp value of zero that marks mass shedding in QE sequences. The Roche-lobe approach treats the deformed relativistic star as a sphere in a Newtonian potential, which the paper itself calls 'only indicative.' So the central third-fate conclusion depends on preferring the less GR-informed estimate over the sequence's own deformation indicator. That is a load-bearing inconsistency, not a minor caveat.\n\nOther concerns are smaller: the ISCO radii come from parabolic fits with no error bars, and only one EOS is used. Those are normal for initial-data studies and not fatal.\n\nThe paper is worth a serious referee because the QE sequences are useful and the question it asks is real, but the current version does not establish the third-fate scenario. A referee should focus on forcing either a reconciliation of the two mass-shedding diagnostics or an explicit statement that the Roche-lobe result is not yet viable. The paper's own recommendation of time evolution is the right one. I would send it to peer review with a request for major revision, and I would not cite the third-fate conclusion as fact until the inconsistency is addressed.","headline":"A useful set of QE sequences for GW200105/GW200115, but the third-fate claim rests on a Newtonian Roche-lobe comparison that contradicts the paper's own κ-based mass-shedding estimate, so the central conclusion is not yet established.","tokens_in":17745,"tokens_out":2174,"would_cite":false,"duration_ms":28748,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.25.D-","04.30.-w","97.60.Jd","97.80.-d"],"model":"deepseek-v4-flash","headline":"Quasi-equilibrium models of GW200115-like binaries show the neutron star overflowing its Roche lobe, and thus starting mass transfer, before the binary reaches the innermost stable circular orbit.","keywords":["black hole–neutron star binaries","quasi-equilibrium sequences","Roche lobe overflow","innermost stable circular orbit","mass transfer","GW200115","FUKA initial data","neutron star tidal disruption"],"falsifier":"Evolve the $q=0.4167$ configuration from its last converged FUKA initial data (separation $d=24.3$ solar masses, the TNTYST equation of state [23], spins −0.2 and +0.2) with a full general-relativistic hydrodynamics code and check whether mass transfer begins before the orbit reaches the ISCO and whether an accretion disk or electromagnetic counterpart forms; a cheaper check is to recompute the Roche condition from the actual deformed stellar surface rather than a sphere.","tokens_in":16625,"feed_emoji":"🕳️","tokens_out":15142,"duration_ms":149323,"temperature":0.7,"pith_summary":"Black hole–neutron star binaries whose masses and spins match the gravitational-wave events GW200105 and GW200115 are modelled as quasi-equilibrium sequences with the FUKA initial-data code, and the location of the innermost stable circular orbit is compared with two independent estimates of where the neutron star begins shedding mass. In most configurations the ISCO lies outside the mass-shedding radius, matching the standard picture that the neutron star plunges without yielding an electromagnetic counterpart. For the most probable GW200115 parameters ($q=0.4167$, black hole 3.6 solar masses, neutron star 1.5 solar masses), however, the neutron star's proper radius at the last converged separation exceeds both the Eggleton formula [30] and a second-post-Newtonian Roche lobe radius [32], while the ISCO sits at a slightly smaller separation. The authors therefore suggest that mass transfer can begin before the ISCO and, if it develops catastrophically, would suppress electromagnetic emission, a third possible fate beyond plunge and tidal disruption that could explain events like GW200115.","feed_headline":"Neutron star overflows its Roche lobe before the ISCO","feed_subtitle":"In a GW200115-like binary, mass transfer may start before plunge — a third fate beyond tidal disruption.","key_machinery":"The argument pivots on comparing the proper circumferential radius of the neutron star extracted from a FUKA quasi-equilibrium solution with two Roche lobe radii: the Eggleton formula [30] and a second-post-Newtonian formula [32], together with the Newtonian effective potential (Eq. 40) that locates the $L_1$ point. The ISCO of each sequence is found as the minimum of the binding energy, fitted as a quadratic in the orbital angular velocity, and the mass-shedding radius is estimated from the Newtonian formula (Eq. 38) with the elongation factor $1/\\kappa$ supplied by a mass-shedding indicator [26]. The $q=0.4167$ sequence is singled out because at its last converged separation the comparison reverses: $R_{NS} > R_{LE}$ and $R_{NS} > R_{LPN}$ while $r_{ISCO}$ is only slightly smaller.","core_discovery":"The central claim is that a neutron star in a black hole–neutron star binary can fill and overflow its Roche lobe at a separation larger than the ISCO, so that mass transfer precedes the final plunge. In the sequence matching GW200115 (mass ratio $q=0.4167$, black hole mass 3.6 solar masses, neutron star mass 1.5 solar masses, dimensionless spins −0.2 and +0.2), the last converged solution at $d=24.3$ solar masses gives a neutron star proper radius $R_{NS}=7.76$ solar masses, larger than the standard Roche lobe radius [30] of 7.44 solar masses and the second-post-Newtonian radius [32] of 6.90 solar masses, while the fitted ISCO is at the smaller separation $r_{ISCO}=23.34$ solar masses. The paper concludes that the neutron star is filling its Roche lobe before reaching the last stable orbit, and that if the ensuing mass transfer is catastrophic—on the model of the runaway instability seen in self-gravitating disks—no electromagnetic counterpart would be expected. This is presented as a third fate of the neutron star, beyond the usual plunge and tidal-disruption outcomes.","pith_inferences":["If a single full evolution validates the Roche-lobe screen, the same cheap comparison could map out large black hole–neutron star parameter spaces without solving for full initial data at every point.","The third fate is likely sensitive to the neutron star equation of state: a softer equation of state makes the star more compact and could push Roche filling inside the ISCO, so the prediction is testable with future events of different inferred compactness.","Because the paper compares a proper circumferential radius with spherical Newtonian and second-post-Newtonian Roche radii, a fully relativistic Roche criterion might shift the critical mass ratio; checking this would clarify how general the third fate is."],"forward_implications":["For GW200115-compatible parameters, a missing electromagnetic counterpart can be explained by catastrophic mass transfer before the ISCO, not only by plunge.","The $q=0.4167$ sequence is a concrete initial-data target for dynamical simulations that could confirm the onset of a runaway-type instability.","Mass transfer before the ISCO changes the mass ratio and angular momentum of the binary, so its subsequent inspiral would follow a different quasi-equilibrium track than the standard plunge picture.","The same Roche-lobe-versus-ISCO comparison can be applied to other black hole–neutron star candidates, for instance GW230529-like lower-mass-gap binaries, to predict whether they should show electromagnetic emission."],"supporting_citations":[{"why":"Supplies the standard mass-shedding radius formula (Eq. 38) and the expectation that these binaries plunge without an electromagnetic counterpart.","marker":"[3]"},{"why":"The FUKA initial-data code used to construct the quasi-equilibrium sequences and to read off the neutron star proper radius.","marker":"[14]"},{"why":"Provides the quasi-equilibrium sequence methodology, the mass-shedding indicator convention, and the fitting ansatz for the binding-energy minimum.","marker":"[16]"},{"why":"Gives the Newtonian effective-potential Roche-lobe analysis (Eq. 40) applied to the quasi-equilibrium configurations.","marker":"[4]"},{"why":"Formula for the Eggleton Roche lobe radius compared against the neutron star radius.","marker":"[30]"},{"why":"Second-post-Newtonian Roche lobe radius formula used in the comparison.","marker":"[32]"},{"why":"Defines the mass-shedding indicator kappa used to set the elongation factor in Eq. (38).","marker":"[26]"},{"why":"The TNTYST nuclear equation of state used for the neutron star matter in the sequences.","marker":"[23]"}],"fun_headline_variants":["Mass transfer starts before ISCO in black hole-neutron star binaries","Neutron star overflows Roche lobe, hinting at a third fate","GW200115: Neutron star fills Roche lobe before plunge","Pre-plunge mass transfer: a new fate for neutron stars","Roche lobe overflow precedes ISCO in BHNS binaries"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The assumption that carries the conclusion is that a Newtonian or second-post-Newtonian Roche-lobe radius, computed by treating the deformed neutron star as a sphere in an effective potential, correctly marks the true general-relativistic mass-shedding threshold; the paper itself labels this comparison 'only indicative' in Section V.","fun_headline_variants_meta":{"raw":{"variants":["Mass transfer starts before ISCO in black hole-neutron star binaries","Neutron star overflows Roche lobe, hinting at a third fate","GW200115: Neutron star fills Roche lobe before plunge","Pre-plunge mass transfer: a new fate for neutron stars","Roche lobe overflow precedes ISCO in BHNS binaries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00106,"raw_usage":{"total_tokens":4514,"prompt_tokens":1082,"completion_tokens":3432,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":698,"completion_tokens_details":{"reasoning_tokens":3342}},"tokens_in":698,"tokens_out":3432,"duration_ms":27505,"temperature":1.0,"reasoning_tokens":3342,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:02:29.540102+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evolve the $q=0.4167$ configuration from its last converged FUKA initial data (separation $d=24.3$ solar masses, the TNTYST equation of state [23], spins −0.2 and +0.2) with a full general-relativistic hydrodynamics code and check whether mass transfer begins before the orbit reaches the ISCO and whether an accretion disk or electromagnetic counterpart forms; a cheaper check is to recompute the Roche condition from the actual deformed stellar surface rather than a sphere.","supporting_citations":[{"cited_title":"From the equation of hy- drostatic equilibrium Eqs","cited_arxiv_id":null,"evidence_quote":"Supplies the standard mass-shedding radius formula (Eq. 38) and the expectation that these binaries plunge without an electromagnetic counterpart."},{"cited_title":"The equations result- ing from the spectral representation are solved by using Newton-Raphson iterative method [15]","cited_arxiv_id":null,"evidence_quote":"The FUKA initial-data code used to construct the quasi-equilibrium sequences and to read off the neutron star proper radius."},{"cited_title":"Nishida, A","cited_arxiv_id":null,"evidence_quote":"Provides the quasi-equilibrium sequence methodology, the mass-shedding indicator convention, and the fitting ansatz for the binding-energy minimum."},{"cited_title":"Ferrari, L","cited_arxiv_id":null,"evidence_quote":"Gives the Newtonian effective-potential Roche-lobe analysis (Eq. 40) applied to the quasi-equilibrium configurations."},{"cited_title":"Togashi, K","cited_arxiv_id":null,"evidence_quote":"Formula for the Eggleton Roche lobe radius compared against the neutron star radius."},{"cited_title":"Taniguchi, T","cited_arxiv_id":null,"evidence_quote":"Second-post-Newtonian Roche lobe radius formula used in the comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the mass-shedding indicator kappa used to set the elongation factor in Eq. (38)."}],"review_version":1}