{"id":"54fe825c-5ba5-4b7b-b35a-990075307d43","arxiv_id":"2607.07668","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":5,"one_line_summary":"Subsolar strange star mergers produce a lower post-merger-to-cutoff GW frequency ratio than neutron star mergers, cleanly separating the two classes across equations of state and mass ratios.","lead":"This paper presents the first numerical-relativity simulations of subsolar-mass binary strange star mergers, comparing them to subsolar neutron star mergers. It identifies a gravitational-wave frequency ratio that cleanly distinguishes the two star types, which matters for interpreting future subsolar-mass gravitational-wave detections.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The f2/fcut separation hinges on the q=1.5 SS case, whose convergence is non-monotonic: finest resolution gives 2.48 but coarser grids give 2.57–2.60, narrowing the gap to the NS branch to as little as ~0.05.","rationale":"The reader correctly identified the q=1.5 SS case as the marginal point and flagged the limited EOS sample as a concern. However, the reader framed the issue primarily as one of EOS representativeness and extrapolation to untested configurations (more extreme mass ratios, spinning binaries). The more immediate and concrete concern is the numerical convergence of the existing data point that defines the boundary: the q=1.5 SS case shows non-monotonic convergence where the finest resolution is the low outlier, and the spread across resolutions (2.48–2.60) brings it within ~0.05 of the NS minimum. This is a problem internal to the existing simulation grid, not just an extrapolation concern. That said, the reader's CONDITIONAL verdict remains appropriate. The paper is a well-executed first exploration with clear physical motivation, convergence checks, and honest discussion of limitations. The non-monotonic convergence of the boundary case reinforces why CONDITIONAL rather than ACCEPT is the right call, but it does not escalate to REJECT — the simulations are sound, the concern is about the precision of a boundary claim, and the paper itself acknowledges the marginality of this case. No code or data release is noted, which further supports keeping the verdict conditional pending independent verification.","tokens_in":17962,"tokens_out":4306,"duration_ms":142097,"concrete_test":"Re-run the SS1 0.4+0.6 binary at a fourth, higher resolution (N=120, Δx_finest ≈ 100 m) and perform Richardson extrapolation on f2 and fcut separately, then on their ratio. If the extrapolated f2/fcut exceeds ~2.55, the gap to the NS minimum (2.65) drops below 0.10 and the 'clean separation' claim weakens substantially. If it remains at or below 2.50, the separation is numerically robust for this boundary case.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of 'clean' separation rests disproportionately on a single extreme case: the q=1.5 SS binary (SS1 0.4+0.6), which yields f2/fcut = 2.48 at the finest resolution, versus the NS minimum of 2.65 (WFF1 0.5+0.5). However, Table S1 reveals that this exact case shows strongly non-monotonic convergence across the three resolutions: f2/fcut = 2.598, 2.567, and 2.484 from coarsest to finest. The finest resolution is the outlier on the low side, not a converged limit. If the true converged value lies closer to the coarser-resolution values (~2.57–2.60), the gap to the NS branch shrinks to ~0.05–0.08, which is comparable to the resolution-induced scatter the paper itself quotes (≲3.7% for f2/fcut). The paper uses the finest-resolution value as the headline number without Richardson extrapolation or a fourth resolution to confirm convergence direction. This is the single data point that determines whether the separation is truly 'clean' (no overlap) or merely 'suggestive' (thin margin with uncertain convergence). The paper acknowledges this case is 'extreme and isolated' and that f2 'rises onto the NS relation,' but does not flag that its own convergence for this case is the most problematic in the grid. The claim of non-overlapping separation is only as secure as this one number.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"This Letter presents the first numerical-relativity simulations of subsolar-mass binary strange star (SS) mergers, systematically comparing them with subsolar binary neutron star (NS) mergers across three equations of state (EOSs) per class, component masses 0.3–0.7 M⊙, and mass ratios q=1, 1.22, and 1.5. The simulations use established NR infrastructure (FUKA initial data, SACRA-K evolution, Z4c formulation, HLLC/HLLE solvers) with three grid resolutions per binary. The central claim is that the ratio f2/fcut of the post-merger peak frequency to the gravitational-wave cutoff frequency cleanly separates SS mergers (2.1–2.5) from NS mergers (2.65–2.97) with no overlap. The paper also reports ejecta properties and discusses electromagnetic counterpart prospects.","tokens_in":18247,"tokens_out":1524,"duration_ms":455868,"significance":"The paper is a genuine first: no prior NR simulations of subsolar-mass binary SS mergers exist. The f2/fcut discriminant is a falsifiable, parameter-free prediction derived from simulation dynamics rather than fitted to target frequencies. The quasi-universal relations are fitted to simulation outputs as functions of the tidal deformability Λ̃ computed from stellar models. The ejecta analysis, including the distinction between neutron-rich NS ejecta and decompressed quark-matter SS ejecta, adds astrophysical utility. The simulation grid spans a reasonable EOS range and multiple mass ratios, and the convergence assessment (half-spreads of 2.6–5% across resolutions) is reported transparently. These are substantial strengths for a Letter.","major_comments":[{"comment":"The central claim of 'clean' non-overlapping separation in f2/fcut rests disproportionately on the single extreme case SS1 0.4+0.6 (q=1.5), which yields f2/fcut = 2.48 at the finest resolution versus the NS minimum of 2.65 (WFF1 0.5+0.5). Table S1 reveals that this exact case shows strongly non-monotonic convergence: f2/fcut = 2.598, 2.567, and 2.484 from coarsest to finest. The finest-resolution value is the outlier on the low side, not a converged limit. If the true converged value lies closer to the coarser-resolution values (~2.57–2.60), the gap to the NS branch shrinks to ~0.05–0.08, comparable to the resolution-induced scatter the paper itself quotes (≲3.7% for f2/fcut). The paper uses the finest-resolution value as the headline number without Richardson extrapolation or a fourth resolution to confirm convergence direction. The authors should either (a) add a fourth resolution for,","section":null},{"comment":"this case to establish convergence, or (b) reframe the claim from 'clean separation with no overlap' to 'suggestive separation with a thin margin whose convergence is uncertain for the most extreme case.' The current phrasing in the abstract ('cleanly separates') and the main text ('do not overlap across the numerical models') overstates what the data support given this single load-bearing data point.","section":null},{"comment":"The quasi-universal relations and the f2/fcut discriminant are established on a grid of three SS EOSs and three NS EOSs. The paper itself notes that the q=1.5 SS case 'rises onto the NS relation' for f2, and the separation gap narrows to ~0.17 at the finest resolution (2.48 vs 2.65). Whether more extreme mass ratios (q>1.5), different SS EOS parameterizations (e.g., color-superconducting gaps, non-MIT-bag models), or rapidly spinning configurations would bridge this gap remains untested. The authors should explicitly acknowledge this limitation in the discussion of the discriminant's robustness, rather than stating it 'is a robust discriminant' without qualification. A brief statement that the claim is conditional on the simulated EOS and mass-ratio range would suffice.","section":null}],"minor_comments":[{"comment":"The abstract states f2/fcut 'cleanly separates the two classes'; given the convergence concern for the q=1.5 SS case, consider softening to 'separates the two classes across the simulated grid' or similar.","section":null},{"comment":"In the description of Fig. 3 (lower panel), the shaded band is stated to be at f2/fcut = 2.5–2.6, but the q=1.5 SS data point at 2.48 falls below this band. Clarify whether the band represents the proposed classification threshold or simply marks the visual gap.","section":null},{"comment":"The ejecta for the unequal-mass SS binaries (q=1.22 and 1.5) approaches the baryonic mass-conservation error level (M_ej < 2×10^{-4} M⊙), as noted in Fig. 4. Table S1 shows the Bernoulli-criterion ejecta for SS1 0.45+0.55 at the finest resolution is 1.50×10^{-2} M⊙, well above this floor, but the geodesic criterion gives 1.45×10^{-2}. The text should briefly comment on the consistency between the two criteria for the cases near the conservation floor.","section":null},{"comment":"The thermal adiabatic indices differ between SS (Γ_th = 4/3) and NS (Γ_th = 1.75). A brief justification for these choices, or a reference to where they are validated, would help readers assess sensitivity.","section":null},{"comment":"Reference [70] (SACRA-K) is listed as 'in preparation' (2026). If the code description is not yet available, the manuscript should provide enough detail in the Supplemental Material for reproducibility, which Table S2 partially addresses.","section":null},{"comment":"The text mentions 'animations of the snapshots and corresponding GWs can be found at [42]' — ensure the linked URL is persistent and accessible at publication.","section":null},{"comment":"Minor typographical issue: in the abstract and main text, 'cutofffrequencyfcut' appears to be missing a space (likely a LaTeX formatting artifact).","section":null}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the q=1.5 SS convergence is well-founded and is the primary reason for the major_revision recommendation. The non-monotonic convergence for this single load-bearing data point is visible in the authors' own Table S1, and they do not flag it. If a fourth resolution or Richardson extrapolation confirms the finest-resolution value, the claim stands; if not, the separation is suggestive rather than clean. The paper is otherwise a solid and timely contribution. The EOS sample limitation is a standard concern for quasi-universal-relation papers and is not by itself disqualifying, but the authors should acknowledge it explicitly."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive report. The two major comments are well-taken and we will revise the manuscript accordingly. Below we address each point.","responses":[{"response":"The referee is correct that the SS1 0.4+0.6 (q=1.5) case is the most marginal data point and that its convergence is non-monotonic. We have re-examined Table S1 carefully. The three-resolution sequence for f2/fcut is 2.598, 2.567, 2.484 (coarse to fine), so the finest value is indeed the lowest and the convergence direction is not established. We acknowledge that if the true converged value lies closer to the coarser-resolution values (~2.57–2.60), the gap to the NS minimum (2.65 at finest resolution for WFF1 0.5+0.5) narrows to ~0.05–0.08, which is comparable to the ~3.7% resolution half-spread we quote for f2/fcut. This is a fair concern. We do not currently have a fourth resolution for this case and cannot perform a reliable Richardson extrapolation given the non-monotonic trend. We will therefore revise the manuscript as follows: (1) We will soften the abstract from 'cleanly separates' to 'separates' and add a qualifier noting that the margin is thinnest for the most asymmetric SS case, whose convergence is not yet fully established. (2) In the main text, we will replace 'do not overlap across the numerical models' with a statement that the two branches are separated across our simulation grid, with the smallest margin (~0.16 at finest resolution, possibly as small as ~0.05 if the coarser-resolution values for the q=1.5 SS case are closer to convergence) occurring for the most asymmetric SS binary. (3) We will explicitly flag the non-monotonic convergence of the SS1 0.4+0.6 case in the text and note it as a limitation. We agree that 'clean separation with no overlap' overstates what the current data support for this single load-bearing data point.","revision_made":"yes","referee_comment":"The central claim of 'clean' non-overlapping separation in f2/fcut rests disproportionately on the single extreme case SS1 0.4+0.6 (q=1.5), which shows strongly non-monotonic convergence. The finest-resolution value is the outlier on the low side. The paper uses the finest-resolution value as the headline number without Richardson extrapolation or a fourth resolution."},{"response":"We agree. The claim of robustness is currently stated without sufficient qualification regarding the explored parameter space. We will add a sentence in the discussion of the discriminant explicitly noting that the separation is established within the simulated range of EOSs (three MIT-bag-family SS models and three nuclear NS models), mass ratios (q ≤ 1.5), and non-spinning configurations, and that its persistence under more extreme mass ratios, alternative SS EOS parameterizations (e.g., color-superconducting gap models or non-MIT-bag constructions), or rapidly spinning progenitors remains untested. We will also adjust the phrase 'is a robust discriminant' to 'is a promising discriminant within the simulated EOS and mass-ratio range.' This is an honest reflection of what the simulations cover.","revision_made":"yes","referee_comment":"The quasi-universal relations and the f2/fcut discriminant are established on a grid of three SS EOSs and three NS EOSs. More extreme mass ratios, different SS EOS parameterizations, or rapidly spinning configurations could bridge the gap. The authors should explicitly acknowledge this limitation."}],"tokens_in":17813,"tokens_out":1340,"duration_ms":108739,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"This is the first numerical-relativity study of subsolar-mass binary strange star mergers, and it systematically compares them against subsolar binary neutron star mergers across three EOSs per class, several masses, and mass ratios up to q=1.5. The physical picture is clear and well-argued: self-bound strange stars stay compact through inspiral, collide harder, drive a stronger shock and radial bounce, and end up with a lower post-merger peak frequency f2 despite being more compact before merger. The quasi-universal relations for fcut, fmerger, and f2 versus tidal deformability are a useful organizing framework, and the ejecta analysis — including honest reporting of the mass-conservation error floor for unequal-mass SS binaries — is solid. The simulation infrastructure (SACRA-K, FUKA initial data, Z4c, three resolutions per binary) is established and the convergence reporting is mostly adequate. The animations linked in Figure 1 are a nice touch. This is real computational work and the authors deserve credit for opening up this parameter space.","headline":"First NR simulations of subsolar-mass binary strange star mergers; the f2/fcut discriminator is promising but the 'clean separation' claim rests on a case with non-monotonic convergence.","tokens_in":18770,"tokens_out":1390,"would_cite":true,"duration_ms":104982,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Ratio of two gravitational-wave frequencies cleanly separates strange stars from neutron stars","keywords":[],"falsifier":"If a simulation with a different strange-star equation of state, a more extreme mass ratio, or rapidly spinning components produced f_2/f_cut above 2.5 — or if a neutron-star model produced it below 2.65 — the clean separation would break.","tokens_in":18025,"feed_emoji":"⭐","tokens_out":1529,"duration_ms":136428,"temperature":0.7,"pith_summary":"This paper presents the first full numerical-relativity simulations of subsolar-mass binary strange star mergers, comparing them head-to-head with binary neutron star mergers across a grid of equations of state, masses, and mass ratios. Strange stars are hypothesized compact objects made of deconfined quark matter, self-bound by the strong interaction rather than by gravity. Their self-bound nature makes them more compact and less tidally deformable than neutron stars of the same mass, and this structural difference drives qualitatively different merger dynamics. The strange star's sharp surface and compactness mean it stays intact until contact, then collides violently with strong shock heating and a large radial bounce; the neutron star, being more extended, develops tidal spiral arms and sheds mass before contact. These differing dynamics leave opposing imprints on two characteristic gravitational-wave frequencies: the cutoff frequency f_cut (where tidal effects accelerate the late inspiral) is higher for strange stars because they are more compact, while the dominant post-merger frequency f_2 is lower for strange stars because the violent bounce lowers the average density of the remnant. The paper shows that within each class, these frequencies follow quasi-universal power-law relations with the tidal deformability. The compounding of the two opposing shifts in the ratio f_2/f_cut produces a clean separation: strange stars cluster at 2.1-2.5, neutron stars at 2.65-2.97, with no overlap across the entire simulated grid. Both classes eject roughly 10^-2 solar masses of material, but of different composition — neutron-rich matter for neutron stars and decompressed quark matter for strange stars — offering a complementary electromagnetic test.","feed_headline":"Two GW frequencies separate strange stars from neutron stars","feed_subtitle":"First simulations of subsolar strange-star mergers reveal a ratio that cleanly distinguishes them from neutron stars — no overlap across the","key_machinery":"The central mechanism is the structural difference between self-bound strange stars (compact, sharp-surfaced, R ∝ M^{1/3}) and gravitationally bound neutron stars (extended, growing radius at low mass). This difference propagates through the merger in opposite directions for two frequencies: it raises f_cut for strange stars (less tidal deformation delays the inspiral breakdown to higher frequency) but lowers f_2 (the violent radial bounce and shock re-expansion lower the remnant density). The ratio f_2/f_cut compounds both shifts, producing the discriminant.","core_discovery":"The ratio f_2/f_cut of the post-merger peak frequency to the gravitational-wave cutoff frequency cleanly separates subsolar-mass binary strange star mergers (2.1-2.5) from binary neutron star mergers (2.65-2.97) with no overlap across the simulated grid of equations of state, masses, and mass ratios. This separation arises because the self-bound, compact strange star reaches a higher cutoff frequency (less tidal deformation during inspiral) but a lower post-merger frequency (violent shock and radial bounce lower the remnant's average density), and the two effects compound in the ratio. The paper establishes this by running the first numerical-relativity simulations of subsolar strange star二进","pith_inferences":["The discriminant's robustness against rapid stellar rotation is untested. If either class can be spun up before merger, centrifugal flattening could alter the effective compactness and tidal deformability enough to shift f_cut and f_2 in ways that might narrow or bridge the gap.","If the f_2/f_cut separation holds for mixed binaries (one strange star, one neutron star), the ratio could also diagnose the composition of individual components, not just homogeneous binaries — though the paper does not simulate this case.","The detectability of f_2 requires high signal-to-noise in the post-merger signal, which for subsolar masses at realistic distances likely demands third-generation detectors. The practical utility of the discriminant may therefore be gated by detector sensitivity rather than physics.","If strange quark matter ejected from an SS merger does not fragment into electromagnetically dark nuggets but instead evaporates into neutron-rich nucleons, the kilonova signatures of the two classes could be more similar than expected, making the gravitational-wave discriminant the primary rather than complementary diagnostic."],"forward_implications":["If a subsolar-mass compact-binary merger is detected by third-generation gravitational-wave detectors, measuring f_2/f_cut would allow a direct classification as either a strange star or neutron star binary, testing the Bodmer-Witten conjecture that strange quark matter is the true ground state of baryonic matter.","The quasi-universal relations between characteristic frequencies and tidal deformability, if they hold beyond the simulated grid, could reduce the parameter space needed for waveform templates in subsolar-mass searches.","The different ejecta composition — neutron-rich matter versus decompressed quark matter — means an electromagnetic counterpart (or its absence) to a subsolar merger would provide an independent test of the strange-star hypothesis, complementary to the gravitational-wave discriminant.","The rapidly rotating remnant with ~10^52 erg of rotational energy could power a synchrotron transient from radio to X-ray, offering a third observational channel to distinguish the two classes."],"fun_headline_variants":["A frequency ratio distinguishes subsolar strange stars from neutron stars","First simulations show strange stars and neutron stars have distinct GW ratios","Subsolar strange stars and neutron stars separated by GW frequency ratio","Strange star or neutron star? A gravitational wave ratio reveals the difference","Simulated subsolar mergers reveal a distinct GW frequency ratio for strange stars"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The clean, non-overlapping separation rests on a grid of three strange-star equations of state and three neutron-star equations of state, with mass ratios up to 1.5 and no stellar rotation. Whether more extreme mass ratios, different equation-of-state parameterizations, or rapidly spinning configurations would bridge the gap between the two classes remains untested. The paper itself notes that the q=1.5 strange-star case already pushes f_2/f_cut to 2.48, approaching the NSs' ","fun_headline_variants_meta":{"raw":{"variants":["A frequency ratio distinguishes subsolar strange stars from neutron stars","First simulations show strange stars and neutron stars have distinct GW ratios","Subsolar strange stars and neutron stars separated by GW frequency ratio","Strange star or neutron star? A gravitational wave ratio reveals the difference","Simulated subsolar mergers reveal a distinct GW frequency ratio for strange stars"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":999,"prompt_tokens":524,"completion_tokens":475,"prompt_tokens_details":null},"tokens_in":524,"tokens_out":475,"duration_ms":29672,"temperature":1.0,"reasoning_tokens":458,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T03:00:01.154839+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If a simulation with a different strange-star equation of state, a more extreme mass ratio, or rapidly spinning components produced f_2/f_cut above 2.5 — or if a neutron-star model produced it below 2.65 — the clean separation would break.","supporting_citations":[],"review_version":1}