Magnetic Braking Model (MBM)

• Magnetic Braking Model (MBM) is a framework for describing angular momentum loss through magnetized stellar winds, incorporating dependencies like stellar rotation, convective turnover, and magnetic topology.
• MBMs are employed to model the spin-down of low-mass stars and the orbital evolution of close binaries, including cataclysmic variables and low-mass X-ray binaries.
• MBMs encompass a range of prescriptions—from classical Skumanich-type laws to saturated and topology-modulated models—each calibrated against observables such as cluster rotation sequences and period gaps in binaries.

Searching arXiv for recent and foundational papers on magnetic braking models in stars and close binaries. Magnetic Braking Model (MBM) denotes a class of prescriptions for angular-momentum loss through magnetized stellar winds, used to model the spin-down of low-mass stars and the secular evolution of close binaries such as cataclysmic variables (CVs), low-mass X-ray binaries (LMXBs), black-hole LMXBs, polars, and AM CVn progenitors. Across the literature, an MBM is typically defined by a torque law $\dot{J}_{\rm MB}$ or $\tau_{\rm MB}$ coupled to stellar structure and orbital evolution, with dependencies that may include rotation rate, Rossby number, convective turnover time, magnetic topology, wind mass-loss rate, and whether the donor retains a radiative core. The term is therefore not restricted to a single formula: it encompasses Skumanich-type laws, saturated prescriptions, topology-modulated laws, double-dynamo formalisms, and empirical boost/disruption models, each calibrated against different observables such as cluster rotation sequences, the CV period gap, the period minimum spike, LMXB mass-transfer rates, or the period distribution of detached eclipsing binaries (Gossage et al., 2020, Sarkar et al., 2024, El-Badry et al., 2022).

1. Conceptual scope and physical basis

In the standard physical picture, a magnetized stellar wind co-rotates with the star out to an Alfvén radius $R_A$, so that the wind extracts specific angular momentum $\sim \Omega R_A^2$. For low-mass stars with outer convection zones, magnetic activity is tied to a rotation–convection dynamo, and magnetic braking acts as the dominant long-term mechanism for spin-down; in synchronized close binaries, the same torque is communicated to the orbit through tides and thereby drives Roche-lobe overflow and orbital evolution (Gossage et al., 2020, Van et al., 2018).

Several distinct physical ingredients recur across MBMs. One is the rotation dependence of the torque, often expressed as a steep Skumanich-like $\Omega^3$ law in the unsaturated regime and a weaker $\Omega$ scaling in the saturated regime (Gossage et al., 2020, Sarkar et al., 2023). A second is the role of convective turnover time $\tau_c$, either through the Rossby number $R_o=P_{\rm rot}/\tau_c$ or directly through power-law factors in the torque (Gossage et al., 2020, Echeveste et al., 2024). A third is the donor’s internal structure: the presence or absence of a radiative core is central in CV models that invoke “disrupted” or sharply weakened braking at the fully convective boundary (Sarkar et al., 2024, Dodon et al., 15 May 2026). A fourth is magnetic geometry: some models weaken the torque for complex multipolar fields or for field-line closure in dead zones (Gossage et al., 2020, Belloni et al., 2019).

This diversity reflects the fact that MBM is both a physical mechanism and a phenomenological framework. Some prescriptions are calibrated to single-star rotation evolution, some to interacting-binary observables, and some attempt unification across both domains. A plausible implication is that the term “MBM” is best understood as a model family defined by how it closes the relation between stellar structure, magnetic activity, wind properties, and angular-momentum loss, rather than by a unique canonical equation.

2. Classical and saturation-based prescriptions

The classical formulation in binary evolution is the Skumanich-type law. In CV and LMXB work this is commonly written in the Rappaport–Verbunt–Joss form, with $\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$, often with $\gamma=3$ or $\tau_{\rm MB}$0 (Belloni et al., 2019, Van et al., 2018, Echeveste et al., 2024). In detached and tidally locked low-mass binaries, this unsaturated scaling predicts strong orbital shrinkage at short periods and a short-period period distribution $\tau_{\rm MB}$1 (El-Badry et al., 2022). Observationally, however, the ZTF detached eclipsing-binary distribution is “basically flat,” which was found to be “strongly inconsistent with predictions of classical MB models based on the Skumanich relation” (El-Badry et al., 2022).

This discrepancy motivated saturation-based MBMs. In these, magnetic activity saturates below a critical Rossby number or above a critical angular velocity, causing the torque to transition from $\tau_{\rm MB}$2 to $\tau_{\rm MB}$3 (Gossage et al., 2020, El-Badry et al., 2022). In the Matt et al. (2015) implementation used in MESA, the unsaturated torque is

$\tau_{\rm MB}$4

while the saturated torque is

$\tau_{\rm MB}$5

with $\tau_{\rm MB}$6, $\tau_{\rm MB}$7, and a calibrated $\tau_{\rm MB}$8 in the Gossage et al. implementation (Gossage et al., 2020). The structural prefactor $\tau_{\rm MB}$9 scales as

$R_A$0

with $R_A$1 (Gossage et al., 2020).

A related saturated view is supported by detached eclipsing binaries: the observed orbital-period distribution is best reproduced by torque laws “that scale roughly as $R_A$2, as opposed to $R_A$3 in the standard Skumanich law” (El-Badry et al., 2022). That work further found “no significant difference between the period distributions of binaries containing fully and partially convective stars,” supporting saturated MB in tidally locked binaries across the convective boundary (El-Badry et al., 2022). This directly challenges CV scenarios that require a complete switch-off of MB at full convection, although it does not by itself exclude reduced braking in mass-transferring systems.

3. Rossby number, convection, and magnetic topology

Many modern MBMs are organized around the Rossby number $R_A$4, because it links rotation and convection in dynamo theory. In stellar-evolution implementations, $R_A$5 may be evaluated from local convective properties, for example

$R_A$6

evaluated near the base of the convection zone (Gossage et al., 2020), or in local and global-envelope forms derived directly from the stellar structure (Dodon et al., 15 May 2026).

The Garraffo et al. (2018) prescription adds explicit magnetic-topology dependence. Its base dipolar torque is

$R_A$7

and the full torque is

$R_A$8

with a simplified topology factor

$R_A$9

and magnetic complexity

$\sim \Omega R_A^2$0

In the calibrated MESA implementation, $\sim \Omega R_A^2$1, $\sim \Omega R_A^2$2, and $\sim \Omega R_A^2$3 (Gossage et al., 2020). At low $\sim \Omega R_A^2$4, $\sim \Omega R_A^2$5 is large, fields are complex, and the torque is suppressed; at intermediate $\sim \Omega R_A^2$6, $\sim \Omega R_A^2$7 and Skumanich-like spin-down is recovered; at high $\sim \Omega R_A^2$8, complexity rises again and braking weakens (Gossage et al., 2020).

This topology-modulated view is conceptually distinct from simple saturation. In the Matt et al. framework, rapid rotators are weakly braked because the torque saturates. In the Garraffo et al. framework, rapid rotators are weakly braked because their fields are magnetically complex (Gossage et al., 2020). In open-cluster comparisons, the former matches the slow sequence well but tends to overproduce rapid rotators, whereas the latter reproduces bimodality more naturally but can over-brake sub-solar masses at $\sim \Omega R_A^2$9 Myr to 1 Gyr (Gossage et al., 2020).

Convective turnover time itself is also uncertain. Several prescriptions use powers of $\Omega^3$0: Matt et al. (2015) scales as $\Omega^3$1, Garraffo et al. (2018) as $\Omega^3$2, and Van et al. (2018) as $\Omega^3$3 (Sarkar et al., 2024). Because empirical and MLT-based $\Omega^3$4 values can differ by factors of several at low masses, MB strengths can shift by factors of several as well (Sarkar et al., 2024). This sensitivity is central in both stellar spin-down and close-binary evolution.

4. Unified and double-dynamo formulations

One major extension of MBM is the double-dynamo (DD) formalism developed to connect CVs, AM CVn stars, and single low-mass stars. In this framework, stars with $\Omega^3$5 can host two distinct $\Omega^3$6 dynamos: a boundary-layer dynamo at the radiative–convective interface and a convective-envelope dynamo in the outer envelope or throughout the whole star if fully convective (Sarkar et al., 2024). The total torque is written schematically as

$\Omega^3$7

where $\Omega^3$8 vanishes when the radiative core disappears, leaving only $\Omega^3$9 in fully convective stars (Sarkar et al., 2024).

This transition is used to explain the CV period gap. Above the gap, partially radiative donors are dominated by the boundary-layer dynamo and evolve from $\Omega$0 hr to $\Omega$1 hr. At $\Omega$2 the donor becomes fully convective, the boundary layer disappears, strong MB ceases abruptly, and the system detaches near $\Omega$3 hr. It later re-establishes contact near $\Omega$4 hr under gravitational radiation plus weaker convective-dynamo MB (Sarkar et al., 2024). Below the gap, the convective-dynamo torque increases with decreasing mass and shifts the period minimum from the GR-only value $\Omega$5 min to $\Omega$6 min, matching the observed spike at $\Omega$7–$\Omega$8 min (Sarkar et al., 2024).

The same DD prescription has been extended to the spin-down of fully convective M dwarfs. In that application, the torque is built to scale as $\Omega$9 or $\tau_c$0 for fast rotators and $\tau_c$1 or $\tau_c$2 for slow rotators, with free parameters $\tau_c$3 and $\tau_c$4; an example fit uses $\tau_c$5, $\tau_c$6 (Sarkar et al., 2024). This is consistent with the dedicated FCMD study, which models fully convective M dwarfs with a torque approximately scaling as $\tau_c$7 at shorter periods before transitioning into a Skumanich-type $\tau_c$8 torque, and includes a parametrized reduction of wind mass loss due to wind entrapment in dead zones (Sarkar et al., 2023). That work finds that the model matches open clusters such as NGC2547, Pleiades, NGC2516, and Praesepe, explains the long spin periods of field stars, and predicts that a spread in spin distribution persists till over 3 Gyr (Sarkar et al., 2023).

A plausible implication is that the DD formalism is one of the few MBMs explicitly designed as a cross-domain model, rather than a prescription calibrated to a single population.

5. Empirical boost, disruption, and reduced-wind models

Another important branch of MBM development treats the torque as saturated, boosted, and disrupted. In the SBD formulation, the base saturated torque is

$\tau_c$9

with $R_o=P_{\rm rot}/\tau_c$0, $R_o=P_{\rm rot}/\tau_c$1, and $R_o=P_{\rm rot}/\tau_c$2 (Dodon et al., 15 May 2026). The SBD extension applies a boost factor $R_o=P_{\rm rot}/\tau_c$3 above the fully convective boundary and a disruption factor $R_o=P_{\rm rot}/\tau_c$4 below it: $R_o=P_{\rm rot}/\tau_c$5 (Dodon et al., 15 May 2026).

A more physical interpretation of these factors is provided by the irradiation–$R_o=P_{\rm rot}/\tau_c$6-based SBD model, or i$R_o=P_{\rm rot}/\tau_c$7SBD. There, $R_o=P_{\rm rot}/\tau_c$8 is computed directly from the donor structure rather than from an empirical mass relation, and irradiation from the accreting white dwarf drives an additional wind that enhances MB (Dodon et al., 15 May 2026). The structure-based $R_o=P_{\rm rot}/\tau_c$9 calculation reveals a robust 2–3× spike as the donor approaches full convection, and because the saturated torque scales as $\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$0, this spike yields a disruption factor

$\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$1

typically $\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$2–9, sufficient to initiate the period gap (Dodon et al., 15 May 2026). The irradiation-driven wind leads to a boost factor

$\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$3

with $\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$4–0.5 and plausible CV parameters giving effective boosts $\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$5 during accreting phases above the gap (Dodon et al., 15 May 2026).

This framework links the empirical $\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$6 and $\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$7 parameters to two explicit mechanisms: irradiation-driven winds for boosting and a structure-driven $\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$8 spike for disruption (Dodon et al., 15 May 2026). The resulting evolutionary tracks are reported to be broadly consistent with the main CV properties, including the period gap, the period minimum, donor inflation, and some high-$\dot{J}_{\rm MB}\propto M R^\gamma \Omega^3$9 nova-like states (Dodon et al., 15 May 2026).

By contrast, in polars a different “reduced MB” picture is favored. There the strong white-dwarf field suppresses the donor’s wind zone by creating additional closed field lines and a second dead zone. The torque is written as

$\gamma=3$0

where $\gamma=3$1 is the fraction of open donor field lines, ranging from $\gamma=3$2 for a non-magnetic WD to $\gamma=3$3 for complete suppression (Belloni et al., 2019). The same factor modifies the critical donor mass for MB disruption and the donor bloating factor (Belloni et al., 2019). Population synthesis showed that this reduced MB model reproduces the filled period gap of polars, their lower mass-transfer rates, and their space densities, providing evidence that WD–donor magnetic coupling suppresses MB in highly magnetic CVs (Belloni et al., 2019).

6. Applications and observational tests across stellar and binary populations

The strongest support for any MBM comes from the observables it reproduces. In single stars, MESA models with the Matt et al. and Garraffo et al. prescriptions reproduce many features of open-cluster rotation sequences, though both tend to over-brake sub-solar masses at ages $\gamma=3$4 Gyr (Gossage et al., 2020). Fully convective M-dwarf modeling with physically motivated field and wind scalings reproduces cluster trajectories and field-star spin periods while tracking spin-dependent wind mass losses, Alfvén radii, and surface magnetic fields (Sarkar et al., 2023).

In detached eclipsing binaries, the nearly flat period distribution from $\gamma=3$5 d down to contact strongly supports saturated MB and argues against classical Skumanich braking in tidally locked low-mass binaries (El-Badry et al., 2022). This result has direct implications for progenitors of CVs and LMXBs, because it suggests longer detached lifetimes and weaker AML than assumed in classical binary-evolution calculations (El-Badry et al., 2022).

In CVs, the choice of MBM “fundamentally determines whether CV systems develop the characteristic period gap” (Zhou et al., 5 Feb 2026). Systematic MESA calculations comparing Skumanich, Matt, RM12, intermediate, and convection-boosted laws found that the intermediate prescription provides optimal consistency with observations of nonmagnetic CVs, reproducing both gap location and donor properties, whereas Matt12 and RM12 maintain weak AML and fail to form a gap (Zhou et al., 5 Feb 2026). Another systematic comparison of standard MB, CARB, $\gamma=3$6-boosted, and SBD models found that CARB and $\gamma=3$7-boosted are too strong for CVs, while SBD can reproduce some features better than the standard model but worsens others, especially by overproducing period bouncers and, in population synthesis, by filling the gap with post-bounce systems (Tang et al., 24 May 2026).

A different empirical line of attack comes from direct period derivatives of CVs. A compilation of 52 systems argued that the “required MBM $\gamma=3$8 predictions all fail,” with 44 per cent of the sample showing positive $\gamma=3$9, negative-$\tau_{\rm MB}$00 systems being more negative than required by MBM by an average factor of 110, chaotic variations in some systems, and unmodeled nova-induced period changes dominating long-term evolution in many novae (Schaefer, 2024). That work concluded that the AML recipe of standard CV MBM is wrong by orders of magnitude (Schaefer, 2024). This is a strong critique of any single, universal CV AML track, though it addresses a different level of modeling than prescriptions aimed at reproducing ensemble properties.

In LMXBs, MBM strongly affects whether systems become persistent or transient, whether they can produce binary millisecond pulsars near $\tau_{\rm MB}$01 d, and how easily they reach the ultracompact regime. Detailed calculations found that the standard Skumanich law underpredicts mass-transfer rates by about an order of magnitude and suffers a fine-tuning problem for UCXB formation (Van et al., 2018, Echeveste et al., 2024). Improved prescriptions include MB2, MB3, and CARB. In the Van et al. parameterization,

$\tau_{\rm MB}$02

with MB0: $\tau_{\rm MB}$03, MB2: $\tau_{\rm MB}$04, and MB3: $\tau_{\rm MB}$05 (Echeveste et al., 2024). MB3 is the strongest and broadens the UCXB progenitor window from 0.2 d in MB0 to 1.69 d for a $\tau_{\rm MB}$06 donor, largely eliminating the fine-tuning problem in that grid (Echeveste et al., 2024). CARB improves over MB0 while remaining less extreme (Echeveste et al., 2024). In earlier work, the “intermediate MB” prescription,

$\tau_{\rm MB}$07

was identified as a practical MBM for LMXBs because it reproduces most persistent systems and many transients while yielding neutron-star masses consistent with observations (Van et al., 2018).

Black-hole LMXBs yield a contrasting constraint. Population synthesis plus detailed evolution showed that relatively low-efficiency MB laws, specifically RM12 and RVJ83, agree better with observed BH masses, companion masses, temperatures, orbital periods, and mean accretion rates than stronger prescriptions, and predict that only a very small fraction ($\tau_{\rm MB}$08) of BH LMXBs become UCXBs, consistent with the absence of observed BH UCXBs (Deng et al., 2024). This suggests that MB may not follow a single unified law across NS and BH systems (Deng et al., 2024).

Finally, orbital-period derivatives provide a direct test in LMXBs as well. Detailed MESA modeling of six NS and three BH systems found that the CARB prescription can reproduce the observed donor masses, orbital periods, and $\tau_{\rm MB}$09 values of four NS LMXBs and one BH LMXB, while the standard MB model struggles to explain rapid orbital shrinkage or expansion (Fan et al., 2024). This supports CARB in some accreting systems even though CARB appears too strong for nonmagnetic CVs (Fan et al., 2024, Tang et al., 24 May 2026).

7. Controversies, limitations, and prospects

A central controversy is whether a universal MBM exists at all. Detached binaries and single-star rotation strongly favor saturated braking with $\tau_{\rm MB}$10 or $\tau_{\rm MB}$11 in the rapid-rotation regime (El-Badry et al., 2022). Nonmagnetic CVs appear to require substantially stronger pre-gap braking to produce the observed period gap under the disrupted-MB scenario (Zhou et al., 5 Feb 2026). NS LMXBs often require stronger-than-Skumanich braking or additional wind/convection boosts to match persistent systems and UCXB formation (Van et al., 2018, Echeveste et al., 2024). BH LMXBs favor weaker, saturated MB (Deng et al., 2024). This suggests either that the effective torque is system-dependent or that additional AML channels beyond MB are being absorbed into different MBM calibrations.

Another controversy concerns the fully convective boundary. Classical CV theory ties the period gap to abrupt MB disruption at full convection. Yet detached eclipsing binaries show no significant difference between binaries containing fully and partially convective stars (El-Badry et al., 2022). SBD-type models attempt to reconcile these facts by allowing MB to persist in fully convective stars but with a strong reduction rather than a complete switch-off (Tang et al., 24 May 2026, Dodon et al., 15 May 2026). The double-dynamo formalism offers another reconciliation by distinguishing the disappearance of an interface dynamo from the survival of a convective dynamo (Sarkar et al., 2024).

A third uncertainty is the treatment of $\tau_{\rm MB}$12. Different prescriptions—empirical mass fits, local MLT evaluations, and global envelope integrals—can shift MB strengths and period-gap properties significantly (Sarkar et al., 2024, Tang et al., 24 May 2026, Dodon et al., 15 May 2026). This is not a minor calibration issue: in some SBD implementations, the sign of the $\tau_{\rm MB}$13 exponent in the saturated regime determines whether a spike in $\tau_{\rm MB}$14 weakens or strengthens MB (Dodon et al., 15 May 2026).

There are also domain-specific limitations. MBMs alone do not solve the lithium problem in rotating stellar models (Gossage et al., 2020). Reduced-MB models for polars depend on simplified dipole geometries and interpolated open-flux fractions (Belloni et al., 2019). Irradiation-driven boost models introduce poorly constrained efficiencies for accretion luminosity, irradiation absorption, and wind launching (Dodon et al., 15 May 2026). CV period-derivative studies argue that nova-induced discontinuities and other unknown mechanisms may dominate secular evolution in many systems (Schaefer, 2024).

Taken together, the literature suggests that MBM is evolving from a single empirical law into a layered framework. At minimum, a contemporary MBM must specify the torque’s dependence on rotation, saturation, convective turnover time, magnetic topology, and structural state, and it must state clearly whether it is intended for isolated stars, detached binaries, interacting binaries, or some unified subset of these. The available results support saturated braking in rapidly rotating detached systems, indicate that additional structural or environmental effects are required in accreting binaries, and leave open whether those effects should be described as modified magnetic braking or as distinct AML channels (El-Badry et al., 2022, Sarkar et al., 2024, Deng et al., 2024).