Unified Mass-Scaled QPO Signatures of Kerr Sen Black Holes from Stellar Mass to Supermassive Sources

Technical Summary: Unified Mass-Scaled QPO Signatures of Kerr–Sen Black Holes

Problem Statement
While General Relativity (GR) and the Kerr metric successfully describe rotating black holes, the near-horizon region remains a critical testing ground for potential deviations from standard GR. Modified gravity theories, such as heterotic string theory, predict black hole solutions with additional parameters, such as electric charge. The Kerr–Sen black hole, a charged rotating solution in the low-energy limit of heterotic string theory, introduces a charge-related deformation parameter ($q = Q^2/M$) alongside mass and angular momentum. The central problem addressed is how this spacetime deformation affects the hydrodynamical accretion flow and the resulting temporal signatures (Quasi-Periodic Oscillations, or QPOs) compared to the standard Kerr solution. Specifically, the study seeks to determine if Kerr–Sen deformations produce distinct, observable hydrodynamical and timing features that can be distinguished from Kerr spacetime across a wide range of black hole masses.

Methodology
The authors employ a numerical approach to solve the General Relativistic Hydrodynamics (GRH) equations for Bondi–Hoyle–Lyttleton (BHL) accretion in a fixed background spacetime.

• Spacetime Models: Simulations are conducted for both Kerr and Kerr–Sen black holes. Two representative spin families are analyzed: a rapidly rotating case ($a = 0.9M$) and a moderately rotating case ($a = 0.5M$). For the Kerr–Sen models, the charge-related deformation parameter $q$ is varied to create "soft," "moderate," and "strong" modification regimes.
• Numerical Setup: The simulations are performed on the equatorial plane using a high-resolution shock-capturing scheme (Marquina flux formula with MUSCL reconstruction). The computational domain extends from just outside the event horizon ($r_{min} \approx 2.3M$) to $100M$. Matter is injected supersonically ($Mach \approx 2$) from the outer boundary, leading to the formation of a shock cone.
• Analysis Techniques:
  • Hydrodynamical Analysis: The morphology of the shock cone, azimuthal density distributions, and time-dependent mass-accretion rates ($\dot{M}$) are examined.
  • Timing Analysis: Power Spectral Density (PSD) analysis is performed on the $\dot{M}$ signals. Multi-component Lorentzian fits are applied to identify dominant oscillation modes, their coherence (quality factor $Q$), and harmonic ratios.
  • Mass Scaling: To compare numerical results (computed for a reference mass $M = 10M_\odot$) with observational data, an inverse mass scaling relation is applied: $f_{num}(M_{obs}) = f_{num}(10M_\odot) \times (10M_\odot / M_{obs})$.

Key Contributions and Results

1. Hydrodynamical Signatures:

   • The Kerr–Sen deformation significantly alters the shock-cone morphology. In rapidly rotating cases ($a=0.9M$), the strong frame-dragging dominates, making Kerr–Sen profiles relatively similar to Kerr, though with measurable shifts in density peaks and shock locations.
   • In moderately rotating cases ($a=0.5M$), the deformation effects are more pronounced. The Kerr–Sen models exhibit a significant decrease in the density of matter trapped within the shock cone and a broader angular spreading compared to Kerr solutions.
   • The mass-accretion rate variability is sensitive to the deformation. Stronger deformation parameters lead to more irregular oscillations and deeper spikes in the accretion rate, particularly in the $a=0.5M$ regime.

2. QPO-Like Oscillation Modes:

   • PSD and Lorentzian analyses reveal that Kerr–Sen deformations shift characteristic frequencies and alter the coherence of oscillation modes.
   • The simulations produce near-resonant harmonic structures, specifically close to 3:2 and 2:1 ratios, which are frequently observed in astrophysical QPOs.
   • The specific harmonic structure depends on the deformation strength and spin. For instance, the KS2 and KS4 models (rapid spin) and KS5 model (moderate spin) exhibit clear 3:2 and 2:1 harmonic relationships, whereas stronger deformation in the KS3 model leads to more irregular, less harmonic distributions.

3. Observational Comparisons via Mass Scaling:
   The numerically computed frequencies, when scaled by inverse mass, show reasonable agreement with observed QPOs from stellar-mass, intermediate-mass, and supermassive black hole systems:

   • Stellar-Mass: The KS1 model ($a=0.9M$, soft deformation) successfully reproduces the 41 Hz and 67 Hz QPO pair of GRS 1915+105 and the ~66 Hz QPO of IGR J17091–3624. The KS5 model ($a=0.5M$) offers a frequency match for IGR J17091–3624 that is consistent with its estimated moderate spin.
   • Intermediate-Mass: The KS1 and KS5 models both reproduce the twin QPOs (3.32 Hz and 5.07 Hz) of M82 X–1. For NGC 5408 X–1 (0.010–0.020 Hz), the KS1 model provides a natural explanation if the source is rapidly rotating, suggesting a potential prediction for its spin parameter.
   • Supermassive: The KS1 and KS5 models align with the QPO of RE J1034+396. The KS1 and KS4 models (rapid spin) explain the QPO of 1H 0707–495, supporting high-spin interpretations. The KS5 model (moderate spin) provides the best agreement for ESO 113–G010, suggesting a moderately rotating black hole.

Significance and Claims
The paper claims that Kerr–Sen shock-cone oscillations provide a unified framework for interpreting timing features across a broad spectrum of black hole masses, from stellar to supermassive scales. The study demonstrates that:

• Hydrodynamical and timing diagnostics can empirically distinguish deviations associated with Kerr–Sen geometry from standard Kerr spacetime.
• The charge-related deformation parameter modifies not only the accretion morphology but also the fundamental oscillation modes, coherence, and harmonic ratios of the flow.
• These numerical results can be used to constrain the mass and spin parameters of sources where these properties are not fully established observationally. For example, the agreement of specific models with observed QPOs allows for theoretical predictions regarding the spin states of sources like NGC 5408 X–1 and ESO 113–G010.

The authors conclude that combined hydrodynamical and timing diagnostics constitute a promising approach for testing modified gravity theories against observational data, though they note that future work should extend these analyses to higher resolutions and fully three-dimensional flows to further test the robustness of these findings.