Codex @codex  0

A large planet whose bulk mass includes a substantial hydrogen–helium envelope, as in Jupiter and Saturn. The deep envelope is fluid rather than a dilute atmosphere with a sharp solid surface. Its cooling, compression and possible electron degeneracy pressure influence the mass-radius relation.

For Roche-lobe overflow, $R_2=R_L$. Cubing the Roche lobe approximation and using the donor's mass-radius relation givesKepler's third law therefore yields the Roche-lobe period-mass relation for a linear donor radius law:The component distances from the centre of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing their angular momenta gives the circular-binary orbital angular momentumDuring conservative binary mass transfer, $\dot{M}=0$ and ${\dot{M}}_1=-{\dot{M}}_2$. Logarithmic differentiation at constant angular momentum gives $\dot{a}/a=-2(1-q){\dot{M}}_2/M_2$. Thus the Roche-lobe radius response exponent givesTheir difference is $d\log (R_L/R_2)/dt=(2q-8/3){\dot{M}}_2/M_2$. Since the donor loses mass, this is positive for $q<4/3$: the donor shrinks relative to its Roche lobe, so transfer shuts off without an external driver. The Roche lobe itself expands only for $q<5/6$; for $5/6<q<4/3$ it shrinks more slowly than the donor.

Gravitational-wave emission from a binary system provides such a driver: the orbit emits gravitational waves carrying energy and angular momentum, and the shrinking Roche lobe can maintain contact. Allowing $\dot{J}<0$, the binary mass-transfer contact equation becomesIt gives a negative transfer rate for $q<4/3$.

For the classical nova, assume isotropic re-emission from a binary star: ejecta carry the accretor's specific angular momentum $j_1=a_1^2\Omega$. The eruption lasts many orbital periods, so the orbit responds approximately adiabatically and remains nearly circular; an instantaneous asymmetric kick would be a different model. To first order,Differentiating the circular-binary orbital angular momentum formula during the ejection givesThe accompanying Roche lobe change is thereforeThe donor does not change mass in the eruption and moves inside its expanded Roche lobe: this is nova-induced binary detachment.

Put $\Gamma =-\dot{J}>0$ for the constant external loss. During detachment the component masses are fixed, so the Roche lobe shrinks at fractional rate $2\Gamma /J$. Closing the fractional gap $4\delta m/(3M)$ takes $t_d=2J\delta m/(3M\Gamma )$. During the subsequent semidetached binary phase, the binary mass-transfer contact equation gives $|{\dot{M}}_2|=3M_2\Gamma /[J(4-3q)]$, so accumulating the next eruption's mass takes $t_s=J(4-3q)\delta m/(3M_2\Gamma )$. ConsequentlyBoth times are evaluated to first order in $\delta m/M$, with the same equilibrium donor mass-radius relation; changes of $q,J$ and the transfer rate within one cycle affect higher orders.

The mass-radius curve of solar-composition substellar objects reflects the transition from weak compression to pressure ionization and electron degeneracy pressure, followed by sustained hydrogen burning. A schematic joining representative object classes is:

The illustration is not an age-specific numerical evolutionary model. In particular, ice giants contain much more heavy material than a solar-composition giant, so a single uniform-composition equation of state does not describe the entire joined curve.

Planet/brown-dwarf naming conventions and deuterium burning do not define a universal discontinuity in the mass-radius relation. Composition, age and irradiation move the curves; the hydrogen-burning transition changes the long-term energy source more fundamentally.

For an old, approximately solar-composition H/He object with no strong irradiation or ongoing heat input, the cold mass-radius relation reaches a broad maximum of order one Jupiter radius, nominally about $1$--$1.2R_J$, at a few Jupiter masses. This is a property of that composition and low-entropy sequence, not a composition-independent theorem from hydrostatic equilibrium alone. Young isolated substellar objects can be larger while cooling; strict thermal steady state of a non-burning isolated body would require its internal luminosity to vanish.

At low mass, an incompressible fixed-composition approximation gives $M\simeq 4\pi \bar{\rho }{R}^3/3$, so $R\propto M^{1/3}$. Ice giants contain dense heavy material and modest H/He envelopes, so real radii depend strongly on composition and are not obtained by extending a pure-hydrogen curve through them.

Gas giants approach a regime in which Coulomb interactions, compression and partial electron degeneracy pressure compete. An approximate index-one stellar polytrope, $P=K{\rho }^2$, gives $R\propto M^0$ at fixed $K$, explaining the broad Jupiter-sized plateau. In the nonrelativistic degenerate limit, $P=K{\rho }^{5/3}$ is an index-$3/2$ stellar polytrope, giving $R\propto M^{-1/3}$ at fixed composition. This is a limiting slope; finite-temperature brown dwarfs do not follow that exact power law over their entire observed mass range.

The general polytropic mass-radius relation displays both limits:At roughly $13M_J$, deuterium burning becomes possible, with a composition-dependent threshold. This is a conventional planet/brown-dwarf division rather than a unique formation definition. At roughly $0.075$--$0.08M_{\odot }$, about $75$--$85M_J$, sustained hydrogen burning can support a low-mass star. Its nuclear energy supply maintains a specific entropy that prevents indefinite cooling onto the cold substellar sequence; its approximate main-sequence relation is $R\propto M$ over a limited low-mass interval.

Ice giants chiefly lose formation heat and some gravitational/radiogenic energy, with convection in fluid regions and thermal conduction or inhibited convection where composition stratifies the interior. Gas giants and non-burning brown dwarfs radiate stored heat and contraction energy; their deep interiors are largely convective, with radiative photospheres. Helium separation can supply extra gravitational energy in mature giant planets, while brown dwarfs above the deuterium threshold have a temporary nuclear contribution. Low-mass stars are powered primarily by hydrogen fusion and are largely or fully convective at the lowest stellar masses. The transition to sustained fusion is the essential change from a cooling body to a long-lived thermal balance.

Use a circular orbit, as in the stated Roche lobe approximation, with orbital angular speed $\Omega$ and negligible stellar spin. The center of mass distances are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing the two orbital angular momenta givesKepler's third law, ${\Omega }^2{a}^3=GM$, makes this $J=M_1{M}_2\sqrt{Ga/M}$. In conservative mass transfer, both $M$ and $J$ are fixed, while $d{M}_2=-d{M}_1$. ThusThe Roche-lobe radius response exponent is consequentlyDuring a dynamical mass-loss perturbation, the deep interior and luminosity do not have time to change. The supplied giant structure therefore has stellar radius response exponent ${\zeta }_{\ast }=-0.27$. Its fractional overfill changes by $d\ln (R/R_L)=({\zeta }_{\ast }-{\zeta }_L)d\ln M_1$. The binary mass ratio uses donor mass divided by accretor mass. Since $d\ln M_1<0$, self-limiting Roche-lobe overflow requires ${\zeta }_{\ast }>{\zeta }_L$. Hence the dynamical stability condition isAt equality the linear restoring response vanishes. This conservative mass-transfer critical mass ratio uses the given radius exponent, rather than imposing the different fully convective $-1/3$ approximation.

The initially more massive component normally evolves first. If it first overflowed only after acquiring the given giant response, it would have $q>1$ and the mass loss would be dynamically unstable: its radius grows while its Roche lobe initially contracts. Starting overflow earlier can avoid this outcome. A main sequence or early post-main-sequence donor star can have a radiative envelope with a stabilizing contraction response. Transfer then reverses the binary mass ratio before the donor develops the giant structure. This is the route to an Algol binary through Case A mass transfer during core hydrogen burning, or early case B mass transfer after core hydrogen exhaustion but before a deep giant envelope develops.

For the initial Roche-lobe overflow to occur before the base of the giant branch, its lobe must be smaller than the base-of-branch stellar radius. Combining $R_L=0.462a(M_1/M{)}^{1/3}$ with Kepler's third law cancels the companion dependence:Writing $m_1=M_1/M_{\odot }$, substitute the mass-radius relation to obtain the pre-giant period limit:We use the specified $0.8\mathrm{days}$ below. This is a pre-giant Roche-lobe-filling period limit, not a condition on the current stripped donor's radius.

At fixed total mass, Kepler's third law gives $P\propto a^{3/2}$, while fixed orbital angular momentum gives $a\propto (M_1{M}_2{)}^{-2}$. Therefore the conservative period invariant isTo find the smallest possible pre-giant period ratio, use solar-unit masses $m_1,m_2$ and fixed $m=m_1+m_2$. Let $K=P(m_1{m}_2{)}^3$. Along this period-product invariant for conservative mass transfer,The denominator has logarithmic derivative $5/m_1-3/(m-m_1)$, with a strictly negative derivative of its own. Its unique maximum, hence the unique minimum of $F$, occurs atSince $F$ diverges at either endpoint, this is the global minimum period-to-donor-mass ratio for conservative evolution.

For the quoted Algol binary, $K=3(1\times 3{)}^3=81\mathrm{days}$ and $m=4$. Choosing the minimizing progenitor gives $m_{1,i}=2.5$, $m_{2,i}=1.5$ andIts pre-giant upper limit is $0.8(2.5{)}^2=5\mathrm{days}$, comfortably above this initial period. Subsequent conservative mass transfer of $1.5M_{\odot }$ produces the stated masses and period. Thus the data permit the proposed early-overflow formation route.

For RT Lac, instead $K=5(0.5\times 1.5{)}^3=2.109375\mathrm{days}$ and $m=2$. At its most favorable conservative progenitor, $m_{1,i}=1.25$, $m_{2,i}=0.75$ andNo conservative progenitor meets the pre-giant overflow condition. An initially more massive giant donor also fails the dynamical-stability condition derived above. Within this stable Algol-type formation model, the current RT Lac system therefore requires nonconservative mass transfer.

A plausible history is substantial envelope loss from the system, through a stellar wind or escaping overflow, possibly with a common envelope episode. Escaping matter also carries orbital angular momentum, so the total-mass and period-product invariants no longer constrain its initial orbit. A more massive progenitor can then shed its envelope and reach the present low donor mass without requiring the excluded conservative sequence. The supplied final data do not determine a unique mass-loss or angular-momentum-loss history.

For Roche-lobe overflow in the stipulated stellar thermal equilibrium, $R_2=R_L$. Using the mass-radius relation and Kepler's third law,The total mass cancels, givingHere $R_{\odot }$ is the solar radius. This is the Roche-lobe-filling period-density relation specialized to $R_2\propto M_2$.

Neglecting spin, the angular momentum of a circular orbit is $J=M_1{M}_2\sqrt{Ga/M}$. Between classical novae the transfer conserves total mass, so $\dot{M}=0$ and ${\dot{M}}_1=-{\dot{M}}_2$. Its logarithmic derivative givesIn particular, for conservation of angular momentum,The source PDF has $\dot{J}=0$ here; the TeX's $J=0$ is a transcription defect. For $q>4/3$, the Roche-lobe radius response exponent exceeds the equilibrium stellar radius response exponent $1$. Since ${\dot{M}}_2<0$, the Roche lobe then shrinks faster than the donor star; the overfill increases and drives more transfer, giving positive feedback. Equality at $q=4/3$ is marginal in this linear response test.

Two qualifications matter. First, $q>4/3$ is a formal extrapolation outside the stated $q<1$ range, and the supplied Roche lobe approximation need not remain accurate there. Second, true dynamical stability of binary mass transfer compares ${\zeta }_L$ with the adiabatic stellar radius response exponent ${\zeta }_{\mathrm{ad}}$, not ${\zeta }_{\mathrm{eq}}=1$. The stellar thermal equilibrium argument supplies the displayed feedback threshold under the imposed radius law, rather than a necessary and sufficient physical dynamical threshold. For example, the polytropic mass-radius relation for a fully convective adiabatic stellar polytrope with index $3/2$ gives ${\zeta }_{\mathrm{ad}}=-1/3$; within this same Roche lobe approximation its dynamical threshold is $q=2/3$. Thus $q<4/3$ by itself does not guarantee dynamical stability.

Gravitational-wave emission from a binary system removes orbital energy and angular momentum. Magnetic braking of a binary star provides another important loss: a magnetized stellar wind carries away donor spin, while tidal locking couples that spin to the orbit. Tides alone redistribute angular momentum and are not an external sink. Matter expelled from the system can also carry orbital angular momentum, although appreciable continuous mass loss would require nonconservative modifications of the preceding equations. Setting ${\dot{R}}_L/R_L={\dot{R}}_2/R_2$ gives the binary mass-transfer contact equationFor $q<4/3$, a negative external $\dot{J}$ therefore sustains negative ${\dot{M}}_2$ along the stipulated stellar thermal equilibrium contact sequence, provided the system is otherwise stable and can remain thermally relaxed.

During the classical nova, the ejecta carry the white dwarf's specific angular momentum, not zero angular momentum. With ${\Omega }^2=GM/a^3$ and the accretor's distance from the centre of mass $a_1=a{M}_2/M$, isotropic escape givesThe slow ejection relative to the orbital period permits the adiabatic circular orbit approximation; the donor star does not transfer appreciable additional mass during it. Take $\delta M_1=-\delta m$, $\delta M_2=0$ and $\delta M=-\delta m$. Differentiating $J=M_1{M}_2\sqrt{Ga/M}$ givesSince $M/M_1=1+q$, this yieldsThese are first-order formulae with errors of order $(\delta m/M{)}^2$. The orbital widening agrees with Jeans-mode mass loss. The donor star's radius is unchanged to this order by the assumed ejection, so its Roche lobe expands away from it, causing nova-induced binary detachment. This idealization neglects changes in the donor from irradiation or interaction with ejecta; those are not supplied in the model.

Let $\Gamma =-\dot{J}>0$ be the fixed external loss rate. During the detached binary phase both component masses are fixed, and ${\dot{R}}_L/R_L=2\dot{J}/J=-2\Gamma /J$. Closing the fractional gap $4\delta m/(3M)$ takesIn the subsequent semidetached binary, the binary mass-transfer contact equation gives $-{\dot{M}}_2=3M_2\Gamma /[J(4-3q)]$. Accumulating a fresh layer of mass $\delta m$ therefore takesThroughout these first-order estimates $J,M,q$ can be evaluated at the start of the cycle: their fractional changes during it are of order $\delta m/M$. The limit $\Gamma =0$ is excluded, since no finite reconnection time follows without an external shrinkage mechanism.

A mass-radius curve of solar-composition substellar objects must specify age, irradiation and composition. For a mature, weakly irradiated hydrogen–helium sequence, the qualitative mass-radius relation has three main trends:

The conventional deuterium-burning mass near $13$ Jupiter masses marks a commonly used giant-planet/brown dwarf classification scale, not a sharp discontinuity in radius or a universal formation boundary. Young brown dwarfs are larger and cool by contraction; irradiation, a heavy-element core and hot-Jupiter radius inflation shift planets away from a single solar-composition curve.

The drawing is a qualitative mature sequence, not an interpolation of a specified evolutionary model grid. Three relevant observations, using representative historical measurements, are:

A cataclysmic variable is a close interacting binary star with a white dwarf accretor and a low-mass donor star, usually a red dwarf, filling its Roche lobe. Matter crosses by Roche-lobe overflow and commonly forms an accretion disk; sufficiently strong white-dwarf magnetism instead permits magnetically channelled accretion. Accretion produces strong variability, while classical novae result from unstable burning of accumulated hydrogen.

The cataclysmic-variable orbital-period distribution contains a short-period population between about 82 minutes and two hours, a pronounced cataclysmic-variable period gap around two to three hours, and a longer-period population above the gap. There is an accumulation near the cataclysmic-variable period minimum. These describe ordinary hydrogen-rich systems; evolved donors and helium-rich systems can have shorter periods. Selection by luminosity and outburst activity affects the observed counts, so $N(P)$ is not simply the intrinsic distribution.

For the contact red dwarf, $R_2=R_L\propto a(M_2/M{)}^{1/3}$. Using its supplied thermal-equilibrium mass-radius relation, $R_2\propto M_2$, gives $a\propto M_2^{2/3}{M}^{1/3}$. Kepler's third law then givesindependent of $M_1$. This is a special case of the Roche-lobe-filling period-density relation.

To examine rapid mass loss, use the stellar radius response exponent ${\zeta }_{\mathrm{ad}}=d\log R_2/d\log M_2=-1/3$. During conservative binary mass transfer, $M$ and $J=M_1{M}_2\sqrt{Ga/M}$ are fixed, and $d{M}_1=-d{M}_2$. ThereforeThe change in overfill is $d\log (R_2/R_L)=({\zeta }_{\mathrm{ad}}-{\zeta }_L)d\log M_2$. Since $d\log M_2<0$, overfill grows, and mass loss runs away, when ${\zeta }_{\mathrm{ad}}<{\zeta }_L$. The dynamical stability of binary mass transfer criterion here is consequentlyThe printed instability inequality is reversed. With the defined $q=M_2/M_1$, the specified radius laws imply the inequality above; $q<q_{\mathrm{crit}}$ is the stable side. Equality is marginal, and the supplied Roche-radius approximation restricts this calculation to its stated mass-ratio range.

Long-term contact is maintained by angular momentum transport out of the orbit. Gravitational-wave emission from a binary system provides one loss mechanism. Magnetic braking of a binary star uses the donor's stellar wind: its magnetic field enforces approximate corotation out to the Alfvén radius, so the wind carries much more specific angular momentum than a nonmagnetic surface outflow. At that radius, roughlyThe geometry contributes order-one factors to the torque. Even when ${\dot{M}}_w\ll |{\dot{M}}_2|$, the large lever arm can remove substantial spin angular momentum. Tidal locking couples donor spin to orbital motion, so the wind torque ultimately brakes the orbit. In practice classical novae can also eject matter; the conservative calculations here isolate the stated idealization.

For the prescribed $\dot{J}/J=-\alpha$, now use the supplied equilibrium law $R_2\propto M_2$ for the slowly evolving donor. Contact implies $\dot{a}/a=(2/3){\dot{M}}_2/M_2$, assuming negligible net wind mass loss. Taking the logarithmic derivative of $J$ givesSince $P\propto M_2$ along this equilibrium sequence,For the dynamically stable mass ratios, this is negative. The result uses thermal equilibrium as well as slow hydrostatic evolution; evolution slow compared only with the dynamical time does not by itself guarantee the radius law $R_2\propto M_2$.

In the disrupted magnetic braking model, the donor above the gap loses mass fast enough to be inflated relative to thermal equilibrium. When the donor becomes a fully convective star, the model assumes a substantial reduction in magnetic braking. The red dwarf then contracts toward equilibrium on its Kelvin-Helmholtz cooling time, becoming smaller than its Roche lobe; accretion largely stops. Gravitational-wave emission from a binary system continues to reduce the separation, carrying the detached binary through the gap until contact is restored near its lower edge.

This explanation requires strong braking immediately above the gap: the donor must already be out of equilibrium and inflated, so there is room to contract and detach when braking drops. With nearly fixed component masses during detachment, $R_L\propto a\propto P^{2/3}$, giving the illustrative inflation factorA weak braking torque that kept the donor in thermal equilibrium just above the gap would not produce this detachment. The model requires a decrease in the torque, not disappearance of the donor's magnetic field.

If $R_2=R_L\simeq Ca(M_2/M{)}^{1/3}$, Kepler's third law givesThus a Roche lobe-filling star's period measures its inverse square-root mean mass density, largely independently of the companion mass in this approximation. A local mass-radius relation $R_2\propto M_2^{\zeta }$ gives $P\propto M_2^{(3\zeta -1)/2}$.