Open Access
Issue
A&A
Volume 698, May 2025
Article Number L22
Number of page(s) 6
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202554828
Published online 16 June 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

The standard theory of cataclysmic variable (CV) evolution mostly ignores the impact that strong white dwarf (WD) magnetic fields may have on evolutionary sequences (e.g., Belloni & Schreiber 2023), and recent works have shown that this might be the reason for some of the disagreements between predictions and observations. Postulating that strong (i.e., ≥1 MG) magnetic fields can be generated in rapidly rotating crystallizing WDs, Schreiber et al. (2021) proposed a scenario that explains several otherwise inexplicable observational facts. According to this model, WDs that crystallize before mass transfer starts generate a strong magnetic field as soon as accretion has spun up the WD. If the emerging magnetic field connects with that of the donor star, synchronization torques transfer spin angular momentum to the orbit, which causes the separation between the two stars to increase and the CV to convert into a detached binary. Angular momentum loss through gravitational radiation and/or reduced magnetic braking (Belloni et al. 2020) then slowly shrinks the orbit until the donor fills its Roche-lobe and the system becomes a magnetic CV.

This new evolutionary sequence explains why virtually all detached close binaries containing a strongly magnetic WD observed to date are close to Roche-lobe filling and typically have periods between 3–5 h. The sequence also explains the existence of the radio pulsing WD binaries AR Sco (Marsh et al. 2016) and SDSS J191213.72-441045.1 (Pelisoli et al. 2023). In agreement with the proposed crystallization dynamo, Bagnulo & Landstreet (2021) found that also in single WDs strong magnetic fields appear late, that is, at WD ages of 2–3 Gyr, roughly when the WDs start to crystallize. In addition, the proposed dynamo seemed to explain a wide range of related phenomena (Belloni et al. 2021, 2024a; Schreiber et al. 2022).

Recently, Schreiber et al. (2023) suggested that strong WD magnetic fields produced by a rotation- and crystallization-driven dynamo could even contribute to solving the problem of the missing period bouncers. The standard theory for the formation and evolution of CVs predicts that ∼60 − 80% of all CVs should already have passed the period minimum (Kolb 1993; Goliasch & Nelson 2015; Belloni et al. 2018), while the fraction derived from observations ranges from a few percent (Pala et al. 2020; Inight et al. 2023a,b) to a quarter (Rodriguez et al. 2025). As the accretion rate decreases after a CV passes the period minimum, the core temperature also decreases and the WD might eventually crystallize. Assuming that this leads to the generation of a strong WD magnetic field, the binary might detach and no longer be detectable as a period bouncer (Schreiber et al. 2023). The success in explaining the otherwise inexplicable observational facts listed above has made the crystallization-driven dynamo a popular idea for the generation of WD magnetic fields.

However, more recent theoretical works on the origin of WD magnetic fields indicate two main reasons why the dynamo idea might not be viable. First, the energy in the convective zone of crystallizing WDs might not be sufficient to generate strong magnetic fields (e.g., Castro-Tapia et al. 2024, and references therein). Second, the diffusion timescale for the field generated in the core to reach the surface might be too long. In other words, even if the dynamo were somehow able to produce strong fields, WD magnetic fields would be expected to appear even later than observed (Camisassa et al. 2024).

As an alternative explanation for the late appearance of WD magnetic fields, Camisassa et al. (2024) suggested that WD magnetic fields are generated in the convective cores of main-sequence stars and diffuse to the surface during the WD cooling phase. The corresponding diffusion timescales could agree with the observed delay in the appearance of the WD magnetic fields. According to this scenario, for a given WD mass, the appearance of strong magnetic fields depends only on the age of a given WD and not fundamentally on the core temperature as in the crystallization dynamo idea.

To keep our study model independent, for this work we used observational constraints for the appearance of strong WD magnetic fields from single WDs and WDs in detached binaries. As WD masses in CVs cluster around ∼0.8 M we used the age at which single WDs appear to be frequently strongly magnetic as the typical age for the emergence of WD magnetic fields in CVs. We then investigated the implications of magnetic field appearance as a simple age effect on CVs by performing a population synthesis and determining the ages of CV WDs at different orbital periods. We found that if strong magnetic fields appear 2–3 Gyr after the WD formed, the fraction of predicted magnetic WDs in CVs and their orbital periods broadly agree with the observations. In addition, many CVs should become magnetic CVs close to the period minimum. This means that the evolutionary sequence proposed by Schreiber et al. (2021) as well as the suggested idea for solving the problem of the missing period bouncers (Schreiber et al. 2023) do not depend on the assumption that the fields are generated by the crystallization-driven dynamo, but instead remain plausible if a fossil field needs 2–3 Gyr to diffuse to the surface.

2. Population synthesis

To investigate the impact of the late appearance of WD magnetic fields on populations of CVs, we carried out CV population synthesis with the Binary Star Evolution (BSE) code (Hurley et al. 2000, 2002), updated by Belloni et al. (2018). This version of the BSE code includes state-of-the-art prescriptions for CV evolution, which allows us to accurately carry out population synthesis.

We adopted the initial binary population described in detail in Belloni et al. (2024b), corresponding to ≈6.15 × 105 zero-age main-sequence binaries. We picked the primary mass from the canonical Kroupa (2001) initial mass function in the range between 1 and 8 M. The orbital period and mass ratio distributions were constructed following Moe & Di Stefano (2017). These distributions are correlated, that is, the distribution of orbital periods depends on the primary mass as well as the binary fraction, while the distributions of eccentricity and mass ratio depend on the orbital period and the primary mass.

For common-envelope evolution, a low efficiency was adopted by assuming that 25% of the dissipated binary orbital energy is used to expel the common envelope, which is consistent with the evidence that cataclysmic variable progenitors experience strong orbital shrinkage during common-envelope evolution (e.g., Zorotovic et al. 2010). The envelope-structure parameter depends on the red giant properties (Claeys et al. 2014) and contributions from recombination energy were assumed negligible.

The subsequent evolution is driven by orbital angular momentum loss due to magnetic braking, consequential angular momentum loss, and gravitational radiation. We assumed the standard prescription for gravitational radiation and for magnetic braking Rappaport et al. (1983), with the normalization factors determined by Zorotovic et al. (2016). For consequential angular momentum loss, we adopted the prescription put forward by Schreiber et al. (2016).

During the Galactic disk life-time, assumed to be ≈10 Gyr, we adopted a constant star formation rate. This implies that the birth-time distribution of the binary main-sequence stars is uniform. The binary main-sequence stars were then evolved with the BSE code from the birth time until the assumed Galactic age of 10 Gyr.

3. Predicted white dwarf ages in CVs

To compare the model prediction with observations we separated the predicted CVs into two groups, according to the type of the donor. If the donor is still burning hydrogen in its interior, we have a pre-bounce CV, while if the donor is a brown dwarf, we have a period bouncer. For the threshold between the two classes we assumed 0.07 M. We furthermore separated pre-bounce CVs into those with donors that still have a radiative core (with periods exceeding ∼3 h and donor masses ≥ 0.2 M) experiencing strong angular momentum loss through magnetic braking and those that have fully convective donors (periods shorter than ∼3 h and donor masses below 0.2 M).

To illustrate the number of strongly magnetic WDs our population synthesis predicts, we show the period distributions of predicted present-day CVs that did not yet reach the period minimum highlighting the fraction of CVs that are older than the assumed threshold for the appearance of strong WD magnetic fields (Fig. 1). Although most pre-bounce systems host young WDs, a significant fraction of 36%, 23%, and 16% of them host WDs older than 2, 2.5, and 3 Gyr, respectively. For systems above the period gap, that is, CVs with orbital periods longer than ∼3 h, the fractions of systems having WDs older than 2, 2.5, 3 Gyr are 12%, 9%, and 6%, respectively. The same fractions increase to 45%, 27%, and 20% for systems below or inside the gap, with orbital periods shorter than 3 h.

thumbnail Fig. 1.

Orbital period distribution of all pre-bounce CVs. We also show the distributions of pre-bounce CVs with WDs older than 2, 2.5, and 3 Gyr. The overall fractions of systems with WDs older than 2, 2.5, and 3 Gyr are 36%, 23%, 16%, respectively. Above the period gap (i.e., Porb ≳ 3 h), the same fractions are 12%, 9%, and 6%, while for systems inside and below the gap, they are 45%, 27%, and 20%.

In Fig. 2 we include predicted period bouncers and show the age distributions of CVs. The distribution of all CVs is dominated by that of period bouncers at ages ≳ 2 Gyr, while at lower ages it is dominated by pre-bounce systems, with peaks at ∼0.25 Gyr (above the gap) and ∼1.75 Gyr (below or inside the gap). The fraction of period bouncers older than 2, 2.5, and 3 Gyr is 99.6, 96.9, and 92.8%, respectively. The absolute number of CVs with WDs exceeding a given age are given in Table 1.

thumbnail Fig. 2.

WD age distribution for different groups of CVs. The distribution of all CVs is dominated by that of period bouncers at ages ≳ 2 − 3 Gyr, while at younger ages it is dominated by pre-bounce systems, with peaks at ∼0.25 Gyr (above the gap) and ∼1.75 Gyr (below or inside the gap).

Table 1.

Numbers of CVs predicted in our population synthesis.

Given that the emergence of the magnetic field may cause a detached phase (generated by the transfer of spin angular momentum to the orbit), which might take several Gyr if the WD magnetic field appears at the period minimum or beyond (see Schreiber et al. 2023 and Appendix A), it is important to evaluate at which period the field appeared for the predicted present-day CVs. The distribution of secondary masses at which the present-day CVs became strongly magnetic (assuming as above three different ages for the appearance of the field) is shown in Fig. 3. Interestingly, for ≈35%, ≈69%, and ≈81% of the current period bouncers the WD reached the age of 2, 2.5, and 3 Gyr when the donor mass was below 0.08 M, that is, close to or after the system reached the period minimum.

thumbnail Fig. 3.

Donor mass distribution of the present-day CVs when their WD ages reached 2, 2.5, and 3 Gyr illustrating at which donor mass a strong magnetic field could have appeared. The peak of the distributions is at 0.076, 0.045, and 0.033 M for age limits of 2, 2.5, and 3 Gyr respectively. For 35%, 69%, and 81% of the present-day period bouncers the magnetic field appeared at donor masses ≤ 0.08 M, that is, after or around the period minimum.

4. Strongly magnetic WDs in CVs

We simulated CV populations and determined the age distribution of the accreting WDs. In what follows we estimate the number of strongly magnetic (≥1 MG) white dwarfs in the simulated sample and compare these predictions with observations.

4.1. Predicting magnetic WDs in CVs

To estimate when a strong magnetic field appears we used the available observational constraints from single WDs and WDs in detached binaries. First, we assumed that WDs that formed through common envelope evolution are all born without a strong magnetic field. While a few very young (i.e., ≤0.6 Gyr) and massive single WDs are strongly magnetic (Bagnulo & Landstreet 2022, their Fig. 2), we assumed that these WDs become strongly magnetic early through a channel that is not available to close binaries (e.g., the merging of two C/O core WDs). This is reasonable because not a single young and strongly magnetic WD has been found in thousands of detached post-common envelope binaries (Liebert et al. 2005; Schreiber et al. 2010, 2021). Instead, all strongly magnetic WDs found among these detached progenitors of CVs are older than ∼2 Gyr (Parsons et al. 2021; van Roestel et al. 2025).

We assumed a fixed WD age for the appearance of strong WD magnetic fields of 2, 2.5, and 3 Gyr. We selected these thresholds because of the ages of magnetic WDs in detached binaries (see above) and because the fraction of strongly magnetic (i.e., ≥1 MG) single WDs with masses typical for CVs (i.e., 0.75–0.95 M Zorotovic et al. 2011; Pala et al. 2022) significantly increases beyond these ages. In the sample analyzed by Bagnulo & Landstreet (2022, their Fig. 2) only one out of 27 (= 3.7 ± 3.6%) WDs younger than 2 Gyr is strongly magnetic, while this fraction increases to 50 ± 17% (4/8) and 60 ± 22% (3/5) for WDs older than 2 and 3 Gyr, respectively. We thus assumed that 30–80% of white dwarfs in CVs become magnetic as soon as they pass the defined age threshold of 2–3 Gyr. The large range of possible fractions (i.e., 30–80%) of strongly magnetic WDs results from the small number of observed WDs with masses typical for CVs (0.75–0.95 M) that are older than 2 Gyr.

Unsurprisingly, based on the above assumptions, our simulations predict a steady increase in the fraction of magnetic WDs moving from CVs above the gap, to CVs in and below the gap, to period bouncers. Taking into account that 30–80% of WDs that pass the age threshold become strongly magnetic, the predicted fraction of strongly magnetic WDs in CVs is ∼2 − 10% above the period gap, 6–36% in and below the gap, and reaches 28–80% for period bouncers. In addition, in 7–65% of all present-day period bouncers the magnetic field is predicted to appear when the donor star mass decreases below 0.08 M. In what follows we relate these findings to observed samples of CVs and discuss implications for the problem of the missing period bouncers.

4.2. Comparison with observations

The small but almost complete 150 pc sample of CVs (Pala et al. 2020) contains 1/6 = 17 ± 14% magnetic systems above the gap, and this number increases for systems in and below the gap to 14/33 = 40 ± 9%. While clearly subject to low number statistics, this is broadly consistent with the 2–10 and 6–36% of CVs with strongly magnetic WDs in our predicted CV population. Only three period bouncers are in the 150 pc sample and none of them contains a magnetic WD. This seems to disagree with our prediction, but magnetic period bouncers might be absent in observed samples as they detach for long periods of time when the field appears close to the period minimum or beyond.

For the larger SDSS CV samples from Inight et al. (2023a,b, 2025) we find the fraction of magnetic CVs to be 16.5 ± 2.6% above the period gap and 24.4 ± 2.1% otherwise. Both values are still in reasonable agreement with the prediction, although the predicted fraction of strongly magnetic WDs above the period gap is slightly smaller than observed (by a factor of ∼1.5). Unfortunately, it remains unclear whether the differences between the larger SDSS sample and the small volume-limited sample is mainly caused by observational biases or by low number statistics. However, we believe that observational biases are more likely the cause of the differences between the two samples as the period gap also looks different for both samples (Schreiber et al. 2024). A large and complete volume-limited sample of CVs is needed to perform more detailed comparisons between observations and theoretical predictions. The predicted and observed fractions are summarized in Table 2.

Table 2.

Observed and predicted fractions of strongly magnetic WDs in CVs.

Given the general agreement, however, we conclude that if the age of the WD is the one and only essential parameter for producing strong magnetic fields, as perhaps indicated by studies of single WDs (Bagnulo & Landstreet 2021, 2022), the predicted number of magnetic CVs roughly agrees with the observations if most WDs in CVs become strongly magnetic at ages of 2–3 Gyr.

4.3. Magnetism and detachment in CVs

As described above, an age-dependent appearance of strong magnetic fields in WDs provides a reasonable agreement between the observed and predicted fractions of magnetic WDs in CVs. If this scenario for the occurrence of magnetic WDs in CVs is correct, a relatively large number of period bouncers should contain strongly magnetic WDs (∼27 − 80%).

As described in the introduction, Schreiber et al. (2021) developed an evolutionary sequence connecting detached binaries containing a magnetic WD (Parsons et al. 2021), WD pulsars (Marsh et al. 2016; Pelisoli et al. 2023), and magnetic CVs. The idea of emerging magnetic fields converting CVs to detached systems has later been extended to period bouncers (Schreiber et al. 2023). However, in the latter work, the appearance of the magnetic field was assumed to depend on the onset of crystallization, which in turn depends on the density, the composition, and the core temperature of the WD. Given the recent criticism on the crystallization dynamo idea (Castro-Tapia et al. 2024; Camisassa et al. 2024) this scenario might not apply.

In Appendix A we present MESA calculations illustrating that the same gigayear-long detachment is predicted if the magnetic field appears close to the period minimum or beyond because the WD age reached the threshold. Figure A.1 shows that the detached phase takes longer when the companion’s mass is lower. From our population synthesis we estimated that in 7–65% of all present-day period bouncers, the magnetic field appeared close to the period minimum or later. If the occurrence of the magnetic field in these systems indeed causes a gigayear-long detached phase, the predicted number of accreting period bouncers is potentially significantly reduced. This represents an important result because it implies that the evolutionary sequence proposed by Schreiber et al. (2021) and also the idea of reducing the predicted number of period bouncers through the appearance of WD magnetic fields do not depend on postulating the existence of a crystallization- and rotation-driven dynamo. It is unclear whether the crystallization- and rotation-driven dynamo can explain the relative numbers as effectively as the purely time dependent appearance of WD magnetic fields. This is because for high accretion rates the WD core is heated on timescales of ∼108 yr, which may prevent accreting WDs from crystallizing, and might melt previously crystallized cores (Schreiber et al. 2023).

The idea of period bouncers detaching because of the appearance of the magnetic field predicts a large number of magnetic WDs with brown dwarf companions that are not filling their Roche-lobe. A few systems that might fit this description have been identified recently through cyclotron (van Roestel et al. 2025) or X-ray (Muñoz-Giraldo et al. 2023; Cunningham et al. 2025) emission. However, the interpretation of these objects as detached period bouncers remains uncertain as the observed low accretion rates might also result from low states, as first observed in EF Eri (Filor et al. 2024).

4.4. Dependence on unconstrained parameters

While the outlined scenario appears plausible, we have to admit that our predictions depend on a number of assumed parameters that are neither observationally nor theoretically well constrained. While the strength of magnetic braking in CVs is relatively well constrained by the measured donor star radii (McAllister et al. 2019) and the onset of the period gap (Knigge et al. 2011; Schreiber et al. 2024), magnetic braking might be significantly stronger than predicted by the standard prescription (Rappaport et al. 1983) for binaries with longer orbital periods.

There is growing evidence that magnetic braking saturates in close binaries (El-Badry et al. 2022; Belloni et al. 2024b), which implies orders of magnitude stronger angular momentum loss for detached binaries with periods between a few hours and a few days (see also Barraza-Jorquera et al. 2025). Stronger angular momentum loss between common envelope evolution and the second mass transfer phase (the CV stage) would lead to generally younger CVs, and thus would change the fractions of 2–3 Gyr CVs predicted by population synthesis. Assuming the same WD age for strong field appearance (2–3 Gyr), fewer systems would be predicted to be magnetic among pre-bounce CVs, while the fraction of period bouncers among all CVs would further increase, as would the number of period bouncers that detach.

The second parameter that could impact the predictions is the assumed age at which the WD magnetic field appears. We simply assumed constant ages between 2 and 3 Gyr when estimating the fraction of strongly magnetic WDs for different evolutionary stages (orbital period ranges). At these ages single WDs with typical CV WD masses are frequently strongly magnetic (Bagnulo & Landstreet 2022). However, given the small number of observed WDs in this parameter range we cannot exclude a younger age for the field appearance as, for example, the diffusion timescale for typical CV WD masses (0.83 M, Zorotovic et al. 2011) seems to be shorter than 2 Gyr (Camisassa et al. 2024). This would have the opposite effect of decreased evolutionary timescales toward the CV phase, that is, this would cause the number of predicted magnetic pre-bounce CVs to increase and the number of detached period bouncers would decrease. The general agreement between observations and predictions would persist.

We also note that the durations of the detached phases calculated in Schreiber et al. (2023) and in the Appendix assume a high efficiency for the accreted angular momentum to spin up the WD, a synchronization timescale (1 Myr), and that all the spin angular momentum is transferred to the orbit. Even for magnetic fields appearing in period bouncers, the detached phase might well be shorter than we predict here. In this case, the number of accreting period bouncers would be reduced by a smaller amount, but the prediction of a relatively large number of magnetic WDs in period bouncers would remain.

Finally, as discussed above, the fraction of magnetic WDs in observed samples of CVs that have not yet passed the period minimum is slightly higher than predicted (see Table 2). Perfect agreement could be reached by assuming a higher fraction of WDs becoming strongly magnetic at the threshold or assuming a lower age limit for the field appearance. However, this would lead to predicting fewer detached period bouncers.

We conclude that assuming that the late appearance of WD magnetic fields can be modeled as a simple age effect leads to predicted fractions of strongly magnetic CVs that are in general agreement but slightly lower than observed. For now, investigating the full parameter space to potentially reach even better agreement would represent a rather futile exercise as several key processes, such as angular momentum loss through magnetic braking, remain poorly understood and representative and large samples of observed single WDs (with typical CV WD masses) and CVs are not yet available.

5. Conclusions

We performed binary population synthesis of CVs and assumed that strong WD magnetic fields appear at WD ages of 2–3 Gyr. Our simulations predict relative numbers of magnetic CVs above and below the period gap that are broadly in agreement with the observations and that a large fraction of period bouncers should contain magnetic WDs. Furthermore, in 7–65% of the present-day CVs the magnetic field appeared when the system had reached or evolved past the period minimum. At this evolutionary stage the appearance of the WD magnetic fields and the transfer of spin angular momentum to the orbit through synchronization torques can cause extended detached phases (several gigayears). We therefore conclude that the evolutionary sequence connecting AR Sco, detached binaries containing a strongly magnetic WDs (Parsons et al. 2021), and magnetic CVs above the gap (Schreiber et al. 2021), and also the idea of reducing the number of predicted accreting period bouncers through detachment (Schreiber et al. 2023), do not depend on assuming that the fields are generated by the crystallization- and rotation-driven dynamo, but instead could be reproduced by simply assuming that strong magnetic fields on the surface of WDs appear as a result of aging.


Acknowledgments

MRS acknowledge financial support from FONDECYT (grant number 1221059). DB acknowledges support from FONDECYT (grant number 3220167) and the São Paulo Research Foundation (FAPESP), Brazil, Process Numbers #2024/03736-2 and #2025/00817-4.

References

  1. Bagnulo, S., & Landstreet, J. D. 2021, MNRAS, 507, 5902 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bagnulo, S., & Landstreet, J. D. 2022, ApJ, 935, L12 [NASA ADS] [CrossRef] [Google Scholar]
  3. Barraza-Jorquera, J. A., Schreiber, M. R., & Belloni, D. 2025, A&A, 696, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Belloni, D., & Schreiber, M. R. 2023, in Handbook of X-ray and Gamma-ray Astrophysics, eds. C. Bambi, & A. Santangelo, 129 [Google Scholar]
  5. Belloni, D., Schreiber, M. R., Zorotovic, M., et al. 2018, MNRAS, 478, 5626 [NASA ADS] [CrossRef] [Google Scholar]
  6. Belloni, D., Schreiber, M. R., Pala, A. F., et al. 2020, MNRAS, 491, 5717 [Google Scholar]
  7. Belloni, D., Schreiber, M. R., Salaris, M., Maccarone, T. J., & Zorotovic, M. 2021, MNRAS, 505, L74 [NASA ADS] [CrossRef] [Google Scholar]
  8. Belloni, D., MikoÅ‚ajewska, J., & Schreiber, M. R. 2024a, A&A, 686, A226 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Belloni, D., Schreiber, M. R., Moe, M., El-Badry, K., & Shen, K. J. 2024b, A&A, 682, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Camisassa, M., Fuentes, J. R., Schreiber, M. R., et al. 2024, A&A, 691, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Castro-Tapia, M., Zhang, S., & Cumming, A. 2024, ApJ, 975, 63 [NASA ADS] [CrossRef] [Google Scholar]
  12. Claeys, J. S. W., Pols, O. R., Izzard, R. G., Vink, J., & Verbunt, F. W. M. 2014, A&A, 563, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Cunningham, T., Caiazzo, I., Sienkiewicz, G., et al. 2025, MNRAS, 540, 633 [Google Scholar]
  14. El-Badry, K., Conroy, C., Fuller, J., et al. 2022, MNRAS, 517, 4916 [NASA ADS] [CrossRef] [Google Scholar]
  15. Filor, L. W., Mori, K., Bridges, G., et al. 2024, ApJ, submitted [arXiv:2412.11273] [Google Scholar]
  16. Goliasch, J., & Nelson, L. 2015, ApJ, 809, 80 [Google Scholar]
  17. Hurley, J. R., Pols, O. R., & Tout, C. A. 2000, MNRAS, 315, 543 [Google Scholar]
  18. Hurley, J. R., Tout, C. A., & Pols, O. R. 2002, MNRAS, 329, 897 [Google Scholar]
  19. Inight, K., Gänsicke, B. T., Breedt, E., et al. 2023a, MNRAS, 524, 4867 [NASA ADS] [CrossRef] [Google Scholar]
  20. Inight, K., Gänsicke, B. T., Schwope, A., et al. 2023b, MNRAS, 525, 3597 [NASA ADS] [CrossRef] [Google Scholar]
  21. Inight, K., Gänsicke, B. T., Schwope, A., et al. 2025, MNRAS, 536, 1057 [Google Scholar]
  22. Jermyn, A. S., Bauer, E. B., Schwab, J., et al. 2023, ApJS, 265, 15 [NASA ADS] [CrossRef] [Google Scholar]
  23. Knigge, C., Baraffe, I., & Patterson, J. 2011, ApJS, 194, 28 [Google Scholar]
  24. Kolb, U. 1993, A&A, 271, 149 [NASA ADS] [Google Scholar]
  25. Kroupa, P. 2001, MNRAS, 322, 231 [NASA ADS] [CrossRef] [Google Scholar]
  26. Liebert, J., Wickramasinghe, D. T., Schmidt, G. D., et al. 2005, AJ, 129, 2376 [NASA ADS] [CrossRef] [Google Scholar]
  27. Marsh, T. R., Gänsicke, B. T., Hümmerich, S., et al. 2016, Nature, 537, 374 [NASA ADS] [CrossRef] [Google Scholar]
  28. McAllister, M., Littlefair, S. P., Parsons, S. G., et al. 2019, MNRAS, 486, 5535 [NASA ADS] [CrossRef] [Google Scholar]
  29. Moe, M., & Di Stefano, R. 2017, ApJS, 230, 15 [Google Scholar]
  30. Muñoz-Giraldo, D., Stelzer, B., de Martino, D., & Schwope, A. 2023, A&A, 676, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Pala, A. F., Gänsicke, B. T., Breedt, E., et al. 2020, MNRAS, 494, 3799 [NASA ADS] [CrossRef] [Google Scholar]
  32. Pala, A. F., Gänsicke, B. T., Belloni, D., et al. 2022, MNRAS, 510, 6110 [NASA ADS] [CrossRef] [Google Scholar]
  33. Parsons, S. G., Gänsicke, B. T., Schreiber, M. R., et al. 2021, MNRAS, 502, 4305 [CrossRef] [Google Scholar]
  34. Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [Google Scholar]
  35. Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4 [Google Scholar]
  36. Paxton, B., Marchant, P., Schwab, J., et al. 2015, ApJS, 220, 15 [Google Scholar]
  37. Paxton, B., Schwab, J., Bauer, E. B., et al. 2018, ApJS, 234, 34 [NASA ADS] [CrossRef] [Google Scholar]
  38. Paxton, B., Smolec, R., Schwab, J., et al. 2019, ApJS, 243, 10 [Google Scholar]
  39. Pelisoli, I., Marsh, T. R., Buckley, D. A. H., et al. 2023, Nat. Astron., 7, 931 [NASA ADS] [CrossRef] [Google Scholar]
  40. Rappaport, S., Verbunt, F., & Joss, P. C. 1983, ApJ, 275, 713 [Google Scholar]
  41. Rodriguez, A. C., El-Badry, K., Suleimanov, V., et al. 2025, PASP, 137, 014201 [NASA ADS] [CrossRef] [Google Scholar]
  42. Schreiber, M. R., Gänsicke, B. T., Rebassa-Mansergas, A., et al. 2010, A&A, 513, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Schreiber, M. R., Zorotovic, M., & Wijnen, T. P. G. 2016, MNRAS, 455, L16 [NASA ADS] [CrossRef] [Google Scholar]
  44. Schreiber, M. R., Belloni, D., Gänsicke, B. T., Parsons, S. G., & Zorotovic, M. 2021, Nat. Astron., 5, 648 [NASA ADS] [CrossRef] [Google Scholar]
  45. Schreiber, M. R., Belloni, D., Zorotovic, M., et al. 2022, MNRAS, 513, 3090 [NASA ADS] [CrossRef] [Google Scholar]
  46. Schreiber, M. R., Belloni, D., & van Roestel, J. 2023, A&A, 679, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Schreiber, M. R., Belloni, D., & Schwope, A. D. 2024, A&A, 682, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. van Roestel, J., Rodriguez, A. C., Szkody, P., et al. 2025, A&A, 696, A242 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Zorotovic, M., Schreiber, M. R., Gänsicke, B. T., & Nebot Gómez-Morán, A. 2010, A&A, 520, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Zorotovic, M., Schreiber, M. R., & Gänsicke, B. T. 2011, A&A, 536, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Zorotovic, M., Schreiber, M. R., Parsons, S. G., et al. 2016, MNRAS, 457, 3867 [Google Scholar]

Appendix A: MESA simulations

To evaluate the plausibility of the suggestion that the number of period bouncers is reduced by the late appearance of magnetic fields as suggested by Schreiber et al. (2023), the duration of the detached phase generated by the appearance of magnetic fields plays an important role. The duration is of the order of several tens of Myr above the gap (Schreiber et al. 2021) but can reach ∼ 1 Gyr if coupled with the binary entering the period gap (e.g., Schreiber et al. 2023, red track in their fig. 1).

To evaluate how long the detached phase may take at shorter periods, we carried out MESA simulations of CVs that reach the age of field appearance at different orbital periods. We used the version r15140 of the MESA code (Paxton et al. 2011, 2013, 2015, 2018, 2019; Jermyn et al. 2023) to calculate CV evolution. Our approach follows closely that by Schreiber et al. (2023)1. The only difference is that we assumed that the magnetic field appears on the WD surface when its age has a certain value, instead of assuming the rotation- and crystallization-driven dynamo. To illustrate the impact of late magnetic field appearance around the period minimum and during period bouncer evolution, we assumed an age of 2 Gyr.

We simulated four systems and started all simulations immediately after the WD formation. We assumed a WD mass of 0.8 M and an initial companion mass of 0.6 M. The initial orbital periods (0.4, 0.5, 0.6, 0.7 d) were chosen so that the WD magnetic fields appear close to the period minimum or during period bouncer evolution.

Figure A.1 shows the evolutionary tracks we obtained and clearly illustrates that the duration of the detached phase increases with decreasing donor mass at which the magnetic field appears. This is simply because the relative importance of the spin angular momentum of the WD is increasing for lower donor masses. The detached phase therefore can become fundamentally important for period bouncers.

thumbnail Fig. A.1.

Donor mass against orbital period during evolution below the period gap. The magnetic field is assumed to appear on the WD surface when it is 2 Gyr old and the timescales for the detached phase triggered by the magnetic field appearance are given in the figure. The lower the donor mass when the magnetic field appears, the stronger the impact of the synchronizing torque, and consequently, the longer the detachment phase. For magnetic fields appearing around the period minimum or during the period bouncer evolution, the detachment phases last ≳1 Gyr.

All Tables

Table 1.

Numbers of CVs predicted in our population synthesis.

Table 2.

Observed and predicted fractions of strongly magnetic WDs in CVs.

All Figures

thumbnail Fig. 1.

Orbital period distribution of all pre-bounce CVs. We also show the distributions of pre-bounce CVs with WDs older than 2, 2.5, and 3 Gyr. The overall fractions of systems with WDs older than 2, 2.5, and 3 Gyr are 36%, 23%, 16%, respectively. Above the period gap (i.e., Porb ≳ 3 h), the same fractions are 12%, 9%, and 6%, while for systems inside and below the gap, they are 45%, 27%, and 20%.

In the text
thumbnail Fig. 2.

WD age distribution for different groups of CVs. The distribution of all CVs is dominated by that of period bouncers at ages ≳ 2 − 3 Gyr, while at younger ages it is dominated by pre-bounce systems, with peaks at ∼0.25 Gyr (above the gap) and ∼1.75 Gyr (below or inside the gap).

In the text
thumbnail Fig. 3.

Donor mass distribution of the present-day CVs when their WD ages reached 2, 2.5, and 3 Gyr illustrating at which donor mass a strong magnetic field could have appeared. The peak of the distributions is at 0.076, 0.045, and 0.033 M for age limits of 2, 2.5, and 3 Gyr respectively. For 35%, 69%, and 81% of the present-day period bouncers the magnetic field appeared at donor masses ≤ 0.08 M, that is, after or around the period minimum.

In the text
thumbnail Fig. A.1.

Donor mass against orbital period during evolution below the period gap. The magnetic field is assumed to appear on the WD surface when it is 2 Gyr old and the timescales for the detached phase triggered by the magnetic field appearance are given in the figure. The lower the donor mass when the magnetic field appears, the stronger the impact of the synchronizing torque, and consequently, the longer the detachment phase. For magnetic fields appearing around the period minimum or during the period bouncer evolution, the detachment phases last ≳1 Gyr.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.