🔭 astrophysics

Pulsar Timing Array Sensitivity to Anisotropy: Empirical Sensitivity Curves, Scaling Relations, and the Multi-Resolution Pixel Basis

This paper quantifies Pulsar Timing Array sensitivity to gravitational wave background anisotropy using Fisher information matrices and simulations to derive empirical scaling relations with array parameters, project future NANOGrav capabilities, and introduce a multi-resolution pixel basis to account for sky-location-dependent sensitivity variations.

Original authors: Taha T. Moursy, Nihan S. Pol, Gabriella Agazie, Nikita Agarwal, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Anjana Ashok, Jeremy G. Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Taha T. Moursy, Nihan S. Pol, Gabriella Agazie, Nikita Agarwal, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Anjana Ashok, Jeremy G. Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Lankeswar Dey, Timothy Dolch, Graham M. Doskoch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Emiko C. Gardiner, Nate Garver-Daniels, Peter A. Gentile, Kyle A. Gersbach, Joseph Glaser, Deborah C. Good, Kayhan Gültekin, Aiden Gundersen, C. J. Harris, Jeffrey S. Hazboun, Ross J. Jennings, Aaron D. Johnson, Megan L. Jones, David L. Kaplan, Anala K. Sreekumar, Luke Zoltan Kelley, Matthew Kerr, Joey S. Key, Nima Laal, Michael T. Lam, William G. Lamb, Bjorn Larsen, T. Joseph W. Lazio, Natalia Lewandowska, Tingting Liu, Duncan R. Lorimer, Jing Luo, Ryan S. Lynch, Chung-Pei Ma, Dustin R. Madison, Ashley Martsen, Cayenne Matt, Alexander McEwen, James W. McKee, Maura A. McLaughlin, Natasha McMann, Bradley W. Meyers, Patrick M. Meyers, Matthew T. Miles, Chiara M. F. Mingarelli, Andrea Mitridate, Cherry Ng, David J. Nice, Shania A. Nichols, Stella Koch Ocker, Daniel J. Oliver, Ken D. Olum, Timothy T. Pennucci, Benetge B. P. Perera, Polina Petrov, Henri A. Radovan, Scott M. Ransom, Paul S. Ray, Joseph D. Romano, Jessie C. Runnoe, Alexander Saffer, Shashwat C. Sardesai, Ann Schmiedekamp, Carl Schmiedekamp, Kai Schmitz, Levi Schult, Brent J. Shapiro-Albert, Xavier Siemens, Joseph Simon, Sophia V. Sosa Fiscella, Ingrid H. Stairs, Daniel R. Stinebring, Kevin Stovall, Abhimanyu Susobhanan, Joseph K. Swiggum, Jacob Taylor, Stephen R. Taylor, Mercedes S. Thompson, Jacob E. Turner, Michele Vallisneri, Rutger van Haasteren, Sarah J. Vigeland, Haley M. Wahl, Kevin P. Wilson, Caitlin A. Witt, David Wright, Olivia Young

Published 2026-08-11

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Taha T. Moursy, Nihan S. Pol, Gabriella Agazie, Nikita Agarwal, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Anjana Ashok, Jeremy G. Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Lankeswar Dey, Timothy Dolch, Graham M. Doskoch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Emiko C. Gardiner, Nate Garver-Daniels, Peter A. Gentile, Kyle A. Gersbach, Joseph Glaser, Deborah C. Good, Kayhan Gültekin, Aiden Gundersen, C. J. Harris, Jeffrey S. Hazboun, Ross J. Jennings, Aaron D. Johnson, Megan L. Jones, David L. Kaplan, Anala K. Sreekumar, Luke Zoltan Kelley, Matthew Kerr, Joey S. Key, Nima Laal, Michael T. Lam, William G. Lamb, Bjorn Larsen, T. Joseph W. Lazio, Natalia Lewandowska, Tingting Liu, Duncan R. Lorimer, Jing Luo, Ryan S. Lynch, Chung-Pei Ma, Dustin R. Madison, Ashley Martsen, Cayenne Matt, Alexander McEwen, James W. McKee, Maura A. McLaughlin, Natasha McMann, Bradley W. Meyers, Patrick M. Meyers, Matthew T. Miles, Chiara M. F. Mingarelli, Andrea Mitridate, Cherry Ng, David J. Nice, Shania A. Nichols, Stella Koch Ocker, Daniel J. Oliver, Ken D. Olum, Timothy T. Pennucci, Benetge B. P. Perera, Polina Petrov, Henri A. Radovan, Scott M. Ransom, Paul S. Ray, Joseph D. Romano, Jessie C. Runnoe, Alexander Saffer, Shashwat C. Sardesai, Ann Schmiedekamp, Carl Schmiedekamp, Kai Schmitz, Levi Schult, Brent J. Shapiro-Albert, Xavier Siemens, Joseph Simon, Sophia V. Sosa Fiscella, Ingrid H. Stairs, Daniel R. Stinebring, Kevin Stovall, Abhimanyu Susobhanan, Joseph K. Swiggum, Jacob Taylor, Stephen R. Taylor, Mercedes S. Thompson, Jacob E. Turner, Michele Vallisneri, Rutger van Haasteren, Sarah J. Vigeland, Haley M. Wahl, Kevin P. Wilson, Caitlin A. Witt, David Wright, Olivia Young

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe is filled with a faint, cosmic hum—a background noise made of ripples in space-time called gravitational waves. These ripples are like the sound of a distant ocean, but instead of water, they are made of the fabric of reality itself. For decades, scientists have been trying to "listen" to this hum using a unique kind of telescope: not a dish that looks at the sky, but a collection of cosmic lighthouses called pulsars. These are dead stars that spin incredibly fast, beaming radio waves toward Earth like a lighthouse beam sweeping across the sea. Because they spin so steadily, they act as the most precise clocks in the universe.

When a gravitational wave passes through the space between Earth and a pulsar, it stretches and squeezes space, causing the light from the pulsar to arrive a tiny bit early or late. By comparing the arrival times of many pulsars across the sky, scientists can try to detect the pattern of these waves. The big question is: Is this cosmic hum coming from everywhere equally (like white noise from a radio), or does it have a specific shape, with louder spots and quieter spots? Finding out if the hum is "lumpy" or "smooth" could tell us if the waves are caused by giant black holes colliding in our neighborhood or by something exotic from the very beginning of the universe.

The Cosmic Detective Work: Mapping the Lumpy Hum

In this paper, the NANOGrav collaboration (a team of scientists using North American pulsars) acts like a group of detectives trying to figure out how good their "ears" are at hearing the shape of this cosmic hum. They didn't just listen to the real sky; they built a massive digital simulation—a virtual universe filled with fake pulsars and fake gravitational waves—to test their tools. Their goal was to create a "sensitivity map," which is basically a report card showing exactly how well they could spot a lumpy gravitational wave background depending on where in the sky it was hiding.

The Tools of the Trade: Pixel Maps and Spherical Harmonics
To understand the shape of the cosmic hum, the scientists used two different ways to slice up the sky, kind of like choosing between a grid map and a globe with latitude and longitude lines.

The Pixel Basis: Imagine the sky is a giant screen covered in a grid of tiny squares (pixels). The scientists checked the power of the waves in each square.
The Spherical Harmonic Basis: Imagine the sky is a beach ball. Instead of squares, they looked at the waves using patterns that wrap around the ball, like the stripes on a soccer ball or the swirls on a marble.

By using these two methods, they could measure how sensitive their array of pulsars was to finding anisotropy (which is just a fancy word for "not being the same in every direction").

What They Found: The Rules of the Game
The team ran thousands of simulations to see how different factors changed their ability to hear the cosmic hum. They discovered some very specific rules about how their "ears" work:

More Pulsars = Much Better Hearing: The most important finding is that adding more pulsars to the array makes a huge difference. They found that sensitivity scales with the number of pulsars ( $N_{psr}$ ) to the power of about 0.8. In everyday terms, if you double the number of pulsars, you don't just double your hearing ability; you get a massive boost, making it much easier to spot the lumpy parts of the universe.
Better Clocks Help, But Not as Much: They also tested what happens if the pulsars' clocks are more accurate (meaning less error in the arrival times, or $\delta t$ ). They found that improving the clock precision helps, but only a tiny bit. The sensitivity scales with the error to the power of -0.08. This means that even if you make the clocks incredibly precise, it won't help you as much as simply adding more pulsars to the team.
The Size of the "Lump" Matters: The scientists found that it is easier to spot a big, wide "lump" in the cosmic hum than a tiny, sharp one. The sensitivity scales with the size of the feature ( $\Delta\Omega$ ) to a power between 1.6 and 2.1. If the lumpy part of the universe is huge, the pulsars can hear it easily. If it's a tiny speck, it's much harder to detect.
Time Alone Isn't the Magic Bullet: A common hope is that if we just keep listening for longer (say, 30 years instead of 16), we will hear everything clearly. The simulations showed that simply waiting longer doesn't help much, unless you are looking at the very lowest frequencies. The biggest gains come from having more pulsars and better sky coverage, not just more time.

The "Dirty" Secret: Why Math Matters
One of the paper's key insights is about how the scientists do their math. When they look at pairs of pulsars, the data from one pair isn't totally independent of another pair; they are "covariant," meaning they are linked. The paper shows that if you ignore these links (a common shortcut), you might think you are twice as sensitive as you actually are, especially at low frequencies. It's like thinking you have super-hearing because you forgot to account for the echo in the room. By doing the full, complex math that includes these links, they get a much more honest (and slightly less optimistic) picture of their sensitivity.

A New Way to Look at the Sky: The Multi-Resolution Map
Finally, the authors introduced a clever new tool called a "multi-resolution pixel basis." Imagine trying to draw a map of the world. If you use a grid where every square is the same size, you might have tiny squares over the ocean (where there are no pulsars) and huge squares over the land (where there are many). This is wasteful.

The new method allows the map to change its zoom level. In areas where the scientists have many pulsars (like the northern sky for NANOGrav), they can use tiny, high-resolution pixels to see fine details. In areas with few pulsars, they use larger, blurry pixels. This lets them balance the uncertainty: they can get a sharp picture where they have the data and a fuzzy picture where they don't, rather than forcing the whole sky to be the same resolution. They tested this with simulations and showed it works, proving that they can trade "sharpness" for "uncertainty" depending on where they are looking.

The Bottom Line
This paper doesn't claim to have found a new black hole or a new type of wave. Instead, it provides the rulebook for how to find them. It tells us that to hear the shape of the universe's gravitational hum, we need to focus on finding more pulsars and spreading them out across the sky, rather than just waiting longer or making our clocks slightly more precise. It also warns us to be careful with our math, ensuring we don't overestimate how well we can hear the cosmic ocean. With these new sensitivity curves and tools, the NANOGrav team is better equipped than ever to map the lumpy, beautiful structure of the gravitational wave background.

Technical Summary: Pulsar Timing Array Sensitivity to Anisotropy

Problem Statement
Pulsar Timing Arrays (PTAs) have recently reported evidence for a gravitational wave background (GWB) consistent with the Hellings-Downs correlation signature. While the origin of this background is debated between an astrophysical population of supermassive black hole binaries (SMBHBs) and cosmological sources, distinguishing between these scenarios is critical. Cosmological backgrounds are generally expected to be statistically isotropic, whereas astrophysical backgrounds may exhibit significant anisotropy due to loud, nearby SMBHBs or correlations with the large-scale distribution of galaxies. To validate these interpretations, PTAs require robust methods to quantify their sensitivity to anisotropy in the GWB. Existing sensitivity curves often assume isotropy or neglect the covariance between pulsar-pair cross-correlations, potentially leading to overestimated sensitivity. Furthermore, the dependence of anisotropy sensitivity on PTA parameters—such as the number of pulsars ( $N_{psr}$ ), timing precision ( $\delta t$ ), observation duration ( $T_{span}$ ), and the angular scale of the anisotropy—requires empirical characterization.

Methodology
The authors utilize a cross-correlation-based Fisher information matrix approach to quantify PTA sensitivity to anisotropy in two bases: the radiometer pixel basis (using HEALPix tessellation) and the spherical harmonic basis. The analysis employs the per-frequency optimal statistic (PFOS) to compute pulsar-pair cross-correlations, their uncertainties, and the full covariance matrix, explicitly including pulsar-pair covariance which is induced by the isotropic null hypothesis.

To derive scaling relations, the authors generate 1,200 simulated datasets based on the NANOGrav 15-year dataset. These simulations vary the number of pulsars ( $N_{psr} \in \{17, 34, 51, 67\}$ ), the root-mean-square timing residual uncertainty ( $\delta t \in \{100, 500, 1000\}$ ns), and the gravitational wave frequency ( $f_{GW}$ ). Intrinsic red noise and an isotropic GWB (consistent with NANOGrav 15-year results) are injected into the simulations. A second set of 50 realistic simulations extends the NANOGrav dataset to a 30-year baseline to forecast future sensitivity growth.

The Fisher matrix is computed analytically using the inverse of the covariance matrix of the cross-correlations. Sensitivity is quantified as the square root of the diagonal elements of the inverse Fisher matrix ( $\sqrt{\text{diag}(M^{-1})}$ ), representing the uncertainty in the recovered power or spherical harmonic coefficients. The authors also introduce a multi-resolution pixel basis using the UNIQ scheme in HEALPix to address the strong dependence of sensitivity on sky location and pulsar density.

Key Contributions and Results

Empirical Scaling Relations:
The study establishes empirical scaling laws for PTA sensitivity to anisotropy (represented by $\sqrt{\text{diag}(M)}$ ):
- Number of Pulsars: Sensitivity scales approximately as $N_{psr}^{0.8 \pm 0.1}$ . This strong dependence is attributed to the use of pulsar-pair cross-correlations.
- Timing Uncertainty: Sensitivity scales as $\delta t^{-0.08 \pm 0.01}$ . The dependence is weaker than that on the number of pulsars.
- Frequency: Sensitivity scales as $f_{GW}^{-0.12 \pm 0.04}$ (for the Fisher matrix elements alone).
- Angular Scale: For spherical harmonics with $8 \le \ell \le 20$ , sensitivity scales with angular scale $\Delta\Omega$ as $\Delta\Omega^{1.58 \pm 0.03}$ . The scaling is shallower at low $\ell$ and depends on the $m$ value (orientation) due to the anisotropic distribution of pulsars in the array.
Impact of Pulsar-Pair Covariance:
The authors demonstrate that neglecting the covariance between pulsar-pair cross-correlations leads to a significant overestimation of sensitivity, particularly at low frequencies where the GWB signal is strongest relative to noise. At the lowest frequencies, sensitivity is overestimated by a factor of approximately two. This effect diminishes at higher frequencies.
Sky Location Dependence:
Sensitivity to anisotropy varies significantly across the sky, with the most sensitive regions (high pulsar density) being 2–3 times more sensitive than the least sensitive regions. This variation is driven by the non-uniform distribution of pulsars in current PTAs.
Forecasting Future Sensitivity:
Projections to a 30-year baseline indicate that increasing the observation duration ( $T_{span}$ ) alone, without increasing the number of pulsars or improving timing precision, yields negligible improvements in average sensitivity to anisotropy. While some growth is observed at the lowest frequencies (where previously short baselines become sufficient to contribute), the overall gain is small compared to the potential gains from increasing $N_{psr}$ or sky coverage.
Multi-Resolution Pixel Basis:
Motivated by the strong sky-location dependence, the authors introduce a multi-resolution pixel basis. This framework allows for varying the resolution of the sky map (and thus the trade-off between localization and uncertainty) independently in different regions. Injection-and-recovery simulations using this basis demonstrate its ability to balance uncertainties across the sky, leveraging regions of high pulsar density for higher resolution.

Significance
The paper provides a rigorous framework for quantifying PTA sensitivity to anisotropy, correcting previous overestimations by including pulsar-pair covariance. The derived scaling relations highlight that increasing the number of pulsars and improving sky coverage (e.g., through the International Pulsar Timing Array) are far more effective strategies for enhancing anisotropy sensitivity than simply extending observation times. The introduction of the multi-resolution pixel basis offers a novel tool to adaptively map the GWB anisotropy, accounting for the heterogeneous nature of current PTA datasets. These results are essential for interpreting future anisotropy searches and distinguishing between astrophysical and cosmological origins of the GWB.

More like this