Recognition: unknown
Exploring the statistical anisotropy of primordial curvature perturbations with pulsar timing arrays
Pith reviewed 2026-05-09 20:54 UTC · model grok-4.3
The pith
A dipole anisotropy in the primordial curvature power spectrum generates dipolar and quadrupolar patterns in scalar-induced gravitational waves observed by pulsar timing arrays.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We demonstrate that the primordial dipole induces both dipolar and quadrupolar anisotropies in the energy density spectrum of scalar-induced gravitational waves (SIGWs), without generating extra polarization modes. Based on this anisotropic spectrum, we derive the corresponding PTA overlap reduction functions (ORFs), which exhibit frequency dependence, with the anisotropies enhanced on small scales. Furthermore, owing to the non-uniform distribution of millisecond pulsars over the sky in current PTA dataset, the ORFs exhibit a morphology that explicitly depends on the preferred direction of the anisotropy. However, our bayesian analysis of the NANOGrav 15-year dataset still yields no sign
What carries the argument
Dipole-type statistical anisotropy in the primordial power spectrum that induces anisotropic energy density in scalar-induced gravitational waves and direction-dependent overlap reduction functions for pulsar timing arrays.
If this is right
- The overlap reduction functions for PTAs become frequency dependent.
- Anisotropies are enhanced on small scales in the SIGW spectrum.
- The specific form of the ORFs depends on the preferred direction of anisotropy and the distribution of pulsars.
- Broader frequency coverage in future PTA observations will allow tighter constraints on the anisotropy amplitude.
Where Pith is reading between the lines
- This method provides a new probe of early universe directional preferences that could complement cosmic microwave background measurements of statistical anisotropy.
- If a signal is detected at higher frequencies, it could distinguish between different models of inflation that break isotropy.
- The current weak limits suggest that combining PTA data with other gravitational wave observatories might be necessary to fully test the phenomenological model.
Load-bearing premise
The analysis relies on a phenomenological dipole parameterization of anisotropy in the primordial power spectrum and assumes that the PTA observational frequencies lie below the SIGW spectral peak where anisotropic contributions are suppressed.
What would settle it
Detection of the predicted frequency-dependent dipolar and quadrupolar anisotropies in future PTA data at frequencies closer to or above the spectral peak, consistent with an amplitude g greater than 0.5, or the absence of any such signal across an extended frequency range would test the model.
Figures
read the original abstract
The recent detection of a stochastic gravitational wave background by pulsar timing arrays has opened a new window in understanding supermassive black hole binaries and in probing the universe at the early time. Recently, pulsar timing array (PTA) collaborations have been further paving the way to probe anisotropies in the stochastic gravitational wave background. This study investigates dipole-type statistical anisotropy in the primordial power spectrum within a phenomenological framework. We demonstrate that the primordial dipole induces both dipolar and quadrupolar anisotropies in the energy density spectrum of scalar-induced gravitational waves (SIGWs), without generating extra polarization modes. Based on this anisotropic spectrum, we derive the corresponding PTA overlap reduction functions (ORFs), which exhibit frequency dependence, with the anisotropies enhanced on small scales. Furthermore, owing to the non-uniform distribution of millisecond pulsars over the sky in current PTA dataset, the ORFs exhibit a morphology that explicitly depends on the preferred direction of the anisotropy. However, our bayesian analysis of the NANOGrav 15-year dataset still yields no significant evidence for a preferred direction and a weak upper limit on anisotropy amplitude $(g\lesssim0.5)$. This result arises because the observational frequency band lies below the spectral peak, where our models predict suppressed anisotropic contributions. This limitation highlights the potential of future PTA observations. Specifically, datasets with broader frequency coverage are expected to tighten constraints on dipole-type anisotropy.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a phenomenological investigation of dipole-type statistical anisotropy in the primordial curvature power spectrum. It demonstrates that this induces both dipolar and quadrupolar anisotropies in the energy density spectrum of scalar-induced gravitational waves (SIGWs) without generating extra polarization modes, derives the corresponding frequency-dependent PTA overlap reduction functions (ORFs) that depend on the preferred direction due to non-uniform pulsar sky distribution, and performs a Bayesian analysis on the NANOGrav 15-year dataset. The analysis finds no significant evidence for a preferred direction and reports a weak upper limit g ≲ 0.5, which the authors attribute to the PTA frequency band lying below the SIGW spectral peak where anisotropic contributions are suppressed. The work highlights the potential for future PTA datasets with broader frequency coverage to improve constraints.
Significance. If the derivations of the anisotropic SIGW spectra and ORFs are correct, this provides a concrete framework for using PTA data to constrain primordial statistical anisotropies, with the frequency dependence and directional morphology of the ORFs offering testable predictions. The application to real NANOGrav data and the explicit discussion of observational limitations represent a strength, though the reported upper limit is conditional on the choice of scalar spectrum parameters.
major comments (2)
- [Abstract and Bayesian analysis section] Abstract and discussion of Bayesian results: The upper limit g ≲ 0.5 and lack of directional evidence are explicitly attributed to the NANOGrav frequency band lying below the SIGW spectral peak, where the kernel averages directions more isotropically and suppresses anisotropic contributions from the dipole term in P_ζ(k, n̂) = P0(k)[1 + g k̂·n̂]. However, the peak location is set by the free phenomenological parameter k* in the scalar spectrum, which is not independently constrained by the NANOGrav spectrum shape or other observables. Shifting k* so the peak enters the observed band would increase the directional modulation in the frequency-dependent ORFs and alter the posterior on g. The analysis should include robustness checks by varying k* or marginalizing over it to support the interpretation of the limit.
- [Section on ORF derivation] Derivation of anisotropic ORFs: The claim that the ORFs exhibit frequency dependence with anisotropies enhanced on small scales is central to the PTA predictions. The manuscript should provide explicit expressions or numerical examples showing how the directional integrals over the anisotropic P_ζ weighted by the SIGW kernel produce the dipolar/quadrupolar terms in Ω_GW(f) and the resulting ORF morphology, particularly for the non-uniform pulsar distribution.
minor comments (1)
- [Abstract] The abstract refers to 'our models predict suppressed anisotropic contributions' without specifying the exact scalar spectrum form or kernel used; a brief equation or reference in the abstract or introduction would improve clarity.
Simulated Author's Rebuttal
We thank the referee for their careful reading of our manuscript and for the constructive comments. We address each major comment below and agree to strengthen the paper with additional robustness checks and explicit derivations as requested.
read point-by-point responses
-
Referee: [Abstract and Bayesian analysis section] Abstract and discussion of Bayesian results: The upper limit g ≲ 0.5 and lack of directional evidence are explicitly attributed to the NANOGrav frequency band lying below the SIGW spectral peak, where the kernel averages directions more isotropically and suppresses anisotropic contributions from the dipole term in P_ζ(k, n̂) = P0(k)[1 + g k̂·n̂]. However, the peak location is set by the free phenomenological parameter k* in the scalar spectrum, which is not independently constrained by the NANOGrav spectrum shape or other observables. Shifting k* so the peak enters the observed band would increase the directional modulation in the frequency-dependent ORFs and alter the posterior on g. The analysis should include robustness checks by varying k* or marginalizing over it to support the interpretation of the limit.
Authors: We agree that k* is a free phenomenological parameter whose value is not independently fixed by the NANOGrav data, and that our interpretation of the weak upper limit on g depends on the assumed location of the SIGW peak. To address this, we will add robustness checks in the revised manuscript by repeating the Bayesian analysis for several representative values of k* (both inside and outside the observed band) and by exploring marginalization over k* where computationally feasible. These results will be reported in an updated results section and discussion, with the abstract revised to qualify the interpretation accordingly. revision: yes
-
Referee: [Section on ORF derivation] Derivation of anisotropic ORFs: The claim that the ORFs exhibit frequency dependence with anisotropies enhanced on small scales is central to the PTA predictions. The manuscript should provide explicit expressions or numerical examples showing how the directional integrals over the anisotropic P_ζ weighted by the SIGW kernel produce the dipolar/quadrupolar terms in Ω_GW(f) and the resulting ORF morphology, particularly for the non-uniform pulsar distribution.
Authors: The manuscript already contains the derivation of the anisotropic Ω_GW(f) via directional integrals of the dipole-modulated P_ζ(k, n̂) against the SIGW kernel, leading to the dipolar and quadrupolar anisotropies that enter the frequency-dependent ORFs. To make this more transparent, we will insert the explicit integral expressions for the anisotropic contributions to Ω_GW(f) and add numerical examples (including plots) that illustrate the frequency dependence, the enhancement of anisotropies at higher frequencies, and the directional morphology induced by the actual NANOGrav pulsar sky distribution. revision: yes
Circularity Check
No significant circularity; derivation and data analysis remain independent
full rationale
The paper introduces a phenomenological parameterization of dipole anisotropy in the primordial spectrum, derives the resulting dipolar/quadrupolar anisotropies in the SIGW energy density and the associated frequency-dependent ORFs via standard kernel integrals, and then performs a separate Bayesian analysis on the external NANOGrav 15-year dataset to place an upper limit on the free parameter g. None of these steps reduce to their inputs by construction: the derivation follows from the chosen ansatz without tautology, the ORFs are computed outputs rather than renamed fits, and the reported limit g≲0.5 is the direct posterior result from fitting the data under the explicitly stated assumption that the PTA band lies below the scalar spectrum peak. The assumption itself is flagged by the authors as a limitation rather than smuggled in as a derived result. No self-citation chain or uniqueness theorem is invoked to force the central claims. This is a standard model-dependent constraint exercise with no load-bearing reduction to the inputs.
Axiom & Free-Parameter Ledger
free parameters (2)
- g =
≲0.5
- preferred direction
axioms (2)
- domain assumption Standard assumptions for the generation of scalar-induced gravitational waves from primordial curvature perturbations.
- domain assumption The stochastic gravitational wave background in the PTA band is dominated by scalar-induced contributions.
Reference graph
Works this paper leans on
-
[1]
G. Agazieet al.(NANOGrav), Astrophys. J. Lett.951(2023), 10.3847/2041-8213/acdac6, arXiv:2306.16213 [astro-ph.HE]
-
[2]
Antoniadiset al.(EPTA, InPTA:), Astron
J. Antoniadiset al.(EPTA, InPTA:), Astron. Astrophys.678, A50 (2023), arXiv:2306.16214 [astro- ph.HE]
-
[3]
D. J. Reardonet al., Astrophys. J. Lett.951, L6 (2023), arXiv:2306.16215 [astro-ph.HE]
work page internal anchor Pith review arXiv 2023
-
[4]
H. Xuet al., Res. Astron. Astrophys.23, 075024 (2023), arXiv:2306.16216 [astro-ph.HE]
work page internal anchor Pith review arXiv 2023
-
[5]
G. Agazieet al.(NANOGrav), Astrophys. J. Lett.952, L37 (2023), arXiv:2306.16220 [astro-ph.HE]
-
[6]
G. Agazieet al.(NANOGrav), Astrophys. J. Lett.951, L50 (2023), arXiv:2306.16222 [astro-ph.HE]
-
[7]
Agazieet al.(NANOGrav), Astrophys
G. Agazieet al.(NANOGrav), Astrophys. J. Lett.956, L3 (2023), arXiv:2306.16221 [astro-ph.HE]
-
[8]
J. Antoniadiset al.(EPTA, InPTA), Astron. Astrophys.685, A94 (2024), arXiv:2306.16227 [astro- ph.CO]
-
[9]
G. Agazieet al., Astrophys. J. Lett.978, L29 (2025), arXiv:2408.10166 [astro-ph.HE]
-
[10]
Agazieet al., (2026), arXiv:2601.09481 [astro-ph.HE]
G. Agazieet al., (2026), arXiv:2601.09481 [astro-ph.HE]
- [11]
-
[12]
Exploring the spectrum of stochastic gravitational-wave anisotropies with pulsar timing arrays,
G. Sato-Polito and M. Kamionkowski, Phys. Rev. D109, 123544 (2024), arXiv:2305.05690 [astro- ph.CO]
- [13]
- [14]
- [15]
- [16]
-
[17]
S. Kuroyanagi, K. Takahashi, N. Yonemaru, and H. Kumamoto, Phys. Rev. D95, 043531 (2017), arXiv:1604.00332 [astro-ph.CO]
- [18]
- [19]
- [20]
- [21]
- [22]
-
[23]
X. Zhang, J.-Z. Zhou, and Z. Chang, Eur. Phys. J. C82, 781 (2022), arXiv:2208.12948 [astro-ph.CO]
-
[24]
C. Chen and A. Ota, Phys. Rev. D106, 063507 (2022), arXiv:2205.07810 [astro-ph.CO]
-
[25]
Y.-T. Kuang, J.-Z. Zhou, D. Wu, and Z. Chang, (2023), arXiv:2309.06676 [astro-ph.CO]
- [26]
- [27]
-
[28]
D. Mukherjee, H. V. Ragavendra, and S. K. Sethi, Phys. Rev. D113, 023533 (2026), arXiv:2506.23798 [astro-ph.CO]
- [29]
-
[30]
J.-P. Li, S. Wang, Z.-C. Zhao, and K. Kohri, (2025), arXiv:2505.16820 [astro-ph.CO]
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [31]
-
[32]
Vilenkin, Phys
A. Vilenkin, Phys. Lett. B107, 47 (1981)
1981
-
[33]
C. J. Hogan and M. J. Rees, Nature311, 109 (1984)
1984
-
[34]
Vachaspati and A
T. Vachaspati and A. Vilenkin, Phys. Rev. D31, 3052 (1985)
1985
-
[35]
T. Damour and A. Vilenkin, Phys. Rev. D71, 063510 (2005), arXiv:hep-th/0410222
- [36]
-
[37]
C. Ringeval and T. Suyama, JCAP12, 027 (2017), arXiv:1709.03845 [astro-ph.CO]
-
[38]
J. Baeza-Ballesteros, E. J. Copeland, D. G. Figueroa, and J. Lizarraga, Phys. Rev. D110, 043522 (2024), arXiv:2308.08456 [astro-ph.CO]
- [39]
-
[40]
Z.-Y. Qiu and Z.-H. Yu, Chin. Phys. C47, 085104 (2023), arXiv:2304.02506 [hep-ph]
- [41]
-
[42]
Gravitational waves from metastable cosmic strings in the delayed scaling scenario,
Y. Hu and K. Kamada, JCAP04, 044 (2025), arXiv:2501.18380 [astro-ph.CO]
-
[43]
C. J. Hogan, Mon. Not. Roy. Astron. Soc.218, 629 (1986)
1986
- [44]
-
[45]
A. Chatrchyan, M. C. D. Marsh, and C. Nikolis, Phys. Rev. Lett.136, 041005 (2026), arXiv:2507.01191 [hep-ph]
-
[46]
K. Inomata, M. Kamionkowski, K. Kasai, and B. Shakya, Phys. Rev. D112, 083523 (2025), arXiv:2412.17912 [astro-ph.CO]
- [47]
- [48]
-
[49]
L. Niemi and T. V. I. Tenkanen, Phys. Rev. D111, 075034 (2025), arXiv:2408.15912 [hep-ph]
-
[50]
S. Jiang, F. P. Huang, and C. S. Li, Phys. Rev. D108, 063508 (2023), arXiv:2305.02218 [hep-ph]
- [51]
- [52]
- [53]
-
[54]
D. Baumann, P. J. Steinhardt, K. Takahashi, and K. Ichiki, Phys. Rev. D76, 084019 (2007), arXiv:hep-th/0703290
-
[55]
J. R. Espinosa, D. Racco, and A. Riotto, JCAP09, 012 (2018), arXiv:1804.07732 [hep-ph]
work page Pith review arXiv 2018
-
[56]
K. Kohri and T. Terada, Phys. Rev. D97, 123532 (2018), arXiv:1804.08577 [gr-qc]
-
[57]
Primordial Black Holes - Perspectives in Gravitational Wave Astronomy -
M. Sasaki, T. Suyama, T. Tanaka, and S. Yokoyama, Class. Quant. Grav.35, 063001 (2018), arXiv:1801.05235 [astro-ph.CO]
work page Pith review arXiv 2018
-
[58]
Dom` enech, Universe7, 398 (2021), arXiv:2109.01398 [gr-qc]
G. Dom` enech, Universe7, 398 (2021), arXiv:2109.01398 [gr-qc]
-
[59]
A topic review on probing primordial black hole dark matter with scalar induced gravitational waves,
C. Yuan and Q.-G. Huang, iScience24, 102860 (2021), arXiv:2103.04739 [astro-ph.GA]
-
[60]
G. Dom` enech, S. Pi, A. Wang, and J. Wang, JCAP08, 054 (2024), arXiv:2402.18965 [astro-ph.CO]
-
[61]
S. Wang, Z.-C. Zhao, and Q.-H. Zhu, Phys. Rev. Res.6, 013207 (2024), arXiv:2307.03095 [astro- ph.CO]
- [62]
- [63]
-
[64]
C. Yuan, Z.-C. Chen, and Q.-G. Huang, Phys. Rev. D100, 8 (2019), arXiv:1906.11549 [astro-ph.CO]
-
[65]
Planck 2018 results. VI. Cosmological parameters
N. Aghanimet al.(Planck), Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO]
work page internal anchor Pith review arXiv 2020
-
[66]
T. Nakama and T. Suyama, Phys. Rev. D92, 121304 (2015), arXiv:1506.05228 [gr-qc]
-
[67]
Afzalet al.(NANOGrav), Astrophys
A. Afzalet al.(NANOGrav), Astrophys. J. Lett.951(2023), 10.3847/2041-8213/acdc91, arXiv:2306.16219 [astro-ph.HE]
- [68]
-
[69]
Induced gravitational waves with kination era for recent pulsar timing array signals,
K. Harigaya, K. Inomata, and T. Terada, Phys. Rev. D108, 123538 (2023), arXiv:2309.00228 [astro- ph.CO]. 22
-
[70]
V. De Luca, A. Kehagias, and A. Riotto, Phys. Rev. D108, 063531 (2023), arXiv:2307.13633 [astro- ph.CO]
-
[71]
S. Choudhury, K. Dey, A. Karde, S. Panda, and M. Sami, Phys. Lett. B856, 138925 (2024), arXiv:2310.11034 [astro-ph.CO]
-
[72]
X. Wang, X.-H. Ma, and Y.-F. Cai, Int. J. Mod. Phys. D34, 2550027 (2025), arXiv:2412.19631 [astro-ph.CO]
-
[73]
S. Choudhury, K. Dey, S. Ganguly, A. Karde, S. K. Singh, and P. Tiwari, Eur. Phys. J. C85, 472 (2025), arXiv:2409.18983 [astro-ph.CO]
- [74]
-
[75]
T. Suyama and J. Yokoyama, Phys. Rev. D84, 083511 (2011), arXiv:1106.5983 [astro-ph.CO]
-
[76]
Y.-F. Cai, X. Tong, D.-G. Wang, and S.-F. Yan, Phys. Rev. Lett.121, 081306 (2018), arXiv:1805.03639 [astro-ph.CO]
work page Pith review arXiv 2018
- [77]
- [78]
-
[79]
On the slope of the curvature power spectrum in non-attractor inflation,
O. ¨Ozsoy and G. Tasinato, JCAP04, 048 (2020), arXiv:1912.01061 [astro-ph.CO]
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.