{"id":"94d95b6c-c505-492f-a272-7e2753f6c300","arxiv_id":"2608.06944","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"By injecting eccentric binary signals into simulated PTA data, the paper shows individual masses are recoverable at high frequencies, while low-frequency signals can be confused with the stochastic background and yield Bayes factors near unity.","lead":"Simulations of pulsar timing array data with eccentric supermassive black hole binaries show that the individual masses of the two black holes become measurable when the signal is strong and at high frequency. They also show that low-frequency eccentric signals can be partially absorbed into the stochastic background model, making them hard to identify.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The mass-measurement claim is built entirely on the authors' own post-leading-order eccentric waveform [9]; because injections and fits share that model, no injection-recovery test can expose waveform systematics, so an independent cross-check is required.","rationale":"I agree with the reader's weakest assumption and consider it load-bearing. The abstract makes a physical claim about measurability of masses, not merely a statement of numerical self-consistency. Section IVA's demonstration depends on the post-leading-order evolution of the binary over tau_alpha; if that evolution is wrong, the likelihood surface can still be sharp but shifted. No independent waveform or code comparison is provided. The proposed test, generating data with an independent eccentric-waveform implementation and recovering with the authors' model, would settle whether the mass measurement is robust or model-dependent. I do not see an internal inconsistency that would justify rejection: the Bayes factor near unity and the free-spectrum results are honestly reported, and the CRN and Hellings-Downs limitations are acknowledged. The absence of released analysis scripts and the qualitative recovery-accuracy statements are secondary reproducibility concerns, not independent grounds for changing the verdict. The appropriate verdict therefore remains CONDITIONAL, as the reader stated.","tokens_in":15900,"tokens_out":6794,"duration_ms":77895,"concrete_test":"Generate the Section IVA high-frequency dataset (Forb = 14.85 nHz, e0 = 0.5, log10 M = 9.2, sigma_pd = 20%) with an independent eccentric-binary waveform implementation, for example the quasi-Keplerian 3PN formalism or the publicly available PTA eccentric waveform of Susobhanan et al. (2020), keeping the same injected parameters, and then run the authors' post-leading-order [9] recovery pipeline on those residuals. If the injected m1 and m2 fall inside the 90% credible intervals of the recovered posteriors, the concern is resolved; if either falls outside, the mass-measurement claim is a self-consistency artifact of using the same model for injection and recovery.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result of Section IVA is that the symmetric mass ratio nu, and hence the individual masses m1 and m2, is constrained by periapsis precession, harmonic structure, and the Earth-term-versus-pulsar-term evolution over the light-travel time tau_alpha. The simulated datasets are generated with the same post-leading-order model [9] that is used for the recovery. This makes the pipeline internally consistent, but it cannot test whether the evolution equations for xi, gamma, x, and e over tau_alpha, or the harmonic amplitudes in Eqs. (4)-(5), are physically faithful. The strongest part of the claim, unbiased individual masses at high frequency when the pulsar term and high-order PN terms are included, is exactly the part that would be biased by an incomplete or approximate dynamical model. The manuscript provides no independent validation of [9]: no comparison with a separately implemented eccentric waveform, no public code release, and no cross-check against established PN results. The limitation subsection discusses CRN and Hellings-Downs approximations and computational cost, but it does not list waveform fidelity as a limitation. Therefore the headline claim is conditional on [9] being correct. If [9] contains an unmodeled 2PN/3PN term, a spin effect, or an integration error over tau_alpha, the recovered masses could fall outside the quoted 90% credible intervals while the posteriors still appear well behaved.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents Bayesian injection-recovery studies for continuous gravitational-wave signals from eccentric supermassive-black-hole binaries in simulated PTA datasets. The signal model is the authors' earlier post-leading-order eccentric waveform [9], and the main scientific claims are (i) individual component masses can be measured when the source is at high orbital frequency, provided both the pulsar term and post-leading-order PN dynamics are included; (ii) dropping the pulsar term or truncating the dynamics to leading order degrades or biases the recovered parameters; and (iii) at low frequencies an eccentric CGW can be partially absorbed by a common red-noise process, giving a bimodal posterior and a Bayes factor of order unity. The paper also reports a GPU-compatible implementation and downsampling tests.","tokens_in":16063,"tokens_out":7894,"duration_ms":80716,"significance":"Within the self-consistency of the adopted waveform, this is a substantial and useful parameter-estimation study. The demonstration that the pulsar term and high-order PN evolution can break the mass-ratio degeneracy is physically well motivated, and the honest reporting of near-unity Bayes factors and the prior-volume penalty is a strength. The CRN-eCGW degeneracy and its impact on model selection are important for real PTA analyses. However, the physical significance of the headline mass-measurement claim is currently conditional on the fidelity of the authors' own waveform model [9], which is not independently validated in this manuscript.","major_comments":[{"comment":"The headline claim that the symmetric mass ratio ν and hence m1 and m2 are measurable is established only within the authors' own post-leading-order eccentric waveform model. The simulated residuals are generated with the same dynamical equations and harmonic amplitudes used for recovery, so the injection-recovery tests demonstrate pipeline self-consistency but not the physical fidelity of the waveform. In particular, the leading-order-model biases shown in Figs. 2-3 are measured against the same higher-order model used to make the injections, not against an independent truth. The limitations subsection (Sec. V.c) does not list waveform fidelity as a limitation. To support the central claim, the authors should either cross-check [9] against an independent eccentric waveform implementation, for example another PN code or established 2PN/3PN results, and examine sensitivity to unmodeled effects such as spins, or explicitly reframe the result as a pipeline demonstration conditional on [9].","section":"Sec. IVA, Eqs. (4)-(5), ref. [9]"},{"comment":"The limitation paragraph states 'we fix the CRN and WN parameters,' but the low-frequency CRN+CGW analysis in Sec. IVB explicitly says 'We infer the parameters of the eccentric SMBHB together with the CRN parameters,' and Fig. 10 presents posterior distributions for log10 A_CRN and γ_CRN. This contradiction is not cosmetic: whether the CRN is fitted or fixed changes the meaning of the bimodal likelihood, the CRN-only comparison, and the reported Bayes factor. Please correct the limitation statement and make clear in Table II which parameters were held fixed in each run.","section":"Sec. V.c vs. Sec. IVB, Figs. 7, 10"}],"minor_comments":[{"comment":"The low-frequency, near-circular injection is given as Forb = 4 nHz and e0 = 0.003 in the text and Fig. 5, but as Forb = 5 nHz and e0 = 0.01 in Table I and in Sec. V. Please harmonize these values.","section":"Table I and Sec. IVA/Conclusion"},{"comment":"The text refers to 'the EPTA DR2new dataset [9]', but reference [9] is the authors' waveform paper; the EPTA data papers should be cited instead.","section":"Sec. IVB"},{"comment":"The Fig. 11 caption describes the free spectrum for the 'high frequency' dataset, but this figure appears in the low-frequency subsection and Fig. 14 is the high-frequency analogue; please check that the captions and table references are not swapped.","section":"Figs. 11 and 14"},{"comment":"There are several typographical and consistency issues that should be cleaned up: 'correponding', 'obtaine', 'Shanon', 'One the other hand', and 'Forb = 15,10nHz' in Sec. V.","section":"Throughout"},{"comment":"The comparison of posteriors with and without the pulsar term would be easier to assess if quantitative summaries, such as credible-interval widths or prior-posterior distances, were reported in addition to the corner plots.","section":"Sec. IVA, Fig. 6"},{"comment":"The abstract says the paper focuses on the 'detection strategy' for individual binaries, while Sec. I states that detectability is not considered and will be the subject of a separate publication; consider aligning the wording.","section":"Abstract and Sec. I"}],"recommendation":"major_revision","confidential_remarks":"The main concern is the absence of any independent validation of the waveform model that underpins the headline claim. This is not a question of internal logic, which is sound, but of physical risk: if the waveform of [9] is incomplete, the recovered masses could be biased while the posteriors still appear well behaved. If the authors cannot provide such validation, they should at minimum soften the claim to 'within our model.' The paper otherwise fits the journal's scope and the work is potentially important for PTA analyses."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is a careful injection-recovery study that delivers a concrete result: for an eccentric SMBHB at high orbital frequency, individual component masses can be measured only when both the pulsar term and post-leading-order PN dynamics are included. That is genuinely new, since circular-orbit CGW analyses do not constrain the mass ratio. Second, the paper documents a low-frequency degeneracy between an eccentric CGW and a common red noise process, with an honest Bayes factor of about 0.7, so the field gets a useful warning about model selection in that regime.\n\nWhat the paper does well: the simulated-data framework is internally consistent, the limitations are stated in the text, and the diagnostic work with free-spectrum fits and Hellinger/JS distances is sensible. The GPU-based implementation is a practical contribution. I credit the authors for reporting the bimodal likelihood rather than claiming a detection.\n\nThe soft spots are real but not fatal for the claims as scoped. The mass-measurement result depends entirely on the authors' own eccentric waveform [9] for both injection and recovery, so no amount of internal consistency can expose systematics in that model. The paper does not cross-check against an independent implementation, and no code or data is released, so the quantitative posteriors cannot be reproduced. The discussion also does not list waveform fidelity as a limitation, which I would have expected. The CRN approximation is acknowledged and the authors are upfront that a Hellings-Downs background could make degeneracies stronger. Minor point: recovery accuracy is described qualitatively; numerical errors on the recovered masses would help.\n\nThe central argument holds for what it claims: given an accurate waveform model, the pipeline recovers the masses. The missing validation of [9] is the load-bearing caveat, and the authors should say so explicitly and ideally provide an independent cross-check or public code.\n\nWho this is for: PTA data analysts and anyone doing CGW parameter estimation. It deserves a serious referee.","headline":"A careful injection-recovery study that credibly demonstrates individual mass measurement for high-frequency eccentric SMBHBs via the pulsar term and PN dynamics, with an honest low-frequency CRN degeneracy; the main caveat is that the waveform model is self-validated only.","tokens_in":16713,"tokens_out":1962,"would_cite":true,"duration_ms":20593,"reading_group":"yes","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.Gb"],"model":"deepseek-v4-flash","headline":"Eccentric supermassive-black-hole binaries at high GW frequencies can have both masses measured by pulsar timing arrays, via precession and harmonics, given the pulsar term and higher-order post-Newtonian dynamics.","keywords":["pulsar timing arrays","eccentric supermassive black hole binaries","continuous gravitational waves","parameter estimation","Bayesian inference","post-Newtonian dynamics","pulsar term","common red noise"],"falsifier":"Inject a signal produced by an independent, independently implemented eccentric waveform code, with the same source parameters and a comparable pulsar term, and recover it with this paper's model; if the recovered $m_1$ and $m_2$ fall outside the $90\\%$ credible interval, the mass-measurement claim is wrong. A real high-frequency eccentric binary detection with individual masses measured independently, for example through electromagnetic observations, would settle the question directly.","tokens_in":15563,"feed_emoji":"🔭","tokens_out":13377,"duration_ms":114012,"temperature":0.7,"pith_summary":"The paper aims to show that pulsar timing arrays can do more than measure the chirp mass of a supermassive black hole binary: for binaries on eccentric orbits detected at high gravitational-wave frequencies, the individual component masses $m_1$ and $m_2$ can be measured. The reason is that the symmetric mass ratio $\\nu=m_1m_2/(m_1+m_2)^2$ leaves visible traces in periapsis precession and in the relative amplitudes of the waveform's many harmonics, and the binary's frequency and eccentricity evolution between the Earth term and the pulsar term breaks the degeneracy that plagues circular templates. These traces are only usable when the model includes both the pulsar term and post-leading-order post-Newtonian corrections; dropping either one leaves the masses unconstrained or biased. A separate result concerns detectability: at low orbital frequencies the eccentric signal's harmonics overlap the common red noise that mimics the stochastic gravitational-wave background, producing a bimodal posterior and a Bayes factor near unity even though the data are informative about the signal parameters. A sympathetic reader would care because this is a concrete route to measuring the masses of individual supermassive black holes, not just their chirp mass, from nanohertz observations.","feed_headline":"Pulsar timing arrays can weigh both black holes in eccentric binaries","feed_subtitle":"Periapsis precession plus the pulsar term pin down both masses of high-frequency eccentric binaries.","key_machinery":"The load-bearing mechanism is the difference between the Earth term and the pulsar term in the timing residuals of equation (1). The pulsar term samples the same binary at an earlier orbital epoch, separated by the light-travel time $\\tau_\\alpha$, so when the binary evolves appreciably over $\\tau_\\alpha$ the data carry frequency and eccentricity evolution information. That information is read through four time-dependent orbital variables: the true anomaly $\\xi(t)$, the periapsis advance $\\gamma(t)$, the dimensionless azimuthal frequency $x(t)=(M\\omega_\\phi)^{2/3}$, and the eccentricity $e(t)$. Periapsis precession and the harmonic amplitudes $a(e,\\xi)$, $b(e,\\xi)$, $c(e,\\xi)$ depend on the symmetric mass ratio $\\nu$ in ways that break the total-mass/mass-ratio degeneracy, but only if the post-leading-order post-Newtonian evolution equations from the authors' earlier model paper are used to compute the pulsar term; at leading order the predicted evolution is too small and biases the recovered masses.","core_discovery":"On the paper's own terms, the central discovery is that the degeneracy between total mass and mass ratio in a continuous-wave search is broken by eccentricity. For a high-frequency ($F_{\\rm orb}\\approx 15$ nHz) binary with initial eccentricity $e_0=0.5$, simulated six-pulsar data with white noise only yield informative posterior distributions for the symmetric mass ratio $\\nu$, and therefore for $m_1$ and $m_2$, when the waveform includes both Earth and pulsar terms and post-leading-order post-Newtonian dynamics. The same data fitted with leading-order dynamics produce broader, biased posteriors, and with pulsar-distance uncertainty $\\sigma_{\\rm pdist}=5\\%$ the injected $m_2$ falls outside the $90\\%$ credible interval. Fitting with the Earth term only leaves all eccentric-CGW parameters unconstrained at high frequency, because the pulsar term is what carries the frequency-evolution information. At low frequency and low eccentricity ($F_{\\rm orb}\\approx 4$ nHz, $e_0\\approx 0.003$), the binary hardly evolves over the pulsar light-travel time, so the Earth-term-only model still detects the signal but cannot constrain mass or distance. The paper also establishes a correlation between the eccentric continuous-wave signal and common red noise at low frequencies: a highly eccentric $3$ nHz binary produces a broad harmonic spectrum that the common red noise can partially absorb, yielding a bimodal likelihood and a Bayes factor $B\\approx 0.7$ against the signal model, even though per-parameter prior-posterior distance metrics show that the data are informative.","pith_inferences":["If pulsar timing arrays detect a handful of high-frequency eccentric binaries, the individual-mass measurements could yield a direct dynamical census of supermassive black hole masses that bypasses scaling relations used in population studies.","The bimodal likelihood and near-unity Bayes factor imply that global model selection may miss low-frequency eccentric continuous-wave candidates; targeted searches with external priors, for example from galaxy catalogues, could break the degeneracy with the background.","The mass-measurement claim is only as strong as the orbital-evolution model used to generate the pulsar term; an independent cross-validation of that model against a differently implemented waveform, or against a real detection with an independent mass measurement, would be the natural next test.","The GPU-accelerated likelihood makes joint background-plus-source searches at full array scale computationally feasible, so the same methodology can be applied to the next generation of pulsar timing data sets."],"forward_implications":["A detected eccentric supermassive black hole binary at tens of nanohertz yields individual masses, not just a chirp mass, turning pulsar timing arrays into instruments for dynamical mass measurement of supermassive black holes.","The pulsar term is not optional at high frequency: dropping it halves the matched-filter SNR from 10.52 to 5.7 for the shown eccentric source and leaves parameters unconstrained, so eccentric-search pipelines must include it.","Using leading-order orbital dynamics for high-frequency sources biases $m_1$ and $m_2$ outside the $90\\%$ credible interval when pulsar distances are known to $5\\%$, so higher-order post-Newtonian evolution is required for unbiased masses.","Tighter pulsar-distance priors, from $20\\%$ to $5\\%$, mainly sharpen sky localisation rather than the binary parameters, identifying distance knowledge as the current limit on source position.","Low-frequency eccentric sources can be partially absorbed by the common red noise; a Bayes factor near unity does not mean the signal is absent, and per-parameter prior-posterior distances are a useful complementary diagnostic."],"supporting_citations":[{"why":"Supplies the eccentric binary waveform model and the post-leading-order post-Newtonian evolution equations on which the mass-measurement claim rests.","marker":"[9]"},{"why":"Supplies the simulated astrophysical supermassive-black-hole-binary population from which the injected source parameters are drawn.","marker":"[8]"},{"why":"Provides the 25-pulsar array configuration, observation epochs, and noise characteristics used to build the realistic simulated datasets.","marker":"[1]"},{"why":"Supplies the GPU-accelerated likelihood implementation used in the Bayesian inference.","marker":"[10]"},{"why":"Supplies the parallel-tempered ensemble Markov-chain Monte Carlo sampler used to draw posterior samples.","marker":"[11]"},{"why":"Selects the six-pulsar subset used in the idealized white-noise dataset for the mass-measurement study.","marker":"[26]"},{"why":"Provides the pulsar selection and noise treatment used in the same idealized dataset.","marker":"[27]"},{"why":"Supplies the product-space method used to compute the Bayes factor between the common-red-noise-only and common-red-noise-plus-signal models.","marker":"[29]"}],"fun_headline_variants":["Eccentric binaries allow PTA to weigh both black holes individually","Pulsar term is key to resolving masses in eccentric SMBH binaries","Eccentricity breaks mass degeneracy in pulsar timing array searches","PTAs can measure both masses of eccentric supermassive black hole binaries"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the post-leading-order eccentric binary evolution model from the authors' earlier paper faithfully describes the real orbital dynamics over the pulsar light-travel time; the mass measurement is extracted from that model's predicted frequency and eccentricity evolution, and the same model both generates and fits the simulated data.","fun_headline_variants_meta":{"raw":{"variants":["Eccentric binaries allow PTA to weigh both black holes individually","Pulsar term is key to resolving masses in eccentric SMBH binaries","Eccentricity breaks mass degeneracy in pulsar timing array searches","PTAs can measure both masses of eccentric supermassive black hole binaries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000485,"raw_usage":{"total_tokens":2448,"prompt_tokens":1056,"completion_tokens":1392,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":672,"completion_tokens_details":{"reasoning_tokens":1313}},"tokens_in":672,"tokens_out":1392,"duration_ms":11083,"temperature":1.0,"reasoning_tokens":1313,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:52:37.229351+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Inject a signal produced by an independent, independently implemented eccentric waveform code, with the same source parameters and a comparable pulsar term, and recover it with this paper's model; if the recovered $m_1$ and $m_2$ fall outside the $90\\%$ credible interval, the mass-measurement claim is wrong. A real high-frequency eccentric binary detection with individual masses measured independently, for example through electromagnetic observations, would settle the question directly.","supporting_citations":[{"cited_title":"curse of dimensionality","cited_arxiv_id":null,"evidence_quote":"Supplies the simulated astrophysical supermassive-black-hole-binary population from which the injected source parameters are drawn."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the 25-pulsar array configuration, observation epochs, and noise characteristics used to build the realistic simulated datasets."},{"cited_title":"Quality over Quantity: Optimizing pulsar timing array analysis for stochastic and continuous gravitational wave signals","cited_arxiv_id":"2211.03201","evidence_quote":"Provides the pulsar selection and noise treatment used in the same idealized dataset."},{"cited_title":"Bartolucci, L","cited_arxiv_id":null,"evidence_quote":"Supplies the product-space method used to compute the Bayes factor between the common-red-noise-only and common-red-noise-plus-signal models."}],"review_version":1}