Observational constraints on the spin/anisotropy of the CCOs of Cassiopeia A, Vela Jr. and G347.3-0.5 and a single surviving continuous gravitational wave candidate

Using Einstein@Home computing resources to analyze LIGO data from observing runs O3a through O4a, this study conducts the most sensitive search to date for continuous gravitational waves from three central compact objects in supernova remnants, setting unprecedented constraints on their ellipticity and crustal anisotropy while identifying a single surviving candidate from G347.3-0.5 that warrants further investigation with upcoming data.

The Gist

The Big Picture: Listening for a Whisper in a Hurricane

Imagine the universe is a massive, chaotic concert hall. Most of the time, it's filled with the deafening roar of black holes smashing together (these are the "bursts" scientists usually look for). But hidden in that noise, there might be a faint, steady hum—a single note played by a spinning neutron star that never stops. This is a Continuous Gravitational Wave (CW).

The problem? This hum is incredibly quiet. It's like trying to hear a single person whispering a secret in the middle of a roaring hurricane. To hear it, you need to listen for a very long time and use a very clever trick to separate the whisper from the wind.

This paper is about a team of scientists (and thousands of volunteers) who tried to find this whisper coming from three specific "concert halls" in our galaxy: Cassiopeia A, Vela Jr., and G347.3. These are the leftovers of massive stars that exploded recently (in cosmic terms).

The Cast of Characters

1. The Targets (The Whisperers):

   • Cassiopeia A (Cas A): A young, violent explosion remnant.
   • Vela Jr.: A very young, nearby remnant.
   • G347.3: An ancient remnant, possibly from a star seen by Chinese astronomers in 393 AD.
   • Why them? These are young neutron stars. Because they are young, they might still be wobbling or spinning fast enough to make a gravitational "hum."

2. The Volunteers (The Super-Computers):

   • The scientists couldn't do this alone. The math required to listen for these whispers is so heavy it would take a supercomputer centuries to crunch.
   • Instead, they used Einstein@Home. This is a project where regular people (like you and me) donate the idle time on their home computers to help crunch the numbers. Think of it as a global choir of thousands of computers singing in unison to find the signal.

The Detective Work: How They Searched

The search was like a multi-stage sieve, getting finer and finer at every step.

• Stage 0 (The Wide Net): They cast a huge net over the data from the LIGO detectors (the ears that hear gravitational waves). They looked at about 45 million potential "whispers" (candidates). This stage was fast but not super precise.
• The Follow-Up (The Magnifying Glass): Most of those 45 million candidates were just noise (static). The team took the top survivors and looked at them with a much sharper lens. They used more data and longer listening times.
  • Analogy: Imagine finding a suspicious footprint in the sand. Stage 0 says, "Hey, look at that footprint!" Stage 1 says, "Let's measure the length." Stage 2 says, "Let's check the soil composition." Stage 3 says, "Let's see if it matches a specific shoe."
• The Result: After all that filtering, 44,999,999 candidates were proven to be just noise. Only one candidate survived the entire process.

The Lone Survivor: The Mystery Candidate

One candidate from the G347.3 search survived all the tests.

• Is it a real signal? Maybe. It's "marginal," meaning it's just barely loud enough to be interesting, but not loud enough to be a confirmed discovery yet.
• The Twist: The scientists checked the data again and again. Sometimes it looks like a signal; sometimes it looks like a glitch. It's like finding a shadow that looks like a person, but when you walk closer, it might just be a coat rack.
• What's next? They need new data (from a later time period) to confirm if this shadow is real. Unfortunately, that new data isn't public yet, so the mystery remains unsolved for now.

The Real Victory: Setting the "Silence" Limits

Even though they didn't find a confirmed gravitational wave, this paper is a massive success. Why? Because they set the strictest rules ever on what these neutron stars can't be doing.

Think of it like a police investigation where they didn't catch the criminal, but they proved the criminal couldn't be wearing a red hat, couldn't be taller than 5 feet, and couldn't be driving a blue car.

Here is what they learned about the "shape" of these stars:

1. The "Bumpy" Star (Ellipticity):

   • If a neutron star has a "mountain" on it (even a tiny one), it wobbles as it spins, creating gravitational waves.
   • The Finding: For stars spinning faster than a certain speed, the "mountains" must be smaller than a grain of sand on a planet the size of Earth. If they were any bigger, we would have heard the hum.

2. The "Wobbly" Star (R-Modes):

   • Imagine a spinning top that starts to wobble violently.
   • The Finding: These stars aren't wobbling nearly as much as some theories predicted. They are much more stable than we thought.

3. The "Cracked" Star (Anisotropy):

   • This is the new discovery in this paper. Neutron stars have a crust (a solid shell). Scientists wondered if this crust is "stiff" in one direction and "squishy" in another (anisotropic). If it is, the star's spin could change its shape slightly.
   • The Finding: They proved that the crust of these stars is likely very uniform. It's not "cracked" or "stiff" in weird directions. If it were, we would have detected a signal. This is the first time anyone has put a limit on this specific property of neutron star crusts.

The Bottom Line

• Did they find a gravitational wave? Not definitively. They found one "maybe," but it needs more proof.
• Did they fail? Absolutely not. They used the power of thousands of volunteers to listen harder than ever before.
• What did they achieve? They proved that the neutron stars in these three supernova remnants are incredibly smooth, stable, and symmetrical. They ruled out many wild theories about how these stars are shaped.

In short: They didn't find the ghost, but they proved the house is empty, and they mapped out exactly where the ghost couldn't be hiding. That is a huge step forward in understanding the universe.

Technical Summary

1. Problem Statement

The paper addresses the search for Continuous Gravitational Waves (CWs) emitted by three young neutron stars located at the centers of supernova remnants (SNRs): Cassiopeia A (Cas A), Vela Jr., and G347.3-0.5.

• Context: Despite a decade of gravitational wave observations, no CWs have been detected. Detecting them requires integrating detector strain data over long periods (months) to accumulate signal power, which is computationally prohibitive for large parameter spaces using fully coherent methods.
• Targets: These specific SNRs are prime targets because their embedded neutron stars (Central Compact Objects or CCOs) are young, potentially possessing large non-axisymmetric deformations (ellipticity) or unstable $r$-modes that could power strong CW emission. Their sky positions are known, allowing for "directed" searches that are more sensitive than all-sky surveys.
• Goal: To perform the deepest and broadest search to date for CWs from these three sources, setting stringent upper limits on gravitational wave amplitude, equatorial ellipticity, $r$-mode saturation, and, for the first time, neutron-star crustal anisotropy.

2. Methodology

The analysis employs a hierarchical, semi-coherent search strategy leveraging distributed computing and advanced follow-up techniques.

Data Sets

• Initial Search: Used data from the first half of the third LIGO observing run (O3a, April–October 2019).
• Follow-up: Used independent data from the second half of O3 (O3b) and the first part of the fourth observing run (O4a, May 2023–January 2024) to vet candidates.
• Preprocessing: Data was cleaned of instrumental lines (calibration lines, mains harmonics) and short-duration glitches.

Search Pipeline

1. Stage 0 (Einstein@Home):

   • Platform: Distributed computing via thousands of volunteers (BOINC architecture).
   • Method: A "stack-slide" semi-coherent search using the Global Correlatation Transform (GCT).
   • Statistic: Candidates were ranked using a line-robust statistic ($\hat{\beta}_{S/GLtL}$) to distinguish signals from stationary noise and instrumental lines.
   • Parameter Space: Covered frequencies $20\text{--}1500$ Hz. The search included non-standard spin-down models (braking index $n$) to account for timing noise and glitches.
   • Scale: Approximately $2 \times 10^{18}$ waveform templates were searched across $\sim 45$ million initial candidates.

2. Hierarchical Follow-up (Stages 1–4):

   • Pixeling: A "pixeling" procedure selected the top 5 candidates from $100 \times 10$ sub-volumes (pixels) in the $(f, \dot{f})$ plane within $0.5$ Hz bands to ensure uniform coverage.
   • Progressive Coherence:
     • Stage 1: Semi-coherent search on O3a (longer coherence time $T_{coh}$).
     • Stage 2: Fully coherent search on O3a.
     • Stage 3: Fully coherent search on combined O3a + O3b.
     • Stage 4: Fully coherent search on O4a data.
   • Growth Ratio ($R_a$): Candidates were vetted based on the growth of their detection statistic relative to the coherent integration time. Candidates failing to show signal-like growth were rejected.

Astrophysical Constraints

The authors derived constraints on physical parameters using the upper limits on strain amplitude ($h_0$):

• Ellipticity ($\epsilon$): Derived from $h_0$ assuming a triaxial rotator.
• $r$-mode Amplitude ($\alpha$): Derived assuming emission from unstable $r$-modes.
• Crustal Anisotropy ($\langle \phi \rangle$): A novel constraint linking ellipticity to the anisotropy of the neutron star crust, assuming a specific spin-down history from birth frequency $\nu_0$ to current frequency $\nu$.

3. Key Contributions

• Deepest Search to Date: This study provides the most sensitive upper limits on CWs from Cas A, Vela Jr., and G347.3, improving upon previous results by factors ranging from 25% to 4.5 times depending on the frequency band and target.
• First Constraints on Crustal Anisotropy: The paper presents the first observational constraints on the degree of anisotropy ($\langle \phi \rangle$) in neutron star crusts.
• Candidate Survival: Only one candidate survived the entire hierarchical pipeline, originating from the G347.3 search.
• Optimization Framework: The study utilizes an optimized allocation of computational resources (J. Ming et al. 2016) to maximize detection probability across different frequency bands and targets.

4. Results

Candidate Status

• Surviving Candidate: One candidate from the G347.3 low-frequency search ($f \approx 31.7$ Hz) survived all four stages.
  • Significance: The candidate has a growth ratio $R_4 = 2.3$ (marginal) and a Gaussian $p$-value of $\approx 0.75\%$. However, after accounting for the multiple trials (13 follow-ups), the false-alarm probability is $\approx 10\%$.
  • Nature: It is not considered a confirmed detection. Phase consistency between O3 and O4a data is weak, and the O4a evidence is lower than O3. Further investigation with newer data (O4b/O4c) is required but currently unavailable to the public.
• Other Targets: All candidates from Cas A and Vela Jr. were rejected during the follow-up stages.

Upper Limits and Astrophysical Constraints

• Gravitational Wave Amplitude ($h_0$):
  • Cas A: Improved depth by 42% (200–500 Hz) and 25% (500–978 Hz) over previous best results. Above 978 Hz, sensitivity improved by a factor of 2.8.
  • Vela Jr.: Improved depth by 37% up to 978 Hz.
  • G347.3: Improved depth by 93% up to 1300 Hz.
• Equatorial Ellipticity ($\epsilon$):
  • For spin periods $< 2$ ms, $\epsilon$ is constrained to be $< 4 \times 10^{-7}$ for all targets.
  • For Vela Jr. (assuming optimistic distance of 200 pc), $\epsilon \lesssim 4 \times 10^{-8}$ at 200 Hz.
• $r$-mode Amplitude ($\alpha$):
  • Excludes $\alpha \gtrsim 4 \times 10^{-5}$ for spin frequencies $> 250$ Hz.
  • For Vela Jr. (200 pc), $\alpha \lesssim 2 \times 10^{-6}$.
• Crustal Anisotropy ($\langle \phi \rangle$):
  • Key Finding: For spin periods between 1.3 ms and 100 ms, crustal anisotropy values larger than $5 \times 10^{-3}$ are excluded.
  • If the anisotropy is as high as $5 \times 10^{-3}$, the neutron stars in these targets must be spinning slower than 10 Hz or faster than 750 Hz (which is unlikely given their ages).
  • Assuming a more conservative anisotropy of $5 \times 10^{-5}$, spin frequencies above 391–580 Hz are excluded.

5. Significance

• Physical Implications: The constraints on crustal anisotropy provide the first observational test of theoretical models regarding the internal structure and "mountain" formation of neutron stars. The results suggest that if these young neutron stars have significant crustal anisotropy, they are likely spinning much slower than the frequencies probed, or the anisotropy is significantly lower than some theoretical maximums.
• Methodological Advancement: The successful application of the pixeling strategy and the hierarchical follow-up using independent data (O3b, O4a) demonstrates a robust framework for reducing false alarms in deep CW searches.
• Future Outlook: The single surviving G347.3 candidate highlights the need for continued monitoring. The authors note that upcoming public data releases (O4b, O4c) are critical for confirming or refuting this candidate, as the current O4a data is insufficient for a definitive conclusion.

In summary, this paper represents a major step forward in the search for continuous gravitational waves, setting the most stringent limits to date on the deformations and internal anisotropy of young neutron stars, while maintaining rigorous statistical standards for candidate vetting.