Internal resonance in non-linear disk oscillations and the amplitude evolution of neutron star kilohertz QPOs

We explore some properties of twin kilohertz quasiperiodic oscillations (QPOs) in a simple toy-model consisting of two oscillation modes coupled by a general nonlinear force. We examine resonant effects by slowly varying the values of the tunable, and nearly commensurable, eigenfrequencies. The behavior of the actual oscillation frequencies and amplitudes during a slow transition through the 3:2 resonance is examined in detail and it is shown that both are significantly affected by the nonlinearities in the governing equations. In particular, the amplitudes of oscillations reflect a resonant exchange of energy between the modes, as a result the initially weaker mode may become dominant after the transition. We note that a qualitatively similar behavior has been recently reported in several neutron star sources by Torok (2008, arXiv:0812.4751), who found that the difference of amplitudes in neutron star twin peak QPOs changes sign as the observed frequency ratio of the QPOs passes through the value 3:2.

💡 Research Summary

The paper presents a theoretical investigation of internal resonance between two non‑linear oscillation modes that are thought to represent the twin kilohertz quasi‑periodic oscillations (QPOs) observed in neutron‑star low‑mass X‑ray binaries. The authors construct a minimal “toy model” in which the two modes have eigenfrequencies ω₁ and ω₂ that are nearly in a 3:2 ratio. A general non‑linear coupling force containing cubic and quartic terms links the modes, allowing energy exchange. By slowly varying the eigenfrequencies so that the system drifts through the exact 3:2 commensurability, they apply a multiple‑scale analysis to derive amplitude and phase evolution equations for the two modes. Numerical integration of these equations reveals two key effects. First, the observable frequencies (the non‑linear corrected frequencies) deviate sharply as the system passes through resonance; instead of a strict frequency locking predicted by linear theory, the frequencies change continuously but with a steep slope near the resonance. Second, a resonant amplitude exchange occurs: the initially dominant mode loses energy to the weaker one, so that after the transition the second mode becomes the stronger. The magnitude of this exchange depends sensitively on the sign and size of the non‑linear coupling coefficients and on the damping rates.

The authors compare their results with the observational findings reported by Torok (2008), who noted that the difference in amplitudes of the twin QPOs reverses sign when the observed frequency ratio crosses 3:2 in several neutron‑star sources. The model reproduces this sign change naturally, suggesting that the observed QPOs may indeed be manifestations of a non‑linear internal resonance in the accretion disc. Moreover, the study shows that the width of the resonant capture region (the range of frequency ratios over which the exchange occurs) is controlled by the non‑linear coefficients and damping, providing a possible explanation for the variability of QPO properties across different sources and epochs.

Overall, the work demonstrates that a simple non‑linear dynamical system can capture essential features of neutron‑star kilohertz QPOs, offering a bridge between phenomenological observations and underlying disc physics. It highlights the importance of non‑linear mode coupling, suggests directions for more realistic magnetohydrodynamic simulations (including external driving, magnetic stresses, and radiation pressure), and motivates further observational tests involving other rational frequency ratios such as 4:3 or 5:4.