Assessing the similarity of continuous gravitational-wave signals to narrow instrumental artifacts

Continuous gravitational-wave signals (CWs) are long-lasting quasi-monochromatic gravitational-wave signals expected to be emitted by rapidly-rotating non-axisymmetric neutron stars. Depending on the rotational frequency and sky location of the source, certain CW signals may behave in a similar manner to narrow-band artifacts present in ground-based interferometric detectors. Part of the detector-characterisation tasks in the current generation of interferometric detectors (Advanced LIGO, Advanced Virgo, and KAGRA) aim at understanding the origin of these narrow artifacts, commonly known as ``spectral lines''. It is expected that similar tasks will continue after the arrival of next-generation detectors (e.g. Einstein Telescope and Cosmic Explorer). Typically, a fraction of the observed lines in a given detector can be associated to one or more instrumental causes; others, however, have an unknown origin. In this work, we assess the similarity of CW signals to spectral lines in order to understand whether a CW may be mistaken for a noise artifact. Albeit astrophysically unlikely, our results do not rule out the possibility of a CW signal being visible in the detector's power spectrum.


Introduction
Continuous gravitational waves (CWs) are long lasting gravitational wave signals (GW) whose detection remains, so far, unattained [1,2].Amongst the expected sources we find rapidly-spinning non-axisymmetric neutron stars (NS) [3], but also other more exotic ones such as evaporating boson clouds formed around spinning black holes [4,5], or planetary-mass compact binary systems [6,7].Their expected frequency lies on the audible band of the GW spectrum, which makes them prime targets to be detected using the current generation of advanced ground-based interferometric detectors (Advanced LIGO [8], Advanced Virgo [9] and KAGRA [10]), as well as third-generation detectors (3G) (Einstein Telescope (ET) [11] and Cosmic Explorer (CE) [12]).
The frequency evolution of a CW signal as emitted by a rapidly-spinning NS can usually be described using a Taylor series [13] where f 0 is the initial CW frequency at a fiducial reference time t ref and f k≥1 are the spindown parameters.The number of terms to include in Eq. ( 1) is generally dependent on the age of the source: older objects tend to spin down more slowly and thus require a lower number of terms to be accurately modelled [14][15][16].Typical searches for unknown CW sources include one or two spindown terms in order to remain computationally affordable [17].
Searches for younger objects, such as newborn neutron stars, try to detect CW emission on a much shorter timescale (hours to days) and tend to be conducted using the general torque equation as discussed in [18][19][20].

arXiv:2403.03027v1 [gr-qc] 5 Mar 2024
Upon arrival to the detector, the frequency of a CW signal is Doppler-modulated due to the detector's motion around the Solar System Barycenter (SSB) where c is the speed of light, n is the sky position of the source in the SSB, and ⃗ v is the velocity of the detector with respect to the SSB.The detector's velocity can be further expanded into two components ⃗ v = ⃗ v o + ⃗ v r , where v o ≈ 10 −4 c is the Earth's orbital speed and v r ≈ 10 −6 c is the Earth's rotational speed.This translates into two time-dependent and sky-location-dependent frequency modulations.The magnitude of the Doppler modulation is proportional to ⃗ v • n.Since the time-averaged orientation of ⃗ v is close to the ecliptic plane, sky positions near the poles (for which ⃗ v • n ≈ 0) are expected to suffer the smallest frequency modulation [21,22].Ground-based interferometric detectors are affected by a wide variety of noise sources that behave in a similar manner to GW signals [23].For example, short noise transients ("glitches") are known to overlap with GWs produced in the coalescense of two compact objects [24][25][26][27].If unaddressed, noise artifacts significantly degrade the sensitivity of a GW search, in the sense that astrophysically interestiong candidates will likely be recovered at a higher false alarm probability (i.e.lower significance level).
CW searches are mainly affected by spectrally narrow and persistent noise artifacts, also known as "lines" due to their line-like appearance in a power spectrum [28].Lines are initially identified by inspecting the power spectral density of a detector, either manually or using a peak-finding algorithm, and then further classified into "instrumental lines" if they can be associated to an instrumental cause, or "unknown lines" if otherwise.Well-known instrumental lines include, for example, lines at multiples of the 60 Hz electrical power frequency in the Advanced LIGO detectors or at multiples of the the resonance frequency of the mirror suspensions (about 500 Hz).Lists of instrumental and unknown lines in the third observing run of the Advanced LIGO detectors (O3) can be found in [29,30], respectively.
Lines tend to cause a high number of outliers in CW searches, as their relatively strong power and persistence with respect to the detector's background PSD causes them to look like excess power for an analysis.This is specially problematic in broad searches for unknown sources, such as all-sky searches, as their CW models tend to be particularly sensitive to lines due to the use of low coherence times [17] 1 .To deal with such an elevated number of candidates, CW searches follow two main approaches: First, CW candidates can be tested using a broad suite of consistency vetoes in order to find an anomalous behavior in their amplitude or frequency evolution [22,[33][34][35][36]. Second, the frequency evolution of a candidate can be cross-checked against a list of instrumental lines in order to understand whether its significance is caused by the crossing of one or several lines.
Moderately-strong CW signals may appear as as narrow features in the power spectrum (see Fig. 1).These features will become more obvious as the sensitivity of the interferometric detectors improves [34].Even for the current generation of detectors, the possibility of a CW signal being visible in the detector has not been ruled out.Vetoing candidates near narrow spectral features, without additional evidence of their terrestrial origin, would result in an as-yet unquantified increase in the false dismissal probability of true CW signals.
Current practices in the LVK collaboration therefore discourage the use of unknown lines lists to veto CW candidates.The presence of a visible artifact in the data without a clear instrumental origin is not considered enough evidence for a candidate to be deemed non-astrophysical.However, No systematic study exists, to date, to back up such a recommendation.
In this paper, we study the similarity of CW signals and spectral lines in the current and next generation of interferometric detectors.Specifically, we construct a statistical criterion to quantify whether a CW signal is "visible" in the power spectrum of an interferometric detector.This criterion will be used to understand if artifacts in detector data can be consistent with astrophysical CW signals and strengthen the recommendation of not using unidentified lines to veto the results of a CW search for unknown sources.The paper is structured as follows: in Sec. 2 we introduce basic statistical notation and a criterion to quantify the "visibility" of a CW signal in a power spectrum.In Sec. 3, we compare the distribution of visible CW signals to an optimistic astrophysical distribution of sources.Conclusions are drawn in Sec. 4.

Characterizing the visibility of continuous-wave signals
In this section, we construct a quantitative criterion to flag a CW signal as "visible" in a power spectrum.A visible signal would be susceptible of being flagged as an unknown line.This criterion will be used in Sec. 3 to understand whether astrophysically-possible CW sources may produce visible signals and understand the severity of using unknown lines in a search veto procedure.
The output of a ground-based GW detector can be described as a time series of additive zero-mean Gaussian noise n and, possibly, a GW signal h CW analyses usually make use of short-time Fourier transforms (SFTs) [37], which are the Fourier transforms of short data segment with a duration T SFT , typically less than a few hours: Here, the superscript α refers to the time at the beginning of an SFT, M is the number of data samples within an SFT, and ∆t = T SFT /M.Frequency resolution is related to SFT duration as δ f = T −1 SFT .Noise within an SFT can be assumed to be white and stationary due to their short duration.As a result, the noise is fully characterized by its single-sided power spectral density (PSD) S n .
The optimal strategy to identify a monochromatic signal in Gaussian noise is to identify local maxima in the frequency-domain spectrum [38].This strategy remains useful for quasi-monochromatic signals as long as the frequency modulations are small compared to the frequency resolution of the dataset and the noise distribution remains stationary.For the case of long-duration narrow-band signals, non-stationarities in the noise distribution can be dealt with by normalizing the power in each SFT according to their PSD where we used the well-known relation between variance and PSD [14] The collection of time-frequency normalized-power values P α ( f ) is usually referred to as a "spectrogram".This spectrogram can then be averaged along time (i.e.SFTs) to reveal the presence of persistent signals: We shall refer to the frequency-dependent quantity P( f ) as "power spectrum".Line identification tasks consists in identifying narrow artifacts in P( f ) and cross-correlating them with an instrumental cause in the detector [28].
The statistical properties of P are well-known.Assuming Gaussian noise, P α is the sum of the squares of two zero-mean unitary-variance Gaussian random variables, thus it follows chi-squared distribution with two degrees of freedom P α ∼ χ 2  2 .P, on the other hand, is the average of N SFT identical and independent random variables, with usually N SFT ∼ O(10 3 − 10 4 ).Thus, due to the central limit theorem, P follows a Gaussian distribution with mean µ = 2 and standard deviation σ = 2, P ∼ Gauss(2, 2).As a result, the presence of artifacts in the data is related to deviations from Gaussianity in the distribution of P.
We will assume throughout this work that data consist of Gaussian noise and a single CW signal.Thus, the visibility of a CW signal is directly related to how much P deviates from Gaussianity.The characterization of a similar method to be applied as a decision criterion on a real-data situation is beyond the scope of this work and left for future work.
Tests for deviations from Gaussianity are commonly referred to as "normality tests" [39].These have customarily been used in a CW searches to identify parameter-space regions contaminated by noise artifacts [17].In order to choose a specific normality test, one needs to understand what kind of deviations from Gaussianity are expected in the data.For instance, Ref. [40] uses the Cramér-von Mises statistic [41], to flag overpopulations of outliers within a few standard deviations of the mean; this is motivated by the fact that CW signals, on the other hand, tend to populate the very far end of a distribution's tail and thus have a negligible contribution to the Cramér-von Mises statistic .
In Fig. 1, we show three power spectra containing CW signals at different sky locations.The behavior displayed by a CW signal is that of a relatively narrow disturbance that significantly shifts the power of the affected frequency bins well into the positive tail of the background noise distribution.As thoroughly discussed in [42], the presence of such outliers can be readily measured by the excess kurtosis where and angle brackets denote ensemble average.The excess kurtosis is constructed so that κ = 0 for a Gaussian variable, since ⟨(X − µ) 4 ⟩ = 3σ 2 for any Gaussian random variable X with mean µ and standard deviation σ.Kurtosis is expected to be positive if the distribution of P contains outliers significantly displaced from the background average, which makes it an appropriate tool to construct a visibility criterion.
In the following subsections, we characterize κ to quantify the visibility of a CW signal in Gaussian noise.These results will then be used in Sec. 3 to assess the visibility of astrophysical CW signals.

Expected kurtosis from a finite Gaussian sample
Several kurtosis estimators have been proposed in the literature [43,44].The distribution of these estimators is generally dependent on the distribution of the underlying samples, and must be characterized in order to cogently quantify deviations from Gaussianity.In this work, we compute the sample kurtosis using the scipy.stats.kurtosisfunction as implemented in SciPy [45].To characterize kurtosis distribution, shown in Fig. 2, we numerically generate 1.6 • 10 5 noise power spectra and compute the sample kurtosis.
Samples are generated using lalpulsar_Makefakedata_v5 [46] to simulate a 1-year Gaussian-noise datastream with a PSD of S n = 1 × 10   T SFT = 1800 s, 7200 s, 18 000 s.Each realization lasts for 1 year and contains 4500 frequency bins.Thresholds corresponding to the 99% quantile are highlighted in gray and collected in Table 1.
T SFT 1800 s 7200 s 18 000 s κ th 0.182 0.185 0.188 Table 1.Kurtosis thresholds corresponding to the 99% quantile of the numerically-generated kurtosis distribution shown in Fig. 2. We note that lower T SFT values tend to cause lower κ th .Since the duration of the observing run is fixed to 1 year, a lower T SFT value causes a higher number of SFTs to be generated.This causes the distribution of the resulting P to be closer to a Gaussian, which in turn diminishes the spread of the kurtosis distribution.
using three different time lengths, namely T SFT = 1800 s, 7200 s, 18 000 s.For each SFT, we computed their power spectrum [Eq.( 7)], which then was used to compute the sample's kurtosis κ.The frequency band ∆ f was adjusted according to T SFT so that all samples contained ∆ f /δ f = 4500 frequency bins; as a result, sample kurtosis is always estimated using the same number of samples and thus the kurtosis distributions for all T SFT values are comparable.We select a kurtosis threshold κ th corresponding to the 99% quantile of the resulting kurtosis distribution for each T SFT , as shown in Table 1.Any Gaussian-noise power spectrum yielding a kurtosis above κ th will be flagged as containing a visible CW signal.We exemplify this visibility criterion in Fig. 3, where we show "visible" CW signal according to our criterion κ > κ th .This criterion will be used in Sec.2.2 to compute the visible fraction of an astrophysicallymotivated population of CW sources.Concretely, we will compute the CW amplitude at which a certain fraction of the population of CW signals returns a kurtosis above κ th .The choice of κ th as the 99% quantile implies that the amplitudes corresponding to small visible fractions (e.g. less than 10% of visible signals) will be slightly overestimated, as 1% of the strictly non-visible signals will pass the threshold due to noise fluctuations.The results here reported, thus, will be conservative in the sense that the risks associated to a CW signal being visible in the power spectrum will not be underestimated.

Kurtosis from CW signals
To characterize the visibility of CW signals in a power spectrum, we simulate a population of CW sources with different amplitudes and determine the fraction of signals above the kurtosis threshold κ th computed in Sec.2.1.We then compute, using a numerical fit, the amplitude at which a representative fraction of signals.This process is akin to sensitivity estimation procedures used in CW searches (e.g.[47] and references therein).
We select 8 representative frequency bands across the frequency range expected for a CW signal, namely [10, 30, 100, 200, 300, 500, 700, 1000] Hz.This choice is motivated by the fact that the bandwidth of the Doppler modulations induced on the CW are proportional to the CW's frequency [Eq.( 2)]: higher frequencies are expected to modulate along broader frequency bands, and thus produce less prominent peaks than lower frequencies.
For simplicity, we set all the spindown parameters to zero.This is consistent with the usual spindown values observed in the known pulsar population [48].Non-zero spindown values would cause a broader modulation of the CW's frequency, which in turn would diminish the prominence of the resulting peak in the power spectrum, reducing the visibility of a signal.This choice is inline with producing a conservative result that does not underestimate the visibility of a CW signal in the power spectrum.
The CW amplitude h 0 is specified in terms of the sensitivity depth D [2,49,50] We select a range of sensitivity depth values from 0.1 to 80 in order to sample the full range of visibility fractions.PSD is fixed to S n = 1 × 10 −46 Hz −1 .Throughout this analysis, we simulate signals as seen in a 1-year observing run of the Advanced LIGO Hanford detector.Since 1) the duration of the observing run is on the order of a year, 2) CW sources are uniformly distributed across the sky, 3) CW signals span a very narrow frequency band, in the sense that the PSD can be considered constant, and 4) amplitudes are expressed relative to PSD by means of the sensitivity depth, the resulting sensitivity depths are valid for any other ground-based interferometric detector, and will be used in Sec.2.2 to estimate the visibility of a possible astrophysical population of CW sources for different detectors.
For each depth and frequency band, we simulate N I = 3000 CW signals in Gaussian noise and compute the fraction of visible signals p V (D) as where κ (n) is the kurtosis of the n-th simulated signal.This quantity describes the fraction of visible signals in a population with a given constant depth.More generally, p V is the probability of drawing a visible signal from the specified population at a given depth.The results are shown in Fig. 4. The specific p V value of interest depends on the application at hand.we fit a sigmoidlike ansatz using the function NonlinearModelFit as implemented in Mathematica [51] to and compute the visible sensitivity depth at p V = 10% and p V = 90%.The resulting sensitivity depths are listed in Table 2 and shown in Fig. 5, and follow the expected trend: higher frequencies produce a broader Doppler modulation, and thus  Table 2. Sensitivity depth D for p V = 10% and p V = 90% as estimated from the results shown in Fig. 4. The uncertainty in each value, which we estimated by propagating the binomial uncertainty on the empirical p V using the covariance matrix of the numerical fit, is ±1 for all the values in the table.Depth values appear to be systematically biased towards higher values as we reduced T SFT .This is because κ th is lower for lower T SFT values, as discussed in Table 1.  2. require a higher amplitude (lower depth) in order to produce a visible peak in the power spectrum.The higher p V , the bigger the fraction of visible signals: as a result, at a given frequency, D tends to decrease (higher amplitudes) as we increase p V .As discussed in Table 2, the results appear to be shifted for different T SFT values.This is a result of the different kurtosis thresholds, as discussed in Table 1.From Fig. 5, we can conclude that low-frequency CW signals tend to be "more visible" than high-frequency signals at a similar amplitude, in the sense that the lower Doppler modulation allows for a more prominent peak to be displayed.
These sensitivity depths should be interpreted as the relative amplitude with respect to the background noise required for a CW signal to be visible in a power spectrum.For a given detector PSD, these values can be converted to the corresponding CW amplitude h 0 using Eq.(11).To understand the risks associated to a visible CW signal, such as mistakenly flagging it as an unknown line or whether visible artifacts are enough to ascribe a nonastrophysical origin, however, we should compare these results to the expected amplitude from a population of astrophysical CW sources.

Implications for astrophysical CW sources
The visible sensitivity depths recorded in Table 2 must be compared to a realistic estimate of the expected amplitude of a CW signal in order to assess the risks of associated to a visible CW signal in the power spectrum.Since no CW has been detected to date, we will follow the conservative approach taken so far and construct an optimistic astrophysical CW source (Sec.3.1) to compare against the obtained results for different detector configurations (Sec.3.2).

An optimistic CW source
The nominal amplitude of a CW signal h 0 emitted by a NS with a quadrupolar deformation parametrized by the equatorial ellipticity ε is given by [13] where f 0 is the CW frequency (twice the rotational frequency in this model), d is the distance from the NS to the detector, and I = 10 38 kg m2 is the canonical moment of inertia of a NS around the spinning axis [52] 2 .The results from this section could also be re-interpreted assuming other other emission mechanisms, such as r-modes or free precession [2,53].
The latests theoretical bounds on the maximum ellipticity sustained by a NS are ε ∼ 10 −6 [54,55].This upper bound increases by a few orders of magnitude if more exotic objects are considered [56].On the lower end, the observed population of millisecond pulsars appears to sustain an ellipticity of about ε ∼ 10 −9 [57].
We construct an optimistic CW source by considering the full range of plausible ellipticities ε ∈ [10 −9 , 10 −6 ] for a NS located at d = 20 pc, which corresponds to the closest location to the detector at which a NS is expected to be found [58].Note that h 0 depends on the ratio ε/d; for example, the expected amplitude of a source with ε = 10 −6 at 2 kpc is equal to that of a source with ε = 10 −8 at 20 pc.In Fig. 6. we show the expected sensitivity depth for optimistic sources in the sensitive band of the Advanced LIGO (O3, O4, and O5 sensitivities) [59][60][61][62], Einstein Telescope [63,64], and Cosmic Explorer [65] detectors.This is computed by using Eq. ( 14) and the available PSD curves in the literature.
The resulting sensitivity depth curves are such that D ∝ √ S n f −2 0 .For a given ellipticity and distance to the source, low frequencies tend to produce high depth values (weak signals) as h 0 ∝ f 2 0 and the PSD of the detector tends ot be very steep.Towards higher frequencies, where shot-noise is dominant, h 0 grows faster than the detector's PSD, resulting in a monotonically decreasing sensitivity depth (increasing signal amplitude).For a given ellipticity and distance to the source, a CW signal is "easier" to detect at a high frequency than at a low frequency.
We note that all-sky CW search upper limits in O3 data [47] rule out this optimistic population.Said upper limits, however, are only valid for deterministic CW signals whose frequency evolution is not affected by stochastic effects, such as glitches [66,67] or spinwandering [68].For a generic CW signal, the results as reported by the SOAP pipeline [69,70] are less constraining.

Visibility of CW signals in a power spectra
In Figs. 7, 8, and 9 we compare the visible sensitivity depths derived in Sec. 2 to the expected sensitivity depths from astrophysical CW signals introduced in Sec.3.1.To ease the discussion, we use the results using T SFT = 1800 s, which are the most conservative.Other T SFT values produce comparable results as the variations are small.Figure 6.Sensitivity depth for an optimistic CW source located at d = 20 pc for different groundbased interferometric detectors.The lower limit of the shaded area corresponds to ε = 10 −6 ; the upper limit corresponds to ε = 10 −9 .The presence of narrow artifacts in the O3 sensitivity depths is due to the use a PSD estimated using real-data [59,60].The shoulder-like plateau at low frequencies for the ET detector corresponds to the transition of the low-frequency interferometer to the high-frequency interferometer [64].
We are interested in identifying frequency bands for which the visible sensitivity depth is higher than the expected depth for an optimistic CW signal.In such bands, CW signals will be visible in the power spectrum.
Overall, we find that visible signals tend to be located in the upper end of the frequency spectrum.This fact seems to contradict the results from Sec. 2.2, where we concluded that high frequencies are less visible than low frequencies due to a broader Doppler modulation.The detector's sensitivity curve, however, degrades very steeply towards low frequencies; towards high frequencies, where shot-noise is dominant, the frequency dependence is much more gentle compared to h 0 's quadratic dependency on f 0 , and thus the resulting sensitivity depth quickly becomes lower enough to be visible.
For high ellipticity values, ε ∈ 10 −7 , 10 −6 , and for the Advanced LIGO detectors in their O3 configuration (Fig. 7), CW signals would start to become visible at about 100 Hz.As we progress into O4 and O5 sensitivities (Fig. 8), this frequency may reduce close to 50 Hz.Typical searches for CW signals from unknown sources survey frequencies within [20,1000] Hz [47].As a result, as we approach the design sensitivity of the Advanced detectors, CW signals are expected to be visible for 95% of the frequency band.For 3G detectors (Fig. 9), high ellipticity signals will be visible across practically the whole frequency band.

Conclusion
We have studied whether an astrophysical CW signal could be visible as a narrow spectral artifact in the power spectrum of a gravitational-wave detector.To do so, we have calibrated the sample kurtosis of a power spectrum as a measure of the visibility of a CW signal.We found that, depending on the signal's frequency, sensitivity depths between D = 5 and D = 60 could become visible in a power spectrum, and thus mimic the behavior of a line.
These results were compared to the expected amplitude of an optimistic astrophysical signal.For high ellipticities ε ∈ 10 −7 , 10 −6 , CW signals become visible from 100 Hz onwards for the Advanced LIGO detectors which amount to 95% of the frequency bands surveyed in a typical all-sky search; for 3G detectors, signals will be visible across the whole frequency band.For low ellipticities ε ∈ 10 −9 , 10 −8 , CW signals would start to become visible at 700 Hz for the O4 detectors, 500 Hz for the O5 detectors, and 100 Hz for 3G detectors.We conclude that it is unsafe to discard CW search outliers near narrow spectral artifacts of unknown origin.Safely discarding outliers near spectral artifacts requires evidence that the artifact is instrumental.
Throughout this study we assumed a population of CW signals with a negligible spindown parameter f 1 ≈ 0 Hz/s.Given the typical resolution of an all-sky CW search, this assumption is consistent with the observed pulsar population [48].The impact of higher spindown values would be to spread the power of a CW signal along neighboring frequency bins, reducing the peak's prominence in a power spectrogram.This is more likely to affect high ellipticity sources, which are expected to have a higher spindown value.The results here presented tend to overestimate the visibility of high ellipticity sources.
Also, in this study we neglected two properties of CW signals, namely their distinctive double-horned shape in a power spectrum [34] and the fact that they would be expected to show up in multiple detectors at once.The former requires the CW signal to be located in a relatively quiet frequency band, as otherwise the expected horns would end up tarnished by instrumental or unknown lines.The latter requires a network of detectors with comparable sensitivities at a given frequency band, which makes it depend heavily on the specific configuration of the detectors.
Our results draw attention to the importance of conducting detector characterization studies on narrow spectral artifacts, and to the non-negligible possibility of missing an       obvious CW signal if such tasks are neglected.Although a significant (and increasing) number of lines have been characterized and mitigated thanks to detector characterization efforts within the LIGO-Virgo-KAGRA collaboration, the origin of a comparable amount still remains, to date, unknown.The importance of rigorous investigation to distinguish between instrumental artifacts and genuine CW signals will only increase as we progress into the era of design sensitivity of the advanced detectors and further beyond into 3G.

Figure 1 .
Figure 1.(Top) Power spectra for the three CW signals at f 0 = 50 Hz and null spindown parameter.(Bottom) Distribution of power spectra and kurtosis for the three same CW signals.As discussed in Sec. 1, the width of the spectral artifact is related to the sky position of the CW source.Sources away from the sky poles (low declination |δ| ≲ 1) tend to produced broader signals than sources closer to the sky poles (high declination |δ| ≈ 1.5),The amplitudes of these signals corresponds to a sensitivity depth of D = 15.for illustration purposes.We refer the reader to Sec. 2.2 and Sec. 3 for a detailed discussion on relevant CW amplitudes.

Figure 2 .
Figure2.Distribution of sample kurtosis using 1.6 • 10 5 Gaussian noise realizations for T SFT = 1800 s, 7200 s, 18 000 s.Each realization lasts for 1 year and contains 4500 frequency bins.Thresholds corresponding to the 99% quantile are highlighted in gray and collected in Table1.

Figure 4 .
Figure 4. Fraction of visible signals p V for a population of all-sky isotropically-oriented signals at different frequencies and sensitivity depths.Each simulated signal is injected in a 1-year data set for the Advanced LIGO H1 detector.Each computed p V value has an associated binomial uncertainty of less than 1%.

Figure 7 .
Figure 7.Expected sensitivity depth for an optimistic source consistent with a NS located at 20 pc as observed by the Advanced LIGO Hanford[59] (top) and Livingston[60] (bottom) detectors during O3.Solid lines correspond to different ellipticity values (ε = 10 −9 is beyond D = 100).Note that these are a discretized version of the shaded areas in Fig.6.Triangles correspond to visible sensitivity depths at p V = 10% and p V = 90% as reported in Table2.
Figure 7.Expected sensitivity depth for an optimistic source consistent with a NS located at 20 pc as observed by the Advanced LIGO Hanford[59] (top) and Livingston[60] (bottom) detectors during O3.Solid lines correspond to different ellipticity values (ε = 10 −9 is beyond D = 100).Note that these are a discretized version of the shaded areas in Fig.6.Triangles correspond to visible sensitivity depths at p V = 10% and p V = 90% as reported in Table2.