Searches for Continuous Gravitational Waves from Supernova Remnants in the first part of the LIGO-Virgo-KAGRA Fourth Observing run

Using data from the first eight months of the LIGO-Virgo-KAGRA O4 run, this study conducted the most sensitive wide-band directed searches to date for continuous gravitational waves from 15 nearby supernova remnants, finding no evidence of signals and establishing the most stringent upper limits on intrinsic strain amplitude, neutron star ellipticity, and $r$-mode amplitudes for these sources.

The Gist

The Cosmic "Hum" Hunt: A Search for Ghostly Neutron Stars

Imagine the universe as a giant, silent ocean. For years, we've been listening for the loud, crashing waves caused by massive collisions—like two black holes smashing together. We've found those! But now, scientists are trying to hear something much quieter: a faint, continuous hum.

This paper is a report from a massive international team (the LIGO-Virgo-KAGRA Collaboration) who spent eight months in 2023–2024 listening for this hum. They weren't looking for a crash; they were looking for a lullaby sung by a specific type of cosmic object: a neutron star.

The Cast of Characters: What is a Neutron Star?

Think of a neutron star as a cosmic diamond. It's the crushed core of a star that exploded. If you took a mountain (like Mount Everest) and squished it down until it fit inside a city block, you'd have something with the density of a neutron star.

These objects spin incredibly fast—sometimes hundreds of times per second. If a neutron star is perfectly round, it spins silently. But if it has a tiny bump on its surface (even if that bump is smaller than an atom), it creates a wobble as it spins. This wobble sends out ripples in space-time called Gravitational Waves.

Because the star spins constantly, these ripples aren't a one-time "crash"; they are a continuous, steady tone, like a violin string being bowed forever. This is the "Continuous Gravitational Wave" (CW) the scientists are hunting.

The Target List: Supernova Remnants

The team focused on 15 specific locations in the sky known as Supernova Remnants.

• The Analogy: Imagine a supernova is a massive firework that just exploded. The "remnant" is the smoke and debris left behind.
• The Theory: Inside the smoke of these fireworks, there should be a new, baby neutron star hiding.
• The Problem: We can't see these stars with telescopes because they don't flash light like normal pulsars. They are "ghosts" hiding in the debris. The scientists hoped to find them by listening for their gravitational hum.

The Detective Work: Five Different Ears

To catch this faint signal, the team didn't just use one method. They used five different "pipelines" (search strategies), like using five different types of microphones to find a whisper in a noisy room.

1. The "Pattern Matcher" (BSD): This looks for specific patterns in the data, like a detective looking for a specific fingerprint.
2. The "Cross-Checker" (PyStoch): This compares the noise from two different detectors (in Washington and Louisiana) to see if they hear the same thing at the same time.
3. The "Horse-Tracker" (Viterbi): This method is smart. It knows that these stars might wobble or change speed slightly over time. It tracks the signal like a horse running a race, adjusting its prediction as the race goes on.
4. The "Double-Track" (Dual-Harmonic): Some stars might sing two notes at once (a fundamental note and a higher harmony). This pipeline listens for both.
5. The "Template Hunter" (Weave): This uses a massive library of pre-calculated "templates" (like a sheet music book) to see if the data matches any known song.

The Results: The Great Silence

After listening to the data for eight months, the answer was: Nope.

They didn't find any gravitational waves.

• The Analogy: It's like standing in a crowded stadium, holding a very sensitive microphone, and listening for a specific person humming a tune. You hear the roar of the crowd (noise), but you don't hear the hum.

Does this mean they failed?
Absolutely not! In science, a "null result" is still a huge discovery. Here's what they actually found:

1. The "Silence is Loud" Discovery: Even though they didn't hear the hum, they proved that if the stars were humming, it would have to be extremely quiet.
2. Setting the Volume Limit: They set a "volume limit" for the universe. They can now say, "If there is a neutron star in that cloud, its gravitational wave volume must be lower than X."
   • For the closest target (Vela Jr.), they proved the star is smoother than a billiard ball. If it had a bump, it would have to be smaller than a single atom.
3. Beating the Competition: Their search was the most sensitive ever done for these specific targets. They improved upon previous searches by a factor of 1.5 to 2, meaning their "microphones" are much better now.

Why Does This Matter?

Even though they didn't find the signal, this paper is a victory lap for physics.

• Ruling Out Bad Ideas: By proving the stars aren't wobbling as much as some theories predicted, they are helping physicists understand how neutron stars are built.
• The Future: They are building better "ears" for the next round of listening. Every time they don't find a signal, they learn more about what not to expect, narrowing down the search for the day they finally hear that cosmic hum.

In short: The universe is still keeping its secrets, but thanks to this paper, we know exactly how quiet those secrets are, and we are getting better at listening.

Technical Summary

1. Problem Statement

Continuous gravitational waves (CWs) are long-duration, nearly monochromatic signals expected from non-axisymmetric, rapidly rotating neutron stars (NSs). While transient gravitational waves from binary mergers have been detected, CWs remain elusive.

• Target: Young to middle-aged core-collapse supernova remnants (SNRs) are prime candidates because they likely host isolated neutron stars (Central Compact Objects or Pulsar Wind Nebulae) that have not yet been identified as radio pulsars.
• Challenge: These signals are orders of magnitude weaker than merger signals. Furthermore, the lack of electromagnetic timing data (spin frequency $f$ and spin-down $\dot{f}$) for these hidden NSs requires "directed" searches over wide parameter spaces, which is computationally expensive.
• Goal: To conduct the most sensitive wide-band directed searches to date for CWs from 15 nearby SNRs using data from the first eight months (O4a) of the LIGO-Virgo-KAGRA (LVK) fourth observing run.

2. Methodology

The analysis utilized data from the LIGO Hanford and Livingston detectors collected between May 24, 2023, and January 16, 2024. The study employed five complementary search pipelines, each with different signal models and semi-coherent strategies to cover various emission scenarios (e.g., triaxial deformation, r-modes, spin wandering).

Target Selection

• 15 SNRs were selected from the Green catalogue and SNRcat, including well-known sources like Cas A, Vela Jr., and SNR 1987A.
• Targets were chosen based on the presence of a Central Compact Object (CCO) or Pulsar Wind Nebula (PWN) within the remnant.
• Updates: SNR 1987A was added (now in the remnant stage), while G354.4+0.0 was removed (reclassified as an HII region).

Search Pipelines

1. Band-Sampled-Data (BSD): A semi-coherent search using the Frequency-Hough transform. It targets 8 sources (including new additions G291.0−0.1 and G330.2+1.0) in the 10–600 Hz band, searching for both spin-down and spin-up scenarios.
2. PyStoch: A cross-correlation-based radiometer method. It analyzed 4 specific targets (Cas A, Vela Jr., SNR 1987A, G347.3−0.5) across a broad band (20–1726 Hz) using folded data. It employs adaptive binning to account for signal frequency drift.
3. Single-Harmonic Viterbi (SHV): Uses a Hidden Markov Model (HMM) to track signals with secular spin-down and stochastic spin wandering. Applied to 14 targets.
4. Dual-Harmonic Viterbi (DHV): An extension of SHV that simultaneously tracks two frequency components ($f_\star$ and $2f_\star$), improving sensitivity for specific emission models. Applied to 7 targets.
5. Weave: A template-based semi-coherent search using the $\mathcal{F}$-statistic. It focused on Cas A, Vela Jr., and G347.3−0.5 in the 20–475 Hz range, assuming power-law spin-down with braking indices $n \in [2, 7]$.

Data Processing

• Data underwent glitch gating and cleaning.
• Candidates surviving initial thresholds were subjected to rigorous follow-up procedures, including vetoes against known instrumental lines, consistency checks between detectors, and coherence-time increases to distinguish astrophysical signals from noise.

3. Key Contributions

• Multi-Pipeline Approach: The simultaneous application of five distinct pipelines provides robust coverage of the parameter space, mitigating the risk of missing signals due to model mismatches.
• O4a Sensitivity: This represents the first directed CW search using LVK O4a data, leveraging improved detector sensitivity (factor of $\sim 1.5$ over O3) and a longer observation time (8 months vs. 6 months in O3a).
• Target Updates: Inclusion of SNR 1987A and the re-evaluation of target ages/distances based on the 2025 SNRcat updates.
• Addressing Previous Candidates: The search specifically investigated a candidate reported by the Einstein@Home project for G347.3−0.5 (Ming et al. 2025). The authors found no compatible candidates with their search configurations, suggesting the previous candidate was likely a noise fluctuation or required different search parameters.

4. Results

• No Detection: No statistically significant continuous gravitational wave signals were detected from any of the 15 targets. All surviving candidates were consistent with instrumental artifacts or noise fluctuations.
• Upper Limits: The study set 95% confidence-level upper limits on the intrinsic strain amplitude ($h_0$).
  • Most Stringent Limit: For the nearby source Vela Jr. (G266.2−1.2), the BSD pipeline achieved a limit of $h_0 \approx 5.94 \times 10^{-26}$ near 250 Hz. The Weave pipeline reached $\approx 3.92 \times 10^{-26}$ near 266 Hz.
  • Comparison: These limits surpass previous O3 results by a factor of roughly 1.5 to 2, depending on the source and frequency.
  • Spin-Down Limits: For several sources (e.g., G65.7+1.2, G189.1+3.0, Vela Jr.), the achieved upper limits surpassed the indirect "spin-down limit" (the maximum possible strain if all energy loss were due to GWs), particularly under optimistic assumptions of young, nearby sources.
• Astrophysical Constraints:
  • Ellipticity ($\epsilon$): Limits were converted to constraints on the neutron star's equatorial ellipticity. For Vela Jr., constraints reached $\epsilon \lesssim 10^{-7}$ at frequencies $>400$ Hz, approaching or exceeding theoretical limits for normal neutron stars.
  • r-mode Amplitude ($\alpha$): Constraints on r-mode oscillation amplitudes reached $\alpha \lesssim 10^{-5}$ (and down to $\sim 10^{-6}$ for SHV at high frequencies), significantly surpassing theoretical predictions for nonlinear saturation mechanisms ($\alpha \sim 10^{-3}$).

5. Significance

• State-of-the-Art Sensitivity: This work establishes the current state-of-the-art for directed CW searches from SNRs, providing the tightest constraints to date on the deformation and energy loss mechanisms of young neutron stars.
• Exclusion of Hidden Populations: By ruling out signals down to specific strain amplitudes, the results constrain the population of hidden, non-pulsing neutron stars in the Galaxy.
• Methodological Benchmark: The successful integration of multiple pipelines and the rigorous follow-up of candidates (including the refutation of a high-profile candidate from another collaboration) demonstrates the maturity of the LVK search infrastructure.
• Future Outlook: The results highlight the potential of future observing runs (O4b, O5) with longer datasets and enhanced sensitivity to probe even deeper into the parameter space, potentially reaching the theoretical limits of neutron star deformability.