Continuous Gravitational Wave Signals Research Articles

Deciphering accretion-driven starquakes in recycled millisecond pulsars using gravitational waves

ABSTRACT Recycled millisecond pulsars are susceptible to starquakes as they are continuously accreting matter from their binary companion. A starquake happens when the rotational frequency of the star crosses its breaking frequency. In this study, we perform a model analysis of an accreting neutron star suffering a starquake. We analyse two models: a spherical star with accreting mountains and a deformed star with accreting mountains. We find that as the star crosses the breaking frequency and suffers a starquake, there is a sudden change in the continuous gravitational wave signal arriving from it. The amplitude of the gravitational wave signal increases suddenly both for the spherical and deformed star. For the spherical star, the accreting matter entirely dictates the amplitude of the gravitational wave. For the deformed star, both the accreting matter and the deformation from spherical symmetry play a significant role in determining the amplitude of the gravitational wave signal. This sudden change in the continuous gravitational wave signal in recycled millisecond pulsars can be a unique signature for such pulsars undergoing a starquake.

The second data release from the European Pulsar Timing Array V. Search for continuous gravitational wave signals

We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favour of the presence of the CGW over a common uncorrelated noise process. In contrast, the frequentist analysis estimates the p-value of the candidate to be $< 1<!PCT!>$, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only $0.7$. Therefore, we cannot conclusively determine the origin of the observed feature, nor can we rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations for the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature.

Assessing the Similarity of Continuous Gravitational-Wave Signals to Narrow Instrumental Artifacts

Continuous gravitational-wave (CW) signals 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 characterization 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 signal 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.

The LISA Data Challenge Radler analysis and time-dependent ultra-compact binary catalogues

Context. Galactic binaries account for the loudest combined continuous gravitational wave signal in the Laser Interferometer Space Antenna (LISA) band, which spans a frequency range of 0.1 mHz–1 Hz. Aims. A superposition of low frequency Galactic and extragalactic signals and instrument noise comprise the LISA data stream. Resolving as many Galactic binary signals as possible and characterising the unresolved Galactic foreground noise after their subtraction from the data are a necessary step towards a global fit solution to the LISA data. Methods. We analysed a simulated gravitational wave time series of tens of millions of ultra-compact Galactic binaries hundreds of thousands of years from merger. This dataset is called the Radler galaxy and is part of the LISA Data Challenges. We used a Markov chain Monte Carlo search pipeline specifically designed to perform a global fit to the Galactic binaries and detector noise. Our analysis was performed for increasingly larger observation times of 1.5, 3, 6 and 12 months. Results. We show that after one year of observing, as many as ten thousand ultra-compact binary signals are individually resolvable. Ultra-compact binary catalogues corresponding to each observation time are presented. The Radler galaxy is a training dataset, with binary parameters for every signal in the data stream included. We compare our derived catalogues to the LISA Data Challenge Radler catalogue to quantify the detection efficiency of the search pipeline. Included in the appendix is a more detailed analysis of two corner cases that provide insight into future improvements to our search pipeline.

Erratum: All-sky search in early O3 LIGO data for continuous gravitational-wave signals from unknown neutron stars in binary systems [Phys. Rev. D 103 , 064017 (2021)

Erratum: All-sky search in early O3 LIGO data for continuous gravitational-wave signals from unknown neutron stars in binary systems [Phys. Rev. D <b>103</b> , 064017 (2021)

Searches for continuous gravitational waves from neutron stars: A twenty-year retrospective

Seven years after the first direct detection of gravitational waves, from the collision of two black holes, the field of gravitational wave astronomy is firmly established. A first detection of continuous gravitational waves from rapidly-spinning neutron stars could be the field’s next big discovery. I review the last twenty years of efforts to detect continuous gravitational waves using the LIGO and Virgo gravitational wave detectors. I summarise the model of a continuous gravitational wave signal, the challenges to finding such signals in noisy data, and the data analysis algorithms that have been developed to address those challenges. I present a quantitative analysis of 297 continuous wave searches from 80 papers, published from 2003 to 2022, and compare their sensitivities and coverage of the signal model parameter space.

Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

ABSTRACT The International Pulsar Timing Array 2nd data release is the combination of data sets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95 per cent upper limits on their amplitude h95. The most sensitive frequency is 10 nHz with h95 = 9.1 × 10−15. We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit.

All-sky search for continuous gravitational waves from isolated neutron stars using Advanced LIGO and Advanced Virgo O3 data

We present results of an all-sky search for continuous gravitational waves which can be produced by spinning neutron stars with an asymmetry around their rotation axis, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. Four different analysis methods are used to search in a gravitational-wave frequency band from 10 to 2048 Hz and a first frequency derivative from −10−8 to 10−9 Hz/s. No statistically significant periodic gravitational-wave signal is observed by any of the four searches. As a result, upper limits on the gravitational-wave strain amplitude h0 are calculated. The best upper limits are obtained in the frequency range of 100 to 200 Hz and they are ∼1.1×10−25 at 95% confidence level. The minimum upper limit of 1.10×10−25 is achieved at a frequency 111.5 Hz. We also place constraints on the rates and abundances of nearby planetary- and asteroid-mass primordial black holes that could give rise to continuous gravitational-wave signals.10 MoreReceived 3 January 2022Accepted 17 October 2022DOI:https://doi.org/10.1103/PhysRevD.106.102008© 2022 American Physical SocietyPhysics Subject Headings (PhySH)Research AreasGravitational wave sourcesPhysical SystemsNeutron stars & pulsarsTechniquesGravitational wave detectionGravitation, Cosmology & Astrophysics

Search for continuous gravitational waves from HESS J1427-608 with a hidden Markov model

We present a search for continuous gravitational wave signals from an unidentified pulsar potentially powering HESS~J1427-608, a spatially unresolved TeV point source detected by the High Energy Stereoscopic System (HESS). The search uses a semi-coherent algorithm, which combines the maximum likelihood $\mathcal{F}$~statistic with a hidden Markov model to efficiently detect and track quasi-monochromatic signals that wander randomly in frequency. It uses data from the second observing run of the Advanced Laser Interferometer Gravitational-Wave Observatory. Multi-wavelength observations of the HESS source are combined with the proprieties of the population of TeV-bright pulsar wind nebulae to constrain the search parameters. We find no evidence of gravitational wave emission from this target. We set upper limits on the characteristic wave strain $h_0^{95\%}$ (for circularly polarized signals) at $95\%$ confidence level in sample sub-bands and interpolate it to estimate the sensitivity in the full band. We find $h_0^{95\%} = 1.3\times 10^{-25}$ near 185~Hz. The implied constraints on the ellipticity and \textit{r}-mode amplitude reach $\epsilon\leq 10^{-5}$ and $\alpha \leq 10^{-3}$ at 200~Hz, respectively.

Quality over quantity: Optimizing pulsar timing array analysis for stochastic and continuous gravitational wave signals

ABSTRACT The search for gravitational waves using Pulsar Timing Arrays (PTAs) is a computationally expensive complex analysis that involves source-specific noise studies. As more pulsars are added to the arrays, this stage of PTA analysis will become increasingly challenging. Therefore, optimizing the number of included pulsars is crucial to reduce the computational burden of data analysis. Here, we present a suite of methods to rank pulsars for use within the scope of PTA analysis. First, we use the maximization of the signal-to-noise ratio as a proxy to select pulsars. With this method, we target the detection of stochastic and continuous gravitational wave signals. Next, we present a ranking that minimizes the coupling between spatial correlation signatures, namely monopolar, dipolar, and Hellings &amp; Downs correlations. Finally, we also explore how to combine these two methods. We test these approaches against mock data using frequentist and Bayesian hypothesis testing. For equal-noise pulsars, we find that an optimal selection leads to an increase in the log-Bayes factor two times steeper than a random selection for the hypothesis test of a gravitational wave background versus a common uncorrelated red noise process. For the same test but for a realistic European PTA (EPTA) data set, a subset of 25 pulsars selected out of 40 can provide a log-likelihood ratio that is 89 % of the total, implying that an optimally selected subset of pulsars can yield results comparable to those obtained from the whole array. We expect these selection methods to play a crucial role in future PTA data combinations.

Impact of signal clusters in wide-band searches for continuous gravitational waves

In this paper we present a study of some relevant steps of the hierarchical frequency-Hough (FH) pipeline, used within the LIGO and Virgo Collaborations for wide-parameter space searches of continuous gravitational waves (CWs) emitted, for instance, by spinning neutron stars (NSs). Because of their weak expected amplitudes, CWs have not been still detected so far. These steps, namely the spectral estimation, the {\it peakmap} construction and the procedure to select candidates in the FH plane, are critical as they contribute to determine the final search sensitivity. Here, we are interested in investigating their behavior in the (presently quite) extreme case of signal clusters, due to many and strong CW sources, emitting gravitational waves (GWs) within a small (i.e. <1 Hz wide) frequency range. This could happen for some kinds of CW sources detectable by next generation detectors, like LISA, Einstein Telescope and Cosmic Explorer. Moreover, this possibility has been recently raised even for current Earth-based detectors, in some scenarios of CW emission from ultralight boson clouds around stellar mass black holes (BHs). We quantitatively evaluate the robustness of the FH analysis procedure, designed to minimize the loss of single CW signals, under the unusual situation of signal clusters. Results depend mainly on how strong in amplitude and dense in frequency the signals are, and on the range of frequency they cover. We show that indeed a small sensitivity loss may happen in presence of a very high mean signal density affecting a frequency range of the order of one Hertz, while when the signal cluster covers a frequency range of one tenth of Hertz, or less, we may actually have a sensitivity gain. Overall, we demonstrate the FH to be robust even in presence of moderate-to-large signal clusters.

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO-Virgo data

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO-Virgo run in the detector frequency band [10,2000] Hz have been used. No significant detection was found and 95% confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about 7.6×10−26 at ≃142 Hz. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass–boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.Received 9 April 2022Accepted 29 June 2022DOI:https://doi.org/10.1103/PhysRevD.106.042003© 2022 American Physical SocietyPhysics Subject Headings (PhySH)Research AreasCosmic rays & astroparticlesDark matterGeneral relativityGravitational wavesGravitation, Cosmology & Astrophysics

A Fast Data Processing Technique for Continuous Gravitational Wave Searches

This article discusses the potential advantages of a data processing technique for continuous gravitational wave signals searches in the data measured by ground-based gravitational wave interferometers. Its main advantage over other techniques is that it does not need to search over the signal’s direction of propagation. Although it is a “coherent method” (i.e., it coherently processes year-long data), it is applied to a data set obtained by multiplying the original time-series with a (half-year) time-shifted copy of it. As a result, the phase modulation due to the interferometer motion around the Sun is automatically canceled in the signal of the synthesized time-series. Although the resulting signal-to-noise ratio is not as high as that of a coherent search, it equals that of current hierarchical methods. In addition, since the signal search is performed over a parameters space of smaller dimensionality, the associated false-alarm probability should be smaller than those characterizing hierarchical methods and result in an improved likelihood of detection.

Search Methods for Continuous Gravitational-Wave Signals from Unknown Sources in the Advanced-Detector Era

Continuous gravitational waves are long-lasting forms of gravitational radiation produced by persistent quadrupolar variations of matter. Standard expected sources for ground-based interferometric detectors are neutron stars presenting non-axisymmetries such as crustal deformations, r-modes or free precession. More exotic sources could include decaying ultralight boson clouds around spinning black holes. A rich suite of data-analysis methods spanning a wide bracket of thresholds between sensitivity and computational efficiency has been developed during the last decades to search for these signals. In this work, we review the current state of searches for continuous gravitational waves using ground-based interferometer data, focusing on searches for unknown sources. These searches typically consist of a main stage followed by several post-processing steps to rule out outliers produced by detector noise. So far, no continuous gravitational wave signal has been confidently detected, although tighter upper limits are placed as detectors and search methods are further developed.

Search for Continuous Gravitational-wave Signals in Pulsar Timing Residuals: A New Scalable Approach with Diffusive Nested Sampling

Abstract Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov Chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution is to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semianalytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling has shown capability in coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including  ( 10 2 ) pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also highly beneficial to include DNest in current toolboxes of PTA analysis.

Electromagnetic signatures of dark photon superradiance

Black hole superradiance is a powerful tool in the search for ultra-light bosons. Constraints on the existence of such particles have been derived from the observation of highly spinning black holes, absence of continuous gravitational-wave signals, and of the associated stochastic background. However, these constraints are only strictly speaking valid in the limit where the boson's interactions can be neglected. In this work we investigate the extent to which the superradiant growth of an ultra-light dark photon can be quenched via scattering processes with ambient electrons. For dark photon masses $m_{\gamma^\prime} \gtrsim 10^{-17}\,{\rm eV}$, and for reasonable values of the ambient electron number density, we find superradiance can be quenched prior to extracting a significant fraction of the black-hole spin. For sufficiently large $m_{\gamma^\prime}$ and small electron number densities, the in-medium suppression of the kinetic mixing can be efficiently removed, and quenching occurs for mixings $\chi_0 \gtrsim \mathcal{O}(10^{-8})$; at low masses, however, in-medium effects strongly inhibit otherwise efficient scattering processes from dissipating energy. Intriguingly, this quenching leads to a time- and energy-oscillating electromagnetic signature, with luminosities potentially extending up to $\sim 10^{57}\,{\rm erg / s}$, suggesting that such events should be detectable with existing telescopes. As a byproduct we also show that superradiance cannot be used to constrain a small mass for the Standard Model photon.

Search for continuous gravitational waves from small-ellipticity sources at low frequencies

We present the results of an all-sky search for continuous gravitational wave signals with frequencies in the 20-500 Hz range from neutron stars with ellipticity of 1e-8. This frequency region is particularly hard to probe because of the quadratic dependence of signal strength on frequency. The search employs the Falcon analysis pipeline on LIGO O2 public data. Compared to previous Falcon analyses the coherence length has been quadrupled, with a corresponding increase in sensitivity. This enables us to search for small ellipticity neutron stars in this low frequency region up to 44 pc away. The frequency derivative range is up to 3e-13 Hz/s easily accommodating sources with ellipticities of 1e-7 and a corresponding factor of 10 increase in reach. New outliers are found, many of which we are unable to associate with any instrumental cause.

Fourier transform of the continuous gravitational wave signal

The direct detection of continuous gravitational waves from pulsars is a much anticipated discovery in the emerging field of multimessenger gravitational wave (GW) astronomy. Because putative pulsar signals are exceedingly weak large amounts of data need to be integrated to achieve desired sensitivity. Contemporary searches use ingenious ad hoc methods to reduce computational complexity. In this paper we provide analytical expressions for the Fourier transform of realistic pulsar signals. This provides description of the manifold of pulsar signals in the Fourier domain, used by many search methods. We analyze the shape of the Fourier transform and provide explicit formulas for location and size of peaks resulting from stationary frequencies. We apply our formulas to analysis of recently identified outlier at 1891.76 Hz.

Continuous Gravitational-Wave Data Analysis with General Purpose Computing on Graphic Processing Units

We present a new approach to searching for Continuous gravitational Waves (CWs) emitted by isolated rotating neutron stars, using the high parallel computing efficiency and computational power of modern Graphic Processing Units (GPUs). Specifically, in this paper the porting of one of the algorithms used to search for CW signals, the so-called FrequencyHough transform, on the TensorFlow framework, is described. The new code has been fully tested and its performance on GPUs has been compared to those in a CPU multicore system of the same class, showing a factor of 10 speed-up. This demonstrates that GPU programming with general purpose libraries (the those of the TensorFlow framework) of a high-level programming language can provide a significant improvement of the performance of data analysis, opening new perspectives on wide-parameter searches for CWs.

A method for detecting continuous gravitational wave signals from an ensemble of known pulsars

This work describes a method to improve the detection efficiency for continuous gravitational waves by combining information from known pulsars, which emission is too low to be individually detectable. We propose a new ensemble statistic t, which combines statistics from individual pulsars, defined through the five-vector method. Using a collection of band sample data files from the LIGO science run O2, we test the method with software injections in order to study the statistical properties of t. We also study the performance of our method by varying the number of injected signals and their amplitude. We propose and test a possible strategy for the application of the ensemble method to a set of real pulsars.