From System Parameters to Flares: Radiation Hydrodynamics of Star-Disk Collisions
Quasi-periodic eruptions (QPEs) — repeating soft X-ray flares in galactic nuclei — are thought to be caused by a star colliding with the accretion disk of a supermassive black hole. Three-dimensional radiation-hydrodynamic simulations of such collisions were performed, varying the surface density and vertical profile of the disk, the star’s speed and radius, and the intersection angle; the star is treated as unperturbed. The morphology of the shock wave and outflows is nearly independent of the star’s speed and disk density, but a faster star increases flare brightness, and a dense disk increases duration. Increasing the star’s radius enhances direct ejecta, making the flare brighter and longer. Disks with strong vertical concentration produce bright, short flares because radiation escapes more easily through the tenuous outer layers. Oblique intersections reduce the asymmetry of the two ejections and lengthen the event. Empirical relations for peak luminosity and duration as functions of parameters are derived. Applying to source GSN 069, the optimal model is a star with a radius about solar on a retrograde orbit, colliding with a dense disk after tidal disruption of another star, having a vertically concentrated profile.
Context
Quasi-periodic eruptions (QPEs)—a recently discovered class of nuclear transients linked to supermassive black holes—appear as bright recurrent X-ray flares with a still-mysterious mechanism. A model in which a star periodically crosses the accretion disk offers a natural explanation for the recurrence, duration, and alternation of bright and dim flares. Connecting local collision parameters—disk density, speed, and impact angle—to observed emission properties is crucial for testing this hypothesis and for future detection of such systems via gravitational waves with space-based interferometers. It has long been clear that matter’s behavior near the Schwarzschild radius requires accounting for both hydrodynamics and radiative transfer.
Methods
The researchers performed a series of 3D local simulations based on smoothed particle hydrodynamics (SPH) with radiation diffusion in the flux-limited diffusion approximation. The star was modeled as a rigid sphere crossing a local patch of the accretion disk at speeds up to 0.1 the speed of light. Five key parameters were varied: star speed, disk surface density, stellar radius, vertical density profile (from uniform to Gaussian with various dispersions), and local impact angle. For each configuration, the dynamics of the shock and ejecta were computed, along with the bolometric emission emerging through the photosphere.
Results
The analysis revealed the following patterns. Increasing the star speed (up to 0.15 the speed of light) boosts peak luminosity (L_peak ∝ v^2) and flare duration but does not alter shock morphology. Denser disks (with surface density Σ_d up to 100 times the fiducial value) prolong flares (Δt ∝ Σ_d^0.5) while only weakly affecting peak luminosity. A larger stellar radius enhances ejecta asymmetry: the forward flow becomes dominant, and its luminosity scales as R_⋆^1.5. Disks with a vertically concentrated profile (Gaussian distribution with dispersion σ) produce brighter and shorter flares due to efficient radiation escape through tenuous outer layers. Decreasing the local impact angle i (from 90° to 30°) reduces the momentum asymmetry of the two ejecta and lengthens the flares; for the reverse flow, L_peak ∝ (sin i)^-0.8. Applying these empirical scalings to the source GSN 069 showed the best match is achieved for a star with radius ~R_⊙ on a retrograde orbit (sin i ~ 0.4) and a dense post-tidal-disruption disk concentrated toward the midplane.
Implications
These results for the first time systematically connect local physical collision conditions with global observational properties of QPEs, paving the way to use these flares to probe accretion disk structure near supermassive black holes. The discovered dependence of brightness asymmetry on impact angle positions QPEs as a potential tool for determining orbital geometry and density distributions in disks, which could be refined by future high-resolution spectroscopy in the X-ray band.
Future development
Future model development will incorporate a realistic stellar structure, disk rotation and shear, and directional luminosities for different viewing angles. This will refine predictions for specific sources and lay the groundwork for interpreting data from upcoming missions like Athena. Moreover, concurrent detection of gravitational waves from such systems by next-generation laser interferometers will allow independent measurement of orbital parameters.
Impact
The results will impact high-energy astrophysics, accretion theory, and the dynamics of relativistic objects, as well as the planning of observing campaigns in the X-ray and gravitational-wave bands.
Next steps
Next steps include running simulations with more realistic equations of state, incorporating general relativistic effects, and developing semi-analytic models for rapid estimation of QPE parameters from observational data.
Key open problems
The study is directly linked to the unsolved problem of fueling active galactic nuclei and variability mechanisms in accretion disks. It also provides insight into the evolution of stellar populations in the vicinity of supermassive black holes and opens a new way to probe strong gravitational fields.
🎯 If such a star really collided with a disk at 0.1c, the released energy would be enough to momentarily outshine an entire galaxy, and the shock wave would generate temperatures of hundreds of millions of degrees.
🎬 The idea of an object periodically plunging into a scorching disk around a black hole recalls scenes from Interstellar, where the protagonists explore Gargantua’s accretion disk, but in reality such collisions produce recurring X-ray flares.
Key numbers
- star speed: 0.05–0.15 c (up to 45,000 km/s)
- peak QPE bolometric luminosity: ~10^41–10^43 erg/s
- flare duration: tens of minutes to several hours
- QPE period: hours to days
- radiation temperature: 100–200 eV (~1–2 million K)