{"id":"9bc295a6-53ee-4678-80e5-932a5976dc7a","arxiv_id":"2608.07095","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper derives analytic anisotropic overlap reduction functions for all six gravitational-wave polarizations and all spherical-harmonic orders with the full pulsar term, recovering the Hellings-Downs curve in the isotropic limit.","lead":"This paper derives closed-form mathematical expressions for the angular correlation response of pulsar timing arrays to anisotropic gravitational-wave backgrounds, covering all six polarization modes. If correct, it removes a standard approximation that fails for scalar longitudinal waves and for closely spaced pulsars, supporting future precision mapping of the gravitational-wave sky.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central deliverable — complete analytic ℓ≤5 ORFs for all six modes — is delegated to an unreviewed GitHub repository, and the manuscript shows no numerical verification of the claimed exact formulas; until the repository is checked against Eq. (16), the headline claim is unverified.","rationale":"The reader's weakest assumption is the response function in Eq. (13). I agree that a wrong response model would invalidate every ORF expression, but that premise is standard in the PTA literature and is less immediately checkable from the manuscript alone. The more acute and testable weakness is that the paper's own central artifact—the exact ℓ≤5 formulas—is not present in the manuscript. The paper explicitly delegates these expressions to an external GitHub repository and does not show a numerical comparison with Eq. (16). The key transformation Eq. (34) also leaves the branch structure of the two-to-one mapping from the sphere to the (x̃,ỹ) ellipse unstated, so a reader cannot confirm that the repository formulas implement the defining integral rather than a related but differently normalized one. The symmetry theorems and limiting forms that recover known isotropic results are useful independent support, and they make the framework plausible, but they do not prove the ℓ≤5 formulas. A direct quadrature comparison settles the question in either direction; if the repository formulas agree with Eq. (16), the central claim is substantially supported, and if they do not, the claim fails. The reader's CONDITIONAL verdict remains the right verdict pending that check, so no change is needed.","tokens_in":28036,"tokens_out":16095,"duration_ms":148811,"concrete_test":"Download the repository cited as Ref. [36]; for each of the six polarizations and for a grid of (ℓ,m) with ℓ≤5 (e.g., T: (3,1),(3,3); V: (4,0),(4,2); B: (3,1); L: (5,1),(5,5)), evaluate the provided closed-form formulas at parameter triples covering au≠av (e.g., au=3, av=10) and γ∈{0.1,0.5,1.5,3.0}; compare with direct adaptive numerical quadrature of Eq. (16) using the same response Eq. (13), aiming for relative agreement below 1e-6 away from γ=0 and γ=π. Also check that the tensor (0,0) term is consistent with Eq. (42) and the Hellings-Downs curve in the appropriate limit. If any discrepancy exceeds numerical tolerance, the claimed exact analytic expressions are not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract and Sec. I claim a universal, approximation-free framework yielding fully analytic anisotropic ORFs to arbitrary order, with complete expressions up to ℓ≤5. In the manuscript, however, the actual closed forms are not printed: Sec. III states “Here we present the analytical forms of all ORFs for ℓ≤5. For the detailed expressions of these formulas, please see Ref. [36],” where [36] is a GitHub repository. Only ℓ≤2 short-wavelength limits and autocorrelations appear in the text. The claimed agreement with numerical integration is asserted but never shown—the figures compare analytic curves against the short-wavelength approximation, not against numerical quadrature. This matters because Eq. (34), the only bridge from the defining angular integral Eq. (16) to the closed forms, is not derived: the Jacobian 1/(sinγ cosθ) requires handling the two-to-one map from the sphere to the (x̃,ỹ) ellipse and the branch of cosθ, and the paper does not state how this is done. If the repository formulas are absent, erroneous, or do not reproduce direct numerical integration of Eq. (16) at the same (au, av, γ), then the central claim of exact analytic expressions up to ℓ≤5 fails as written. The conditional verdict is therefore appropriate until this check is performed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a universal analytical framework for the overlap reduction functions (ORFs) of a stochastic gravitational-wave background with pulsar timing arrays, covering all six polarization modes and arbitrary spherical-harmonic order. The authors choose a coordinate system with the two pulsars symmetric about the z-axis, introduce the variables x̃ and ỹ as projections of the sky direction onto the two pulsar directions, and state that the ORF integral on the sphere transforms into an integral over an elliptical domain with a Jacobian 1/(sin γ cos θ). They further state that this integral can be evaluated analytically. The paper presents four symmetry theorems for the anisotropic ORFs, short-wavelength limits for ℓ ≤ 2, autocorrelation limits, and plots of the claimed analytical results for tensor, vector, scalar-breathing, and scalar-longitudinal modes. The complete analytical expressions up to ℓ ≤ 5, however, are not printed in the manuscript; the reader is referred to a GitHub repository [36]. The paper repeatedly asserts agreement with numerical integration, but no numerical quadrature comparison is shown; the figures compare the analytical curves only with the short-wavelength approximation.","tokens_in":28325,"tokens_out":3290,"duration_ms":32700,"significance":"If the claimed formulas are correct, the paper would supply exact, closed-form angular correlation kernels for anisotropic PTA searches, removing the need for numerical integration over the sky in each pair evaluation. The symmetry results, especially the vanishing of ORFs for ℓ + m odd and the exchange symmetry under pulsar exchange, are plausible, non-trivial, and could be useful in reducing computational cost and in code validation. The paper also correctly emphasizes the failure of the short-wavelength approximation for the scalar longitudinal mode at small angular separations. However, the central deliverable — the complete ℓ ≤ 5 analytical expressions — is absent from the manuscript, and the claimed agreement with numerical integration is not demonstrated. As written, the headline result is therefore unverified, and the paper cannot be accepted in its current form.","major_comments":[{"comment":"The paper states: \"Here we present the analytical forms of all ORFs for ℓ≤5. For the detailed expressions of these formulas, please see Ref. [36],\" where Ref. [36] is a GitHub repository. This means the central claim of the paper — complete analytic expressions up to ℓ ≤ 5 for all six modes — is not actually contained in the manuscript. A journal publication must be self-contained for its main result; the full formulas need to be included in the paper or in a published appendix/supplement, with the repository serving only as auxiliary code. As it stands, the claimed result cannot be checked from the paper itself.","section":"Sec. III, before Sec. III A"},{"comment":"The crucial step \"This integral can be evaluated analytically\" is asserted without derivation. Equation (34) involves a change of variables from the sphere to the elliptical domain E, with Jacobian 1/(sin γ cos θ). The transformation from (θ, φ) to (x̃, ỹ) is two-to-one, since the sphere is covered twice when cos θ changes sign, and the branch of cos θ in the Jacobian is not specified. The text does not explain how the integration domain is resolved, how the two-to-one map is handled, or how the singular behavior at cos θ = 0 is treated. This missing derivation is load-bearing because all closed-form formulas in the paper and repository rest on the analytic evaluation of Eq. (34). The authors should provide a step-by-step derivation of the reduction, including the domain decomposition and branch handling.","section":"Sec. III, Eq. (34)"},{"comment":"The manuscript repeatedly claims that the analytical expressions are \"in excellent agreement with numerical integration,\" but no comparison with numerical quadrature of the defining integral Eq. (16) is shown. All plotted comparisons are between the claimed analytical curves and the short-wavelength approximation, which is a different object. This is a particular concern because the formulas are delegated to a repository and the analytic integration step is not demonstrated. The authors should add explicit residual plots or tables comparing their closed-form expressions with direct numerical integration of Eq. (16) over the same (a_u, a_v, γ) values, for each polarization and for a representative set of (ℓ, m) up to ℓ = 5. Without such a check, the central claim of exact analytical expressions is not supported by the evidence in the paper.","section":"Secs. III A–D, Figs. 2–13"}],"minor_comments":[{"comment":"The notation R_ij(a_u, a_v, x̃, ỹ) ⊙ Ỹ_ℓm(x̃, ỹ) is introduced without definition. Please define the \"⊙\" operation explicitly, including its normalization and how it relates to the integral in Eq. (16).","section":"Sec. III, Eq. (34)"},{"comment":"The response function in Eq. (13) is taken from Ref. [31] with the choice r2 = 0, but the text later refers to ⃗r2 = L2 r̂2 and L2. Please clarify the notation and the sense in which Earth's position can be set to zero, since the phase factor e^{i 2π f n̂·r2/c} is retained.","section":"Sec. II B, Eq. (13)"},{"comment":"The subscript \"sw\" is used for short-wavelength expressions without a definition at first use. Please define the short-wavelength approximation and state clearly where it is being applied.","section":"Sec. III"},{"comment":"Equation (4) writes δ^2(k̂, k̂′) while Eq. (5) writes δ^(2)(k̂ − k̂′); these notations are inconsistent. Please use a single convention for the two-dimensional delta function on the sphere.","section":"Sec. II A, Eqs. (4) and (5)"},{"comment":"Reference [36] is incomplete and appears as \"(????)\". If the GitHub repository is to be cited, please provide a full citation with a version, a DOI, or an archival link (for example, Zenodo) so that the exact version of the formulas can be identified.","section":"Bibliography"},{"comment":"Several figure captions contain malformed text such as \"au = av = 10, = 0,m = 0\" and \"sw,m = 1\", with missing symbols and inconsistent spacing. These should be corrected throughout Figs. 2–13.","section":"Figure captions"},{"comment":"In the discussion of vector modes, the paper states that for the full expression \"larger a_v is needed for better agreement with the short-wavelength approximation.\" This is qualitative; please provide a quantitative convergence criterion, for example the value of a_v at which the difference falls below a stated tolerance.","section":"Sec. III B"}],"recommendation":"major_revision","confidential_remarks":"The paper would benefit from a clear statement of its novelty relative to the existing analytical isotropic results, such as Hu et al. (2022) and the spherical-harmonic expansions of Gair et al. (2014). In particular, the phrase \"for the first time\" in the introduction should be made precise. If the full formulas are placed in an appendix and validated against numerical integration, the paper could become a useful reference for anisotropic PTA analyses; in its current form, the main result is external to the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the combination: anisotropic sky distribution, full pulsar term, and all six polarization modes in one analytical ORF framework. Hu et al. did the isotropic full-pulsar-term case, and Anil Kumar & Kamionkowski and Inomata et al. did the anisotropic case but dropped the pulsar term. This paper fills that gap, and the symmetry theorems in Sec. II C look correct and genuinely useful — the vanishing of ℓ+m odd, the exchange symmetry in au and av, and the Wigner-D rotation property all check out in the appendix. The printed limiting forms recover the Hellings–Downs curve and match known short-wavelength results for tensor, vector, breathing, and longitudinal modes, which is real evidence the framework is right.\n\nThe soft spot is big and it is the whole ballgame: the actual ℓ≤5 formulas are not in the paper. Sec. III says \"For the detailed expressions of these formulas, please see Ref. [36],\" and Ref. [36] is a GitHub repository. That means the central deliverable is unreviewed and unprintable in the form submitted. A referee cannot verify the claim without auditing an external repo that is not fixed in time. The derivation step in Eq. (34) is also asserted rather than shown — the Jacobian 1/(sinγ cosθ) requires handling the two-to-one map from the sphere to the (x̃,ỹ) ellipse and the branch of cosθ, and the paper never states how that was done. On top of that, the claimed agreement with numerical integration is not demonstrated anywhere: the figures compare the analytic curves against the short-wavelength approximation, not against direct quadrature of Eq. (16). The abstract says the framework is \"free of approximations\" and the expressions are \"fully analytical\" to arbitrary order, but the only order actually presented is ℓ≤2 in the short-wavelength limit plus autocorrelations.\n\nNone of this makes the result wrong. The derivation is self-contained up to the hard integral, there is no fitting to data, and the limiting cases are consistent with external benchmarks. I would bet the repository formulas are correct. But as written, the paper's main claim is unverified.\n\nSend it to peer review — this deserves referee time. The referee should check the GitHub repo against direct numerical integration of Eq. (16) at representative (au, av, γ), and the formula list should be required to appear in the paper or in a supplement that is part of the review. If the formulas check out, this is a solid methodological advance for the PTA community.","headline":"Useful if the GitHub formulas hold up, but the paper as written delegates its main result to an unreviewed repository and shows no direct numerical check.","tokens_in":28836,"tokens_out":1857,"would_cite":false,"duration_ms":19245,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives exact closed-form overlap reduction functions for all six gravitational-wave polarization modes and arbitrary spherical-harmonic order, with no short-wavelength approximation, and shows their monopole component…","keywords":["pulsar timing arrays","stochastic gravitational-wave background","overlap reduction functions","anisotropy","polarization modes","Hellings-Downs curve","scalar longitudinal mode","analytic expressions"],"falsifier":"Evaluate the paper's closed form for a specific low-order ORF, say $\\Gamma^{T,0,0}(a_u,a_v,\\gamma)$ with $a_u=a_v=10$ at $\\gamma=1$ rad, and compare it with direct numerical quadrature of Eq. (18) using the response function of Eq. (13) to machine precision; agreement at the level of $10^{-12}$ confirms the reduction, while any systematic discrepancy—or any nonzero recovered value for an $\\ell+m$ odd multipole—shows the analytic mapping of the integral is incorrect.","tokens_in":27829,"feed_emoji":"📡","tokens_out":11647,"duration_ms":94078,"temperature":0.7,"pith_summary":"This paper claims a universal analytical framework for the overlap reduction functions (ORFs) that pulsar timing arrays use to search for an anisotropic stochastic gravitational-wave background. ORFs are the angular-correlation kernels telling how much correlated timing signal any pair of pulsars sees from a given sky distribution of gravitational waves; previously they were available analytically only under the short-wavelength approximation or for isotropic skies. The authors derive closed-form expressions, valid for arbitrary multipole order $\\ell m$ and for all six polarization modes (tensor, vector, breathing, longitudinal), that keep the full frequency and pulsar-distance dependence, and they show the monopole component recovers the Hellings-Downs curve. The same derivation removes the divergences that plagued vector and scalar-longitudinal modes in the short-wavelength approximation. If correct, this turns a numerical bottleneck into instant evaluation, which matters for full-sky anisotropy maps, polarization separation, and tests of modified gravity with current and next-generation pulsar arrays.","feed_headline":"Exact PTA correlation kernels for all six GW modes","feed_subtitle":"Closed-form overlap reduction functions replace numerical integration and tame the modes that used to diverge.","key_machinery":"The central object is the anisotropic overlap reduction function $\\Gamma^{A,\\ell m}_{ij}(a_u,a_v,\\gamma)$, the projection in Eq. (18) of the product of two single-pulsar Doppler response functions onto the sky harmonic $\\tilde{Y}_{\\ell m}$. Three moves working together carry the argument: a change of integration variables from the sky angles to the projected cosines $\\tilde{x}=\\hat{n}\\cdot\\hat{u}$ and $\\tilde{y}=\\hat{n}\\cdot\\hat{v}$, which maps the unit sphere to an elliptical domain and turns the response product into exponentials times low-order polynomials in $\\tilde{x}$ and $\\tilde{y}$; a set of symmetry theorems (Wigner-$D$ rotation behaviour, exact vanishing when $\\ell+m$ is odd, exchange symmetry in the two distance parameters $a_u$ and $a_v$, and the $\\ell,-m$ reflection identity) that cut the number of independent integrals by roughly three quarters; and evaluation of the surviving integrals in terms of elementary functions plus the exponential integral $\\mathrm{E}_1$ and the sine and cosine integrals $\\mathrm{Si}$ and $\\mathrm{Ci}$.","core_discovery":"On the paper's own terms, the core discovery is that the anisotropic overlap-reduction-function integral—$\\Gamma^{A,\\ell m}_{ij}(a_u,a_v,\\gamma)$, the spherical-harmonic projection of the two-pulsar response product—can be evaluated in closed form for every polarization $A$ and every harmonic $(\\ell,m)$, with no short-wavelength approximation. By changing variables from the sky angles to the two dot products $\\tilde{x}=\\hat{n}\\cdot\\hat{u}$ and $\\tilde{y}=\\hat{n}\\cdot\\hat{v}$, the authors convert the spherical integral into one over an elliptical domain in which the response functions reduce to exponentials times low-order polynomials in $\\tilde{x}$ and $\\tilde{y}$. Using four symmetry theorems—the Wigner-$D$ rotation rule, the vanishing of all $\\ell+m$ odd terms, the exchange symmetry between the two pulsar-distance parameters, and the reflection relation for $m\\leftrightarrow -m$—they reduce the surviving integrals to expressions in elementary functions and the special functions $\\mathrm{E}_1$, $\\mathrm{Si}$, $\\mathrm{Ci}$. The framework's $(0,0)$ component reproduces the Hellings-Downs curve for tensor modes and the known isotropic results for vector, breathing, and longitudinal modes; unlike previous treatments, the vector and scalar-longitudinal ORFs are finite everywhere, including the coincident-pulsar limit. Explicit formulas are given through $\\ell\\leq 5$, with the algorithm extendable to arbitrary order.","pith_inferences":["Beyond the paper: the same change-of-variables and symmetry machinery should transfer to other correlated detectors—space-based interferometers or astrometric surveys—whose ORF integrals share the same spherical-harmonic skeleton, though the paper demonstrates it only for pulsar pairs.","Beyond the paper: the exact imaginary parts of the ORFs for unequal pulsar distances carry information that real-valued Hellings-Downs analyses ignore; a frequency-dependent statistic built on the full complex ORF may detect anisotropies at lower signal-to-noise ratio, a gain the paper does not quantify.","Beyond the paper: the parity rule that $\\ell+m$ odd gives zero offers a cheap null test—a real PTA pipeline that recovers a statistically significant $\\ell+m$ odd coefficient would be flagging a systematics or model error rather than a gravitational-wave signal."],"forward_implications":["Any pulsar pair's ORF at any multipole and any polarization becomes a single closed-form evaluation, so anisotropic sky mapping and polarization-separation pipelines no longer need per-pair numerical quadrature over oscillatory integrands.","Existing isotropic analyses are exactly recovered: the $(0,0)$ tensor component gives the Hellings-Downs curve, and the $(0,0)$ components for the other modes reduce to the published isotropic ORFs.","Because $\\ell+m$ odd coefficients vanish identically and $m<0$ coefficients are fixed by symmetry, roughly three quarters of the naive harmonic grid is forced to zero or related, shrinking the parameter space of anisotropy searches.","The formerly divergent vector and scalar-longitudinal ORFs become finite and frequency-dependent; the model predicts the longitudinal autocorrelation grows linearly with frequency and the vector autocorrelation logarithmically, which is testable once pulsar-distance information is included."],"supporting_citations":[{"why":"Defines the Hellings-Downs curve, the isotropic tensor ORF that the framework's (0,0) component must reproduce.","marker":"[7]"},{"why":"Provides earlier isotropic ORFs for tensor, vector, and scalar modes whose short-wavelength limits the new expressions recover.","marker":"[22]"},{"why":"Identified the divergence of the scalar-longitudinal ORF under the short-wavelength approximation, the problem the full pulsar-term treatment resolves.","marker":"[23]"},{"why":"Introduced the spherical-harmonic expansion of PTA anisotropy that the paper's Γ^{ℓm} framework builds on.","marker":"[25]"},{"why":"Gave complete analytical isotropic ORFs for all six modes without the short-wavelength approximation, which this work extends to anisotropic multipoles.","marker":"[26]"},{"why":"Established a spherical-harmonic ORF framework that neglects the pulsar term, the approximation the paper removes.","marker":"[27]"},{"why":"Supplies the Doppler-frequency response function in Eq. (13), the load-bearing starting point for every ORF integral in the paper.","marker":"[31]"},{"why":"Hosts the detailed analytic expressions for ℓ≤5 that the print version only summarizes.","marker":"[36]"}],"fun_headline_variants":["Exact ORFs for all six GW modes in PTA","Closed-form anisotropy ORFs for PTAs","PTA ORFs: exact analytic for all polarizations","All-sky PTA correlation functions, exact","Taming PTA ORFs: exact across all six modes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the standard Doppler-frequency response model of Eq. (13)—with its pulsar-term phase and the choice $r_2=0$—so if real timing residuals respond differently, every ORF expression inherits the error.","fun_headline_variants_meta":{"raw":{"variants":["Exact ORFs for all six GW modes in PTA","Closed-form anisotropy ORFs for PTAs","PTA ORFs: exact analytic for all polarizations","All-sky PTA correlation functions, exact","Taming PTA ORFs: exact across all six modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000554,"raw_usage":{"total_tokens":2715,"prompt_tokens":1093,"completion_tokens":1622,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":1544}},"tokens_in":709,"tokens_out":1622,"duration_ms":31356,"temperature":1.0,"reasoning_tokens":1544,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:55:42.710344+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the paper's closed form for a specific low-order ORF, say $\\Gamma^{T,0,0}(a_u,a_v,\\gamma)$ with $a_u=a_v=10$ at $\\gamma=1$ rad, and compare it with direct numerical quadrature of Eq. (18) using the response function of Eq. (13) to machine precision; agreement at the level of $10^{-12}$ confirms the reduction, while any systematic discrepancy—or any nonzero recovered value for an $\\ell+m$ odd multipole—shows the analytic mapping of the integral is incorrect.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides earlier isotropic ORFs for tensor, vector, and scalar modes whose short-wavelength limits the new expressions recover."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identified the divergence of the scalar-longitudinal ORF under the short-wavelength approximation, the problem the full pulsar-term treatment resolves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Doppler-frequency response function in Eq. (13), the load-bearing starting point for every ORF integral in the paper."},{"cited_title":"Zhou and J","cited_arxiv_id":null,"evidence_quote":"Hosts the detailed analytic expressions for ℓ≤5 that the print version only summarizes."}],"review_version":1}