{"id":"799799a3-8fc3-4227-96e0-1c955e618501","arxiv_id":"2505.05900","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A Monte Carlo population forecast predicts 6 to 22 spinning neutron star components in LISA-detectable double neutron stars will be detectable by Cosmic Explorer, with a moment-of-inertia accuracy near 8%.","lead":"This paper predicts how many spinning neutron stars inside double-neutron-star binaries could be spotted by the future Cosmic Explorer detector, after LISA finds the binaries first. It estimates 6 to 22 detectable neutron stars and a moment-of-inertia measurement accuracy of about 8 percent.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dual-line and I3 claims rest on the unverified waveform expansion in Eqs. (4)-(5); if the O(gamma*kappa) and O(kappa^2) amplitude coefficients are wrong, the inversion (21) and the ~8% I3 accuracy (23) do not follow.","rationale":"The manuscript is a competent population forecast and the authors disclose many caveats. The 6/22 detection numbers for the h(2a) component are driven by the leading-order ellipticity term, whose amplitude is the standard continuous-wave expression and is therefore not the main risk. The genuinely new claims -- the dual-line inversion of kappa and gamma and the ~8% measurement of I3 -- depend on the relative amplitudes of h(2a), h(2b), and h(2c) in Eq. (20), which are extracted from an unverified expansion in Eqs. (4)-(5). Because the paper provides no derivation and no numerical check of those amplitude factors, a reader cannot distinguish a correct extension from an algebraic mistake. The reader's weakest_assumption identifies exactly this point, and I agree. The appropriate verdict remains CONDITIONAL: the forecast is reasonable conditional on the waveform model, but the model needs independent verification. A secondary reproducibility issue (Table I omits the epsilon sampling needed for kappa in Eq. (14)) also supports the conditional verdict, but the waveform model is the single most load-bearing item.","tokens_in":16826,"tokens_out":14436,"duration_ms":131294,"concrete_test":"Perform an independent derivation of the waveform to O(gamma^2, kappa^2) for a rigid triaxial star under spin-orbit precession, e.g., by numerical integration of the quadrupole formula for the same rigid-body dynamics, and compare the Fourier amplitudes at the six frequencies Omega_r +/- Omega_p, Omega_r + 3 Omega_p, 2 Omega_r, and 2(Omega_r +/- Omega_p) with Eqs. (4)-(5). If any amplitude coefficient (e.g., 16, 64, gamma^2 + 64 kappa^2) differs by more than a few percent, the inversion formulas (21) and the I3 error budget (23) are invalid. Also verify that no same-order terms are missing from h(2a) or h(2c).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central novelty is the dual-line measurement of NS structural parameters, and every quantitative statement in Sec. III C (the inversion formulas Eq. (21) and the moment-of-inertia error budget Eq. (23)) is an algebraic consequence of the amplitude factors in Eq. (20). Those amplitude factors are taken from the waveform expansion in Eqs. (4)-(5), which the paper presents without derivation. The expansion treats kappa and gamma as independent small parameters and keeps terms to O(gamma^2, kappa^2) including O(gamma*kappa), but no check is shown that the retained terms are complete or correctly normalized. In particular, the frequencies and amplitudes of h(2b) and h(2c) (Eqs. (4e)-(4f)) determine whether the combinations h2b0 - h2c0 and h2a0 + 4h2c0 appearing in Eq. (21) isolate gamma and kappa as claimed. A missing or miscomputed term of the same order (e.g., an O(gamma^2) correction to h(2a), or an O(gamma*kappa) term in h(2c)) would bias kappa and gamma and propagate directly into the I3 error estimate. The paper references Refs. [17,52,59] for the underlying formalism, but the reduction to the independent-kappa,gamma case is not shown. Because Eq. (23) and Fig. 5 feed the headline '~8%' accuracy, this unverified expansion is the load-bearing element. A secondary ambiguity is that Table I does not list the equatorial ellipticity epsilon that enters kappa through Eq. (14), so the kappa prior used in the population is not fully specified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies Galactic double neutron star (DNS) systems as dual-line gravitational-wave sources, combining low-frequency binary inspiral signals detectable by LISA with high-frequency continuous waves from individual spinning, precessing, triaxially deformed neutron stars detectable by Cosmic Explorer. The authors extend a previously published waveform expansion to second order in the structural parameters κ (nonsphericity) and γ (wobble angle), treating them as independent, and use a LISA-detectable DNS population from Wagg et al. (2022) to synthesize a spinning-neutron-star population under uniform and log-uniform priors. They report detectable counts and SNRs for the h(1a), h(2a), h(2b), and h(2c) components, derive algebraic inversion formulas for κ and γ from the h(2x) amplitudes, and estimate a median relative moment-of-inertia accuracy of about 8%.","tokens_in":17167,"tokens_out":6337,"duration_ms":63665,"significance":"If the waveform expansion is correct, the paper supplies a concrete, testable pathway for using dual-line DNS detections to constrain neutron-star internal structure, and it provides an explicit comparison of how prior choices affect detectability forecasts. The inversion formulas in Eq. (21) are transparent and internally consistent, and the Monte Carlo framework is clearly described. However, the central quantitative claims—the detection counts, the SNR distributions, and the 8% moment-of-inertia accuracy—depend directly on waveform terms that are asserted without derivation, and on a population-selection rule that is stated ambiguously. The paper would be valuable after these load-bearing issues are resolved, but in its current form the headline numbers are not yet fully supported.","major_comments":[{"comment":"The waveform expansion is the load-bearing element of the paper, yet it is presented without derivation. The inversion formulas in Eq. (21) and the error budget in Eq. (23) are algebraic consequences of the amplitude coefficients in Eq. (20), which in turn come from the expansions in Eqs. (4)-(5). A missing, mis-normalized, or incomplete term of the same order—for example, an O(γ^2) correction to h(2a), an O(γκ) term in h(2c), or a different O(κ^2) coefficient in h(2b)—would directly bias the inferred κ, γ, and I3. Please provide a derivation of the independent-κ,γ expansion starting from Refs. [17,52,59], or at least a detailed consistency check that verifies completeness to O(γ^2,κ^2), including the limit ε→0, the limit κ~O(γ^2), and order-by-order comparison with Ref. [17].","section":"Sec. II B, Eqs. (4)-(5) and Sec. III C, Eq. (20)"},{"comment":"The selection rule for the detectable-source counts is ambiguous. The sentence \"we exclude cases where the combined SNR of all waveform components exceeds 7\" appears to remove the loudest cases, while the following clause suggests that the intended exclusion is of cases where only the combined SNR is above threshold but no individual waveform component is detectable. The counts in Table II and the moment-of-inertia accuracy sample in Fig. 5 depend directly on this choice, so the exact criterion used in the Monte Carlo simulation needs to be stated precisely.","section":"Sec. III B, selection rule"},{"comment":"Table I lists the moment of inertia I3, the oblateness ǫ, and the wobble angle γ, but it does not list the equatorial ellipticity ε, even though Eq. (14) defines κ in terms of both ε and ǫ. Without specifying how ε is sampled—its range and whether uniform or log-uniform sampling is applied—the κ prior is not fully defined and the population simulation is not reproducible. Please add ε to Table I and state its sampling distribution explicitly.","section":"Sec. II C 2, Eq. (14), and Table I"},{"comment":"The sample over which the ~8% moment-of-inertia accuracy is computed is not defined. The inversion in Eq. (21) requires the amplitudes h2a0, h2b0, and h2c0, so systems with only one detectable h(2x) component cannot contribute to this measurement. The paper should state whether the median reported in Fig. 5 is conditioned on all three components being detectable, give the size of that subset, and clarify how the selection rule in Sec. III B affects this subset. As written, the headline accuracy is not tied to a well-defined population.","section":"Sec. III C and Fig. 5"}],"minor_comments":[{"comment":"The phrase \"In the limit as Ωpre→0\" is confusing because the displayed frequencies still contain Ωp, the free-precession frequency; please state explicitly that this limit sets α=0 while leaving Ωp unchanged.","section":"Sec. II B, Eq. (6)"},{"comment":"Reference [56] lists duplicated author names (\"T. Wagg, T. Wagg, K. Breivik, K. Breivik\"); the entry should be corrected.","section":"References"},{"comment":"The histograms in Fig. 3 are labeled \"Probability\" for counts obtained from 1000 Monte Carlo simulations; please specify the binning and use a consistent label such as \"frequency\" or \"probability\".","section":"Fig. 3"},{"comment":"The statement that h(1b) and h(1c) are \"too weak to exceed the detection thresholds\" would be more informative if accompanied by their median SNR or an upper bound, so the reader can assess how far below threshold they fall.","section":"Sec. III B and Table II"}],"recommendation":"major_revision","confidential_remarks":"The main concern is verifiability rather than novelty: the dual-line DNS idea is timely and the population comparison is useful, but the paper's central quantitative claims rest on an unverified waveform expansion and an ambiguously defined selection rule. An appendix deriving or validating Eqs. (4)-(5) would resolve the most serious issue. I see no grounds for rejecting the manuscript outright, provided these points are addressed carefully."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Honest take: this is a competent scenario forecast, not a breakthrough. The genuinely new pieces are the O(γ², κ²) extension of the spinning-NS waveform (including the O(γκ) and O(κ²) terms) and the integration of the Wagg et al. LISA DNS catalog with a population of spinning NS components. The Monte Carlo work is transparent: 1000 iterations, uniform vs log-uniform priors, detector antenna patterns, and the main caveats are openly stated (35 systems is optimistic, only ~60% of DNSs are reliably identifiable, conservative models give ~8). The spin-period/eccentricity fit to the Dewi et al. data is a sensible way to assign spins.\n\nThe soft spot is exactly where the reader and stress-test put it. Eqs. (4)-(5) present the extended waveform without derivation. Since Eq. (21) inverts those amplitudes to solve for γ and κ, and Eq. (23) propagates them into a moment-of-inertia accuracy, any missing or miscomputed term at the same order changes the conclusions. The references [17,52,59] point to the formalism, but the specific reduction to independent κ,γ is not shown. This is load-bearing and fixable: the authors should either supply the derivation or verify the coefficients numerically against a known code.\n\nTwo smaller issues. The selection rule in Sec. III B is worded confusingly — 'exclude cases where the combined SNR exceeds 7' cannot be what they mean, given the counts in Table II; I assume they meant exclude cases where no individual component reaches threshold. And Table I omits the equatorial ellipticity ε that enters κ through Eq. (14), so the κ prior is under-specified. Both are easy fixes. Also, the abstract's ~8% I3 accuracy applies only to the subpopulation where both the h(2a) and h(2c) lines are detected; the joint count is smaller than the individual counts in Table II, so the claim needs that qualifier.\n\nNo circularity or fatal error. The algebra in Eqs. (20)-(21) is internally consistent, and the paper is honest about the caveats. The citation pattern looks fine — self-citations here are to prior work that this paper directly extends.\n\nWho it's for: CW search strategists and anyone interested in a concrete list of targets for LISA + Cosmic Explorer. The population part is useful even if the waveform details are pending. But I would not cite the waveform extension until the derivation is out.\n\nRecommendation: send to peer review. It deserves referee time, with the request that the derivation of Eqs. (4)-(5) be added or verified, and the selection rule and ε prior clarified. If the expansion holds up, this is a solid forecast; if not, the numbers shift. A serious referee can sort that out.","headline":"A useful dual-line DNS forecast whose central waveform terms are asserted rather than derived; the population synthesis is transparent, but the headline 8% accuracy needs its caveats restored.","tokens_in":17764,"tokens_out":3947,"would_cite":false,"duration_ms":37766,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.30.-w","95.55.Ym","97.60.Jd"],"model":"deepseek-v4-flash","headline":"Galactic double neutron stars can act as dual-line gravitational-wave sources, and joint LISA–Cosmic Explorer observations could measure the neutron star moment of inertia to about 8%.","keywords":["double neutron stars","gravitational waves","LISA","Cosmic Explorer","continuous gravitational waves","neutron star moment of inertia","precessing triaxial neutron star","population synthesis"],"falsifier":"Direct a 4-year Cosmic Explorer search at the sky position, orbital frequency, and eccentricity of the loudest LISA-resolved Galactic double neutron star from the same population model; if no h(2a) line appears above an SNR of 7 with the amplitude predicted by Eq. (20a), the optimistic count of 6–22 and the waveform amplitudes behind it are falsified for that source.","tokens_in":16558,"feed_emoji":"🌌","tokens_out":8928,"duration_ms":83208,"temperature":0.7,"pith_summary":"Double neutron stars in the Milky Way can emit gravitational waves on two very different frequency scales at once: the binary orbit radiates in the millihertz band, where LISA can see it, and each triaxially deformed, precessing neutron star radiates near its spin frequency and its second harmonic, where next-generation ground detectors like Cosmic Explorer are sensitive. The paper extends the waveform model for such spinning neutron stars to second order in the two structural parameters that control nonsphericity and wobble, then plugs in a simulated population of 35 LISA-detectable DNS systems. It predicts that a 4-year Cosmic Explorer search directed at these LISA sources will recover 6 spinning components under log-uniform parameter sampling and 22 under uniform sampling, with median signal-to-noise ratios around 20–30. For systems detected in both bands, the three amplitudes near twice the spin frequency can be inverted to solve for the neutron star's nonsphericity and wobble angle, and the moment of inertia can be inferred to roughly 8% relative accuracy. If right, this gives a direct route from gravitational-wave data to neutron star interior physics.","feed_headline":"Dual-line gravitational waves may measure neutron-star inertia to 8%","feed_subtitle":"A LISA+Cosmic Explorer search could recover 6–22 spinning neutron stars and pin down their moment of inertia.","key_machinery":"The load-bearing object is the six-component waveform of a triaxially deformed neutron star undergoing geodetic precession in a binary, written as h+ and h× sums oscillating at Ωr+Ωp, Ωr−Ωp, Ωr+3Ωp, 2Ωr, 2(Ωr+Ωp), and 2(Ωr−Ωp), with amplitude factors h2a0, h2b0, and h2c0 at the 2Ωr family. The crucial identities are the inversion formulas $\\gamma = 2\\sqrt{(h_{2b0}-h_{2c0})h_{2c0}}/(h_{2a0}+4h_{2c0})$ and $\\kappa = h_{2c0}/[4(h_{2a0}+4h_{2c0})]$, plus the moment-of-inertia error propagation that combines the SNRs of h2a, h2c, and the LISA DNS SNR. These formulas convert line amplitudes into structural parameters, which is what lets a dual-line detection measure I3.","core_discovery":"The paper's central claim is that a resolved double neutron star binary is simultaneously a millihertz gravitational-wave source and a source of two high-frequency line families, and that observing both bands with LISA and Cosmic Explorer turns the high-frequency amplitudes into a measurement of neutron-star structure. Concretely, the paper forecasts that among the ~35 DNS systems LISA would resolve in 4 years under the optimistic population model, 6 (log-uniform sampling) or 22 (uniform sampling) spinning neutron star components would also be detected by Cosmic Explorer, with median SNR between roughly 20 and 30. It further claims that once the h(2a), h(2b), and h(2c) lines are measured, the nonsphericity κ and wobble angle γ can be algebraically inverted from amplitude ratios, and with the LISA-derived distance the moment of inertia I3 follows with relative accuracy ~8%. The 8% figure is presented as insensitive to the two sampling choices; detection counts are not.","pith_inferences":["Editorial extension: the same amplitude inversion could be applied to any precessing triaxial neutron star with two detected sidebands, even without LISA, so the method generalizes to isolated pulsars if h(2b) and h(2c) are ever seen.","Editorial extension: the SNR-only error budget ignores covariances among κ, γ, and I3; a full Bayesian recovery on simulated data would likely widen the 8% error, so the quoted accuracy is an idealized lower bound.","Editorial extension: LISA and TianQin both cover the millihertz band, so the orbit-finder role in the dual-line strategy could be played by either space mission, and the same population forecast could be tested against future LISA mock data challenges."],"forward_implications":["A 4-year Cosmic Explorer directed search pointed at the ~35 LISA DNS positions should yield a few to a few dozen h(2a) detections, with central values 22 (uniform) and 6 (log-uniform), so the two sampling priors can be empirically distinguished by counting.","The h(1b) and h(1c) sidebands are predicted to be below threshold for both sampling schemes, meaning the usable dual-line information sits in the 2Ωr family and the h(1a) line.","When h(2a) and h(2c) are both measured, the relative moment-of-inertia error is centered at 8% with a 16%–84% spread from 4% to 13%, independent of the sampling choice.","If multiple dual-line DNS systems are detected, a Bayes factor on the number counts can decide whether NS structural parameters follow a log-uniform or uniform distribution.","Under the more conservative LISA model with only 8 DNS detections, the dual-line counts shrink to roughly a quarter of the optimistic values, so the 6–22 forecast scales with the LISA catalog size."],"supporting_citations":[{"why":"Supplies the 35 LISA-detectable DNS systems and their mass, distance, frequency, and eccentricity distributions used to seed the spinning-NS population.","marker":"[48]"},{"why":"Provides the base spinning-NS waveform with spin-orbit coupling that this work extends to O(γ^2,κ^2).","marker":"[17]"},{"why":"Defines the precessing triaxial NS waveform and the SNR integral used for Cosmic Explorer detectability.","marker":"[51]"},{"why":"Introduces the independent treatment of κ and γ that motivates the higher-order expansion.","marker":"[52]"},{"why":"Provides the simulated pulsar spin period versus orbital eccentricity relation fitted as Eq. (7).","marker":"[61]"},{"why":"Gives LISA's fractional distance error ΔD/D = 1/ρDNS used in the I3 error budget.","marker":"[50]"},{"why":"Supplies the relation Ωp ≃ ǫΩr and the claim that ǫ is measurable to much better than 10^-4, used to set the I3 accuracy.","marker":"[59]"},{"why":"Motivates the minimum neutron-star ellipticity floor ε ≈ 10^-9 adopted in the population priors.","marker":"[70]"}],"fun_headline_variants":["Dual-line GWs pin neutron-star inertia to 8%","LISA+CE: 6-22 spinning NSs, inertia to 8%","Neutron-star inertia to 8% via dual-line GW","Dual-line detection: NS structure at 8% accuracy","GW dual lines: inertia precision 8% for NSs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The forecast rests on treating the neutron star's nonsphericity κ and wobble angle γ as independent and expanding the waveform to second order in both; if the new O(γκ) and O($κ^{2}$) terms in the polarization amplitudes are wrong, the inversion formulas for γ, κ, and I3 fail and the detection counts shift.","fun_headline_variants_meta":{"raw":{"variants":["Dual-line GWs pin neutron-star inertia to 8%","LISA+CE: 6-22 spinning NSs, inertia to 8%","Neutron-star inertia to 8% via dual-line GW","Dual-line detection: NS structure at 8% accuracy","GW dual lines: inertia precision 8% for NSs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000225,"raw_usage":{"total_tokens":1469,"prompt_tokens":953,"completion_tokens":516,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":422}},"tokens_in":569,"tokens_out":516,"duration_ms":5834,"temperature":1.0,"reasoning_tokens":422,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:54:29.099183+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Direct a 4-year Cosmic Explorer search at the sky position, orbital frequency, and eccentricity of the loudest LISA-resolved Galactic double neutron star from the same population model; if no h(2a) line appears above an SNR of 7 with the amplitude predicted by Eq. (20a), the optimistic count of 6–22 and the waveform amplitudes behind it are falsified for that source.","supporting_citations":[{"cited_title":"Double Compact Objects as Low-frequency Gravitational Wave Sources","cited_arxiv_id":"0811.1602","evidence_quote":"Defines the precessing triaxial NS waveform and the SNR integral used for Cosmic Explorer detectability."},{"cited_title":"Gravitational-wave radiation from double compact objects with eLISA in the Galaxy","cited_arxiv_id":"1404.3848","evidence_quote":"Introduces the independent treatment of κ and γ that motivates the higher-order expansion."}],"review_version":1}