Disc-oscillation resonance and neutron star QPOs: 3:2 epicyclic orbital model

The high-frequency quasi-periodic oscillations (HF QPOs) that appear in the X-ray fluxes of low-mass X-ray binaries remain an unexplained phenomenon. Among other ideas, it has been suggested that a non-linear resonance between two oscillation modes in an accretion disc orbiting either a black hole or a neutron star plays a role in exciting the observed modulation. Several possible resonances have been discussed. A particular model assumes resonances in which the disc-oscillation modes have the eigenfrequencies equal to the radial and vertical epicyclic frequencies of geodesic orbital motion. This model has been discussed for black hole microquasar sources as well as for a group of neutron star sources. Assuming several neutron (strange) star equations of state and Hartle-Thorne geometry of rotating stars, we briefly compare the frequencies expected from the model to those observed. Our comparison implies that the inferred neutron star radius “RNS” is larger than the related radius of the marginally stable circular orbit “rms” for nuclear matter equations of state and spin frequencies up to 800Hz. For the same range of spin and a strange star (MIT) equation of state, the inferrred radius RNS is roughly equal to rms. The Paczynski modulation mechanism considered within the model requires that RNS < rms. However, we find this condition to be fulfilled only for the strange matter equation of state, masses below one solar mass, and spin frequencies above 800Hz. This result most likely falsifies the postulation of the neutron star 3:2 resonant eigenfrequencies being equal to the frequencies of geodesic radial and vertical epicyclic modes. We suggest that the 3:2 epicyclic modes could stay among the possible choices only if a fairly non-geodesic accretion flow is assumed, or if a different modulation mechanism operates.

💡 Research Summary

The paper tackles one of the most persistent puzzles in high‑energy astrophysics: the origin of the high‑frequency quasi‑periodic oscillations (HF QPOs) observed in the X‑ray flux of low‑mass X‑ray binaries (LMXBs). While HF QPOs in black‑hole systems tend to appear at nearly constant frequencies that often form a 3:2 ratio, neutron‑star (NS) sources display twin QPO peaks whose frequencies drift over time but still cluster around a 3:2 ratio in many cases. A widely discussed explanation is a non‑linear resonance between two disc‑oscillation modes whose eigenfrequencies are identified with the radial (ν_r) and vertical (ν_θ) epicyclic frequencies of geodesic motion in the strong‑gravity field of the compact object. This “3:2 epicyclic resonance” model predicts that at a specific resonance radius r₍₃:₂₎ the ratio ν_θ/ν_r = 3/2, and that the observed lower and upper QPO frequencies (ν_L, ν_U) correspond directly to ν_r and ν_θ (or to small corrections thereof).

The authors examine whether this model can simultaneously reproduce the observed QPO frequencies and satisfy the Paczyński modulation mechanism, which requires the neutron‑star radius R_NS to lie inside the marginally stable circular orbit radius r_ms (i.e. R_NS < r_ms). The Paczyński mechanism envisions an accretion disc that oscillates across an equipotential surface; when the disc surface crosses a critical potential, a rapid increase in mass inflow onto a hot spot on the stellar surface produces the observed X‑ray modulation.

To test the model, the paper proceeds in several steps:

1. Resonance Scenarios – Two possibilities are considered: (a) the observed frequencies are essentially the eigenfrequencies, with variations caused by a small shift Δr of the resonance radius; (b) the eigenfrequencies are fixed, and the observed drift is due to non‑linear corrections Δν_L, Δν_U. Prior work shows that scenario (a) does not reproduce the observed ν_U–ν_L correlation, so the authors focus on (b).

2. Mass Estimate in the Schwarzschild Limit – Ignoring rotation, the authors set ν_L≈600 Hz and ν_U≈900 Hz (the typical 3:2 pair) and solve the Schwarzschild epicyclic formulas. This yields a neutron‑star mass of roughly 1 M_⊙, a value first noted by Bursa (2004).

3. Equation‑of‑State (EoS) Modelling – Using the Hartle–Thorne formalism for rotating stars, they compute mass‑radius curves for a broad suite of nuclear‑matter EoS (various Skyrme parametrizations, DBHF, APR, FPS, BBB2, GLENDNH3) and for a strange‑matter MIT bag model (B=10¹⁴ g cm⁻³, α_c=0.15). For M≈1 M_⊙, all nuclear‑matter EoS give radii larger than the marginally stable orbit radius r_ms, violating the Paczyński condition. The MIT bag model yields R_NS≈r_ms, but only for masses below 1 M_⊙ and spin frequencies exceeding ~800 Hz.

4. Implications for the 3:2 Model – Because the Paczyński condition is essential for the modulation mechanism assumed in the model, the authors conclude that the pure geodesic 3:2 epicyclic resonance is incompatible with most realistic neutron‑star configurations. Only a narrow parameter space (low‑mass strange stars, very rapid rotation) satisfies the condition, making the model highly implausible for the majority of observed sources.

5. Possible Resolutions – The paper suggests two ways the model might be salvaged: (i) the disc flow could be significantly non‑geodesic (e.g., pressure, magnetic fields, radiation pressure) which can shift epicyclic frequencies by up to ~15 %, potentially allowing the resonance to occur at radii where R_NS < r_ms; (ii) an alternative modulation mechanism (e.g., relativistic lensing, Doppler boosting, magnetic hot‑spot oscillations) could replace the Paczyński mechanism, removing the stringent radius constraint.

Overall, the study provides a rigorous, quantitative test of the 3:2 epicyclic resonance hypothesis for neutron‑star HF QPOs. By combining relativistic orbital dynamics, realistic equations of state, and the Paczyński modulation requirement, it demonstrates that the simplest geodesic resonance model cannot simultaneously satisfy observational frequency ratios and the necessary geometric conditions for most neutron‑star systems. The work therefore pushes the community toward more sophisticated disc‑flow models or alternative QPO generation mechanisms, and highlights the importance of precise mass, radius, and spin measurements in constraining strong‑gravity astrophysics.