The Evolution of Cataclysmic Variables as Revealed by their Donor Stars

We reconstruct the evolutionary path followed by cataclysmic variables (CVs) from the observed mass-radius relationship of their donor stars. Along the way, we update the semi-empirical CV donor sequence of Knigge (2006) and present a comprehensive review of the link between CV evolution and donor physics. After calibrating state-of-the art stellar models for use in the CV setting, we fit self-consistent theoretical evolution sequences to the observed donor masses and radii. In the standard model of CV evolution, AML below the period gap is assumed to be driven solely by gravitational radiation (GR), while AML above the gap is usually described by a magnetic braking prescription due to Rappaport, Verbunt & Joss (1983). We find that simple scaled versions of these recipes match the data quite well. However, the optimal scaling factors turn out to be f_GR = 2.47 +/- 0.22 below the gap and f_MB = 0.66 +/- 0.05 above. The implications and applications of our results include: (1) The revised evolution sequence yields correct locations for the CV minimum period and the upper edge of the period gap; the standard sequence does not. (2) A comparison of predicted and observed WD temperatures suggests an even higher value for f_GR, but this is sensitive to the assumed WD mass. (3) The absolute donor magnitudes predicted by our sequences can be used to set firm lower limits on the distances toward CVs. (4) Both standard and revised sequences predict that short-period CVs should be susceptible to dwarf nova (DN) eruptions, consistent with observations. However, both sequences also predict that the DNe fraction among long-period CVs should decline with P_orb. Observations suggest the opposite behaviour. (5) The ratio of long-period CVs to short-period, pre-bounce CV is about 3x higher for the revised sequence than the standard one. This may resolve a long-standing problem in CV evolution. [abridged]

💡 Research Summary

This paper presents a comprehensive reconstruction of cataclysmic variable (CV) evolution by exploiting the observed mass‑radius (M‑R) relationship of the donor (secondary) stars. The authors begin by reviewing the standard CV evolution framework, in which angular‑momentum loss (AML) above the orbital‑period gap (~3 h) is driven by magnetic braking (MB) and below the gap (~2 h) by gravitational radiation (GR). While this paradigm successfully explains the existence of the period gap and the minimum period, several long‑standing discrepancies remain: the predicted minimum period (~65 min) is shorter than the observed spike at ~82 min, and population synthesis models over‑predict the dominance of short‑period systems.

The core methodology is to treat donor radii as a secular tracer of the long‑term mass‑transfer rate (Ṁ). The authors argue that other common Ṁ diagnostics—white‑dwarf (WD) temperatures, accretion luminosities, and outburst statistics—are sensitive to short‑term fluctuations and therefore unreliable for probing secular evolution. In contrast, donor radii integrate the mass‑loss history over gigayear timescales, making them ideal for constraining AML.

A detailed analysis of non‑mass‑loss contributions to donor inflation is performed. Tidal and rotational deformation of Roche‑lobe‑filling stars inflates radii by ~4.5 % (below the gap) and ~7.9 % (above). Comparison with measurements of isolated low‑mass stars suggests an additional systematic offset of ~1.5 % for fully convective donors and ~4.9 % for partially radiative donors. Irradiation‑induced bloating is found to be minor, comparable to or smaller than the deformation effects.

Using state‑of‑the‑art stellar evolution models (calibrated against the above offsets), the authors fit self‑consistent CV evolution tracks to a compiled sample of donor masses and radii. They adopt simple scaling factors for the AML prescriptions: AML = f_GR·GR below the gap and AML = f_MB·MB above, where GR is the standard gravitational‑radiation loss and MB follows the Rappaport, Verbunt & Joss (1983) formulation. The best‑fit values are f_GR = 2.47 ± 0.22 and f_MB = 0.66 ± 0.05. Thus, the data demand roughly 2½ times stronger AML than pure GR below the gap, while magnetic braking must be reduced to about two‑thirds of the classic prescription above the gap.

The revised evolutionary tracks reproduce the observed locations of the period minimum (P_min ≈ 82 min) and the upper and lower edges of the period gap, whereas the unscaled standard model fails to do so. Predicted donor spectral types, near‑infrared absolute magnitudes, and WD effective temperatures are compared with observations. The near‑IR donor magnitudes provide robust lower limits on distances, useful for systems lacking parallaxes.

Implications are explored in several domains. First, the donor‑based distance method is demonstrated. Second, the authors examine dwarf‑nova (DN) stability by comparing secular Ṁ from the tracks with the critical Ṁ for disk instability. Both the standard and revised models predict that short‑period CVs should be DN‑prone, consistent with observations, but they also predict a declining DN fraction for long‑period CVs, contrary to the observed increase. The paper discusses possible resolutions, including differences between secular and instantaneous Ṁ, and the role of additional heating mechanisms in the disk.

Third, population synthesis using the revised tracks yields a long‑period‑to‑short‑period pre‑bounce CV ratio about three times larger than the standard model, alleviating the long‑standing “missing long‑period CV” problem. The authors also assess the sensitivity of their results to uncertainties in non‑mass‑loss bloating, sample biases, and the physical plausibility of enhanced AML below the gap.

In the discussion, the revised AML scaling is compared with recent empirical studies (e.g., Littlefair et al. 2008; Sirotkin & Kim 2010). The authors argue that while enhanced AML below the gap may appear at odds with simple magnetic‑braking theories, it could arise from additional mechanisms such as consequential angular‑momentum loss, circumbinary disks, or weak residual magnetic braking in fully convective stars.

In summary, by calibrating donor radii against modern stellar models and accounting for deformation and irradiation effects, the paper provides a robust, donor‑based reconstruction of CV secular evolution. The resulting scaled AML prescriptions reconcile several key observational constraints—period minimum, period gap, donor spectral types, and population ratios—while highlighting remaining tensions (e.g., dwarf‑nova fractions) that point to avenues for future theoretical and observational work.