Quasi-Periodic Oscillation

What a QPO is — and what it is not

Look at the X-ray light curve of a bright neutron-star or stellar-mass black-hole binary and you almost never see a clean sine wave. Fold the data into a power spectrum, however — a Fourier transform of the count rate squared, plotted against frequency — and structure jumps out. There is a broad continuum of red noise rolling off toward high frequencies, and superimposed on that continuum are one or several narrow but finite-width peaks. Those peaks are quasi-periodic oscillations: the X-rays prefer a particular timescale, but the phase wanders. A QPO is a clock, but a wobbly clock.

The quantitative definition is the quality factor

Q = ν_0 / Δν_FWHM

where ν_0 is the centroid frequency and Δν_FWHM is the full-width at half-maximum of the Lorentzian profile fitted to the peak. Observed QPOs span Q ≈ 2 (the conventional lower limit; below it the peak is called a "noise component") to Q ≈ 50. By contrast, a young rotation-powered radio pulsar maintains Q ≳ 10⁹ over decades; the millisecond X-ray pulsations of SAX J1808-3658 reach Q ~ 10⁶. The factor 10⁵ gap between QPOs and true pulsations is what people mean by "quasi-periodic" — the clock is real but it is reset by something on a timescale of tens to thousands of cycles.

QPOs live in the X-ray power spectrum specifically because the X-rays come from the innermost, hottest parts of the accretion flow, where the relevant orbital timescales are milliseconds for stellar-mass objects. Optical and UV modulations exist too — they are slower and noisier because they probe larger radii — but the diagnostic value of a QPO is that it is locked to a region where general relativity dominates.

Discovery: van der Klis and the EXOSAT era

In 1985 Michiel van der Klis and collaborators were analysing observations of the Z-source neutron-star binary Sco X-1 with EXOSAT. The power spectrum showed a peak that drifted in frequency around 6–20 Hz as the source moved through its hardness-intensity diagram. They published it under the careful title "Intensity-dependent quasi-periodic oscillations in the X-ray flux of GX 5-1" (Nature 316, 225, 1985); within a year the phenomenon had been confirmed in Sco X-1 itself and in several other LMXBs. Before that paper the existing framework distinguished only coherent pulsations (millisecond X-ray pulsars) from broadband noise; the quasi-periodic regime was new.

QPOs remained an exotic class until the Rossi X-ray Timing Explorer (RXTE) launched in December 1995 with a 6500 cm² Proportional Counter Array (PCA) optimised for microsecond timing of bright sources. RXTE's first big discovery was the kilohertz twin-peak QPOs in neutron-star LMXBs (Strohmayer et al. 1996 on 4U 1728-34, then a flood). Within a few years it had built the modern empirical taxonomy of black-hole LFQPOs and HFQPOs that remain in use today.

The taxonomy: types and frequency ranges

QPOs in stellar-mass black-hole X-ray binaries are sorted by frequency and by their relationship to the source's spectral-timing state. The classification scheme of Casella, Belloni & Stella (2005) is the standard reference.

Class	Frequency	Q	rms	Spectral state	Likely mechanism
Type A	6 – 8 Hz	~3	~2 %	Soft / intermediate	Disk + corona instability
Type B	5 – 6 Hz	~6	~4 %	Soft-intermediate, jet ejection	Compact-jet-related
Type C	0.1 – 10 Hz	7 – 12	3 – 16 %	Hard / hard-intermediate	Lense-Thirring precession
HFQPO (upper)	~150 – 450 Hz	~5	~1 %	Soft / intermediate	Resonance / p-mode
HFQPO (lower)	~100 – 300 Hz	~5	~1 %	Soft / intermediate	Resonance / radial epicyclic
NS kHz pair	~200 – 1300 Hz	2 – 200	1 – 20 %	Atoll, banana branch	Keplerian + sonic point
NS hectohertz	100 – 250 Hz	~1 – 10	~3 %	Various	Trapped disk mode

Two empirical patterns drive the entire field. First, neutron-star kHz QPOs come in pairs whose frequency separation Δν ≈ 250–350 Hz tracks (with some sources showing Δν ≈ ν_spin / 2). Second, black-hole HFQPOs, when they appear in pairs, lock into a 3:2 ratio: 300/450 Hz in GRO J1655-40, 41/67 Hz in GRS 1915+105, 165/240 Hz in XTE J1550-564, 92/184 Hz (later revised) in H1743-322. The HFQPO frequencies are stable: they do not slide with luminosity the way type-C QPOs do, which is one of the reasons people take them as a marker of an underlying geometric resonance rather than a fluid timescale that depends on Ṁ.

Where the frequencies come from

For a particle in a circular orbit around a Schwarzschild black hole at radius r, the Keplerian orbital frequency is

ν_K(r) = (1 / 2π) (GM / r³)^(1/2)
       ≈ 32 kHz × (M/M☉)⁻¹ × (r/r_g)⁻³ᐟ²

where r_g = GM/c² is the gravitational radius. At the Schwarzschild ISCO, r = 6 r_g, this becomes

ν_K(ISCO) ≈ 220 Hz × (10 M☉ / M)

For a Kerr black hole the ISCO moves inward with prograde spin, and the corresponding ν_K(ISCO) rises: at a* = 0.9 it is roughly 800 Hz × (10 M☉/M). The observed HFQPO frequencies in stellar-mass black holes are comfortably inside this window. For a 1.4 M☉ neutron star the same scaling gives ν_K(ISCO) ≈ 1500 Hz, neatly matching the kilohertz pair range.

For a non-circular orbit, two additional frequencies become relevant: the radial epicyclic frequency κ at which a slightly elliptical orbit oscillates in and out, and the vertical epicyclic frequency Ω_θ at which a slightly inclined orbit oscillates up and down through the equatorial plane. In a Newtonian potential κ = Ω_K and Ω_θ = Ω_K, so these are degenerate with the Keplerian frequency. In Kerr spacetime they split:

κ²    = ν_K² (1 − 6/x + 8 a* x^(-3/2) − 3 a*² /x² )
Ω_θ²  = ν_K² (1 − 4 a* x^(-3/2) + 3 a*²/x²)

where x = r/r_g. The radial epicyclic frequency vanishes at the ISCO and falls below ν_K everywhere inside; this is the relativistic origin of orbital instability inside r_ISCO. The vertical epicyclic frequency drops to zero at r = 0 only in extreme Kerr; everywhere relevant Ω_θ < Ω_K. These three frequencies — ν_K, κ, Ω_θ — and their differences ν_K − κ (the periastron-precession frequency) and ν_K − Ω_θ (the nodal-precession or Lense-Thirring frequency) provide the toolkit from which every QPO model is built.

Competing models — what produces the peak

No single model has yet displaced the others. The major contenders for the high-frequency QPOs are:

• Relativistic precession model (Stella & Vietri 1998). Three observed frequencies — the lower HFQPO, the upper HFQPO, and the type-C LFQPO — are identified with the radial epicyclic frequency κ, the Keplerian frequency ν_K, and the nodal (Lense-Thirring) precession ν_K − Ω_θ at one and the same radius. The radius is allowed to vary; the three frequencies are constrained to evolve together. The model has been tested by Motta et al. (2014) on simultaneous detections in GRO J1655-40 and reproduces mass-spin combinations consistent with optical measurements.
• Parametric resonance (Abramowicz & Kluźniak 2001). A non-linear coupling between κ and Ω_θ in Kerr spacetime resonantly amplifies oscillations at radii where their ratio is a small integer, with the 3:2 commensurability the leading mode. Predicts that HFQPO pairs should lock at 3:2, which is what is observed.
• Trapped p-mode oscillations (Wagoner 1999; Kato). Pressure waves on a thin disk can be trapped in cavities defined by the Kerr radial-epicyclic frequency, giving "diskoseismology" eigenmodes whose frequencies depend on (M, a). Predicts a discrete spectrum of QPO frequencies but has struggled with the rms amplitudes observed.
• Precession of a hot inner flow (Ingram, Done & Fragile 2009). The truncated thin-disk-plus-hot-inner-flow geometry of black-hole hard states has a misaligned inner Comptonising torus that precesses bodily about the spin axis at the Lense-Thirring frequency at its mean truncation radius. As the truncation radius shrinks during a hard-to-soft transition, the precession frequency rises — and the LFQPO frequency rises with it. This model has had clear success at reproducing type-C QPO phenomenology, including the inclination-dependent rms and the energy-dependent phase lags.
• Corrugation / magnetic instabilities. Magnetically-driven precession-tearing modes (Dexter & Blaes 2014; Tchekhovskoy) provide an alternative source of nodal-frequency variability that does not require a misaligned disk.

Currently the precessing-inner-flow picture is the strongest framework for type-C LFQPOs and the relativistic-precession + 3:2 resonance combination is the strongest framework for the HFQPOs, but a unified theory remains a research goal rather than a textbook result.

Worked example: deducing mass and spin from a 3:2 HFQPO pair

Take GRO J1655-40, which RXTE detected with a stable HFQPO pair at 300 Hz and 450 Hz. Its companion is an F-type subgiant whose radial-velocity curve gives a dynamical black-hole mass of M = 5.4 ± 0.3 M☉ (Beer & Podsiadlowski 2002).

In the parametric-resonance model the two peaks are identified with the radial epicyclic frequency κ and the vertical epicyclic frequency Ω_θ at the resonance radius where κ : Ω_θ = 2 : 3. For a Kerr black hole the location of that radius depends only on the dimensionless spin a*, and at the same radius Ω_θ is given by the closed-form Kerr expression above. Setting Ω_θ = 450 Hz × 2π and M = 5.4 M☉, one solves numerically and finds

a* ≈ 0.93   (parametric-resonance HFQPO fit)

The relativistic-precession model on the same data, when fitted simultaneously with the simultaneously-observed 17 Hz LFQPO, gives a* ≈ 0.97. Iron-line X-ray reflection spectroscopy of the same source gives a* = 0.94 ± 0.01 (García et al. 2018). The three independent techniques converge on a near-extremally-spinning black hole. The agreement is by no means automatic; in some systems the methods give different answers and the literature treats those as open puzzles.

GRS 1915+105, the other workhorse 3:2 system, gives 41/67 Hz at M = 12.4 M☉, implying a* ≳ 0.98 in the same framework. The convergence of multiple methods on near-maximal spin for GRS 1915+105 is one of the cleaner observational results in stellar-mass black hole spin measurement.

Famous QPO systems

Source	Type	QPO frequencies	Mass	Significance
Sco X-1	NS LMXB (Z)	~6–20 Hz LFQPO; 800/1100 Hz pair	1.4 M☉ NS	First QPO discovery, 1985; prototype kHz twin peaks
GRS 1915+105	BH HMXB-like	41/67 Hz HFQPO; complex LFQPOs	12.4 M☉	Canonical 3:2; quasi-maximal spin a* ≳ 0.98
GRO J1655-40	BH LMXB	300/450 Hz HFQPO; 17 Hz type-C	5.4 M☉	Best simultaneous HFQPO + LFQPO trio
XTE J1550-564	BH LMXB	165/240 Hz HFQPO; type A/B/C LFQPOs	9.1 M☉	Defined the LFQPO type-A/B/C taxonomy
H1743-322	BH LMXB	~165/240 Hz HFQPO	~8 M☉	Inclination-dependent type-C rms benchmark
4U 1608-52	NS LMXB (atoll)	~470/870 Hz kHz pair	1.4 M☉ NS	RXTE-era kHz QPO laboratory
4U 1636-53	NS LMXB (atoll)	~830/1230 Hz kHz pair	1.4 M☉ NS	Cleanest twin-peak coherence Q ≳ 200
RE J1034+396	AGN (NLS1)	3730 s QPO	~10⁶ M☉	First clean AGN QPO (Gierliński 2008)
1ES 1927+654	Changing-look AGN	~18 min QPO appearing 2022	~10⁶ M☉	QPO turned on alongside coronal collapse
ASASSN-14li	TDE	7.65 mHz QPO	~10⁶ M☉	QPO in a transient TDE disk

The instruments that did the work

• EXOSAT (1983–1986). Discovered QPOs (van der Klis 1985). The Medium Energy gas counter array gave the first power spectra capable of resolving Lorentzian peaks against shot-noise continuum.
• RXTE / PCA (1995–2012). The workhorse. 6500 cm² effective area at 2–10 keV, μs timing, near-continuous coverage of bright XRBs for 16 years. The PCA archive at HEASARC still drives the bulk of QPO physics papers.
• NICER (2017–). Soft X-ray (0.2–12 keV) timing payload on the ISS with similar effective area to PCA at 1–3 keV and superior energy resolution. Has refined LFQPO phase-lag spectroscopy for sources like MAXI J1820+070 and 4U 1543-47.
• NuSTAR (2012–). Hard X-ray (3–79 keV) focusing optics; resolves the HFQPO spectrum into the disk-reflection and Comptonised components, sharpening the geometric inference.
• IXPE (2021–) & XRISM (2023–). X-ray polarimetry and high-resolution spectroscopy respectively; both are starting to test QPO geometric models by detecting precession-modulated polarisation angles and energy-resolved reverberation.

Why QPOs are worth the effort

Astrophysical clocks are rare and precious. Coherent radio pulsars give us neutron-star equations of state and Shapiro-delay tests of general relativity, but they tell us nothing about accretion. Black holes have no coherent signal: their spectra are continuum, their event horizons swallow any embedded periodicity. QPOs are the only natural clocks we have inside the innermost ~10 r_g of a stellar-mass black hole, the precise region where strong-field GR effects dominate. Every other technique for measuring black-hole spin — continuum-fitting, iron-line reflection, gravitational-wave ringdown — is a spectral measurement, and each has its own systematics. QPOs are timing measurements, with completely different systematic vulnerabilities. The convergence (or otherwise) of timing- and spectral-derived spins is a check that the underlying picture of GR-dominated accretion is internally consistent.

The pattern of QPOs also encodes the geometry of the inner flow. The fact that type-C LFQPO rms amplitude is larger in high-inclination systems is direct evidence that the modulation is geometric (projected area of a precessing structure) rather than intrinsic (a luminosity oscillation). Combined with X-ray polarimetry from IXPE, this gives an empirical handle on the misalignment between the inner accretion flow and the black-hole spin axis — a key parameter in jet launching, super-Eddington outflows, and the orbital alignment history of the binary itself.

Common pitfalls

• Confusing a QPO with a coherent pulsation. A Q = 5 Lorentzian peak is not a pulsation; its phase drifts. Folding the light curve at the centroid frequency will show a pulse that smears within a few hundred cycles. Treating it as coherent gives spuriously precise spin-parameter constraints.
• Reading 3:2 off noisy data. Many HFQPO claims rest on detections at 3–5σ; the "3:2" frequencies are sometimes set by where the searcher looks. Robust HFQPO pairs need contemporaneous detection of both peaks at well-separated luminosities to rule out chance coincidence.
• Forgetting the deadtime correction. RXTE PCA's microsecond deadtime distorts power spectra at high frequencies; modern reanalyses fit the deadtime model jointly with the QPO. Naive reanalysis of archival data overestimates Q and rms.
• Mixing fluid- and orbital-timescale QPOs. Atoll-source hectohertz QPOs and Z-source horizontal-branch oscillations probably originate in disk pressure modes, not ISCO orbits. Plugging their frequencies into ν_K(ISCO) formulas to derive a "mass" gives nonsense.
• Assuming all type-C LFQPOs are Lense-Thirring. The model is the leading candidate but is not yet proven; alternative magnetic-precession and radiative-instability mechanisms also produce QPOs in the 0.1–10 Hz band. Be wary of model-dependent spin parameters that assume Lense-Thirring without an independent check.