Search for Continuous Gravitational Waves from Scorpius X-1 in LIGO O2 Data

We present the results of a search in LIGO O2 public data for continuous gravitational waves from the neutron star in the low-mass X-ray binary Scorpius X-1. We search for signals with $\approx$ constant frequency in the range 40-180 Hz. Thanks to the efficiency of our search pipeline we can use a long coherence time and achieve unprecedented sensitivity, significantly improving on existing results. This is the first search that has been able to probe gravitational wave amplitudes that could balance the accretion torque at the neutron star radius. Our search excludes emission at this level between 67.5 Hz and 131.5 Hz, for an inclination angle $44^\circ \pm 6^\circ$ derived from radio observations (Fomalont et al. 2001), and assuming that the spin axis is perpendicular to the orbital plane. If the torque arm is $\approx $ 26 km -- a conservative estimate of the \alfven\ radius -- our results are more constraining than the indirect limit across the band. This allows us to exclude certain mass-radius combinations and to place upper limits on the strength of the star's magnetic field. We also correct a mistake that appears in the literature in the equation that gives the gravitational wave amplitude at the torque balance (Abbott et al. 2017b, 2019a) and we re-interpret the associated latest LIGO/Virgo results in light of this.


INTRODUCTION
Fast spinning neutron stars are promising sources of continuous gravitational waves in the frequency range 20 Hz -2 kHz. The emission is typically generated by a non-axisymmetry in the star with respect to its rotation axis. The simplest example is the presence of an equatorial ellipticity that deforms the star into a triaxial ellipsoid rotating around the principal moment of inertia axis (Jaranowski et al. 1998).
The strength of the gravitational wave signal is proportional to the ellipticity of the star. The maximum ellipticity that a neutron star could support before breaking has been estimated to lie in the 10 −7 − 10 −5 range for neutron stars made of normal matter and a few orders of magnitude higher for exotic matter (Horowitz & Kadau 2009;Johnson-McDaniel & Owen 2013;Baiko & Chugunov 2018;Gittins et al. 2020). The minimum ellipticity is harder to estimate: we expect some ellipticity due to magnetic deformation, but the precise value depends strongly on the assumed magnetic field strength and configuration (see for example Haskell et al. 2008;Mastrano et al. 2011;Suvorov et al. 2016). Woan et al. (2018) have argued for a minimum ellipticity ∼ 10 −9 based on the spin-down of millisecond pulsars (due to either magnetic field effects or some other source of ellipticity such as crustal deformation).
For accreting neutron stars, the accretion process provides a potential additional source of asymmetry, par-ticularly if accreting material is channeled unevenly onto the surface by the star's magnetic field. This can lead to thermal and compositional gradients in the crust that generate a crustal 'mountain' (Bildsten 1998;Ushomirsky et al. 2000;Haskell et al. 2006;Singh et al. 2020). Accretion-induced deformation of the star's magnetic field might also result in asymmetries (Melatos & Payne 2005;Vigelius & Melatos 2009). Accretion could also drive the excitation of some kind of internal oscillation that results in gravitational wave emission (Andersson et al. 1999;Haskell 2015). Uncertainty about the accretion process and the stellar response makes it hard to compute firm estimates for the expected size of the resulting ellipticities, but they could be large enough for the resulting gravitational wave emission to be detectable with the current generation of detectors (Lasky 2015).
What effect might such a gravitational wave torque have on an accreting neutron star? It has long been noted (Papaloizou & Pringle 1978;Wagoner 1984) that neutron stars in low mass X-ray binaries, in spite of having accreted matter for millions of years, spin well below the maximum possible spin frequency (Cook et al. 1994;Haensel et al. 2009), with the fastest accreting neutron star spinning at 620 Hz (Hartman et al. 2003;Patruno & Watts 2012;Watts 2012;Patruno et al. 2017). Since gravitational wave torques scale with a high power of the frequency, as the spin rate increases, they naturally provide a mechanism that kicks-in more strongly than other mechanisms, preventing further spin-up. This has led to the idea of torque balance, where gravitational wave and accretion torques reach equilibrium, preventing further spin-up and ensuring continuous gravitational wave emission (Bildsten 1998). Indeed Gittins & Andersson (2019) have shown that a synthetic population of neutron stars evolved without the gravitational wave torque contribution does not produce the observed spin distribution.
The accretion torque on a neutron star having mass M is where G is the gravitational constant, r m is the torque arm andṀ the accretion rate. The correct value to use for r m is not known a priori, but is typically assumed to be either the neutron star radius R or the radius at which the star's magnetic field starts to disrupt the accretion flow. The maximum accretion luminosity is GMṀ R . If some fraction X of this is radiated away by an X-ray flux F X observed at a distance d, then The gravitational wave intrinsic amplitude h 0 at a distance d, for a gravitational wave signal at twice the spin frequency of the star (which is the case if the ellipticity is caused by a magnetic or crustal mountain) and balancing the accretion torque, is where f GW is the gravitational wave frequency. Substituting Eq. 1 and 2 in Eq. 3 one finds h torq.bal.
We note that Eq. 15 in Abbott et al. (2019a) and Eq. 10 in Abbott et al. (2017b) are incorrect and yield the correct numerical value only if r m = R. In those papers such mistake propagates to the Alfvén radius torque balance amplitude curve of Fig. 5 (yellow curve in Abbott et al. 2019a), which is over-estimated. This in turn makes it look like the constrained inclination angle upper limits from that search (for ι = ι orb ≈ 44 • ) probe the Alfvén radius torque balance limit, when in fact they do not.
For ease of notation we define and re-write Eq. 4 as Scorpius X-1 (Sco X-1) is the brightest persistent Xray source after the Sun and hence, given the scaling of gravitational wave amplitude with X-ray flux, it is a particularly promising continuous wave source. The flux value of 3.9 × 10 −7 erg/cm 2 /s used in Eq. 4 is the long-term average X-ray luminosity of Sco X-1 measured from Earth (see Watts et al. 2008, for details of how this value was derived 1 ). This value yields a torque balance h 0 well within the reach of searches for continuous waves from known pulsars (Abbott et al. 2019b(Abbott et al. , 2020Nieder et al. 2019Nieder et al. , 2020.
Many searches have targeted continuous gravitational wave emission from Sco X-1 (Meadors et al. 2017;Abbott et al. 2017aAbbott et al. ,c,d, 2019a, only since 2017), but none have yet been sensitive enough to probe the torque balance amplitudes of Eq. 4. This is because in contrast to the known pulsars targeted in Abbott et al. (2019bAbbott et al. ( , 2020; Nieder et al. (2019Nieder et al. ( , 2020, the rotation frequency and frequency derivative of the Sco X-1-neutron star, as well as some binary parameters, are unknown. This means that a broad range of waveforms must be tested against the data, and this degrades the attainable sensitivity, through the increased trials factor. Another aspect that makes the Sco X-1 signal search challenging is its computational cost: as illustrated in Watts et al. (2008) our ignorance of the system parameters results in a parameter space so broad that the most sensitive search method, a coherent matched filter over the entire observation time, is computationally prohibitive. This is a frequent predicament in searches for continuous gravitational waves and the standard solution is to adopt semi-coherent search methods, where one trades sensitivity in favour of computational efficiency (Messenger et al. 2015;Dergachev & Papa 2019).
In semi-coherent searches the observation time is partitioned in segments spanning the same duration. If data from several instruments is used, these partitions are ≈ coincident in time. The most important quantity is the duration of such partitions, T coh . The larger T coh is, the more sensitive and the more computationally expensive the search is going to be.
We use for this search a cross-correlation method (Whelan et al. 2015, and references therein). Thanks to the much improved computational efficiency of our new search (Meadors et al. 2018), we are able to use a significantly longer T coh than ever used before and reach unprecedented levels of sensitivity. In particular for the first time a search is sensitive to signals at the torque balance limit at both the stellar radius and for reasonable estimates of the magnetospheric radius.
We use LIGO O2 open data from the Hanford and Livingston detectors (LIGO 2019; Vallisneri et al. 2015) between GPS time 1167984930 (January 2016) and GPS time 1187733514 (August 2016). Overall we have 5090 hours of data, 2496 from Livingston and 2594 from Hanford.
We search for a nearly monochromatic signal from the neutron star in Sco X-1 -below we qualify this assumption further. At the detector the signal appears frequency-modulated due to the relative motion between the star and the detector, and amplitude-modulated due to the sensitivity-response of the detectors, which depends on the line-of-sight direction and hence for a fixed source changes with time. If all the source parameters were known, the gravitational waveform at the detector would also be known, and the search would be a perfectly matched filter, like those carried out for known pulsars. This is not the case.
The low-mass X-ray binary (LMXB) Sco X-1 consists of a 1.4 +1.4 −0.5 M neutron star and a 0.7 +0.8 −0.3 M companion star (Wang et al. 2018, 95% confidence intervals). No accretion-powered pulsations or thermonuclear burst oscillations have been so far detected from the neutron star, so its spin frequency is unknown. The orbital parameters projected semi-major axis, a sin i, time of ascending nodes, T asc , and orbital period, P orb , are constrained within ranges larger than our search resolution on those parameters, so these need to be explicitly searched (Wang et al. 2018).
The search parameters are given in Table 1. We search for gravitational wave signal frequencies between 40 Hz and 180 Hz. The computational cost per unit frequency interval is smaller at lower frequencies, so concentrating computational resources in the lower frequency range makes for the highest return in sensitivity. In fact this is the frequency range in which we can match the torque balance limit, even with an unrestricted prior on the star's inclination angle.
We do not explicitly search over frequency derivatives, reflecting the assumption that the system is close to equilibrium. With our search set-up we have measured an average loss in SNR at the 15% level for gravitational wave first frequency derivative |ḟ GW | 2 × 10 −13 Hz/s. This sets the scale for the maximum rate of change of the spin frequency that would not affect our ability to detect a signal, at |ḟ spin | 1 × 10 −13 Hz/s. We recall that for crustal mountains f GW = 2f spin .
The orbital parameter ranges are taken from Table 2 of Wang et al. (2018). T asc is propagated to 1178556229 GPSs which is ≈ the weighted middle of the LIGO data observation span. We note that this is 206 epochs after the T asc in (Abbott et al. 2019a). Following Eq. 5 of a Time of ascension has been propagated to May 11 16:43:31 UTC 2017, close to the weighted-middle of the gravitational wave data, in order to make the metric approximately diagonal (Whelan et al. (2015)). The relation between Tasc and the epoch of inferior conjunction of the companion star T0 presented in Wang et al. (2018) is Tasc = T0 − P orb /4 (Messenger et al. (2015)) . Galloway et al. (2014), we expand the uncertainty associated with T asc to 139 seconds and then consider the 3σ confidence interval.
The grid spacings in every dimension dλ are chosen so that the loss due to signal-template mismatch is at the m = 25% level. The spacings are estimated based on the metric g λλ as dλ = m g λλ . Expressions for the metric can be found in Whelan et al. (2015). This approach results in an overestimate of the actual mismatch (Allen 2019), and in fact we measure an overall average SNR loss of 16%. The grid spacings are given in Table 1.
Our search employs a fixed T coh 19 hrs, which is a factor of 4.5 (10) longer than the longest (shortest) baseline used by Abbott et al. (2017e). This choice, enabled by the efficiency of our code (see Section 1) is the reason for the higher sensitivity of our search.
We consider search results with detection statistic values above the expected Gaussian-noise fluctuations. Since the number of searched waveforms increases with frequency, noise fluctuations alone produce higher fluctuations at higher frequencies. For this reason our threshold for candidate consideration is not constant but rather increases with frequency.
We find over 97 million results above the threshold. As often happens in this type of search, these results are not uniformly distributed in frequency but tend to come in groups, with the elements of each group having similar signal frequency, and due to the same root cause. We cluster these together and examine each cluster. We find 32 such groups, which we will refer to as "outlier clusters". 23 of them are associated with known spectral contaminations (Covas et al. 2018). The rest of the clusters are discarded based on cross-checking the multi-detector detection statistics with the singledetector statistics: When an outlier is due to a disturbance in one of the detectors, the single-detector statistics will often be larger than the multi-detector one. On the contrary a signal produces a higher value when the data from both detectors is used. Most of these discarded clusters also present a range of signal frequencies with enhanced values of the detection statistic, that is much larger than it would be for a signal. The complete list of outlier clusters is given in Table 2. 3. RESULTS

Upper limits on GW amplitude
As no significant candidate is found, we set upper limits at the 95% confidence level, on the gravitational wave intrinsic amplitude h 95% 0 at the detectors, in half-Hz bands. The upper limits are determined by adding fake signals with a fixed amplitude h 0 to the data, and by measuring the detection efficiency, C(h 0 ). The detection criterion is determined by the value of the detection statistic of the most significant result in the band. The procedure is repeated for various values of h 0 and a sigmoid fit is used to determine the value corresponding to 95% confidence: C(h 95% 0 ) = .95 (Fesik & Papa 2020). Two sets of upper limits are derived, reflecting two assumptions: 1) an arbitrary value of the inclination angle, with cos ι uniformly distributed −1 ≤ cos ι ≤ 1, and 2) the inclination angle ι is equal to the orbital inclination angle and hence drawn from a Gaussian distribution with mean 44 • and standard deviation 6 • (Fomalont et al. 2001;Wang et al. 2018). The latter scenario is equivalent to assuming that the spin axis of the neutron star is perpendicular to the orbital plane. The ι = 44 • ±6 • is a more favourable inclination than average for coupling to the gravitational wave detector (Jaranowski et al. 1998, see for instance Eq.s 21 and 22 ) and the resulting upper limits are a factor ≈ 1.7 smaller than those for arbitrary orientation.
Both upper limits are plotted in Fig. 1 and provided in machine-readable format in Zhang et al. (2020, and Suppl. Mat.). For comparison Fig. 1 also shows upper limits from a previous cross-correlation search on O1 LIGO data (Abbott et al. 2017d) and from a recent Viterbi algorithm search on the same O2 data that we use (Abbott et al. 2019a). The most sensitive of the two searches is the Abbott et al. (2017d) search. It employed variable coherence lengths, with longer T coh in the low frequency range, which explains why at lower frequency it is comparatively more sensitive than at higher frequency. The Abbott et al. (2019a) search is less sensitive than a cross-correlation search but is more robust to de- Figure 1. 95% confidence upper limits on the intrinsic gravitational wave amplitude in half-Hz bands. We assume orbital inclination at 44 • ± 6 • (lower black points) and arbitrary inclination (upper black points). The lower dashed and solid curves are the torque-balance upper limits, based on estimates of the mass accretion rate and assuming the accretion torque to be at the neutron star radius (lower solid red curve) or at the Alfvén radius (upper dashed red curve). For comparison we show the upper limits from previous results (the three fainter upper curves) and draw the incorrect torque balance Alfvén radius upper limit that was reported (dash-dot line). (1, 1) ξ 1/4 X 3/7 , and this multiplicative factor is shown here.
viations of the signal waveform from the assumed model (Suvorova et al. 2016a,b). In particular the method of Abbott et al. (2019a) is robust with respect to loss of phase coherence in the signal.
One of the ways in which the signal could lose phase coherence with respect to the template waveforms of the search is through spin-wandering. This is a nondeterministic "jitter" in the spin of star, caused, for instance, by small changes in the mass accretion rate. The resulting frequency variation depends on the accretion torque, hence on the spin frequency of the star, its moment of inertia, the ratio between the torque arm and the co-rotation radius and the mass accretion rate.
Based on RXTE/ASM observations of Sco X-1, Mukherjee et al. (2018) have explored different systemparameter combinations and the gravitational wave fre-quency changes that may accumulate over different observation periods, due to spin wandering. Their results indicate that, in our frequency range, the maximum frequency change during an observation time of 2 × 10 7 s (our observation time) is less than 2 µHz (our frequency resolution) for the vast majority of the simulated systems. This means that the sensitivity of this search should not be impacted by spin-wandering effects.
Our results improve on existing ones by more than a factor ≈ 1.8. This is an extremely large sensitivity improvement in a large parameter space search like this one. For instance consider that in a broad all-sky search on O2 data, Abbott et al. (2019c) improve over the most sensitive results on O1 data (Dergachev & Papa 2019) by a factor of ≈ 1.1.

Interpretation in terms of torque balance model
Our results are also remarkable in absolute terms because they probe gravitational wave amplitudes that could support emission at the torque balance level. It is the first time that this milestone is reached.
From Eq.s 4 or 6 we see that the torque balance gravitational wave amplitude depends on the torque arm and it is smallest at the star surface. If this minimum torque balance amplitude is larger than our h 95% 0 upper limits it means that our search should have detected a signal; The fact that it has not, means that we can exclude such mass-radius combination: The lower panel of Figure 4 shows the mass-radius regions excluded by the ι ≈ 44 • gravitational wave upper limits for f GW = 117.5 Hz and X = 1. If the torque arm is larger than the star radius, the torque balance amplitude increases, and our gravitational wave upper limits constrain the magnetic field strength of the mass-radius combinations not excluded by 7. We illustrate this point in the next paragraphs.
We take the torque arm to be at the magnetospheric radius r m = max(ξr A , R), with 0.5 ≤ ξ ≤ 1 and r A the Alfvén radius where B 1 is the normalised polar magnetic field strength, defined in Eq. 5. We note that in the gravitational wave literature the Alfvén radius has often been placed at 35 km, corresponding toṀ = 10 −8 M /yr, or X = 0.3 in Eq. 2. The Eddington limit is aṫ M = 2 × 10 −8 M /yr, for a fiducial 1.4 M and 10 km radius neutron star.
By combining Eq. 4 and Eq. 8 we find the torquebalance amplitude when r m ≥ R: the last equation simply reflecting the condition r m ≥ R. We note that h torq.bal. 0 (ξ, X) = h torq.bal.

0
(1, 1) ξ 1/4 X 3/7 and this factor is plotted in Fig. 2 to aid evaluate how the torque balance amplitude changes under different assumptions for torque arm r m and the accretion luminosity.
When this torque balance amplitude is larger than our h 95% 0 upper limits it means that our search should have detected a signal; The fact that it has not, means that we can exclude the associated mass-radius-magnetic field strength combinations: This translates, for every mass-radius, into an upper limit on the magnetic field strength. Figure 4 shows the magnetic field upper limits from Eq. 10 from the ι ≈ 44 • gravitational wave upper limits for f GW = 117.5 Hz, X = 1 and ξ = 1, for different equations of state. The upper limits for different gravitational wave frequencies can be easily derived from the gravitational wave upper limit values using Eq. 10. For the specific example shown in Figure 4, provided that the field is higher than ∼ 2 × 10 8 G, the torque balance limit can be matched for all of the considered equations of state, but magnetic fields above ∼ 6 × 10 8 G can be ruled out.
At Zhang et al. (2020, and Suppl. Mat.) we provide plots like the one of Figure 4 for gravitational wave frequencies in the searched range, at 2 Hz intervals.
The gravitational wave upper limits marginalised over all possible inclination angles lead to less stringent constraints on the physical parameters of the neutron star: the torque balance amplitude with torque arm at the neutron star surface is smaller than our upper limits for all equations of state, so no mass-radius combination can be ruled out. Torque-balance amplitudes larger than our upper limits can only be obtained for larger torque arms corresponding to magnetic field strengths 10 9 G, which are higher than those expected from LMXBs.
If the gravitational wave signal is due to a triaxial ellipsoid rotating around a principal moment of inertia axis I, say along theẑ axis, the gravitational wave intrinsic amplitude h 0 is proportional to the ellipticity ε We convert the h 95% 0 upper limits into ellipticity upper limits with Eq. 11, with d = 2.8kpc and a fiducial value of I = 10 38 kg m 2 . We also derive the ellipticity required for torque balance under the two previous assumptions on the lever arm. All these quantities are plotted in Fig. 3, as a function of the gravitational wave signal frequency.
Above ∼ 90 Hz the values of the ellipticity that we are exploring are a few ×10 −5 and smaller. Deformations which are this large may be sustained by a neutron star crust (Johnson-McDaniel & Owen 2013), although very recent work suggests that the maximum deformations may be smaller (Gittins et al. 2020).

DISCUSSION
This search has placed upper limits on stable GW emission that are tighter than the level predicted by torque balance models for Sco X-1, for ι ∼ 44 • . This conclusion is robust to spin wandering at the level expected for this source. If the accretion torque is applied at the neutron star surface, the GW frequency range for which the torque balance limit is beaten is between 67.5-131.5 Hz, for a 1.4 M and 10 km radius fiducial star. If on the other hand the torque is applied at a magnetospheric radius at 25.6 km (see Eq.s 5 and 8), then the range for which the limit is beaten is the entire searched range, 40-180 Hz, for the fiducial star, as shown in Fig. 1.
If we consider a wider range of masses and radii, consistent with our current best understanding of viable equation of state models, we are able to place constraints on mass-radius-magnetic field strength combinations: • independently of the magnetic field value, our results exclude certain mass-radius combinations. Our tightest limits come for spin periods of ∼ 21 ms (f GW ∼ 96 Hz, at twice the spin frequency) with a narrow range of allowed masses extending only between 1.9-2.2 M and magnetic fields larger than ∼ 3×10 8 G being ruled out for all considered equations of state.
• if the magnetic field is larger than 0.19d 1 ξ −7/4 (F 1 /X) 1/2 (R 3 1 M 1 ) −1/4 (r m > R) we can place upper limits on the magnetic field strength. The upper limit on the magnetic field is highest for the highest mass in the range. Stiffer equations of state have a smaller range of masses (and magnetic field strengths) for which balance can still be possible at the level of our upper limits, than softer equations of state. We find that the field must be smaller than (4 − 10) × 10 8 G, depending on frequency (but excluding too disturbed frequency ranges, e.g. 60 Hz), for all equations of state models considered.
It is the first time that constraints on the magnetic field, mass and radius are obtained through continuous wave observations. This is interesting because the magnetic field is in general very poorly constrained and because observations like these probe mass-radius and magnetic fields through an entirely different mechanism than gravitational wave binary inspiral signals (see e.g. Abbott et al. 2019d;Capano et al. 2020).
If the spin of Sco X-1 is such that it is in the range where the limit is beaten (half the GW frequency for mountain models) and torque balance applies, this means that GW emission is not strong enough to balance the assumed accretion torque. This implies that the accretion torque must be less strong than predicted by the models presented in this paper, which could happen if, for example, strong radiation pressure modifies the structure of the inner disk (Andersson et al. 2005) or due to the effect of winds (Parfrey et al. 2016).
The result also puts limits on the size of thermal/compositional crustal or magnetic mountains in Sco X-1. Limits can also be placed on internal oscillation amplitude for models where that is the mechanism that provides the GW torque (for a different range of spin frequencies since the relationship between spin and GW frequency is different for mode models). . We have considered masses in the range 1 to 3 M , radii between 8 and 13 km and we have dropped any equation of state with a maximum mass lower than 1.9 M , consistent with observations (Antoniadis et al. 2013;Cromartie et al. 2019) and with estimates from dense matter theory and experiment, (see for example Hebeler et al. 2013;Kurkela et al. 2014). The lowest value of B for each curve corresponds to B(rm = R). Lower panel: mass-radius relations for the equations of state considered above (solid lines) and for star configurations of constant compactness (dashed lines). The line that delimits the shaded region shows the mass-radius combinations that satisfy Eq. 7, i.e. that are consistent with our upper limits when rm = R. Below the shaded region the torque balance gravitational wave amplitude with rm = R, is larger than our upper limits, so these configurations are excluded by our null results. Above the shaded region the torque balance gravitational wave amplitude with rm = R, is smaller than our upper limits, so these configurations cannot be excluded if rm = R. If, however, rm > R, i.e. a magnetic field above ≈ 2 × 10 8 G, the corresponding torque balance becomes larger than our upper limits and this allows us to constrain the magnetic field (as shown in the top panel).
An alternative, of course, is that the spin of Sco X-1 is outside the range searched, perhaps at higher frequencies more in line with the spin rates measured for the accretion-powered millisecond pulsars and thermonuclear burst oscillation sources (Patruno & Watts 2012;Watts 2012).
At 1 kHz the torque balance upper limit for the "fiducial star" is between 2.6-3.3 ×10 −26 , for r m ∈ [10 − 25.6] km. This is about a factor of 5-10 lower than a signal that we could detect at that frequency with a search like this one -the variation depending on the torque arm and on the inclination angle, and assuming that one could actually perform this search at such high frequencies. If all the parameters of Sco X-1 were known, a search on the same O2 data as used here could probably detect signals at 2.6-3.3 ×10 −26 . We are however quite far from having a complete ephemeris for Sco X-1. The next best thing would be to know the rotation frequency of the neutron star. The reason is that the torque balance amplitude decreases with frequency (so the sensitivity requirement increases, to match the torque balance limit), and the sensitivity of the searches decreases with frequency due to the shot noise in the detectors and to the increased template resolution per-Hz searched. These factors make it difficult to search very broad frequency bands. If it were possible to iden-tify the spin frequency, for example via the detection of weak or intermittent pulsations (a major goal for future large-area X-ray telescopes Watts et al. 2019;Ray et al. 2019), we might be able to carry out a search like this one, that could begin to probe the torque balance limit when the noise level at 1 kHz reaches its design value of ∼ 5.5 × 10 −24 1/ √ Hz (Abbott et al. 2018) and with ∼ two years of data.

ACKNOWLEDGMENTS
The computation of the work was run on the ATLAS computing cluster at AEI Hannover AEI (2017)