A positive period derivative in the quasi-periodic eruptions of ZTF19acnskyy

This paper reports the first direct measurement of a positive period derivative in the quasi-periodic eruptions of ZTF19acnskyy, a finding that challenges standard extreme mass-ratio inspiral models and remains unexplained by current theoretical scenarios.

The Gist

The Cosmic Mystery: A Clock That's Getting Slower

Imagine you are watching a cosmic lighthouse in a distant galaxy. Every few days, it flashes a brilliant burst of X-ray light. For a long time, astronomers thought these flashes were like a perfect metronome, ticking away at a steady rhythm.

But recently, astronomers discovered something weird about a specific lighthouse named "Ansky" (scientifically known as ZTF19acnskyy). The flashes aren't just ticking; they are slowing down.

In fact, the time between the flashes is getting longer by about 2.5 minutes every single day. Over the course of a year, the gap between flashes has stretched from about 9.5 days to nearly 10 days. This is a huge deal because, in the universe, things usually speed up as they get closer to a black hole (like a figure skater pulling in their arms), not slow down.

The Main Characters

• The Black Hole (The Monster): At the center of this galaxy sits a Supermassive Black Hole (SMBH). It's millions of times heavier than our Sun.
• The Companion (The Dancer): Orbiting this monster is a smaller object, likely a star. It's like a tiny dancer spinning around a giant.
• The Eruptions (The Fireworks): Every time the dancer gets close to the monster, it creates a massive explosion of light (a "Quasi-Periodic Eruption" or QPE).

The Big Question: Why is the Dance Slowing Down?

Usually, when a star orbits a black hole, it loses energy and spirals inward, making the orbit faster and faster. This paper asks: Why is Ansky's star spiraling outward, making the orbit slower and slower?

The authors tested four main theories, like detectives trying to solve a crime:

1. The "Leaking Balloon" Theory (Mass Transfer)

The Idea: Imagine the star is a balloon that is slowly leaking air every time it gets close to the black hole. As it loses mass, the physics of the orbit changes, causing it to drift outward.
The Problem: To explain how much the orbit is slowing down, the star would have to lose about 10% of its total weight every year. That's like a human losing 150 pounds in a single year and still being able to run a marathon. While possible, it's extreme, and it's unclear if the star could survive that much weight loss without falling apart.

2. The "Rocket Kick" Theory (Tidal Kicks)

The Idea: When the star gets close to the black hole, the black hole's gravity might rip a chunk of the star off. If that chunk flies off in one direction, the remaining star gets a "kick" in the opposite direction (like a rocket firing its thrusters), pushing it into a wider, slower orbit.
The Problem: To get the star to slow down this much, it would need to be kicked incredibly hard. But if it were kicked that hard, the star would likely be destroyed or the explosions would look very different than what we see. It's like trying to push a car with a feather; the force just isn't there.

3. The "Spinning Top" Theory (Relativistic Precession)

The Idea: Maybe the orbit isn't actually changing speed. Instead, imagine the orbit is shaped like a slightly squashed circle (an ellipse) that is slowly rotating, like a spinning top wobbling. Because the orbit is wobbling, the light from the star has to travel a slightly different distance to reach us each time, making it look like the timing is changing.
The Problem: To create the effect we see, the wobbling would have to happen over a period of 11 to 27 years. But the math says the wobble should happen much faster. It's like trying to explain a fast-spinning fan by saying it's actually a slow-moving clock hand. The numbers just don't add up.

4. The "Hidden Neighbor" Theory (Binary Black Holes)

The Idea: Maybe the black hole isn't alone. Maybe there is a second, hidden black hole nearby. As the first black hole orbits the second one, it moves toward and away from us, changing the time it takes for the light to reach Earth (like the Doppler effect with a siren).
The Problem: To cause the delay we see, the second black hole would have to be incredibly massive and moving very fast. The math shows this scenario would require a setup that is physically unlikely given what we know about the galaxy.

The Conclusion: We Are Still Stumped

The authors ran the numbers on all these theories, and none of them fit perfectly.

• The "Leaking Balloon" requires the star to lose too much mass too fast.
• The "Rocket Kick" requires forces that would destroy the star.
• The "Spinning Top" and "Hidden Neighbor" theories require physics that don't quite match the observations.

What does this mean?
It means our current understanding of how stars and black holes interact is incomplete. The universe is throwing us a curveball. The star in Ansky is doing something we haven't seen before.

The authors conclude that we need to keep watching this cosmic lighthouse. As the star continues its slow dance, it might reveal a new physical law or a new type of interaction between stars and black holes that we haven't even imagined yet.

In short: We found a cosmic clock that is ticking slower every day. We have four guesses for why, but none of them are the right answer yet. The mystery is still unsolved!

Technical Summary

1. Problem Statement

Quasi-periodic eruptions (QPEs) are recurring soft X-ray transients originating from the nuclei of nearby galaxies, typically modeled as Extreme Mass-Ratio Inspirals (EMRIs) where a stellar-mass companion interacts with a Supermassive Black Hole (SMBH). Standard EMRI models predict that gravitational wave emission and tidal dissipation should cause the orbital period to decrease over time (inspiral).

The source ZTF19acnskyy ("Ansky") presents a unique anomaly. While previous observations (2024) showed recurrence times of 4.5–7 days, continued monitoring in 2025–2026 revealed a smoothly increasing period. The central problem addressed by this paper is the measurement and physical interpretation of this positive period derivative ($\dot{P}$), which contradicts standard inspiral expectations and challenges existing QPE models.

2. Methodology

The authors conducted a comprehensive timing analysis of Ansky using multi-mission X-ray data:

• Observations: Data were collected from NICER (263 observations, Jan–Jun 2025), Swift/XRT, and XMM-Newton (July 2025).
• Data Reduction:
  • NICER data were analyzed using the SCORPEON spectral model over 0.25–10 keV (night) and 0.38–10 keV (day).
  • Swift fluxes were converted to cgs units using a specific blackbody conversion factor.
  • XMM-Newton data were reduced using SAS v21.0.0, extracting light curves in the 0.3–2 keV band.
• Timing Analysis:
  • Peak arrival times ($t_{peak}$) for 28 bursts (23 observed in 2025–2026, including 19 consecutive) were fitted using a four-parameter exponential rise-and-decay model.
  • A systematic uncertainty of 0.1 days was added in quadrature to MCMC errors to account for model inaccuracies in describing peak fluxes.
  • O-C (Observed minus Calculated) Diagrams: The authors compared observed peak times against calculated times assuming a fixed period. They tested four distinct timing models:
    1. Constant period derivative ($\dot{P}$).
    2. Long-period sinusoidal oscillation.
    3. $\dot{P}$ plus short-period sinusoidal modulation.
    4. Hierarchical oscillations (long + short period).

3. Key Results

• Measurement of $\dot{P}$: The data strongly favor a model with a constant positive period derivative.
  • Period ($P_0$): $\approx 9.50 \pm 0.02$ days.
  • Derivative ($\dot{P}$): $(1.7 \pm 0.02) \times 10^{-2}$ days per day.
  • This implies the period increases by roughly 1.5 hours every month.
• Model Comparison:
  • Models including a constant $\dot{P}$ (Models 1 and 3) provided the best statistical fit ($\chi^2_{\nu} \approx 2.9$).
  • Pure sinusoidal models (Models 2 and 4) could mathematically mimic the trend over the ~30-epoch baseline but require physically extreme parameters (oscillation periods of 11–27 years and amplitudes of 1000–4500 days) to fit the data.
  • A secondary short-period oscillation ($P_s \approx 155$ days, amplitude $\approx 0.8$ days) was detected in Models 3 and 4, though its physical reality is uncertain given the limited baseline (~1.2 cycles).
• Backward Extrapolation: Extrapolating the 2025–2026 trend backward suggests the 4.5–7 day recurrence times observed in 2024 correspond to epochs roughly 20 cycles prior, providing a consistent evolutionary picture.

4. Physical Interpretations and Constraints

The authors evaluated several physical mechanisms to explain the positive $\dot{P}$, finding significant challenges for each:

• Mass Transfer (EMRI):
  • Mechanism: Impulsive mass loss at pericenter drives orbital expansion.
  • Constraint: To achieve the observed $\dot{P}$, the star must lose $\sim 10^{-3}$ of its mass per orbit, totaling $\sim 10\%$ of its mass within a year. While energetically plausible, this raises questions about the star's structural stability and whether it can sustain constant luminosity eruptions after such rapid mass loss. Additionally, the viscous timescale at pericenter ($\sim 58$ days) conflicts with the observed short burst durations ($\sim 3$ days).
• Velocity Kicks (Tidal Interactions):
  • Mechanism: Asymmetric mass loss or core reformation during partial TDEs imparts a velocity kick, increasing orbital energy.
  • Constraint: Hydrodynamical simulations suggest kicks large enough to produce the observed $\dot{P}$ require deep plunges ($\beta \gtrsim 0.62$). However, such deep plunges typically lead to total disruption or rapid orbital decay, making it difficult to reconcile with the $\sim 30$ consecutive, stable-luminosity eruptions observed.
• General Relativistic Precession:
  • Mechanism: Light travel-time delays due to apsidal (Schwarzschild) or nodal (Lense-Thirring) precession create apparent period changes.
  • Constraint: Calculations show that standard apsidal precession produces $\dot{P}$ values orders of magnitude too small. While high eccentricity ($e \gtrsim 0.97$) could theoretically match the amplitude, it would result in precession timescales much shorter than the observed stability of $\dot{P}$ over a year.
• Hierarchical SMBH Binary:
  • Mechanism: Reflex motion of the SMBH in a binary system causes light travel-time delays.
  • Constraint: The required binary parameters to produce the observed $\dot{P}$ amplitude and long timescale are inconsistent with the system's mass constraints and viewing geometry.
• Disk Instability Models:
  • Mechanism: Thermal/viscous instabilities in an accretion disk driven by a slowly increasing accretion rate ($\dot{m}$).
  • Constraint: This model can qualitatively reproduce the period increase if $\dot{m}$ crosses a critical threshold. However, the observed quiescent X-ray luminosity implies an accretion rate lower than the critical threshold required for the instability, creating a tension that requires further investigation.

5. Significance and Conclusion

• First Direct Measurement: This paper reports the first direct measurement of a positive period derivative in a QPE source, challenging the prevailing EMRI inspiral paradigm.
• Model Failure: No single existing physical model (EMRI mass transfer, tidal kicks, relativistic precession, or disk instability) provides a complete and consistent explanation for the magnitude, stability, and duration of the observed $\dot{P}$.
• Future Outlook: The authors conclude that the "extreme mass-ratio outspiral" (or a combination of mechanisms) may be at play. The constant nature of $\dot{P}$ over a year rules out short-term oscillations but leaves open the possibility of very long-term sinusoidal modulations. Continued monitoring is essential to detect potential changes in $\dot{P}$ (e.g., a sign change predicted by precession models) or luminosity evolution (predicted by mass-loss models), which will be crucial for distinguishing between these competing theories.

This work highlights the complexity of QPE phenomenology and suggests that current theoretical frameworks may need significant extension to account for sources like Ansky.