Quasi-Periodic Oscillations

Point a fast X-ray timer at an accreting neutron star or black hole and the light never sits still: it flickers on every timescale from seconds to a millisecond. Fourier-transform that flicker and, riding on a broad red-noise hump, one or more sharp-ish peaks emerge — quasi-periodic oscillations, QPOs. They are the closest thing astronomy has to a stethoscope pressed against the innermost few gravitational radii of a compact object — a region no telescope can resolve for a stellar-mass source. The peaks are quasi-periodic because the signal is not a clean sinusoid but a stochastically driven, damped oscillator — turbulence re-excites a natural frequency whose phase slowly wanders — so the power piles into a Lorentzian band, not a delta-function spike. What makes QPOs precious is which timescale they sample. The orbital frequency at radius $r$ around mass $M$ is $\nu_\phi=(1/2\pi)\sqrt{GM/r^3}$ , climbing to roughly a kilohertz just outside a neutron star and a few hundred hertz near a stellar black hole's innermost stable orbit. Observed QPO frequencies land in these ranges — direct evidence that we are watching gas whirl a handful of Schwarzschild radii from the horizon. There, general relativity splits the single orbital frequency into three fundamentals — azimuthal, radial, and vertical — whose small differences appear as slow precessions. Read those against a strong-field model and you can, in principle, weigh the black hole, clock its spin, and constrain the equation of state of matter denser than a nucleus. This article builds the dynamics from the Kerr metric, defines the observables rigorously, and works the headline predictions to real numbers.

Key takeaways

Quasi-periodic oscillations (QPOs) are broad, statistically significant peaks in the Fourier power spectrum of the X-ray flux from accreting neutron stars, black holes, and some active galactic nuclei.
Their frequencies (microhertz to ~1300 Hz) match dynamical timescales at a few to tens of gravitational radii, making QPOs a direct probe of matter in strong-field gravity.
They are classified by frequency: low-frequency QPOs (~0.1–30 Hz, usually type-C) and high-frequency QPOs (~40–450 Hz in black holes; twin kilohertz QPOs up to ~1300 Hz in neutron stars).
The quality factor $Q=\nu_0/\Delta\nu$ (centroid over FWHM) separates a coherent QPO ( $Q\gtrsim2$ ) from noise; the sharpest kHz QPOs reach $Q\sim100$ .
Leading models tie the peaks to relativistic epicyclic frequencies near the ISCO — relativistic precession, epicyclic 3:2 resonance, and Lense–Thirring precession of an inner torus.
Because every orbital frequency scales as $c^3/GM\propto1/M$ , HFQPOs encode the mass and spin — a route to testing general relativity and constraining the equation of state.
Key idea: QPO frequencies are set by the fundamental orbital and epicyclic frequencies of test particles in strong-field gravity near a compact object; measuring them constrains mass, spin, and — for neutron stars — the equation of state, because general relativity splits the single Newtonian orbital frequency into three that differ near the ISCO.

The core insight

QPO frequencies are set by the fundamental orbital and epicyclic frequencies of test particles in strong-field gravity near a compact object; measuring them constrains the mass, spin, and — for neutron stars — the equation of state, because in general relativity the single Newtonian orbital frequency splits into three that differ near the ISCO.

Where this lives today

QPOs sit at a live frontier where accretion physics, general relativity, and dense-matter physics intersect. The Rossi X-ray Timing Explorer (1995–2012) built the entire kHz-QPO and type-C catalogue; its successor for soft X-ray timing, NASA's NICER on the ISS (2017–), pairs microsecond timing with the neutron-star mass–radius programme that constrains the equation of state and feeds the Tolman–Oppenheimer–Volkoff curve. The relativistic-precession weighing of GRO J1655−40 to $M\approx5.3\,M_\odot$ , $a\approx0.29$ showcases QPO-based mass/spin measurement, and type-C QPOs now trace the inner-disc truncation radius in real time as black-hole transients cycle between hard and soft states. The open problem is stark: after four decades no single model — relativistic precession, epicyclic 3:2 resonance, Lense–Thirring precession of a tilted torus — is universally accepted, and different sources may favour different mechanisms. The proposed large-area timing missions eXTP and STROBE-X are designed to phase-resolve individual QPO cycles and measure the frequency tracks precisely enough to discriminate among the models — while candidate QPOs in active galactic nuclei test whether the same $1/M$ dynamics scale up over six orders of magnitude in mass.

The physical picture

Follow a single blob of gas orbiting deep in the accretion disc. In Newtonian gravity its orbit is a closed ellipse: it returns to the same radius, height, and azimuth once per revolution — one frequency. Strong-field gravity breaks that threefold degeneracy. Relativity softens the radial "spring", so the blob's radial oscillation cycles slightly slower than its revolution; the leftover $\nu_\phi-\nu_r$ is a slow advance of the orbit's periastron — the effect that precesses Mercury, here enormously amplified. Frame dragging softens the vertical spring too, so $\nu_\phi-\nu_\theta$ is a slow precession of the whole orbital plane — the nodal, or Lense–Thirring, precession that only a spinning spacetime produces. QPOs are the observable heartbeat of these motions. Why quasi-periodic? Because the disc is turbulent, and turbulence is a random driver. Picture pushing a child on a swing at random moments: it still oscillates near its natural frequency, yet its amplitude and phase drift. Fourier power collects in a band whose width is set by the damping — a Lorentzian bump, not a spike. The quality factor $Q=\nu_0/\Delta\nu$ counts how many cycles stay coherent; $Q$ of a few means "a QPO", $Q\to\infty$ a genuine clock. The first figure fuses the precessing orbit and the power spectrum it produces.

Figure

A blob orbiting near a compact object executes three distinct motions — azimuthal ν_φ, radial ν_r, and vertical ν_θ — whose small differences appear as slow orbital precession; the turbulent drive smears each into a Lorentzian power-spectral peak of width Δν, so the quality factor Q = ν₀ divided by Δν measures coherence

The frequencies scale inversely with mass: everything on the right of the Kerr formulas carries a factor $c^3/GM\propto1/M$ , so doubling the mass halves every orbital frequency at fixed radius in gravitational radii. That scaling explains the observed hierarchy — supermassive black holes in AGN at microhertz to millihertz, stellar black holes ( $5$ – $20\,M_\odot$ ) at tens to hundreds of hertz, neutron stars ( $1.4$ – $2\,M_\odot$ ) up to $\sim1300$ Hz. A measured HFQPO, fed into a spin-dependent model, therefore weighs the object: QPOs are a strong-gravity spectrometer. The second figure plots this frequency–radius diagram — the backbone of the relativistic precession model — for a real neutron-star mass.

Figure

Relativistic precession model for a 2 M⊙ non-spinning neutron star: the upper kHz QPO is identified with the orbital frequency ν_φ and the lower with the periastron frequency ν_per = ν_φ − ν_r, both plotted versus orbital radius in gravitational radii — the two peaks merge at the ISCO where r = 6 r_g, ν_r → 0, and ν_φ reaches its ceiling of about 1100 Hz

The core equations

QPOs live in the Fourier domain, so we begin with the power spectral density (PSD) of the light curve $x(t)$ : the modulus-squared of its discrete Fourier transform, normalized so integrating it returns the fractional variance — letting a QPO's strength be quoted as a percent rms.

P(\nu)=\frac{2}{\langle x\rangle^2\,N^2}\,\lvert X(\nu)\rvert^2 \qquad (1)

where: $P(\nu)$ — power spectral density; $X(\nu)$ — discrete Fourier transform of the light curve; $\langle x\rangle$ — mean count rate; $N$ — number of time bins. Reading: with this (rms) normalization $\mathrm{rms}^2=\int P(\nu)\,\mathrm{d}\nu$ , and a flat white-noise floor (level 2 in Leahy units) is subtracted before any feature is fitted. A QPO is fitted as a Lorentzian — the power spectrum of a driven damped oscillator — with centroid $\nu_0$ and width $\Delta$ whose ratio is the coherence.

P(\nu)=\frac{r^2\,\Delta/(2\pi)}{(\nu-\nu_0)^2+(\Delta/2)^2},\qquad Q\equiv\frac{\nu_0}{\Delta} \qquad (2)

where: $\nu_0$ — centroid frequency; $\Delta$ — full width at half maximum; $r$ — integrated fractional rms amplitude; $Q$ — quality factor. Reading: $Q\gtrsim2$ marks a QPO rather than noise; kHz QPOs reach $Q\sim100$ , but never infinity, which would signal a coherent pulsation. For a circular equatorial geodesic in the Kerr metric, the azimuthal (orbital) frequency has a compact form.

\nu_\phi=\frac{1}{2\pi}\frac{c^3}{GM}\frac{1}{\left(\dfrac{r}{r_g}\right)^{3/2}+a} \qquad (3)

where: $\nu_\phi$ — azimuthal orbital frequency; $r_g=GM/c^2$ — gravitational radius; $a=Jc/GM^2\in[-1,1]$ — dimensionless spin; $M,J$ — mass and angular momentum. Reading: $a\to0$ gives the Schwarzschild value $\nu_\phi=(c^3/2\pi GM)(r/r_g)^{-3/2}$ , exactly Kepler's $\tfrac{1}{2\pi}\sqrt{GM/r^3}$ , with the $1/M$ scaling manifest. The radial epicyclic frequency — small in-and-out oscillations about the circular orbit — is suppressed below $\nu_\phi$ and vanishes at the ISCO.

\nu_r^2=\nu_\phi^2\left(1-6\,\frac{r_g}{r}+8a\,\frac{r_g^{3/2}}{r^{3/2}}-3a^2\,\frac{r_g^2}{r^2}\right) \qquad (4)

where: $\nu_r$ — radial epicyclic frequency; $\nu_\phi$ — orbital frequency at that radius. Reading: the $-6\,r_g/r$ term drives $\nu_r^2\to0$ at $r=6r_g$ — the marginally stable orbit that caps a source's highest coherent frequency. The vertical epicyclic frequency equals $\nu_\phi$ in Schwarzschild but is dragged below it once the object spins.

\nu_\theta^2=\nu_\phi^2\left(1-4a\,\frac{r_g^{3/2}}{r^{3/2}}+3a^2\,\frac{r_g^2}{r^2}\right) \qquad (5)

where: $\nu_\theta$ — vertical epicyclic frequency. Reading: at $a=0$ , $\nu_\theta=\nu_\phi$ and the plane does not precess; any $\nu_\theta<\nu_\phi$ is a pure frame-dragging signal isolating the spin. The two differences between the fundamentals are the precession frequencies most models predict.

\nu_{\rm per}=\nu_\phi-\nu_r,\qquad \nu_{\rm nod}=\nu_\phi-\nu_\theta=\nu_{\rm LT} \qquad (6)

where: $\nu_{\rm per}$ — periastron precession; $\nu_{\rm nod}=\nu_{\rm LT}$ — nodal / Lense–Thirring precession. Reading: both vanish in the Newtonian limit where the three fundamentals coincide. Two empirical relations anchor the phenomenology: in neutron stars the twin kHz QPOs separate by roughly the spin (a beat clue), and in black holes the high-frequency pairs cluster near a 3:2 ratio (a resonance clue).

\Delta\nu=\nu_{\rm upper}-\nu_{\rm lower}\approx\nu_{\rm spin},\qquad \frac{\nu_{\rm upper}}{\nu_{\rm lower}}\approx\frac{3}{2} \qquad (7)

where: $\nu_{\rm upper},\nu_{\rm lower}$ — upper and lower QPO centroids; $\nu_{\rm spin}$ — neutron-star spin frequency. Reading: " $\approx$ " does real work — the separation drifts and is generally not exactly the spin, and the ratio is near but not locked to 3:2, which any successful model must explain.

Derivation AThe innermost stable circular orbit and the maximum QPO frequency

The goal is to find the smallest radius at which matter can hold a stable circular orbit around a non-spinning compact object, and the orbital frequency there — because that radius, the innermost stable circular orbit (ISCO), sets a hard ceiling on how fast any orbital-motion QPO can be. We build the relativistic effective potential from the Schwarzschild geodesic, impose marginal stability, solve for the radius, and convert to a frequency we can test against observed high-frequency QPOs.

Step 1. For a timelike geodesic in Schwarzschild spacetime the conserved energy $\tilde E$ and angular momentum $\tilde L$ per unit mass reduce radial motion to $(\mathrm{d}r/\mathrm{d}\tau)^2=\tilde E^2-V_{\rm eff}(r)$ , with $V_{\rm eff}=\left(1-\tfrac{2GM}{c^2r}\right)\left(1+\tfrac{\tilde L^2}{c^2r^2}\right)$ . Why: a particle rolling in $V_{\rm eff}$ makes "stable orbit" a statement about the potential's shape.

Step 2. A circular orbit sits at $V_{\rm eff}'(r)=0$ ; it is stable if $V_{\rm eff}''>0$ . The innermost stable orbit is the marginal case, so we solve $V_{\rm eff}'=0$ and $V_{\rm eff}''=0$ together. Why: the second derivative's sign is the dividing line — inside it a small inward nudge no longer restores the orbit.

Step 3. The two conditions are met only at $r_{\rm ISCO}=6GM/c^2=6\,r_g$ , independent of $\tilde L$ — equivalently where $\nu_r^2\propto(1-6r_g/r)$ hits zero. Why: the ISCO and the zero of $\nu_r$ coincide because both express the same loss of radial stability.

Step 4. Insert $r=6r_g$ into $\nu_\phi$ to get the ceiling, $\nu_{\rm ISCO}=(c^3/2\pi GM)\,6^{-3/2}\approx2200\,(M_\odot/M)$ Hz — near $220$ Hz for a $10\,M_\odot$ black hole, $1100$ Hz for a $2\,M_\odot$ neutron star, exactly the HFQPO and kHz-QPO bands. Why: this number, with mass as the only input, explains why black-hole HFQPOs stall in the hundreds of hertz while lighter neutron stars reach a kilohertz.

Step 5. Prograde spin shrinks the ISCO toward $r_g$ as $a\to1$ , lifting the ceiling several-fold — so an unusually high HFQPO is evidence of rapid rotation. Why: the ISCO's spin dependence turns a frequency ceiling into a spin measurement. *Note: $r_{\rm ISCO}=6GM/c^2$ ; evaluating $c^3/(2\pi GM)\cdot6^{-3/2}$ gives $2198\,(M_\odot/M)$ Hz, reproducing $\sim2.2$ kHz per solar mass; and $\nu_r\to0$ at $6r_g$ is consistent with equation (4).*

Go deeper · Derivation B · 1 minThe relativistic precession model of the twin kilohertz QPOs

Derivation BThe relativistic precession model of the twin kilohertz QPOs

The goal is to show how a single characteristic radius in the inner disc can generate two simultaneous kHz QPOs and one low-frequency QPO, and to turn the observed frequency–frequency track into a measurement of mass and spin. In the relativistic precession model of Stella and Vietri each observed frequency is one of the three fundamentals at the same radius, so as that radius moves with accretion rate the frequencies sweep out a predicted, fittable curve.

Step 1. Assign the upper kHz QPO to the orbital frequency, $\nu_{\rm upper}=\nu_\phi(r)$ ; the lower to periastron precession, $\nu_{\rm lower}=\nu_\phi(r)-\nu_r(r)$ ; the type-C low-frequency QPO to nodal precession, $\nu_{\rm LF}=\nu_\phi(r)-\nu_\theta(r)$ . Why: three observables become three functions of one radius, so the model makes falsifiable predictions rather than fitting each peak alone.

Step 2. Eliminate the radius between the two kHz frequencies using equation (4). With $x=r_g/r$ , the pair obeys the closed relation below, so $(\nu_{\rm lower},\nu_{\rm upper})$ traces a one-parameter curve as $r$ varies.

\nu_{\rm lower}=\nu_{\rm upper}\left[1-\sqrt{1-6x+8a\,x^{3/2}-3a^2x^2}\,\right],\qquad x=\frac{r_g}{r}

Why: the track's scale fixes $M$ and its curvature fixes $a$ , so fitting it measures both.

Step 3. In the ISCO limit ( $r\to6r_g$ , $a=0$ ), $\nu_r\to0$ , so $\nu_{\rm lower}\to\nu_\phi=\nu_{\rm upper}$ : the two peaks converge and merge. This turnover, seen in several neutron-star sources, locates the ISCO and bounds the mass. Why: a merge of two independent frequencies is a clean prediction only the ISCO can produce.

Step 4. Numerically, for a $2\,M_\odot$ , $a=0$ neutron star $\nu_\phi$ runs from $\approx1100$ Hz at the ISCO to $\approx870$ Hz at $7r_g$ and $\approx510$ Hz at $10r_g$ , bracketing observed upper kHz QPOs; $\nu_{\rm per}$ falls from that ceiling toward zero at large radius. Why: the model lands on the observed kHz band with only the mass as input. *Note: $\nu_\phi=1099$ Hz at $6r_g$ for $2\,M_\odot$ ; at large $r$ , $\nu_r\to\nu_\phi$ so $\nu_{\rm per}\to0$ (the Newtonian limit of no precession); and $\nu_r\to0$ at the ISCO reproduces the observed twin-peak convergence.*

Go deeper · Derivation C · 1 minLense–Thirring precession as the low-frequency QPO

Derivation CLense–Thirring precession as the low-frequency QPO

The goal is to show that the nodal precession of equation (6), in the weak-field slow-rotation limit, collapses to the classical Lense–Thirring frame-dragging rate — closing the loop between type-C low-frequency QPOs and spin, and yielding a testable low-versus-high-frequency correlation. We expand the vertical epicyclic frequency to first order in spin, subtract it from the orbital frequency, and rewrite in terms of the angular momentum $J$ .

Step 1. Expand equation (5) to first order in $a$ : $\nu_\theta\approx\nu_\phi\left[1-2a\left(r_g/r\right)^{3/2}\right]$ . Why: type-C QPOs arise in modestly spinning systems far enough out that the leading spin term suffices.

Step 2. Subtract from $\nu_\phi$ : $\nu_{\rm nod}=\nu_\phi-\nu_\theta\approx2a\,\nu_\phi\left(r_g/r\right)^{3/2}$ , then substitute $\nu_\phi$ and $a=Jc/GM^2$ to reach the boxed frame-dragging form.

\nu_{\rm LT}=\nu_\phi-\nu_\theta\;\xrightarrow{\ \text{weak field}\ }\;\frac{2a}{2\pi}\frac{c^3}{GM}\left(\frac{r_g}{r}\right)^{3}=\frac{GJ}{\pi c^2 r^3},\qquad \Omega_{\rm LT}=2\pi\nu_{\rm LT}=\frac{2GJ}{c^2 r^3}

Why: recovering the textbook rate $\Omega_{\rm LT}=2GJ/c^2r^3$ shows the GR nodal frequency and classical frame-dragging are one physics.

Step 3. The payoff: $\nu_{\rm LT}\propto r^{-3}$ while $\nu_\phi\propto r^{-3/2}$ , so eliminating the radius gives $\nu_{\rm LF}\propto\nu_{\rm HF}^2$ — a quadratic correlation seen in real sources. Why: a parameter-free scaling law between two measured frequencies can confirm or kill a model. *Note: $GJ/(\pi c^2 r^3)$ equals $\Omega_{\rm LT}/2\pi$ — the standard 1918 frame-dragging rate — and evaluates to $6.46$ Hz at $r=10r_g$ , $a=0.2$ , $2\,M_\odot$ , confirming it matches the exact weak-field nodal frequency.*

Limiting cases

The three fundamentals encode different physics in different regimes; their limits are the fastest way to see what each QPO measures.

Case	Result	Meaning
Newtonian / large radius	$\nu_\phi=\nu_r=\nu_\theta=\dfrac{1}{2\pi}\sqrt{\dfrac{GM}{r^3}}$	The three coincide (closed Keplerian ellipse): no precession, no QPO splitting.
Approaching the ISCO (Schwarzschild)	$\nu_r\to0$ as $r\to r_{\rm ISCO}=6\,GM/c^2$	The radial epicyclic frequency vanishes as the orbit turns marginally stable, while $\nu_\phi=\nu_\theta$ stay finite; caps HFQPO frequencies.
Kerr, spin $a$	$\nu_\theta<\nu_\phi$ (prograde)	Spin breaks equatorial symmetry, so $\nu_\theta$ drops below $\nu_\phi$ ; the difference is the Lense–Thirring precession that carries $a$ .
Neutron star vs black hole	$r_{\rm ISCO}\gtrless R_{\rm NS}$	If the ISCO lies outside the surface it sets the top orbital frequency; if the surface is larger, the surface orbit does — an equation-of-state test.
Rapid prograde spin, $a\to1$	$r_{\rm ISCO}\to r_g,\ \nu_{\rm ISCO}$ several × higher	Frame dragging pulls the last stable orbit inward, raising the ceiling — an unusually fast QPO implies near-maximal spin.

Worked examples

1 · The frequency ceiling of a neutron star (fully worked)

A $2.0\,M_\odot$ neutron star of radius $R=11$ km accretes from a companion. (a) Where is its ISCO, and does it lie outside the star? (b) What is the maximum orbital (upper kHz QPO) frequency? (c) An observer reports an upper kHz QPO at $1220$ Hz — is that consistent with orbital motion outside the ISCO?

Step 1. Gravitational radius: $r_g=GM/c^2=1.477\ \mathrm{km}\times(M/M_\odot)=\mathbf{2.95\ km}$ , so the ISCO is at $r_{\rm ISCO}=6r_g=\mathbf{17.7\ km}$ — comfortably outside the $11$ km surface. The star behaves like a black hole for orbital dynamics: the ISCO sets the innermost orbit.

Step 2. Maximum orbital frequency at the ISCO: $\nu_{\rm ISCO}=2198\,(M_\odot/M)\ \mathrm{Hz}=\mathbf{1099\ Hz}$ . No stable circular orbit, and no orbital-motion QPO tied to one, can exceed this.

Step 3. The reported $1220$ Hz exceeds the $1099$ Hz ceiling for a non-spinning $2.0\,M_\odot$ star. Something must give: a lower mass (a $1.8\,M_\odot$ star lifts the ceiling to $1221$ Hz), prograde spin, or a QPO that is not the ISCO orbital frequency. At $2.0\,M_\odot$ , $a=0$ it is impossible, so the data push toward lower mass or finite spin — exactly how kHz QPOs constrain the equation of state.

2 · Weighing a black hole from a 3:2 HFQPO pair

The microquasar GRO J1655−40 shows HFQPOs at $300$ Hz and $450$ Hz — a clean 3:2 pair. Use the ISCO ceiling to bound its mass, and compare with the full relativistic-precession fit.

Step 1. The upper HFQPO, $450$ Hz, cannot exceed the ISCO orbital frequency. For $a=0$ , $\nu_{\rm ISCO}=2198\,(M_\odot/M)$ Hz, so $450\le2198/M$ gives $M\le\mathbf{4.9\,M_\odot}$ — a bound with no modelling.

Step 2. Prograde spin raises the ceiling, so a spinning hole with the same frequencies can be more massive. Fitting all three frequencies — the two HFQPOs plus the $\sim17$ Hz type-C LFQPO — in the RPM gives three equations in $(r,M,a)$ .

Step 3. The published solution (Motta and collaborators) yields $M\approx5.3\,M_\odot$ , $a\approx0.29$ , consistent with the dynamical mass from optical spectroscopy. The QPOs alone weigh the black hole to $\sim5\,M_\odot$ and clock a moderate spin, independent of continuum-fitting or iron-line methods — and Step 1's $4.9\,M_\odot$ bound shows why it had to be a few solar masses.

3 · A beat-frequency spin estimate

A neutron-star binary shows twin kHz QPOs at $920$ Hz and $615$ Hz. What spin does a naive beat-frequency reading imply, and what is the catch?

Step 1. The separation is $\Delta\nu=920-615=\mathbf{305\ Hz}$ , so the beat guess sets $\nu_{\rm spin}\approx305$ Hz — a spin period near $3.3$ ms, in the recycled-millisecond-pulsar range.

Step 2. The catch: measured separations track half the spin in some fast rotators and drift by tens of hertz as the frequencies climb. A $305$ Hz separation strongly hints at a $\sim300$ Hz spin, but confirmation needs an independent measurement — a burst oscillation or pulsation — because the separation is not a rigid clock.

Physics of the portrait

QPOs are an observational science, and no one shaped it more than the pioneer of X-ray timing — who first recognized a QPO for what it was and later found the kilohertz oscillations that reached the strong-field regime.

Origins · 1 minMichiel van der Klis and how this was discovered

Michiel van der Klis

b. 1953, Dutch astronomer · X-ray timing and the discovery of QPOs

Michiel van der Klis turned the fast X-ray variability of accreting compact objects from noise into a quantitative probe of strong gravity. Working with EXOSAT data, in 1985 he and collaborators identified the first quasi-periodic oscillation — a broad, drifting peak in the power spectrum of the neutron-star binary GX 5−1 — and interpreted it as a beat between a Keplerian frequency and the star's spin rather than dismissing it as instrumental. That reframing created a subfield. The decisive moment came in 1996, when the Rossi X-ray Timing Explorer, with its large area and microsecond timing, revealed kilohertz QPOs in neutron-star binaries at frequencies matching orbital motion a few gravitational radii out — the closest look yet at matter near the ISCO. His 2006 review in Compact Stellar X-ray Sources codified the taxonomy and the statistical machinery — Leahy and rms normalizations, Lorentzian fitting, the quality factor $Q$ — used throughout this article. His conviction that timing resolves what imaging cannot remains the field's animating idea: the fastest clocks in the sky are also the deepest probes of curved spacetime.

History

1918

— Lense and Thirring predict that a rotating mass drags spacetime, precessing nearby orbital planes

1970

— Uhuru, the first dedicated X-ray satellite, launches and opens systematic X-ray binary timing

1985

— van der Klis and colleagues discover the first QPO, in the neutron-star binary GX 5−1, with EXOSAT

1996

— the Rossi X-ray Timing Explorer discovers kilohertz QPOs in neutron-star binaries

1998

— Stella and Vietri propose the relativistic precession model of the twin kHz QPOs

2001

— Abramowicz and Kluzniak propose the epicyclic 3:2 resonance model for black-hole HFQPOs

Going further

Texts. Carroll and Ostlie, An Introduction to Modern Astrophysics, 2nd ed. (accretion discs, compact-object binaries) · Dodelson, Modern Cosmology (relativistic dynamics background) · Mo, van den Bosch and White, Galaxy Formation and Evolution (accretion and central engines) · Binney and Tremaine, Galactic Dynamics, 2nd ed. (epicyclic frequencies and orbits in axisymmetric potentials) · van der Klis, in Compact Stellar X-ray Sources (2006), Ch. 2 (the definitive QPO review) · Stella and Vietri, ApJ Letters 492, L59 (1998) (relativistic precession model) · Abramowicz and Kluzniak, A&A 374, L19 (2001) (epicyclic resonance model). Key numerical benchmarks: $r_{\rm ISCO}=6GM/c^2$ and $\nu_{\rm ISCO}=(c^3/2\pi GM)6^{-3/2}=2198\,(M_\odot/M)$ Hz follow from the Schwarzschild effective potential; the Kerr epicyclic formulas reduce correctly in the $a\to0$ and large- $r$ limits; the weak-field Lense–Thirring rate is $\nu_{\rm LT}=GJ/(\pi c^2r^3)$ , equal to $\Omega_{\rm LT}/2\pi$ with $\Omega_{\rm LT}=2GJ/c^2r^3$ ; the twin-QPO orbital frequencies run $\nu_\phi=1099$ Hz at $6r_g$ and $872$ Hz at $7r_g$ for $2\,M_\odot$ ; and the ISCO ceiling gives $M\le4.9\,M_\odot$ from a $450$ Hz upper HFQPO, consistent with the full RPM fit $M\approx5.3\,M_\odot$ , $a\approx0.29$ for GRO J1655−40. Frequency scaling $\propto1/M$ and the near-3:2 ratio are standard observational facts.

Prerequisites. Neutron Stars · Neutron Star Structure and Interior · Equation of State of Dense Matter · Millisecond and Recycled Pulsars. Leads to. Tolman–Oppenheimer–Volkoff Equation · Magnetars.

Ask the tutor: Ask me to invert a set of observed QPO frequencies into a black-hole mass and spin, or to show why the radial epicyclic frequency must vanish exactly at the innermost stable circular orbit.

What is the difference between a QPO and a coherent pulsation?

A coherent pulsation — from a spinning, magnetized neutron star — is a sharp spectral line: its phase is stable over billions of cycles, so you can build a timing model and count every rotation. A QPO is a broad peak: the oscillation is re-excited by turbulence and its phase wanders, so it stays coherent only for $Q=\nu_0/\Delta\nu$ cycles. So you can fold data on a pulsation to sharpen it indefinitely, but a QPO cannot be folded coherently — there is no fixed phase to fold on.

Why do neutron-star QPOs reach higher frequencies than black-hole QPOs?

Because orbital frequency scales as $c^3/GM\propto1/M$ , with the maximum near the ISCO at $\nu_{\rm ISCO}\approx2200\,(M_\odot/M)$ Hz. A $\sim10\,M_\odot$ black hole therefore caps near $220$ Hz (higher if spinning), whereas a $1.4$ – $2\,M_\odot$ neutron star is five to seven times lighter and reaches over a kilohertz. It is a pure mass effect: the same strong-field dynamics at a smaller mass runs faster, which is also why supermassive black holes should show analog QPOs at microhertz to millihertz.

How can three QPO frequencies measure both mass and spin?

The relativistic precession model identifies the three simultaneous frequencies — orbital $\nu_\phi$ , periastron precession $\nu_\phi-\nu_r$ , nodal precession $\nu_\phi-\nu_\theta$ — as the three fundamentals at one radius, giving three equations in three unknowns $(r,M,a)$ . The mass sets the frequency scale ( $c^3/GM$ ), the spin distorts the relative spacing (the $a$ -terms in $\nu_r$ and $\nu_\theta$ ), so the two separate when all three peaks are seen together — as for GRO J1655−40.

Why isn't the twin-kHz-QPO separation exactly the spin frequency?

The original beat-frequency idea predicted $\Delta\nu=\nu_{\rm upper}-\nu_{\rm lower}=\nu_{\rm spin}$ , imagining the lower QPO as a beat between a Keplerian frequency and the star's rotation. Real data disagree: the separation is close to the spin (or, in fast rotators, half it) but drifts by tens of hertz as the frequencies climb. That drift is fatal to a rigid beat model and favours the precession and resonance pictures, in which the separation is $\nu_r(r)$ or a resonance detuning that varies with radius — an approximate resonance, not an identity.

What does the 3:2 frequency ratio in black-hole HFQPOs tell us?

High-frequency QPO pairs in black-hole binaries cluster near a 3:2 ratio (for example $450$ and $300$ Hz in GRO J1655−40). A commensurate ratio is the fingerprint of a nonlinear resonance: two modes — plausibly the vertical and radial epicyclic frequencies — lock together where their frequencies are in a small-integer ratio, pumping energy into a coherent oscillation. In the Abramowicz–Kluzniak model this resonance selects a preferred radius and frequency pair, so the ratio recurs across sources even as the absolute frequencies differ with mass. It is not exactly 3:2 — the resonance has finite width — itself a clue to the driving and damping.

Are QPOs seen outside X-ray binaries?

Yes, and testing that is a major goal. If the same strong-gravity dynamics govern all accreting black holes, supermassive holes in active galactic nuclei should show QPOs lower by the mass ratio — microhertz to millihertz for $10^6$ – $10^9\,M_\odot$ . Candidate QPOs are reported in Seyferts such as RE J1034+396, and quasi-periodic eruptions may be related. Confirming them would extend the $1/M$ scaling across six orders of magnitude in mass — though the long timescales make the statistics hard, since one cycle can take hours to days.