The Potential Impact of Noise Correlation in Next-generation Gravitational Wave Detectors

Abstract

Building upon the statistical formulation for parameter estimation in the presence of correlated noise proposed by Cireddu et al., we present the initial study to incorporate the effects of correlated noise into the analyses of various detector designs’ performance. We consider a two L-shaped detector configuration located in the European Union, and compare the expectation of parameter estimation between the non-colocated and a hypothetical colocated configurations. In our study, we posit the existence of low-frequency correlated noise within the $5\text{ Hz}$ to $10\text{ Hz}$ range for the colocated detector configuration, with a varying degree of correlation. In this specific detector setup, our observations indicate an enhancement in the precision of intrinsic parameter measurements as the degree of correlation increases. This trend suggests that higher degrees of noise correlation may beneficially influence the accuracy of parameter estimation. In particular, when the noise is highly correlated, the uncertainty on chirp mass decreases by up to $30\%$. The absence of an inter-European baseline does hinder the estimation of the extrinsic parameters. However, given a realistic global network with the additional detector located in the United States, the uncertainty of extrinsic parameters is significantly reduced. This reduction is further amplified as the degree of noise correlation increases. When the degree of noise correlation exceeds a certain level, the colocated configuration outperforms the non-colocated configuration. For instance, when the degree of correlation is high, the colocated configuration decreases the 90% credible area of sky location by up to $10\%$ compared to the non-colocated configuration. We conclude that the impact of noise correlation is not trivial and can potentially alter both the quantitative and qualitative outcomes in detector performance. We therefore recommend the inclusion of noise correlation for a comprehensive assessment of the design of third-generation gravitational wave detectors.

I Introduction

Since the first detection of gravitational-wave (GW) in 2015 [1], almost a hundred GW signals from binary mergers have been detected [2]. Such detections have significantly impacted various aspects of astronomy and fundamental physics [3, 4, 5]. The next-generation GW detectors promise to unlock mysteries in astrophysics, fundamental physics, and cosmology [6, 7, 8, 9, 10]. One such proposal is the Einstein Telescope (ET) [11], a third-generation ground-based GW detector.

One of the proposals for the ET consists of a colocated detector network. Such a setting gives rise to a novel challenge, namely, correlated detector noise. Due to the proximity of the colocated detectors [12], the noise in each of the detectors is expected to have a non-negligible correlation. The correlation is likely coming from seismic perturbations, Newtonian noise [13, 14] and magnetic fluctuation [15, 16]. While the literature recognizes the profound implications of correlated noise, existing research predominantly addresses its influence on the search of isotropic stocahstic GW background [13, 15, 14, 16]. A recent study [14] reports that the seismic correlations are significant over several hundreds of meters to a few kilometers in the frequency range 0.01 Hz to 40 Hz. Notably, the impact of correlated noise on GW transients has not been factored into the most recent extensive evaluation of various ET designs [17]. This omission highlights an important gap in our understanding and underscores the necessity for comprehensive studies that integrate the effects of correlated noise into the evaluation of detector configurations, thereby ensuring more accurate and robust design comparisons for the ET.

Conventional GW parameter estimation [18, 19] assumes uncorrelated detector noise; therefore, it is unclear how the presence of correlated noise will impact the parameter estimation. Recent work by Cireddu et al. [20] introduces a statistical formulation to account for correlated noise, laying down a foundation for investigating its impact on parameter estimation.

In this manuscript, we present a statistical framework to integrate the effects of correlated noise in the evaluation of different GW detector designs. Our approach involves numerical analyses aimed at understanding how correlated noise influences the accuracy of parameter estimation for GW transients. A colocated two L-shaped configuration is often used as a proxy for a triangular configuration in the literature. However, in the prescence of correlated noise, this is no longer generally the case. Nevertheless, to isolate the effects of correlated noise from benefits specifically due to detector geometry, such as null stream, our study focus exclusively on two L-shaped detectors under two distinct scenarios: one in which the detectors are colocated and another where they are non-colocated. We assess the impact on parameter estimation uncertainty through Fisher information matrix (FIM) analysis to make a global assessment of a population of sources. Remarkably, our findings reveal that, within the particular detector configuration examined, the presence of noise correlation leads to enhanced precision in parameter estimation compared to a comparable, yet non-colocated, detector network.

Hence, our research challenges the prevailing notion that noise correlation invariably compromises the scientific output of a detector network. Instead, it furnishes compelling evidence of the non-trivial impacts of noise correlation. This insight not only broadens our understanding of the interplay between detector configuration and noise characteristics but also underscores the necessity of re-evaluating design strategies in light of these findings.

II Review of Correlated Noise Statistics

In this section, we review the likelihood formulation in the presence of correlated noise presented in Ref. [20] in the context of two detectors.

The data of the two detectors is denoted as the time series $\boldsymbol{x}_{j}$, where $j=1,2$. These time series are sampled at regular intervals of $\Delta t$. We express the time series of multiple detectors as a matrix $\boldsymbol{X}$, where each row $j$ corresponds to the time series of the $j$-th detector. To compactly represent the spatial-temporal correlations in the noise data, we vectorize the matrix $\boldsymbol{X}$ into a single vector $\boldsymbol{x}$ as follows

$$\boldsymbol{x}:=\textrm{vec}(\boldsymbol{X}^{T})=\begin{bmatrix}\boldsymbol{x}_{1}\\ \boldsymbol{x}_{2}\end{bmatrix}.$$ (1)

$\tilde{\boldsymbol{x}}$ is understood to be the vectorization of the Fourier transform of the individual time series as follows

$$\tilde{\boldsymbol{x}}=\begin{bmatrix}\tilde{\boldsymbol{x}}_{1}\\ \tilde{\boldsymbol{x}}_{2}\end{bmatrix}\,.$$ (2)

The strain data is denoted as $\boldsymbol{d}$. The GW signal is denoted as $\boldsymbol{s}(\boldsymbol{\theta})$ given the source parameters $\boldsymbol{\theta}\in\mathbb{R}^{D}$.

Assume the noise follows the stationary Gaussian distribution, the likelihood function can be shown to be

$$p(\boldsymbol{d}|\boldsymbol{\theta})=\frac{1}{\sqrt{\det(\pi\boldsymbol{S}_{n}/(2\Delta f)})}\exp\left(-\frac{1}{2}\braket{\boldsymbol{d}-\boldsymbol{s}(\boldsymbol{\theta}),\boldsymbol{d}-\boldsymbol{s}(\boldsymbol{\theta})}\right)$$ (3)

where $\boldsymbol{S}_{n}$ is the spectral matrix defined as follows

$$\boldsymbol{S}_{n}=\begin{bmatrix}\boldsymbol{S}_{n}^{11}&\boldsymbol{S}_{n}^{12}\\ \boldsymbol{S}_{n}^{21}&\boldsymbol{S}_{n}^{22}\\ \end{bmatrix}$$ (4)

where

$$\boldsymbol{S}_{n}^{\ell m}[j,k]=2\Delta f\delta_{jk}\mathbb{E}\left[\tilde{n}_{\ell}[j]\tilde{n}_{m}^{*}[k]\right]$$ (5)

with $\delta_{jk}$ being the Kronecker delta function. $\boldsymbol{S}_{n}^{\ell\ell}$ is the one-sided power spectral density (PSD) of the noise in the $\ell$-th detector, and $\boldsymbol{S}_{n}^{\ell m}$ is the cross spectral density (CSD) of the noise in the $\ell$-th and the $m$-th detectors. One should be reminded that our definition of CSD differs from the conventional definition by a complex conjugate. $\braket{\boldsymbol{x},\boldsymbol{y}}$ is the generalized noise-weighed inner product defined as

$$\braket{\boldsymbol{x},\boldsymbol{y}}=4\Delta f\Re\left(\sum_{\ell,m=1}^{2}\sum_{k=k_{\textrm{low}}}^{k_{\textrm{high}}}(\boldsymbol{S}_{n}^{-1})^{\ell m}[k,k]\tilde{x}_{\ell}^{*}[k]\tilde{y}_{m}[k]\right),$$ (6)

where $\Delta f$ is the frequency resolution, and the inner product evaluated between $f_{\rm low}=\Delta fk_{\rm low}$ and $f_{\rm high}=\Delta fk_{\rm high}$, where $k_{\rm low}>0$ and $k_{\rm high}<\lfloor N/2\rfloor-1$, where $N$ is the length of the time series.

III Reduced Spread of the noise distribution in the presence of correlated noise

While the noise model in Eq. (3) accounts for the spatial-temporal correlation in the strain data of the ET, the impact of such correlation on parameter estimation has not yet been studied. Two scenarios are considered to detail the effects of correlated noise: 1) Detectors are colocated, experiencing correlated noise, e.g. seismic noise. The corresponding noise spectral matrix is denoted as

$$\boldsymbol{S}_{n}^{\textrm{corr}}=\begin{bmatrix}\boldsymbol{S}_{n}^{11}&\boldsymbol{S}_{n}^{12}\\ \boldsymbol{S}_{n}^{21}&\boldsymbol{S}_{n}^{22}\\ \end{bmatrix}\,.$$ (7)

2) Detectors are non-colocated, with the extent of separation not leading to a complete alteration in environmental noise sources. Nevertheless, they still have a similar noise properties, leading to the same noise covariance matrix as scenario 1, except for the vanishing non-diagonal elements. The corresponding noise covariance matrix is given by

$$\boldsymbol{S}_{n}^{\textrm{uncorr}}=\begin{bmatrix}\boldsymbol{S}_{n}^{11}&\boldsymbol{0}\\ \boldsymbol{0}&\boldsymbol{S}_{n}^{22}\\ \end{bmatrix}\,.$$ (8)

At this point, we can already obtain insights into the impact of correlated noise on the measurement process. The spread of the Gaussian distribution is proportional to the square root of the determinant of the covariance/spectral matrix. From Fischer’s inequality [21] we can deduce the following inequality

$$\sqrt{\det{\boldsymbol{S}_{n}^{\textrm{corr}}}}\leq\sqrt{\det{\boldsymbol{S}_{n}^{\textrm{uncorr}}}},$$ (9)

from which one can conclude that the presence of correlated noise (i.e. non-zero off-diagonal blocks in Eq. (7)) actually reduces of the spread of the likelihood function in Eq. (3). The propagation of the reduced spread to the precision of paramater estimation is however not trivial. We emphasize that one cannot immediately conclude from the reduced spread that it would result in a more precise parameter estimation since a more careful analysis would reveal that the specific structure of the signal gradient is another important factor. It is however tempting to speculate that the impacts of noise correlation is not always negative.

IV Impact of correlated noise on parameter estimation

To compare the performance, we assess how the uncertainty of parameter estimation changes between the two detector configurations when observing the same source. This comparison is characterized by the ratio of the square root of the determinant of the covariance matrix of the posterior distributions, averaged over noise realizations:

$$r_{\textrm{uncorr}}^{\textrm{corr}}:=\frac{\mathbb{E}_{\boldsymbol{d}^{\rm corr}|\boldsymbol{\theta_{\rm true}}}\left[\sqrt{\det\left(\boldsymbol{\Sigma}_{\boldsymbol{\theta}|\boldsymbol{S}_{n}^{\textrm{corr}}}\right)}\right]}{\mathbb{E}_{\boldsymbol{d}^{\rm uncorr}|\boldsymbol{\theta_{\rm true}}}\left[\sqrt{\det\left(\boldsymbol{\Sigma}_{\boldsymbol{\theta}|\boldsymbol{S}_{n}^{\textrm{uncorr}}}\right)}\right]}=\frac{\int\sqrt{\det\left(\boldsymbol{\Sigma}_{\boldsymbol{\theta}|\boldsymbol{S}_{n}^{\textrm{corr}}}\right)}p(\boldsymbol{d}^{\textrm{corr}}|\boldsymbol{\theta}_{\textrm{true}},\boldsymbol{S}_{n}^{\textrm{corr}})\textrm{d}\boldsymbol{d}^{\textrm{corr}}}{\int\sqrt{\det\left(\boldsymbol{\Sigma}_{\boldsymbol{\theta}|\boldsymbol{S}_{n}^{\textrm{uncorr}}}\right)}p(\boldsymbol{d}^{\textrm{uncorr}}|\boldsymbol{\theta}_{\textrm{true}},\boldsymbol{S}_{n}^{\textrm{uncorr}})\textrm{d}\boldsymbol{d}^{\textrm{uncorr}}}$$ (10)

where $\boldsymbol{\Sigma}_{\boldsymbol{\theta}|\boldsymbol{S}_{n}^{\textrm{corr}/\textrm{uncorr}}}$ is the posterior covariance of the model parameters $\boldsymbol{\theta}$ given a noise realization with the noise spectral matrix $\boldsymbol{S}_{n}^{\textrm{corr}/\textrm{uncorr}}$, and $\boldsymbol{\theta}_{\textrm{true}}$ is the true parameters.

Evaluating the quantity $r_{\textrm{uncorr}}^{\textrm{corr}}$ is challenging for general setups. Instead, we turn to the FIM analysis. The inverse of the FIM is known to characterize the covariance of the posterior distribution in the high SNR limit [22]. For a deterministic signal $\boldsymbol{s}(\boldsymbol{\theta})$, the FIM is given by

$$\mathcal{I}_{jk}^{\textrm{corr/uncorr}}(\boldsymbol{\theta})=\left.\braket{\partial_{\theta_{j}}\boldsymbol{s}(\boldsymbol{\theta}),\partial_{\theta_{k}}\boldsymbol{s}(\boldsymbol{\theta})}\right|_{\boldsymbol{S}_{n}=\boldsymbol{S}_{n}^{\textrm{corr/uncorr}}}$$ (11)

where $\partial_{\theta_{j}}\boldsymbol{s}(\boldsymbol{\theta})$ is the derivative of the signal with respect to the $j$-th model parameter. In the high SNR limit, the ratio $r_{\textrm{uncorr}}^{\textrm{corr}}$ is then reduced to:

$$\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}=\frac{1/\sqrt{\det(\boldsymbol{\mathcal{I}}^{\textrm{corr}}(\boldsymbol{\theta}_{\rm true}))}}{1/\sqrt{\det(\boldsymbol{\mathcal{I}}^{\textrm{uncorr}}(\boldsymbol{\theta}_{\rm true}))}}\,,$$ (12)

$\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}<1$ is indicating an improvement in parameter estimation precision, and vice versa. The mathematical properties of $\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}$ is further studied in the Appendix.

IV.1 Simulation

We consider two scenarios. In the first one, we simulate a two L-shaped detectors located in the European Union (EU) with a $45^{\circ}$ relative offset and identical PSDs, the ET-D sensitivity $\boldsymbol{S}_{n}^{\textrm{ET-D}}$ [23]. Within this scenario, we compare two configurations, i) the two detectors are non-colocated, with one placed at Limburg, and one at Sardinia, featuring uncorrelated noise, denoted by the spectral matrix $\boldsymbol{S}_{n}^{\textrm{uncorr}}$, while for ii) the two detectors are colocated at Limburg, featuring correlated noise, represented by the spectral matrix $\boldsymbol{S}_{n}^{\textrm{corr}}$. The orientation, elevation, and arm length are the same across the two setups. The locations and orientations of Ref. [24] are used for the EU detectors.

In the second scenario, we consider a simulated global network with an L-shaped detector located at Hanford, referred as the US detector in the text, with the 40 km sensitivity of the Cosmic Explorer $\boldsymbol{S}_{n}^{\textrm{CE-40km}}$ [25, 26, 27].

In the simulation, we use the frequency cutoff to be $[5,1024)$ Hz. In the absence of a faithful correlated (seismic) noise model, we assume the noise components under $10$ Hz are correlated with a correlation coefficient of $\alpha$ . Consequently, the CSD for the EU detectors are taken to be $\alpha S_{n}^{\textrm{ET-D}}(f)$ for $f\leq 10\textrm{ Hz}$ and $0$ otherwise.

Therefore, the spectral matrices for a sole EU network is given by

$$\boldsymbol{S}_{n}^{\textrm{uncorr}} =\begin{bmatrix}\boldsymbol{S}_{n}^{\textrm{ET-D}}&\boldsymbol{0}\\ \boldsymbol{0}&\boldsymbol{S}_{n}^{\textrm{ET-D}}\\ \end{bmatrix},$$ (13)

$$\boldsymbol{S}_{n}^{\textrm{corr}} =\begin{bmatrix}\boldsymbol{S}_{n}^{\textrm{ET-D}}&\boldsymbol{\alpha}\boldsymbol{S}_{n}^{\textrm{ET-D}}\\ \boldsymbol{\alpha}\boldsymbol{S}_{n}^{\textrm{ET-D}}&\boldsymbol{S}_{n}^{\textrm{ET-D}}\\ \end{bmatrix},$$

and for a global network with the inclusion of a US detector, the matrices are given by

$$\boldsymbol{S}_{n}^{\textrm{uncorr}} =\begin{bmatrix}\boldsymbol{S}_{n}^{\textrm{ET-D}}&\boldsymbol{0}&\boldsymbol{0}\\ \boldsymbol{0}&\boldsymbol{S}_{n}^{\textrm{ET-D}}&\boldsymbol{0}\\ \boldsymbol{0}&\boldsymbol{0}&\boldsymbol{S}_{n}^{\textrm{CE-40km}}\end{bmatrix},$$ (14)

$$\boldsymbol{S}_{n}^{\textrm{corr}} =\begin{bmatrix}\boldsymbol{S}_{n}^{\textrm{ET-D}}&\boldsymbol{\alpha}\boldsymbol{S}_{n}^{\textrm{ET-D}}&\boldsymbol{0}\\ \boldsymbol{\alpha}\boldsymbol{S}_{n}^{\textrm{ET-D}}&\boldsymbol{S}_{n}^{\textrm{ET-D}}&\boldsymbol{0}\\ \boldsymbol{0}&\boldsymbol{0}&\boldsymbol{S}_{n}^{\textrm{CE-40km}}\end{bmatrix},$$

where $\alpha_{jk}=\alpha\delta_{jk}H(10{\text{ Hz}}-k\Delta f)$, with $H(x)$ being the Heaviside step function. Given the current uncertainties surrounding the level of noise correlation, we have chosen to simulate a range of correlation coefficients, $\alpha$, from $0$ $0.9$ in increments of $0.1$. Recent measurement using seismometers at different sites and separations showed that the coherence can range between $\sim 0.01$ and $\sim 0.5$ at the $90^{\rm{th}}$ percentile [14]. The coherence cannot be trivially mapped to $\alpha$ since it depends on the exact detector layout. Nevertheless, it gives us a rough estimate of $\alpha$ in the range between $\sim 0.005$ and $\sim 0.25$, which is estimated by dividing the coherence by $2$ assuming only the end test mass is subject to noise correlation.

A catalog of $200$ aligned spin binary black hole (BBH) merger transients using the waveform approximant IMRPhenomXPHM [28] is generated. The simulated catalog’s source parameter distributions are outlined in Table 1. Subsequently, we perform FIM analyses on this catalog under two different setups.

Our comparative analysis focuses on evaluating the standard deviation of individual source parameters and the spread of the posterior distributions between the two configurations under the two scenarios. It is essential to note that the FIM only approximates the posterior primarily for high signal-to-noise ratio (SNR) signals. Thus, we set a fiducial SNR cutoff, ensuring signals possess an SNR above $50$. The SNR for this cutoff is determined by the non-colocated configuration.

Parameter	Distribution	Range
Primary mass $m_{1}$	Uniform	($5$, $100$) ${\rm M}_{\odot}$
Secondary mass $m_{2}$	Uniform	($5$, $100$) ${\rm M}_{\odot}$
Primary aligned spin $\chi_{1}$	Uniform	($0$, $0.99$)
Secondary aligned spin $\chi_{2}$	Uniform	($0$, $0.99$)
Inclination $\theta_{\rm JN}$	Sine	($0$, $\pi$)
Polarization angle $\psi$	Uniform	($0$, $\pi$)
Right ascension $\alpha$	Uniform	($0$, $2\pi$)
Declination $\delta$	Cosine	(-$\pi/2$, $\pi/2$)
Luminosity distance $d_{\rm L}$	Uniform in	($100$, $10000$) Mpc
	source frame

Table 1: Distribution of the source parameters in the simulation. The definition of the parameters can be referred to Ref. [29]. Cosmology from Ref. [30] is used for the luminosity distance $d_{\rm L}$ prior.

IV.2 Results

Our findings reveal an improvement in measurement precision of the parameters for the vast majority of the injections when comparing the hypothetical colocated-correlated configuration with the non-colocated configuration.

	EU-only network	EU-US network
$\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}$ (all)	$0.75_{-0.66}^{+19.72}$	$0.20_{-0.06}^{+0.13}$
$\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}$ (intrinsic)	$0.38_{-0.17}^{+0.41}$	$0.65_{-0.16}^{+0.18}$
$\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}$ (extrinsic)	$1.50_{-1.26}^{+44.17}$	$0.29_{-0.07}^{+0.10}$

Table 2: Median and $1\sigma$ uncertainty range of $\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}$, assessed across all, intrinsic-only, and extrinsic-only parameters comparing the colocated configuration with a correlation coefficient of $0.9$ and the non-colocated configuration. $\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}<1$ is indicating an improvement in parameter estimation precision, and vice versa. Notably, the median values for intrinsic parameters consistently remain below 1, indicating increased measurement precision. In the absence of a US detector, the $\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}$ values for extrinsic parameters surpass $1$. However, with a broader network with a US detector, a marked reduction below $1$ suggests superior localization precision with colocated EU detectors in a global detector network.

In Figure 1, we present a comparison of the performance variations between colocated and non-colocated European detector configurations as influenced by the correlation coefficient of detector noise. We show the results of the chirp mass as an example of the intrinsic parameters, and the $90\%$ credible area of sky localization as an example of the extrinsic parameters. In the left panel, we observe a decline in the ratio of the standard deviations for the chirp mass $\mathcal{M}_{c}$ as the correlation coefficient increases. This trend is consistent across both EU-only and EU-US detector networks, implying an improvement in the precision of parameter estimation with greater noise correlation. The right panel presents the variation of the ratio of the $90\%$ credible area for sky localization. Here, the sole EU network exhibits poorer performance in the colocated configuration compared to the non-colocated configuration, as anticipated due to detector separation in the non-colocated case. Remarkably, the inclusion of a US detector markedly improves performance, with the ratio nearing 1 and further reducing to 0.9 as the noise correlation reaches 0.9. This showcases that the improvement from the intrinsic parameters outweighs the additional information gained from an inter-European baseline for sky localization when an intercontinental baseline is present. The distribution of the ratio for other parameters follows the same trend, and is presented in the Appendix.

To summarize the results of the FIM analysis, Table 2 presents the median and the $1\sigma$ interval of the ratio of posterior distributions’ spreads $\hat{r}^{\rm corr}_{\rm uncorr}$ in the extreme scenario when the noise is highly correlated ($\alpha=0.9$). For both cases, the $\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}$ for intrinsic parameters are consistently lower than $1$, signifying an improvement in parameter estimation precision. In the absence of the US detector, the $\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}$ for extrinsic parameters is larger than $1$. However, with a broader network with the US detector, a marked reduction below $1$ suggests improved precision in extrinsic parameters with the colocated configuration in a global detector network.

Figure 1: This figure illustrates the variation in the medium of the ratio of the standard deviations of parameters between colocated and non-colocated EU detector configurations as the correlation coefficient of detector noise increases. The figure on the left displays the changes in the ratio for the chirp mass, $\mathcal{M}_{c}$. Both in the EU-only and EU-US detector networks, the ratio decreases with an increasing correlation coefficient, suggesting enhanced precision in parameter estimation. The right panel presents the variation in the ratio of the $90\%$ credible area for sky localization. For the EU network alone, the colocated configuration underperforms relative to the non-colocated configuration. However, with the integration of a US detector, the ratio approaches 1 and further declines to 0.9 as the noise correlation increases to 0.9, indicating improved localization accuracy with higher noise correlation.

Our results highlight the significant, yet positive, impact of noise correlation on parameter estimation accuracy within specific colocated GW detector configurations. This reveals the non-trivial impacts of noise correlation on the precision of parameter estimation. This underscores the important need to consider noise correlation in the design of future GW observatories, revealing its potential to improve detection capabilities.

V Concluding Remarks

Based on the statistical foundation presented in Cireddu et al. [20], this work represents an initial endeavor to integrate the effects of correlated noise into the analysis of various detector designs’ performance. We compared the precision of parameter estimation between a hypothetical colocated two L-shaped, correlated detector with different levels of noise correlation with a non-colocated, uncorrelated one. Our findings reveal a marked discrepancy in detector performance upon accounting for correlated noise impacts. Specifically, within the considered colocated detector configuration, rather than a hindrance, the presence of correlated noise actually boosts our ability to extract science.

For an EU-US network setting in the highly correlated ($\alpha\approx 0.9$) scenario, the colocated-correlated detectors offers a $10-20\%$ improvement across parameters and a $\sim 10\%$ improvement for the sky localization, as compared to the non-colocated detectors.

A similar study comparing different detector configurations was done in Ref. [8]. Our results and Ref. [8], although partially consistent, cannot be compared directly with each other because of the different configurations considered. However, our investigations affirm that noise correlation significantly influences detector design optimization, potentially altering both quantitative and qualitative outcomes.

In conclusion, despite its importance, the influence of correlated noise remains overlooked in current discussions on future GW detector design. We recommend a thorough consideration of noise correlation in future research endeavors, emphasizing its potential to reshape our understanding and approach to detector design optimization in the pursuit of advancing GW astronomy.

VI Acknowledgement

We thank Rico K. L. Lo for insightful discussions. P.T.H.P is supported by the research program of the Netherlands Organization for Scientific Research (NWO). M.W. is supported by the Research Foundation - Flanders (FWO) through grant no. 11POK24N.

VII Software

The analysis and simulations presented in this work were conducted using a combination of established software packages and custom scripts. The figures are produced using Matplotlib [31, 32] and seaborn [33]. Numerical operations and array manipulations are performed using NumPy [34] and SciPy [35]. To evaluate the Fisher information matrix, gwbench[36] is used with extension to analyze scenarios involving correlated noise.

VIII Data availability

All our data will be available in zenodo.org after publication.

GW

gravitational-wave

ET

Einstein Telescope

SNR

signal-to-noise ratio

LF

low-frequency

HF

high-frequency

ASD

amplitude spectral density

PSD

power spectral density

CSD

cross spectral density

BBH

binary black hole

FIM

Fisher information matrix

PE

parameter estimation

CE

Cosmic Explorer

Appendix A Determinant of the Inverse Fisher Information Matrix

The FIM is defined as

$$\mathcal{I}_{jk}(\boldsymbol{\theta})=\mathbb{E}_{\boldsymbol{d}|\boldsymbol{\theta}}\left[(\partial_{\theta_{j}}\log p(\boldsymbol{d}|\boldsymbol{\theta}))(\partial_{\theta_{k}}\log p(\boldsymbol{d}|\boldsymbol{\theta}))\right]\,.$$ (15)

With the time domain likelihood function defined in Ref. [20], one can show that

$$\boldsymbol{\mathcal{I}}(\boldsymbol{\theta})=(\nabla\boldsymbol{s}(\boldsymbol{\theta}))^{T}\boldsymbol{\Sigma}_{n}^{-1}\nabla\boldsymbol{s}(\boldsymbol{\theta})$$ (16)

where $\nabla\boldsymbol{s}(\boldsymbol{\theta})$ is the Jacobian matrix of $\boldsymbol{s}(\boldsymbol{\theta})$ and $\boldsymbol{\Sigma}_{n}$ is the noise covariance matrix. Perform a singular value decomposition on $\nabla\boldsymbol{s}(\boldsymbol{\theta})$, we have

$$\nabla\boldsymbol{s}(\boldsymbol{\theta})=\boldsymbol{U}(\boldsymbol{\theta})\boldsymbol{S}(\boldsymbol{\theta})(\boldsymbol{V}(\boldsymbol{\theta}))^{T}$$ (17)

where $\boldsymbol{U}(\boldsymbol{\theta})\in\mathbb{R}^{MN\times MN}$ and $\boldsymbol{V}(\boldsymbol{\theta})\in\mathbb{R}^{D\times D}$ are orthogonal matrices, and $\boldsymbol{S}(\boldsymbol{\theta})\in\mathbb{R}^{MN\times D}$ is a rectangular matrix with the singular values in the principal diagonal. $M$ is the number of detectors, $N$ is the number of time bins, and $D$ is the dimensionality of the parameter space $\boldsymbol{\Theta}$. Substitute the decomposition into $\hat{r}_{\textrm{uncorr}}^{\textrm{corr}}$ defined in Eq. (12), we have

$$(\hat{r}_{\textrm{uncorr}}^{\textrm{corr}})^{2}$$ (18)

$$=\left.\frac{\det\left((\nabla_{\boldsymbol{\theta}}\boldsymbol{s}(\boldsymbol{\theta}))^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{uncorr}})^{-1}(\nabla_{\boldsymbol{\theta}}\boldsymbol{s}(\boldsymbol{\theta}))\right)}{\det\left((\nabla_{\boldsymbol{\theta}}\boldsymbol{s}(\boldsymbol{\theta}))^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{corr}})^{-1}(\nabla_{\boldsymbol{\theta}}\boldsymbol{s}(\boldsymbol{\theta}))\right)}\right|_{\boldsymbol{\theta}=\boldsymbol{\theta}_{\textrm{true}}}$$ (19)

$$=\left.\frac{\det\left(\boldsymbol{V}(\boldsymbol{\theta})\boldsymbol{S}(\boldsymbol{\theta})^{T}\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{uncorr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\boldsymbol{S}(\boldsymbol{\theta})\boldsymbol{V}(\boldsymbol{\theta})^{T}\right)}{\det\left(\boldsymbol{V}(\boldsymbol{\theta})\boldsymbol{S}(\boldsymbol{\theta})^{T}\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{corr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\boldsymbol{S}(\boldsymbol{\theta})\boldsymbol{V}(\boldsymbol{\theta})^{T}\right)}\right|_{\boldsymbol{\theta}=\boldsymbol{\theta}_{\textrm{true}}}$$ (20)

$$=\left.\frac{\left(\det\boldsymbol{V}(\boldsymbol{\theta})\right)^{2}\det\left(\boldsymbol{S}(\boldsymbol{\theta})^{T}\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{uncorr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\boldsymbol{S}(\boldsymbol{\theta})\right)}{\left(\det\boldsymbol{V}(\boldsymbol{\theta})\right)^{2}\det\left(\boldsymbol{S}(\boldsymbol{\theta})^{T}\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{corr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\boldsymbol{S}(\boldsymbol{\theta})\right)}\right|_{\boldsymbol{\theta}=\boldsymbol{\theta}_{\textrm{true}}}$$ (21)

$$=\left.\frac{\det\left(\hat{\boldsymbol{S}}(\boldsymbol{\theta})^{T}\left[\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{uncorr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\right]_{\alpha,\alpha}\hat{\boldsymbol{S}}(\boldsymbol{\theta})\right)}{\det\left(\hat{\boldsymbol{S}}(\boldsymbol{\theta})^{T}\left[\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{corr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\right]_{\alpha,\alpha}\hat{\boldsymbol{S}}(\boldsymbol{\theta})\right)}\right|_{\boldsymbol{\theta}=\boldsymbol{\theta}_{\textrm{true}}}$$ (22)

$$=\left.\frac{\left(\det\hat{\boldsymbol{S}}(\boldsymbol{\theta})\right)^{2}\det\left(\left[\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{uncorr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\right]_{\alpha,\alpha}\right)}{\left(\det\hat{\boldsymbol{S}}(\boldsymbol{\theta})\right)^{2}\det\left(\left[\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{corr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\right]_{\alpha,\alpha}\right)}\right|_{\boldsymbol{\theta}=\boldsymbol{\theta}_{\textrm{true}}}$$ (23)

$$=\left.\frac{\det\left(\left[\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{uncorr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\right]_{\alpha,\alpha}\right)}{\det\left(\left[\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{corr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\right]_{\alpha,\alpha}\right)}\right|_{\boldsymbol{\theta}=\boldsymbol{\theta}_{\textrm{true}}}$$ (24)

$$=\left.\frac{1/\det\left(\left[\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{corr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\right]_{\alpha,\alpha}\right)}{1/\det\left(\left[\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{uncorr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\right]_{\alpha,\alpha}\right)}\right|_{\boldsymbol{\theta}=\boldsymbol{\theta}_{\textrm{true}}}$$ (25)

where $\alpha=\{1,2,...,D\}$ denotes the index set, and $\left[{}\cdot{}\right]_{\alpha,\alpha}$ represents the submatrix with the index set $\alpha$. It is noteworthy that the spread of the posterior covariance matrix is inversely proportional to $\sqrt{\det\left(\left[\boldsymbol{U}(\boldsymbol{\theta})^{T}(\boldsymbol{\Sigma}_{n}^{\textrm{corr/uncorr}})^{-1}\boldsymbol{U}(\boldsymbol{\theta})\right]_{\alpha,\alpha}\right)}$. From Eq. (21) to Eq. (22), we replace $\boldsymbol{S}(\boldsymbol{\theta})$ with $\hat{\boldsymbol{S}}(\boldsymbol{\theta})$ defined as a square matrix with the singular values as the diagonal entries.

The geometrical interpretation is clear: the uncertainty in parameter estimation is inversely proportional to the spread of the subspace spanned by the $D$ principal vectors in the inverse noise covariance matrix. These vectors represent the directions in which the signal varies the most. A larger subspace spread corresponds to reduced uncertainty in estimating the model parameters.

In Eq. (25), one shall notice that the singular values are cancelled, implying that the ratio is independent of the strength of the signal. Moreover, the uncertainty is determined by the subspace spanned by the signal gradient vectors. Provided that the subspaces $\boldsymbol{U}$ are the same in the numerator and the denominator, the ratio of the uncertainty is largely dependent on the spread of the noise covariance matrix. Thus, Eq. (25) suggests that the reduced spread of the noise distribution is likely to lead to a smaller uncertainty in parameter estimation of GW transients.

Appendix B Distribution of the ratio of the standard deviations of all parameters

In Figure 2, we present the distributions of the ratio of the standard deviations of the individual parameters comparing the colocated configuration with different correlation coefficients and the non-colocated configuration. The left-hand side of the violin plots showcases outcomes from Europe-only network, while the right-hand side illustrates findings from an Europe-US network.

In the Europe-only network, as the correlation coefficient increases, the measurement uncertainty of the intrinsic parameters decreases. In particular, when the noise is highly correlated ($\alpha=0.9$), the uncertainty of chirp mass $\mathcal{M}_{c}$ decreases by $\sim 30\%$, with a $\sim 10\%$ decrease for the symmetric mass ratio $\eta$ and the component aligned spins $\chi_{1z,2z}$. For the extrinsic parameters, as expected, the non-colocated network is offering better estimation. Yet, it would be more realistic to consider a global detector network.

With the inclusion of the US detector at Hanford, as shown on the right-hand side of the plots, the colocated configuration improves the measurement precision of all parameters compared to the non-colocated configuration. The degree of improvement increases with the correlation coefficient. Notably, when the noise is highly correlated ($\alpha=0.9$), the uncertainties for most of the parameters decrease by $\sim 10-20\%$, and the area of the $90\%$ credible interval of the sky location $\Omega_{90\%}$ decreases by $\sim 50\%$. In Table 3 and Table 4, the median and $1\sigma$ intervals of the same distributions in the highly correlated scenario ($\alpha=0.9$) are shown.

Figure 2: The distribution of the ratio of the standard deviations of the individual parameters of the posterior distributions, with different correlation coefficients, including the chirp mass $\mathcal{M}_{c}$, symmetric mass ratio $\eta$, aligned spins $\chi_{1,2z}$, coalescence phase $\phi_{c}$, coalescence time $t_{c}$, polarization angle $\psi$, inclination angle $\iota$, and luminosity distance $d_{L}$, comparing the colocated configuration and the non-colocated configuration of the EU detectors. LHS of the plots presents the results without including the US detector. RHS of the plots presents the results including the US detector. The ratio of the $90\%$ credible area of the sky localization $\Omega_{90\%}$ is presented in the last column.

EU-only
Detector network	$\mathcal{M}_{c}$	$\eta$	$\chi_{1z}$	$\chi_{2z}$	$\phi_{c}$	$t_{c}$	$\psi$	$\iota$	$d_{L}$	$\Omega_{90\%}$
colocated-correlated	$0.68_{-0.10}^{+0.22}$	$0.82_{-0.19}^{+0.15}$	$0.77_{-0.19}^{+0.19}$	$0.77_{-0.18}^{+0.20}$	$1.22_{-0.41}^{+1.52}$	$1.41_{-0.62}^{+8.73}$	$1.63_{-0.75}^{+5.58}$	$0.92_{-0.12}^{+0.34}$	$1.16_{-0.37}^{+3.54}$	$7.65_{-6.68}^{+337.72}$
colocated-uncorrelated	$1.00_{-0.06}^{+0.05}$	$1.00_{-0.06}^{+0.06}$	$1.01_{-0.07}^{+0.05}$	$1.00_{-0.06}^{+0.07}$	$1.37_{-0.41}^{+1.73}$	$1.55_{-0.64}^{+10.06}$	$1.81_{-0.81}^{+6.54}$	$1.04_{-0.09}^{+0.37}$	$1.34_{-0.43}^{+4.30}$	$9.19_{-7.91}^{+378.12}$
non-colocated	1	1	1	1	1	1	1	1	1	1

Table 3: The median and the 1$\sigma$ interval of the ratio of the standard deviations of the individual parameters in different detector networks compared to the non-colocated configuration, without the presence of the US detector. The colocated-uncorrelated EU network acts as a control to compare with colocated-correlated EU network with a correlation coefficient $\alpha$ of $0.9$ to examine the impact of the presence of correlated noise.

EU-US
Detector network	$\mathcal{M}_{c}$	$\eta$	$\chi_{1z}$	$\chi_{2z}$	$\phi_{c}$	$t_{c}$	$\psi$	$\iota$	$d_{L}$	$\Omega_{90\%}$
colocated-correlated	$0.85_{-0.11}^{+0.11}$	$0.87_{-0.13}^{+0.10}$	$0.86_{-0.14}^{+0.10}$	$0.86_{-0.13}^{+0.11}$	$0.82_{-0.20}^{+0.11}$	$0.93_{-0.13}^{+0.13}$	$0.85_{-0.14}^{+0.13}$	$0.89_{-0.12}^{+0.07}$	$0.89_{-0.13}^{+0.07}$	$0.90_{-0.21}^{+0.23}$
colocated-uncorrelated	$1.00_{-0.03}^{+0.03}$	$0.99_{-0.02}^{+0.03}$	$1.00_{-0.03}^{+0.03}$	$1.00_{-0.03}^{+0.03}$	$0.95_{-0.08}^{+0.08}$	$1.03_{-0.06}^{+0.10}$	$0.99_{-0.12}^{+0.14}$	$0.98_{-0.08}^{+0.05}$	$0.98_{-0.09}^{0.06}$	$1.08_{-0.10}^{+0.13}$
non-colocated	1	1	1	1	1	1	1	1	1	1

Table 4: The median and the 1$\sigma$ interval of the ratio of the standard deviations of the individual parameters in different detector networks compared to the non-colocated configuration, with the presence of the US detector. The colocated-uncorrelated EU network acts as a control to compare with colocated-correlated EU network with a correlation coefficient $\alpha$ of $0.9$

to examine the impact of the presence of correlated noise.