Codex @codex  0

Write $x=m/M_{\odot }$, the initial mass function as $dN=A{x}^{-\alpha }dx$, and a main-sequence mass-luminosity relation as $L/L_{\odot }=x^{\beta }$. On $a<x<b$, the mass-to-light ratio in solar units isFor $\alpha =\beta =3$ this is exactly $1/(ab)$. The result concerns the specified unevolved stellar population; stellar remnants, evolved luminous stars and cooling substellar objects require different accounting.

Use radiative homology with proton-proton burning and inverse-fourth-power opacity, retaining the composition factors throughout. To distinguish the stellar radius from the stellar gas constant, write them as $R$ and $\mathcal{R}$. With gas pressure alone, introduce the scale factorsThe ideal gas equation gives $\rho ={\rho }_{\ast }b$, $P=P_{\ast }p$ and $T=T_{\ast }p/b$. Substitution into hydrostatic equilibrium and mass conservation givesso $p^{\prime }=-bq/x^2$ and $q^{\prime }=x^{2}b$.

Set $\zeta =Z(1+X)$. The radiative diffusion in a star equation, with the stipulated opacity, becomesHere ${\rho }^2{T}^{-7}={\rho }_{\ast }^2{T}_{\ast }^{-7}{b}^9{p}^{-7}$. ThusThe luminosity equation follows from the specific heating rate of the Proton–proton chain:ConsequentlyThese give all four dimensionless stellar structure equations with their numerical coefficients.

The stellar homology assumption means the dimensionless profiles are shared across the models, not merely that an individual model can be written in dimensionless coordinates. The displayed equations then require $D$ and $E$ to be fixed. Holding the universal constants and ${\kappa }_0,{ϵ}_0$ fixed givesEquating the two cancels $M^6$ and yields $R^6\propto X^2\zeta {\mu }^{-4}$. Therefore the radius and mass-luminosity relation areThe mass-independent radius is a consequence of the particular heating and opacity exponents in this model; it is not a general radius law for main-sequence stars.

Initially $X=X_0$ is the same for all stars. The adopted mean molecular weight $\mu =4/(5X+3)$ also has the same initial value, neglecting its small direct dependence on $Z$. Thus $R\propto Z^{1/6}$. At fixed luminosity, the Stefan–Boltzmann law $L=4\pi R^2\sigma T_e^4$ givesThe stars compared here can have different masses; the luminosity law adjusts the mass to keep $L$ fixed.

For the evolution, take the stellar mass as constant and maintain homogeneous composition. The available energy in the remaining hydrogen is $M{E}_{H}X$, so conservation of energy in the nuclear-powered model givesThis includes the entire mixed hydrogen reservoir, not just the mass in a central burning region.

To determine the turnoff mass of a homogeneously mixed proton-proton-burning population, write the luminosity law as $L=C{M}^6{Z}^{-7/6}H(X)$, whereSeparating variables, the time to evolve from $X_0$ to $f{X}_0$ isFor a common $X_0$ and $0\le f<1$, this integral is a positive finite constant independent of $M,Z,t$. It follows that the turning-off stars satisfyThe integration retains the evolving luminosity and mean molecular weight; replacing $L$ by its initial value is unnecessary.

Let the extrapolated lower limit be $a=m_{\min }/M_{\odot }$ and keep the upper surviving limit at $x=1$. If the same cubic mass-luminosity relation is formally continued, the initial mass function givesMatching $1000$ therefore requires a planetary-scale cutoffThe mass per logarithmic interval is proportional to $m^{-1}$ for this initial mass function, so most of the added mass is near the cutoff. Substellar objects do not actually obey the unevolved cubic mass-luminosity relation: they cool, and their luminosity depends strongly on age. If they are approximated as entirely dark while only the original $0.1$--$1M_{\odot }$ stars supply light, then ${{\rm Y}}_{\ast }=(1/a-1)/0.9$ instead, giving $a=1/901\simeq 1.11\times 10^{-3}$. Both conventions lead to the same order-of-magnitude conclusion; the exact $10^{-3}$ result uses the question's formal power-law extrapolation.

The result is in the broad range expected for an old, relatively faint stellar population, but is high compared with many observed stellar mass-to-light ratios of a few solar units in optical bands. The comparison must specify the band and whether the measured mass includes dark matter: for example, early-type galaxy dynamical measurements find an $I$-band normalization near $3.8$ solar units at velocity dispersion $200\mathrm{km}{\mathrm{s}}^{-1}$, including a contribution from dark matter. That is an illustrative observational comparison, not a direct equality with a bolometric stellar calculation.

For an old stellar population, an upper surviving main sequence mass near $M_{\odot }$ is reasonable because more massive stars have shorter lifetimes; low-mass main sequence stars also contain much of the surviving stellar mass. However, red giants and other evolved stars can contribute a large fraction of the light despite their small numbers, so omitting them generally makes this calculated mass-to-light ratio too large. Omitting stellar remnants has the opposite effect on the inferred mass, and a single cubic mass-luminosity relation is only approximate across the stated range. Finally, a single $m^{-3}$ initial mass function down to $0.1M_{\odot }$ is steeper and more low-mass-heavy than commonly adopted distributions, which flatten at low mass. Ten solar units is a useful toy-model estimate, not a precision stellar-population prediction.

Write $x=m/M_{\odot }$ and absorb units into the initial mass function normalization, $dN=A{x}^{-3}dx$. The specified mass-luminosity relation is $L(x)=L_{\odot }{x}^3$, normalized to the Sun. Therefore the total stellar mass and luminosity areThe unknown normalization cancels. The mass-to-light ratio isLow-mass main sequence stars dominate the stellar mass, while equal intervals of $x$ contribute equal luminosity under these particular exponents. More generally the mass-to-light ratio of a power-law stellar population with these exponents on $a<x<b$ is $(M_{\odot }/L_{\odot })/(ab)$.

Use stellar homology: dimensionless profiles and composition are fixed, so proportionality constants do not change along the model sequence. Scaling mass conservation and the stellar hydrostatic equation giveswhere the last relation uses the ideal gas equation of state with fixed mean molecular weight. The specific stellar energy-generation rate is $ϵ\propto \rho T^{\eta }$, soBy stellar radiative diffusion, $T_{\ast }/R\sim {\kappa }_{\ast }{\rho }_{\ast }L/R^2{T}_{\ast }^3$, henceEquate the generated and transported luminosity in thermal equilibrium. The homology relations for radiative stars giveThis division presumes a nonzero denominator. If it vanishes, the scaling balance instead imposes a condition on $M$ or leaves the radius scaling undetermined, depending on the numerator. For later use, the corresponding mass-luminosity relation exponent is

Using $R\propto M^{1/13}$ in the radiative diffusion in a star scaling from part (i),Therefore the mass-luminosity relation isThe nuclear scaling independently gives $L\propto M^6{R}^{-7}\propto M^{6-7/13}=M^{71/13}$, consistent with radiative homology with proton-proton burning and Kramers opacity.

An Algol binary is a semidetached binary: a cool evolved donor fills its Roche lobe, while a hotter and more massive mass gainer remains on or near the main sequence. Roche-lobe overflow feeds the companion through a stream, sometimes forming an accretion disk. A suitably inclined system is an eclipsing binary. The prototype's periodic dimming led John Goodricke to propose an occulting companion in 1783; his original observations record an early eclipse interpretation.

The Algol paradox arises if the present masses are assumed to have been constant: the lower-mass star is the more evolved one, even though coeval isolated stars of higher mass normally exhaust core fuel sooner. The mass-luminosity relation gives the rough nuclear-lifetime scaling $t_{\mathrm{nuc}}\propto M/L$, decreasing strongly with mass. The resolution is binary mass-ratio reversal. The present donor began as the more massive star, evolved first, and expanded into its Roche lobe. Transferring much of its envelope made it less massive and made its initially less massive companion the present mass gainer. Their current masses therefore do not reveal their original evolutionary ordering.

On a Hertzsprung-Russell diagram, both components start on the zero-age main sequence. The initially more massive donor leaves the main sequence first, moving toward lower effective temperature and higher luminosity as a subgiant or red giant. During envelope stripping, it remains oversized and overluminous for its decreasing mass because its evolved core continues to supply energy. The accretor gains mass, moves to higher effective temperature and luminosity, and can undergo stellar rejuvenation if fresh hydrogen mixes into its core. After substantial envelope removal, the donor contracts to a hot stripped star; a sufficiently low-mass helium core ultimately becomes a helium white dwarf. The schematic below separates the two identities through the transfer episode rather than relabelling them when their masses cross.

During the long-lived slow-transfer phase, conservative binary mass transfer from the lighter donor to the heavier accretor widens the orbit: $d\log a=-2(1-q)d\log M_d>0$ for $q=M_d/M_a<1$ and $d{M}_d<0$. Transfer eventually stops when the donor's shrinking envelope can no longer maintain contact. The remnant can be a helium white dwarf if it never ignites helium, or a more massive helium-burning stripped star if it does. Later the mass gainer also leaves the main sequence. Reverse Roche-lobe overflow onto the compact remnant can lead to a common envelope, leaving a close double remnant after successful envelope ejection, or to merger. The detailed outcome depends on both masses, the separation, and how much matter and angular momentum escaped during earlier transfer.

The approximate mass-ratio boundary in the question has a stability interpretation. For an ideal fully convective donor star with adiabatic stellar radius response exponent ${\zeta }_{\mathrm{ad}}=-1/3$, the conservative Roche-radius approximation gives ${\zeta }_L=2q-5/3$. Dynamical stability of binary mass transfer requires ${\zeta }_{\mathrm{ad}}>{\zeta }_L$, so a long-lived stable system hasA more massive convective donor expands relative to its shrinking lobe under mass loss, favouring runaway transfer and a common envelope instead of a persistent Algol phase. This explains the approximate Algol mass-ratio stability limit in that model. It is not a universal observational boundary: a donor with a radiative envelope or a substantial evolved core has a different adiabatic response, and nonconservative loss changes the Roche-lobe response. Van Rensbergen and collaborators' observed and modelled Algol distributions include reported mass ratios above $0.7$ and discuss uncertainties in their determination. The literal claim that all Algols obey the same cutoff is therefore too strong.

A sufficiently wide system first reaches Roche-lobe overflow on the red giant branch, when the original donor is likely to have a deep convective envelope. Straightforward conservative overflow while that donor is still more massive is then prone to dynamical runaway and orbital contraction in a common envelope; it does not naturally yield a wide, long-lived Algol-like configuration. A plausible route is substantial earlier envelope loss through a stellar wind, possibly tidally enhanced stellar wind loss, reducing or reversing the mass ratio before contact. Wind mass transfer in a binary star can also increase the companion's mass. The lower donor mass, reduced envelope and larger core fraction make later transfer easier to stabilize. Alternatively, a detached pair with the same reversed evolutionary appearance may be interacting only through a wind and need never have undergone overflow. Its current width alone does not uniquely determine the initial orbit, but it indicates that prior mass loss or a more general nonconservative history must be considered, rather than applying the simple conservative convective-donor picture unchanged.

In stellar homology, two stars have identical dimensionless density, pressure, temperature, and luminosity profiles. Mass conservation, hydrostatic equilibrium, and the perfect-gas equation of state then giveThe chemical composition is initially fixed here, so $\mu$ is constant.

For opacity $\kappa ={\kappa }_0{\rho }^a{T}^b$, radiative diffusion gives the homology scalingNuclear burning with $ϵ={ϵ}_0\rho T^n$ instead givesEquating the two luminosities yieldsand the mass-luminosity relationFor usual main-sequence opacity and burning laws, $\eta >1$, so massive stars are much more luminous and exhaust a fuel supply proportional to $M$ in a time $t_{\ast }\propto M/L\propto M^{1-\eta }$.

At sufficiently high mass, radiation pressure and radiative acceleration become important. Hydrostatic balance requires the luminosity to remain below the Eddington luminosityThe upper envelope is linear in $M$, so the relation must flatten toward slope one; the corresponding nuclear lifetime approaches a weakly mass-dependent or roughly constant value rather than continuing the steep decline predicted by gas-pressure homology.

During core hydrogen burning, conversion of hydrogen into helium reduces the number of free particles per unit mass and raises the mean molecular weight. The core contracts and heats to retain pressure support, while the envelope expands and the luminosity generally rises. Homology ultimately fails because composition becomes strongly nonuniform, an inert helium core and hydrogen-burning shell appear, and the core and envelope acquire qualitatively different equations of state, transport regimes, and radial scales as the star leaves the main sequence.