alphaFeedback

WorksPhysics and Astronomy

Hellings and Downs correlation of an arbitrary set of pulsars

AI contributions
Not recorded for this paper
View details

Abstract

Pulsar timing arrays (PTAs) detect gravitational waves (GWs) via the correlations they induce in the arrival times of pulses from different pulsars. We assume that the GWs are described by a Gaussian ensemble. The mean correlation h2μu(γ)h^2 μ_{\rm u}(γ) as a function of the angle γγ between the directions to two pulsars was predicted by Hellings and Downs (HD) in 1983. The variance σtot2(γ)σ^2_{\rm tot}(γ) in this correlation was recently calculated by Allen[11] for a single noise-free pulsar pair at angle γγ, which shows that after averaging over many pairs, the variance reduces to an intrinsic cosmic variance σcos2(γ)σ^2_{\rm cos}(γ). Here, we extend this to an arbitraryarbitrary set of pulsars at specific sky locations, with pulsar pairs binned by γγ. We derive the linear combination of pulsar-pair correlations which is the optimal estimator of the HD correlation for each bin, illustrating our methods with plots of the expected range of variation away from the HD curve, for the sets of pulsars monitored by three active PTA collaborations. We compute the variance of and the covariance between these binned estimates, and show that these reduce to the cosmic variance and covariance s(γ,γ′)s(γ,γ') respectively, in the many-pulsar limit. The likely fluctuations away from the HD curve μu(γ)μ_{\rm u}(γ) are strongly correlated/anticorrelated in the three angular regions where μu(γ)μ_{\rm u}(γ) is successively positive, negative, and positive. We also construct the optimal estimator of the squared strain h2h^2. When there are very many pulsar pairs, this determines h2h^2 with arbitrary precision because PTAs probe an infinite set of GW modes. To assess observed deviations away from the HD curve, we characterize several χ2χ^2 goodness-of-fit statistics. We also show how pulsar noise and measurement noise can be included.

Introduction

A pulsar timing array (PTA) is a galactic-scale gravitational-wave (GW) detector. It searches for low-frequency (nanohertz) GWs by precisely monitoring the arrival times of pulses from a set of pulsars [1]. Gravitational waves (e.g., from inspiraling supermassive black-hole binaries in the centers of merging galaxies) influence the pulse arrival times in a way that is correlated between different pulsars. The mean correlation between a pair of pulsars depends upon the angular separation γ\gamma between the lines of sight to each member of the pair, as seen from Earth [2].

The mean correlation was calculated by Hellings and Downs [3] for a unit amplitude, isotropic and unpolarized GW background. It has the simple analytic form [4]

μu(γ)=13−16(1−cos⁡γ2)+1−cos⁡γ2ln⁡(1−cos⁡γ2),(1)\mu_u(\gamma)=\frac{1}{3}-\frac{1}{6}\left(\frac{1-\cos\gamma}{2}\right)+\frac{1-\cos\gamma}{2}\ln\left(\frac{1-\cos\gamma}{2}\right), \tag*{(1)}

where γ\gamma is the angular separation between a pair of pulsars, and “u” means “unpolarized.” Observation of a correlation proportional to this Hellings and Downs curve is the “smoking gun” signature that a PTA has detected GWs [5].

Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. Open access publication funded by the Max Planck Society.

Several groups are searching for such correlations. These groups report strong statistical evidence for fluctuations in the individual pulsar arrival times that share the same “red” spectrum [6, 7, 8, 9] that a GW background is expected to produce [10]. However, there is currently little evidence for the Hellings and Downs angular dependence, which makes it difficult to claim that a GW background is responsible.

Could the lack of evidence for the Hellings and Downs angular dependence be a statistical fluctuation? To answer this, it is important to understand what variations away from the Hellings and Downs predicted mean might be expected. The size of such fluctuations is quantified by the variance of the Hellings and Downs correlation.

Total variance, pulsar variance, and cosmic variance

Recent work [11] calculates this variance for several GW source models, neglecting all sources of noise. The most important model contains NN unpolarized point sources, uniformly distributed (statistically) in space, radiating GWs at the same frequency but with independent random phases. For large numbers of sources, this creates stationary “confusion noise.” It is a Gaussian stochastic process [12, 13, 14] provided that the amplitude of the individual sources vanishes as N→∞N \to\infty.

The (total) variance σtot2\sigma_{\mathrm{tot}}^2 computed in [11] is a sum of “pulsar variance” and “cosmic variance” σcos2\sigma_{\mathrm{cos}}^2. Both types of variance have been observed in simulations based on synthetic catalogs of sources [15]. Pulsar variance arises because different pairs of pulsars separated by the same angle γ\gamma have correlations that differ from the average, in a way that depends (unpredictably) upon their sky positions. Reference [11] computes the total variance (there denoted by σ2\sigma^2 without a subscript) for a single randomly selected pulsar pair separated by angle γ\gamma. In contrast, the cosmic variance is the variance of the correlation after the correlation has been averaged over all possible locations/orientations of the pair. This corresponds to employing an infinite number of pulsar pairs separated by angle γ\gamma, uniformly distributed about the sky, and is called pulsar averaging.

In the confusion-noise model (Sec. III A in [11]), cosmic variance arises because, even after pulsar averaging, the correlation depends upon the relative phases of the GW sources. Each realization of the Universe has different phases, and thus exhibits different pulsar-averaged correlations. The cosmic variance is the amount by which the pulsar-averaged correlation curve is expected to differ from the Hellings and Downs prediction in an ideal world containing an infinite number of noise-free pulsars. Unlike the pulsar variance, it cannot be reduced: the cosmic variance is a fundamental limit to the precision with which the Hellings and Downs predicted mean might be observed at a particular angle γ\gamma.

Figure 1 illustrates the total and cosmic variance, as offsets away from the mean μu\mu_u of the Hellings and Downs correlation. The total variance at angle γ\gamma is the uncertainty associated with the determination of the Hellings and Downs correlation when a single (randomly selected) pulsar pair at angle γ\gamma is used to estimate that correlation. The cosmic variance at angle γ\gamma is the uncertainty that remains when an infinite number of uniformly distributed pulsar pairs at angle γ\gamma are used to do the estimation. In this paper, we study the transition between these two limits, when a finite set of pulsars at specific sky locations are used to estimate the correlation. This reflects observational reality because PTA pulsars are nonuniformly distributed on the sky, so the pairs formed from them have no separation angles in common.

Mean Hellings and Downs correlation and its variance

Figure 1. The mean Hellings and Downs correlation μu\mu_u as a function of the angular separation γ\gamma between a pair of pulsars. Also plotted are μu±σtot\mu_u \pm\sigma_{\mathrm{tot}} and μu±σcos\mu_u \pm\sigma_{\mathrm{cos}}, where σtot2\sigma_{\mathrm{tot}}^2 is the total variance for a single pulsar pair and σcos2\sigma_{\mathrm{cos}}^2 is the cosmic variance. (For this plot, we have set R4/h4=1/2\mathcal{R}^4/h^4 = 1/2 and h2=1h^2 = 1; see text for details of the GW source model.) Note that although the correct term is “standard deviation,” we sometimes call σ\sigma the “variance”.

The plots in Figure 1 assume that the GWs arise from the incoherent sum of many weak sources, giving rise, via the central-limit theorem, to a Gaussian ensemble [12, 13, 14]. For such sources, the scaling relation between the (squared) mean and the variance is described in Appendix B; these plots take R4/h4=1/2\mathcal{R}^4/h^4 = 1/2 and h2=1h^2 = 1. This corresponds to the large source-number limit of the “narrow-band” discrete confusion-noise model (Sec. III A in [11]), where source frequency is assumed to be commensurate with the inverse observation time.

Variance of the Hellings and Downs correlation for an arbitrary set of pulsars

We assume that the reader is familiar with Ref. [11], and extend that analysis in five ways:

(i) First, we define the Hellings and Downs correlation [3] for many pulsar pairs. We have found only one sensible way to generalize the standard “single-pulsar-pair” definition. We partition the angular range γ∈[0,180∘]\gamma\in[0,180^\circ] into nonoverlapping intervals called bins. Then, we define the correlation associated with a bin [γlow,γhigh][\gamma_{\mathrm{low}},\gamma_{\mathrm{high}}] to be an average of the Hellings and Downs correlations ρab\rho_{ab} for the npairsn_{\mathrm{pairs}} pulsar pairs that have γlow<γab<γhigh\gamma_{\mathrm{low}} < \gamma_{ab} < \gamma_{\mathrm{high}}, where the subscripts aa and bb label pulsars, and γab\gamma_{ab} is their angular separation on the sky. We say that these pulsar pairs “fall in” or “lie in” or “are in” the bin. If there is only a single pulsar pair in the bin, then our definition reduces to the standard Hellings and Downs definition.

What type of average should we employ? As shown in Figure 2, a uniform average, where the correlation of each pulsar pair gets the same weight, 1/npairs1/n_{\mathrm{pairs}}, is not the best choice. We define the Hellings and Downs correlation for a given bin to be the weighted sum of the pulsar-pair correlations ρab\rho_{ab} in that bin which is (1) unbiased and (2) minimizes the expected variance.

Schematic illustration of pulsar pairs in an angular-separation bin

Figure 2. The bin at angular separation γ\gamma contains npairs=3n_{\mathrm{pairs}}=3 pulsar pairs with γ∼γ12∼γ38∼γ47\gamma\sim\gamma_{12} \sim\gamma_{38} \sim\gamma_{47}. The naive correlation estimator ρnaive=(ρ12+ρ38+ρ47)/3\rho_{\mathrm{naive} }=(\rho_{12}+\rho_{38}+\rho_{47})/3 weights them uniformly, but is sub-optimal. The optimal estimator ρopt=0.3ρ12+0.3ρ38+0.4ρ47\rho_{\mathrm{opt}}=0.3\rho_{12}+0.3\rho_{38}+0.4\rho_{47} gives more weight to the correlation ρ47\rho_{47}. This is because the 12 and 38 pairs are close on the sky, so ρ12\rho_{12} and ρ38\rho_{38} give nearly redundant estimates. The optimal estimator defines the Hellings and Downs correlation of that bin. Note: for large numbers of pairs uniformly distributed about the sky, rotational symmetry implies that all pairs have the same weight.

In this paper, that weighted average (for a particular bin) is called the “optimal estimator.” In several places, we also consider a special case, which we call the “narrow-bin limit.” This is the (mathematical) limit in which the width (in γ\gamma) of an angular-separation bin vanishes, while the number of pulsar pairs in that bin remains constant or grows.

Conditions (1) and (2) lead immediately to a simple formula for the weights. The weights depend upon the sky directions to all of the pulsars which contribute to a particular bin, and are easily computed for any specific set of directions. In the limit of low pulsar noise, the resulting optimal estimator tells an observer (in a universe described by the Gaussian ensemble) the best way to combine measured pulsar-pair correlations to estimate the Hellings and Downs correlation in a particular bin.

As a byproduct, we also obtain simple formulas for (a) the variance of the optimal estimator for a given angular separation bin, and (b) the covariance of the optimal estimators for any pair of bins.

By definition these are the variance and covariance of the Hellings and Downs correlation. In particular, the variance, which by condition (2) is as small as possible, quantifies the uncertainty in the Hellings and Downs correlation for any particular bin. The variance depends upon the sky directions to all of the pulsars that contribute to that bin.

To illustrate this, Figure 3 shows the variance obtained by using the 88 pulsars currently employed by three active PTAs, after dividing the 88×87/2=382888 \times87/2=3828 possible distinct pulsar pairs among 30 evenly spaced 6∘6^\circ bins. If the Gaussian ensemble is a good description of the GWs in our Universe and the correlation measurements are noise-free, then this shows the expected deviations away from the Hellings and Downs mean. Thus, it is an upper limit on the ability to which that set of 88 pulsars could be expected to recover the Hellings and Downs curve. Similar plots for the three individual PTAs (employing only “their” pulsars) are given in Sec. V B.

Predicted GW-induced variance as a function of pulsar separation angle

Figure 3. The dots (lines are to guide the eye) show the predicted GW-induced variance σopt2\sigma_{\mathrm{opt}}^2 for the current set of 88 IPTA pulsars (assumed free of noise) after optimally combining timing-residual correlations in 30 ×\times 6∘6^\circ angular-separation bins. This reduces the variance far below the single pulsar pair (total) variance σtot2\sigma_{\mathrm{tot}}^2, bringing it close to the cosmic variance σcos2\sigma_{\mathrm{cos}}^2. The dots indicate the scale of the expected fluctuations away from the Hellings and Downs curve, for a Gaussian ensemble of GW background with the binary-inspiral spectrum described in Appendix B. This plot has the same scale and assumptions as the IPTA plot in Figure 9, i.e., α=1\alpha=1 for timing residuals and h22≈0.5622h_2^2 \approx0.5622. See Sec. V B and Appendix § 24 for more details.

(ii) Second, we extend [11] by showing how the variance and covariance of the expected Hellings and Downs correlation decrease as pulsar pairs are added. For example, the variance of the optimal estimator for narrow bins drops from the total variance σtot2(γ)\sigma_{\mathrm{tot}}^2(\gamma) for a single pulsar pair to the cosmic variance σcos2(γ)\sigma_{\mathrm{cos}}^2(\gamma) for an infinite number of pairs (distributed uniformly on the sky). This transition is illustrated later in Figure 6.

Variance of the optimal estimator and minimum number of pulsar pairs as a function of pulsar separation angle

Figure 6. Left panel: the variance of the optimal estimator as the number of pulsar pairs at each angular separation γ\gamma is increased. Shown in the plot, from top to bottom, is σopt(γ)\sigma_{\mathrm{opt}}(\gamma) calculated for 1, 2, 5, 10, 25, 50, 100, and 200 pairs. (For this plot we have set κ4=1/2\kappa^{4}=1/2.) Right panel: the minimum number of pulsar pairs needed for σopt2(γ)\sigma_{\mathrm{opt}}^{2}(\gamma) to reach a value that is 1/e≈37%1/e\approx37\% larger than the cosmic variance. For this number of pulsar pairs, σopt2(γ)≈(1+1/e)σcos2(γ)\sigma_{\mathrm{opt}}^{2}(\gamma)\approx(1+1/e)\sigma_{\mathrm{cos}}^{2}(\gamma).

Our approach to deriving the cosmic variance is different than that given in Sec. IV of [11], and provides a useful

alternative interpretation. Because different pulsar-pair correlations are correlated with one another, as pulsar pairs are added to a bin, the variance does not approach zero. The reason is straightforward: once there are a sufficient number of pulsar pairs in the bin, adding more pairs to that bin provides negligible additional information about the GW-induced correlations. (Section III E demonstrates this behavior with a simple example.)

We also derive an elegant analytic expression for the covariance of the optimal estimators in narrow angular separation bins labeled by jj and kk, in the limit where each contains an infinite number of pulsar pairs. This “cosmic covariance matrix” is not given in [11], nor have we found it elsewhere in the literature. We give it both in “position-space” form and in harmonic space, where it is a diagonal sum of products of Legendre polynomials, dominated by the quadrupole plus a few additional multipoles. The cosmic covariance has a surprisingly simple structure, showing that the expected fluctuations away from the Hellings and Downs mean correlation are strongly correlated/anticorrelated in the three γ\gamma-angle regions where the Hellings and Downs curve is successively positive, negative, and positive.

(iii) Third, we address the question: “What is the strength of the GW stochastic background in the PTA band?” For simplicity, we assume that the frequency-dependent shape of the GW spectrum is known (e.g., ∝∣f∣−7/3\propto|f|^{-7/3} for binary-inspiral sources). The goal is then to estimate the overall scale, which is the squared GW strain h2h^2.

We construct optimal estimators of h2h^2 that are unbiased and have minimum variance. These estimators are linear combinations of the measured pulsar-pair correlations ρab\rho_{ab}. We consider three different possible choices of the correlation set: (a) auto+cross: use all possible pairs of pulsars; (b) cross only: use all pulsar pairs where the two pulsars differ; and (c) auto only: use all pulsar pairs where the two pulsars are the same. We then derive the variance (squared uncertainty) in these estimators of h2h^2. As the number of pulsars increases, these uncertainties get smaller. Their scaling behavior, as a function of the number of pulsars, depends upon which correlation set is used.

We show that if cross-correlations are included in the squared-strain estimator, then the variance tends to zero as Npul→0N_{\rm pul} \to0. Initially, we were surprised by this: in the same limit, the Hellings and Downs correlation has a nonzero cosmic variance, as do local estimates of the squared strain, such as hμνhμνh_{\mu\nu}h^{\mu\nu}. We show that this remarkable behavior arises because a pulsar’s response to GWs allows a PTA to probe infinitely many modes of the GW field.

(iv) Fourth, to enable straightforward tests for consistency with the predictions of the Gaussian ensemble, we construct a set of χ2\chi^2 goodness-of-fit statistics. These can be used to assess if observational measurements of the Hellings and Downs correlation are consistent with the expected deviations away from h2μu(γ)h^2\mu_u(\gamma), based on our predictions of the variance and covariance for the Gaussian ensemble. In defining these χ2\chi^2 statistics, there are three possible choices for which pulsar-pair correlations to employ (as explained above in the context of estimating h2h^2) and two possible choices for how the overall scale of the GW strain is determined. By taking all possible combinations, we arrive at six different statistics.

To set the scale of the GW strain, there are two options: (1) assume that h2h^2 has been obtained via other methods/data, or (2) obtain h2h^2 from the goodness-of-fit statistic, via a process we call “projection,” which minimizes χ2\chi^2. The justification for this choice of name, and the six different options, are discussed in detail in Sec. VIII. Note that the projection process leads to the same h2h^2 estimators as discussed in item (iii) above.

The χ2\chi^2 statistics are a “sum of squares”; we denote them with the same symbol as used for the traditional χ2\chi^2 test, whose distribution describes the sum of squares of normal Gaussian variables. However, the distributions of our χ2\chi^2 statistics and the standard χ2\chi^2 statistics differ in an important way: the random variables whose squares are summed are not Gaussian, even though we assume that the GW stochastic background is described by the Gaussian ensemble. The χ2\chi^2 statistics that we define have the same mean values as the standard distribution, but they have much larger variances because the GW background creates non-Gaussian fluctuations in the pulsar-pair correlations.

(v) Our final extension to [11], which assumes noise-free data, is to also consider the effects of pulsar and measurement noise. While most of our paper considers ideal data, noise is easily included. We do this starting in Sec. IX, assuming that the noise is uncorrelated between different pulsars. The noise modifies the pulsar-pair correlations and thus adds additional terms to the pulsar-pair expectation value, covariance matrix, and higher-order moments. This means that our entire approach and methodology can still be used in the presence of noise.

The noise contributions to the pulsar-pair correlations are fundamentally different than those which are induced by the GWs. For the GW-induced correlations, as the number of pulsar pairs increases, the variance in the Hellings and Downs correlation approaches the (nonzero) cosmic variance. In contrast, since we assume that the different pulsars have independent noise, the contribution of their noise to the variance of the Hellings and Downs correlation vanishes as the number of pulsar pairs increases. Thus, there are reasons to hope that during the coming decades, as PTAs add additional pulsars and improve their measurements, the limit of the cosmic variance can be achieved in practice.

Outline

Here, we give a brief outline of the paper, and introduce the most important notation. Note that our methods and results apply both to pulsar redshifts (set α=0\alpha= 0) and to pulsar timing residuals (set α=1\alpha= 1). See Appendix A for details.

We begin in § 2 by defining the time-averaged correlation ρab\rho_{ab} between the redshifts of two pulsars aa and bb, where the subscript denotes the pulsar pair abab. For the Gaussian ensemble, in § 2 A we calculate the first moment ⟨ρab⟩\langle\rho_{ab}\rangle of the correlation. For two distinct pulsars (meaning a≠ba \ne b) this is proportional to the Hellings and Downs curve μu(γab)\mu_u(\gamma_{ab}) given in (1), where γab\gamma_{ab} denotes the angle between the directions to the two pulsars. To gracefully accommodate the case a=ba=b, for which the “pulsar term” doubles the correlation, we also define the closely related quantity μab\mu_{ab}, cf. (5). In § 2 B, we use this to compute the “second moment” ⟨ρabρcd⟩\langle\rho_{ab}\rho_{cd}\rangle. From the moments, we construct the covariance matrix Cab;cd≡⟨ρabρcd⟩−⟨ρab⟩⟨ρcd⟩C_{ab;cd}\equiv\langle\rho_{ab}\rho_{cd}\rangle-\langle\rho_{ab}\rangle\langle\rho_{cd}\rangle, where cdcd denotes a second pulsar pair constructed from pulsars cc and dd. Third and fourth moments of ρab\rho_{ab} are computed in § 2 C and Appendix § 12.6.

In this paper, angle brackets ⟨Q⟩\langle Q\rangle denote the average of a function or functional QQ. This average is always normalized so that ⟨1⟩=1\langle1\rangle=1. If not explicitly stated otherwise, then this is an average over the Gaussian ensemble of GW realizations [12, 13, 14]. However, in some places we (state that we) use it differently. For example, in (185) it denotes an average over a Gaussian ensemble of noise sources and GW sources. In (50) it denotes an average over pulsar pairs at fixed separation angles (indicated by the subscript cd∈γkcd\in\gamma_k on the rightmost angle bracket).

The covariance matrix C\mathbf{C} is the most important quantity in our analysis, and is used throughout. It is a square, real, positive-definite symmetric matrix. We denote it by C\mathbf{C} when referring to the matrix as a whole, and use indices Cab;cdC_{ab;cd} to indicate a particular entry of the matrix. In § 2 D, we derive two identities satisfied by the entries Cab;cdC_{ab;cd}, which are used later to obtain the cosmic variance limit.

The way that C\mathbf{C} is laid out, i.e., the indexing of its rows and columns, is discussed in § 2 E. Ordering pulsar pairs abab by their angular separation γab\gamma_{ab} puts the covariance matrix C\mathbf{C} into “block-matrix form”, cf. (15). The different blocks correspond to different angular bins, which we label with integers jj and kk. The submatrix CjkC_{jk} (we use the term “block” or “matrix block”) of the covariance matrix has rows (columns) corresponding to all pulsar pairs abab (cdcd) that lie in the jjth (kkth) angular bin.

Pulsar pairs abab and cdcd that both lie in the same angular-separation bin jj have particular importance. The covariance between such pairs is described by blocks which lie along the diagonal of C\mathbf{C}. To simplify the formulas and calculations for these, we often drop the indices, denoting that matrix (block) by C=CjjC=C_{jj}.

In Sec. III, we turn our attention to defining the Hellings and Downs correlation for an arbitrary collection of pulsars. This is defined as the optimal estimator of correlation in a particular angular-separation bin, which is a weighted average of pulsar-pair correlations in that bin: ρopt=∑abwabρab\rho_{\mathrm{opt}}=\sum_{ab}w_{ab}\rho_{ab}. In Sec. III A, we derive the weights wabw_{ab} by requiring that the estimator be unbiased and have the minimum possible variance. In Sec. III B, we calculate the covariance BjkB_{jk} between the optimal estimates in two angular separation bins jj and kk. (Later, this plays an important role in § 11, where we define χ2\chi^2 goodness-of-fit tests. These tests assess observed correlations, to see if they are consistent with the likely fluctuations away from the Hellings and Downs mean.) In Sec. III C and Sec. III D, we examine the narrow-bin limit of the optimal estimator, and show how the variance of the optimal estimator simplifies if the vector 1\mathbf{1} containing all ones is an eigenvector of the corresponding block CC of the full covariance matrix C\mathbf{C}. We conclude this section with a simple pedagogical example: Sec. III E shows how an infinite collection of measurements can have a nonzero variance.

References

  1. [1]R. S. Foster and D. C. Backer, Constructing a pulsar timing array, Astrophys. J. 361, 300 (1990).
  2. [2]When we write “Earth” we actually mean the solar system barycenter (SSB), which is located near the surface of the Sun, about 8 minutes from Earth. Our choice of language is justified because this distance is five orders of magnitude smaller than the wavelength of the nanoher tz GWs of interest here (years).
  3. [3]R. W. Hellings and G. S. Downs, Upper limits on the isotropic gravitational radiation background from pulsar timing analysis, Astrophys. J. 265, L39 (1983).
  4. [4]We have chosen to normalize μᵤ(γ) so that μᵤ(0) = 1/3, consistent with [3,11]. Another common normalization choice is μᵤ(0) = 1/2. Many of our equations contain an explicit factor of μᵤ(0), allowing any desired normalization to be used.
  5. [5]S. R. Taylor, L. Lentati, S. Babak, P. Brem, J. R. Gair, A. Sesana, and A. Vecchio, All correlations must die: Assessing the significance of a stochastic gravitational-wave background in pulsar timing arrays, Phys. Rev. D 95, 042002 (2017).DOI
  6. [6]Z. Arzoumanian, P. T. Baker, H. Blumer, B. Bécsy, A. Brazier, P. R. Brook, S. Burke-Spolaor, S. Chatterjee, S. Chen, J. M. Cordes et al. (NANOGrav Collaboration), The NANOGrav 12.5 yr data set: Search for an isotropic stochastic gravitational-wave background, Astrophys. J. Lett. 905, L34 (2020).DOI
  7. [7]B. Goncharov, R. M. Shannon, D. J. Reardon, G. Hobbs, A. Zic, M. Bailes, M. Curyło, S. Dai, M. Kerr, M. E. Lower et al., On the evidence for a common-spectrum process in the search for the nanohertz gravitational-wave background with the Parkes Pulsar Timing Array, Astrophys. J. Lett. 917, L19 (2021).DOI
  8. [8]S. Chen, R. N. Caballero, Y.-J. Guo, A. Chalumeau, K. Liu, G. Shaifullah, K. J. Lee, S. Babak, G. Desvignes, A. Parthasarathy et al., Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: Inferences in the stochastic gravitational-wave background search, Mon. Not. R. Astron. Soc. 508, 4970 (2021).DOI
  9. [9]J. Antoniadis, Z. Arzoumanian, S. Babak, M. Bailes, A.-S. B. Nielsen, P. T. Baker, C. G. Bassa, B. Bécsy, A. Berthereau, M. Bonetti et al., The International Pulsar Timing Array second data release: Search for an isotropic gravitational wave background, Mon. Not. R. Astron. Soc. 510, 4873 (2022).DOI
  10. [10]J. D. Romano, J. S. Hazboun, X. Siemens, and A. M. Archibald, Common-spectrum process versus cross-correlation for gravitational-wave searches using pulsar timing arrays, Phys. Rev. D 103, 063027 (2021).DOI
  11. [11]B. Allen, Variance of the Hellings-Downs correlation, Phys. Rev. D 107, 043018 (2023).
  12. [12]B. Allen, The stochastic gravity-wave background: Sources and detection, in Relativistic Gravitation and Gravitational Radiation, Cambridge Contemporary Astrophysics Vol. 26, edited by J.-A. Marck and J.-P. Lasota (Cambridge University Press, Cambridge, England 1997), pp. 373–417, 10.48550/arXiv.gr-qc/9604033.arxiv.org/abs/gr-qc/9604033
  13. [13]B. Allen and J. D. Romano, Detecting a stochastic background of gravitational radiation: Signal processing strategies and sensitivities, Phys. Rev. D 59, 102001 (1999).DOI
  14. [14]J. D. Romano and N. J. Cornish, Detection methods for stochastic gravitational-wave backgrounds: A unified treatment, Living Rev. Relativity 20, 2 (2017).DOI
  15. [15]B. Bécsy, N. J. Cornish, and L. Z. Kelley, Exploring realistic nanohertz gravitational-wave backgrounds, Astrophys. J. 941, 119 (2022).DOI
  16. [16]LIGO’s use of a Fabry-Perot cavity means that the effective arm length is hundreds of times larger than the actual arm length. This effective arm length is longer than the shortest GW wavelength which LIGO can detect. Remarkably, a simple approximation that combines the long-wavelength Michelson response with the single-pole approximation to the Fabry-Perot transfer function gives results that agree with the exact response to within a few percent [43].
  17. [17]A more general definition of the correlation ρ_ab is possible. Rather than multiplying the redshifts in the time domain, we could first convolve them with an arbitrary “filter function,” before integration. This corresponds to filtering the pulsar redshifts in the frequency domain. Thus, using Z̃(f) to denote the Fourier transform of Z(t), we could define ρ_ab = ∫ df K(f)Z̃_a(f)Z̃_b*(f), where K(f) = K*(−f) is a particular frequency-domain filter. An appropriate (real, non-negative) choice of K, determined by the spectrum of the GW background and the noise spectral densities of the pulsars, would decrease the fractional uncertainties of estimators associated with “optimal combinations” of ρ_ab. These correspond to larger signal-to-noise ratios (SNRs) as explained, e.g., in [13] in the context of cross-correlation searches for GW backgrounds.
  18. [18]Here, “different distances” means differing by more than a few characteristic GW wavelengths, so tens of years.
  19. [19]L. Isserlis, On a formula for the product-moment coefficient of any order of a normal frequency distribution in any number of variables, Biometrika 12, 134 (1918). Note that Isserlis’s theorem is often called “Wick’s theorem” in the physics community, although Wick’s work was three decades later.
  20. [20]Note that there are two definitions of the sinc(x) function which are widely used. Ours is the “math” definition which has its first positive root at x = π. The other is the “signal-processing” definition sin(πx)/πx, used for example in NUMPY, which has its first positive root at x = 1.
  21. [21]Note that the angular bin widths dγ will have different sizes due to the sin γ factor in |d(cos γ)| = sin γ dγ. For an approximately uniform distribution of pulsars on the sky, dγ will be correspondingly larger for γ near 0° and 180°.
  22. [22]We also allow for singular (i.e., positive-semidefinite) covariance matrices. Although such matrices have zero determinant, one can construct pseudoinverse matrices to invert them; see Sec. VIC for an example.
  23. [23]M. Schmelling, Averaging correlated data, Phys. Scr. 51, 676 (1995).DOI
  24. [24]A. Lazzarini, S. Bose, P. Fritschel, M. McHugh, T. Regimbau, K. Reilly, J. D. Romano, J. T. Whelan, S. Whitcomb, and B. F. Whiting, Optimal combination of signals from colocated gravitational wave interferometers for use in searches for a stochastic background, Phys. Rev. D 70, 062001 (2004).DOI
  25. [25]This can be easily proven by requiring that vᵀCv ≥ 0, with equality holding if and only if v = 0.
  26. [26]R. M. Wald, General Relativity (Chicago University Press, Chicago, 1984), 10.7208/chicago/9780226870373.001.0001.
  27. [27]It is also easy to show that the average values of the entries, where the average is done column by column, are column-independent. In fact, with the assumptions we have made, each matrix is characterized by the average value of its entries; the average value along any row or any column is the same.
  28. [28]L. Lewin, Dilogarithms and Associated Functions (Macdonald, London, 1958).
  29. [29]L. Lewin, Polylogarithms and Associated Functions (Elsevier Science Limited, New York, 1981).
  30. [30]D. Zagier, The dilogarithm function, in Frontiers in Number Theory, Physics, and Geometry II: On Conformal Field Theories, Discrete Groups and Renormalization, edited by P. Cartier, P. Moussa, B. Julia, and P. Vanhove (Springer, Berlin Heidelberg, 2007), pp. 3–65, 10.1007/978-3-540-30308-4_1.DOI
  31. [31]The dilogarithm is a fascinating function. According to the mathematician Don Zagier [30], “Almost all of its appearances in mathematics, and almost all the formulas relating to it, have something of the fantastical in them, as if this function alone among all others possessed a sense of humor.” Li₂(z) is conventionally defined by analytic continuation of the series (4.12) outside the unit circle, with a branch cut along the positive real axis from z = 1 to ∞. Note that Li₂(z) may be written in terms of Spence’s function, but beware that it is not the same as the logarithmic integral function li(x).
  32. [32]J. Gair, J. D. Romano, S. Taylor, and C. M. F. Mingarelli, Mapping gravitational-wave backgrounds using methods from CMB analysis: Application to pulsar timing arrays, Phys. Rev. D 90, 082001 (2014).DOI
  33. [33]E. Roebber and G. Holder, Harmonic space analysis of pulsar timing array redshift maps, Astrophys. J. 835, 21 (2017).DOI
  34. [34]W. Qin, K. K. Boddy, M. Kamionkowski, and L. Dai, Pulsar-timing arrays, astrometry, and gravitational waves, Phys. Rev. D 99, 063002 (2019).DOI
  35. [35]R. C. Bernardo and K.-W. Ng, Pulsar and cosmic variances of pulsar timing-array correlation measurements of the stochastic gravitational wave background, J. Cosmol. Astropart. Phys. 11 (2022) 046.
  36. [36]It may be possible to find the exact inverse of s_jk when N_bins is finite in terms of discrete orthogonal polynomials [74].
  37. [37]Our list of pulsars differs slightly from those that appear in the PTA collaborations’ recent data release [58–61]. For example, compared to our list, [58] removed J0614-3329 and added J2234+0944 and J1012+5307. These changes do not significantly modify our figures or results.
  38. [38]C. Cutler (private communication). Cutler suggested the idea of considering alternative detector types to understand this, analogous to employing a short dipole antenna to measure the squared amplitude of an electromagnetic field near Earth.
  39. [39]An equivalent of the “pulsar term” could be added to our imaginary one-arm LIGO detector. A second detector, formed with two mirrors, would be placed along the line of sight, near the pulsar. Its strain measurement would be relayed back to Earth at the speed of light via a radio or optical link, and then subtracted from the strain measured at the detector near Earth. This would appear as an additional factor in Eq. (6.1), identical to the term in square brackets in Eq. (6.2).
  40. [40]In [11], the pulsar term is controlled by a variable χ, which in this paper is used as a statistic. So, in this paper we use the variable η as a pulsar-term controlling switch.
  41. [41]X. Siemens (private communication).
  42. [42]M. Anholm, S. Ballmer, J. D. E. Creighton, L. R. Price, and X. Siemens, Optimal strategies for gravitational wave stochastic background searches in pulsar timing data, Phys. Rev. D 79, 084030 (2009).DOI
  43. [43]M. Rakhmanov, J. D. Romano, and J. T. Whelan, High-frequency corrections to the detector response and their effect on searches for gravitational waves, Classical Quantum Gravity 25, 184017 (2008).
  44. [44]This is true even in the absence of a pulsar term, as can be seen by setting η = 0 in (6.2).
  45. [45]E. H. Moore, On the reciprocal of the general algebraic matrix, in The Fourteenth Western Meeting of the American Mathematical Society, edited by A. Dresden. Moore’s unpublished paper was read at the 14’th Western Meeting of the American Mathematical Society. A. Dresden gives a summary of the results in the report, Bull. Am. Math. Soc. 26, 384 (1920).
  46. [46]E. H. Moore and R. W. Barnard (assisting), General Analysis (The American Philosophical Society, Philadelphia, 1935), pp. 197–209, https://hdl.handle.net/2027/wu.89010281279.
  47. [47]R. Penrose, A generalized inverse for matrices, Math. Proc. Cambridge Philos. Soc. 51, 406 (1955).DOI
  48. [48]O. M. Baksalary and G. Trenkler, The Moore-Penrose inverse: A hundred years on a frontline of physics research, Eur. Phys. J. H 46, 9 (2021).DOI
  49. [49]A. Das and W. S. Geisler, A method to integrate and classify normal distributions, J. Vis. 21, 1 (2021).DOI
  50. [50]S. J. Chamberlin, J. D. E. Creighton, X. Siemens, P. Demorest, J. Ellis, L. R. Price, and J. D. Romano, Time-domain implementation of the optimal cross-correlation statistic for stochastic gravitational-wave background searches in pulsar timing data, Phys. Rev. D 91, 044048 (2015).DOI
  51. [51]S. J. Vigeland, K. Islo, S. R. Taylor, and J. A. Ellis, Noise-marginalized optimal statistic: A robust hybrid frequentist-Bayesian statistic for the stochastic gravitational-wave background in pulsar timing arrays, Phys. Rev. D 98, 044003 (2018).DOI
  52. [52]Note that C⁻¹ab ≡ C⁻¹a,b is not obtained from C⁻¹ac,bd by setting c → a and d → b. Rather, it is obtained by inverting the Npul × Npul square matrix Cab ≡ Caa,bb, indexed by a = 1, …, Npul and b = 1, …, Npul.
  53. [53]A. Messiah, Quantum Mechanics, Dover Books on Physics Vol. 2 (Dover Publications, New York, 1999), Appendix C. I, “Clebsch-Gordon (C.-G.) Coefficients and ‘3j’ Symbols,” pp. 1054–1060.
  54. [54]E. W. Weisstein, Wigner 3j-Symbol, from MathWorld–A Wolfram Web Resource, https://mathworld.wolfram.com/Wigner3j-Symbol.html.
  55. [55]This argument is given more rigorously in Appendix D.
  56. [56]Although the weights wᵃᵇ are unity, the χ² statistic still takes into account the correlations among different pulsar pairs illustrated in Fig. 2. These enter χ² via the (inverse) covariance matrix.
  57. [57]The fractional uncertainty depends upon the ratio of ℏ/ℏ due to the “non-Gaussian” Ξ term in (8.31), which is proportional to ℏ⁸/ℏ⁸. In table entries and figure plots, we assume the binary inspiral GW model for which ℏ⁸/ℏ⁸ = 0.4741 (see Appendix F and Table III).
  58. [58]G. Agazie, A. Anumarlapudi, A. M. Archibald, Z. Arzoumanian, P. T. Baker, B. Bécsy, L. Blecha, A. Brazier, P. R. Brook, S. Burke-Spolaor et al., The NANO-Grav 15 yr data set: Evidence for a gravitational-wave background, Astrophys. J. Lett. 951, L8 (2023).DOI
  59. [59]J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A.-S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti et al., The second data release from the European Pulsar Timing Array. III. Search for gravitational-wave signals, Astron. Astrophys. (2023).DOI
  60. [60]D. J. Reardon, A. Zic, R. M. Shannon, G. B. Hobbs, M. Bailes, V. Di Marco, A. Kapur, A. F. Rogers, E. Thrane, J. Askew et al., Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array, Astrophys. J. Lett. 951, L6 (2023).DOI
  61. [61]H. Xu, S. Chen, Y. Guo, J. Jiang, B. Wang, J. Xu, Z. Xue, R. N. Caballero, J. Yuan, Y. Xu et al., Searching for the nano-hertz stochastic gravitational wave background with the Chinese Pulsar Timing Array data release I, Res. Astron. Astrophys. 23, 075024 (2023).DOI
  62. [62]P. M. Meyers, M. Vallisneri, and S. R. Taylor (private communication). These modifications were subsequently incorporated into the Enterprise software infrastructure and used in [58,59].
  63. [63]G. J. Feldman and R. D. Cousins, Unified approach to the classical statistical analysis of small signals, Phys. Rev. D 57, 3873 (1998).DOI
  64. [64]The variance (3.11) of the optimal Hellings and Downs correlation estimator does depend on ℏ⁴, which is not known a priori. So calculating the uncertainty in the estimate ρopt,j is more involved. Given a model for the spectral shape of the source, one can apply the “one-sigma, self-consistent” uncertainty approach described at the end of Sec. VII D to bound h² (and thus ℏ²). Similar constructions also hold in the presence of noise, as described Sec. IX D.
  65. [65]R. Blandford, R. Narayan, and R. W. Romani, Arrival-time analysis for a millisecond pulsar, J. Astrophys. Astron. 5, 369 (1984).DOI
  66. [66]R. W. Hellings, Pulsar timing and gravitational waves, in Relativistic Gravitational Experiments in Space: Proceedings of NASA Workshop, Annapolis, Maryland, 1988, NASA Conference Publication Vol. 3046 (1989), pp. 93–97, https://articles.adsabs.harvard.edu/pdf/1989NASCP3046...93H.
  67. [67]J. S. Hazboun, J. D. Romano, and T. L. Smith, Realistic sensitivity curves for pulsar timing arrays, Phys. Rev. D 100, 104028 (2019).DOI
  68. [68]Since PTA pulsar observations have a maximum cadence of about once per week, the upper limits of the integrals in (B2) should not extend to infinity, so the ratio h⁴/ℏ⁴ also depends upon this upper cutoff. However, the integrand is dominated by the low-frequency contributions, so this dependence is weak: Extending the upper cutoff to infinity has no significant effect.
  69. [69]Another simple example is to place a bin at the first zero of the Hellings and Downs curve, near γ = 49° with the corresponding choice μ_bin = 0. Then ⟨p_opt⟩ = 0 and σ²_opt = 0. This “beats” the cosmic variance, since the cosmic variance is always greater than zero.
  70. [70]The Santa Claus Principle holds that “there is no Santa Claus,” meaning that one does not get free gifts. Conversely, since here we are giving something up (angular resolution), we must gain something in exchange (higher SNR).
  71. [71]The cumulants κ_n of a random variable X are the Taylor series expansion coefficients in K(t) ≡ ln⟨e^{tX}⟩ = ∑_{n=0}^∞ κ_nt^n/n!.
  72. [72]Efficient evaluation can be done using PYTHON broadcasting and the function NUMPY.EINSUM().
  73. [73]These inversion transformations were exploited in [11], for example see the arguments following Eq. (G7) of that reference.
  74. [74]C. P. Neuman and D. I. Schonbach, Discrete (Legendre) orthogonal polynomials—A survey, Int. J. Numer. Methods Eng. 8, 743 (1974).DOI

Paper details

vlnc:W_4dfaysf7
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Physics and Astronomy

Not recorded for this paper · View statement →

Looking for related papers…
Contents