--- title: Magnetic Braking Model (MBM) url: https://www.emergentmind.com/topics/magnetic-braking-model-mbm type: topic --- # Magnetic Braking Model (MBM) 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 [2011.02470] [2401.04090] [2208.05488]. ## 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 [2011.02470] [1812.08239]. 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 [2011.02470] [2309.06893]. 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 [2011.02470] [2404.16185]. 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 [2401.04090] [2605.15897]. A fourth is magnetic geometry: some models weaken the torque for complex multipolar fields or for field-line closure in dead zones [2011.02470] [1910.06333]. 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 \(4\) [1910.06333] [1812.08239] [2404.16185]. In detached and tidally locked low-mass binaries, this unsaturated scaling predicts strong orbital shrinkage at short periods and a short-period period distribution \({\rm d}N/{\rm d}P_{\rm orb}\propto P_{\rm orb}^{7/3}\) [2208.05488]. 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” [2208.05488]. 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 \(\dot{J}\propto \Omega^3\) to \(\dot{J}\propto \Omega\) [2011.02470] [2208.05488]. In the Matt et al. (2015) implementation used in MESA, the unsaturated torque is \[ \dot{J} = -\mathcal{T}_0 \left(\frac{\tau_c}{\tau_{c,\odot}}\right)^p \left(\frac{\Omega}{\Omega_\odot}\right)^{p+1}, \] while the saturated torque is \[ \dot{J} = -\mathcal{T}_0\,\chi^p \left(\frac{\Omega}{\Omega_\odot}\right), \] with \(\chi=14\), \(m=0.22\), and a calibrated \(p=2.6\) in the Gossage et al. implementation [2011.02470]. The structural prefactor \(\mathcal{T}_0\) scales as \[ \mathcal{T}_0 = K \left(\frac{R}{R_\odot}\right)^{3.1}\left(\frac{M}{M_\odot}\right)^{0.5}\gamma^{-2m}, \] with \(K=1.4\times10^{30}\,{\rm erg}\) [2011.02470]. 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 \(\dot{J}\propto P_{\rm orb}^{-1}\), as opposed to \(\dot{J} \propto P_{\rm orb}^{-3}\) in the standard Skumanich law” [2208.05488]. 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 [2208.05488]. 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_o\), because it links rotation and convection in dynamo theory. In stellar-evolution implementations, \(\tau_c\) may be evaluated from local convective properties, for example \[ \tau_c(r)=\alpha_{\rm MLT}\frac{H_P(r)}{v_c(r)}, \] evaluated near the base of the convection zone [2011.02470], or in local and global-envelope forms derived directly from the stellar structure [2605.15897]. The Garraffo et al. (2018) prescription adds explicit magnetic-topology dependence. Its base dipolar torque is \[ \dot{J}_{\rm dip}=c\,\Omega^3\tau_c, \] and the full torque is \[ \dot{J}=\dot{J}_{\rm dip}\,Q_J(n), \] with a simplified topology factor \[ Q_J(n)=4.05\,e^{-1.4n}, \] and magnetic complexity \[ n(R_o)=\frac{a}{R_o}+1+bR_o. \] In the calibrated MESA implementation, \(a=0.03\), \(b=0.5\), and \(c=3\times10^{41}\,{\rm g\,cm^2}\) [2011.02470]. At low \(R_o\), \(n\) is large, fields are complex, and the torque is suppressed; at intermediate \(R_o\), \(n\approx 1\) and Skumanich-like spin-down is recovered; at high \(R_o\), complexity rises again and braking weakens [2011.02470]. 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 [2011.02470]. 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 \(t\gtrsim 600\) Myr to 1 Gyr [2011.02470]. Convective turnover time itself is also uncertain. Several prescriptions use powers of \(\tau_c\): Matt et al. (2015) scales as \(\tau_c^2\), Garraffo et al. (2018) as \(\tau_c\), and Van et al. (2018) as \(\tau_c^{8/3}\) [2401.04090]. Because empirical and MLT-based \(\tau_c\) values can differ by factors of several at low masses, MB strengths can shift by factors of several as well [2401.04090]. 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 \(M_\ast\lesssim 1.4\,M_\odot\) can host two distinct \(\alpha-\Omega\) 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 [2401.04090]. The total torque is written schematically as \[ \dot{J}_{\rm MB}=\dot{J}_{\rm BL}+\dot{J}_{\rm CE}, \] where \(\dot{J}_{\rm BL}\) vanishes when the radiative core disappears, leaving only \(\dot{J}_{\rm CE}\) in fully convective stars [2401.04090]. 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 \(\sim 6\) hr to \(\sim 3\) hr. At \(M_2\approx 0.25\,M_\odot\) the donor becomes fully convective, the boundary layer disappears, strong MB ceases abruptly, and the system detaches near \(P_{\rm orb}\simeq 3.18\pm0.04\) hr. It later re-establishes contact near \(P_{\rm orb}\simeq 2.15\pm0.03\) hr under gravitational radiation plus weaker convective-dynamo MB [2401.04090]. Below the gap, the convective-dynamo torque increases with decreasing mass and shifts the period minimum from the GR-only value \(\sim 65\) min to \(\sim 82\) min, matching the observed spike at \(80\)–\(86\) min [2401.04090]. 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 \(\dot{J}_{\rm MB}\propto \Omega\) or \(P_{\rm spin}^{-1}\) for fast rotators and \(\dot{J}_{\rm MB}\propto \Omega^3\) or \(P_{\rm spin}^{-3}\) for slow rotators, with free parameters \(f_1\) and \(f_2\); an example fit uses \(f_1=1.5\), \(f_2=2\) [2401.04090]. This is consistent with the dedicated FCMD study, which models fully convective M dwarfs with a torque approximately scaling as \(P_{\rm spin}^{-1}\) at shorter periods before transitioning into a Skumanich-type \(P_{\rm spin}^{-3}\) torque, and includes a parametrized reduction of wind mass loss due to wind entrapment in dead zones [2309.06893]. 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 [2309.06893]. 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 \[ \dot{J}_{\rm MB,SAT} = -C\left(\frac{R}{R_\odot}\right)^{1/2}\left(\frac{M}{M_\odot}\right)^{-1/2} \begin{cases} \left(\dfrac{P_{\rm rot}}{1\,{\rm d}}\right)^{-3}, & P_{\rm rot}\ge P_{\rm sat},\\[0.8em] \left(\dfrac{P_{\rm sat}}{1\,{\rm d}}\right)^{-2}\left(\dfrac{P_{\rm rot}}{1\,{\rm d}}\right)^{-1}, & P_{\rm rot}