Quark mixing arises because the up-type and down-type quark Yukawa matrices need not be diagonal in the same weak interaction basis. In the Standard Model the three left-handed quark doublets transform as under , whereas and transform as and . With a Higgs doublet and its conjugate , the gauge-invariant Yukawa interaction is
The generation indices are matrix-multiplied. Electroweak symmetry breaking, , produces and .
Each complex fermion mass matrix admits a biunitary diagonalization,
where the diagonal entries are real nonnegative quark masses. Define and , and similarly for right-handed fields. The kinetic terms retain their canonical form because the transformations are unitary matrices. However, the weak charged current becomes
This unitary matrix is the Cabibbo-Kobayashi-Maskawa matrix. Its entries multiply transitions between up-type and down-type mass eigenstates. The minimal Standard Model has no analogous right-handed charged weak current.
Quark mixing does not introduce tree-level flavor-changing neutral currents. Photon, gluon and weak neutral current couplings are proportional to the identity within each fixed-charge quark sector, and . The single Higgs doublet also has diagonal neutral scalar field couplings after the same transformations that diagonalize the masses. At higher loop orders, flavor changes are possible. For example, an amplitude contains . Unitarity of the CKM matrix sets the sum of coefficients to zero, cancelling any mass-independent part. This is the GIM mechanism; mass splittings leave a suppressed, generally nonzero amplitude.
For generations a unitary matrix has real parameters: mixing angles and phases. Rephasing the up-type and down-type quark fields removes phases; a common baryon number phase leaves every matrix entry unchanged. Hence the physical numbers are
Two generations therefore give a real Cabibbo rotation,
with strangeness-changing charged-current transitions suppressed by . Three generations have three mixing angles and one irreducible weak CP phase. This counting assumes nondegenerate masses; degeneracies allow additional basis rotations.
Writing and , a standard exact parametrization is
Complex entries alone do not establish CP violation, because quark quantum field rephasings can change their phases. For nondegenerate quark masses the weak mixing sector conserves CP precisely when a basis makes real. A useful invariant test is the Jarlskog invariant,
A nonzero requires all three mixing angles and a nontrivial phase. If two same-charge quark masses coincide, the additional basis freedom removes the corresponding physical CP effect. The possible QCD theta angle is a separate CP parameter, not part of this weak-mixing count.
The Wolfenstein parametrization displays the hierarchy of quark mixing economically. Put , , and . Expanding the exact matrix gives
The observed pattern is predominantly diagonal, with - mixing largest and -, - mixing successively suppressed. The expansion implies . The truncated matrix must not be treated as an exactly unitary matrix.
Column orthogonality gives
The three complex numbers form a CKM unitarity triangle, with area . Dividing by puts its vertices at , and , where . The barred parameters include higher-order corrections to the Wolfenstein parameters. Comparing independent measurements of its sides and angles tests both quark mixing and the absence of additional flavor dynamics.
Semileptonic decay rates determine magnitudes of CKM entries through factors once hadronic matrix elements are controlled. Neutral kaon and neutral B-meson mixing constrain products of entries through weak box scattering amplitudes. CP asymmetries probe their invariant relative phases, with kaon mixing providing indirect CP violation and decay amplitudes providing direct CP violation. These processes complement tests of the unitarity constraints on rows and columns. The minimal three-generation Standard Model describes weak quark mixing by a unitary matrix with three angles and one physical CP phase; it does not predict the Yukawa couplings or their hierarchy.
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 and . With gas pressure alone, introduce the scale factors
The ideal gas equation gives , and . Substitution into hydrostatic equilibrium and mass conservation gives
so and .
Set . The radiative diffusion in a star equation, with the stipulated opacity, becomes
Here . Thus
The luminosity equation follows from the specific heating rate of the Proton–proton chain:
Consequently
These 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 and to be fixed. Holding the universal constants and fixed gives
Equating the two cancels and yields . Therefore the radius and mass-luminosity relation are
The 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 is the same for all stars. The adopted mean molecular weight also has the same initial value, neglecting its small direct dependence on . Thus . At fixed luminosity, the Stefan–Boltzmann law gives
The stars compared here can have different masses; the luminosity law adjusts the mass to keep fixed.
For the evolution, take the stellar mass as constant and maintain homogeneous composition. The available energy in the remaining hydrogen is , so conservation of energy in the nuclear-powered model gives
This 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 , where
Separating variables, the time to evolve from to is
For a common and , this integral is a positive finite constant independent of . It follows that the turning-off stars satisfy
The integration retains the evolving luminosity and mean molecular weight; replacing by its initial value is unnecessary.
For the Schwarzschild criterion, displace a small gas parcel upwards. Assume its composition is unchanged, it remains in pressure balance with its surroundings, and it moves quickly enough that heat exchange is negligible. The parcel obeys the adiabatic equation of state . Combining this with the ideal gas equation gives , so its adiabatic temperature gradient is .
Let the surrounding temperature gradient be . For an upward displacement , the parcel's temperature relative to its new surroundings is
If , this is negative: the parcel is cooler and, at equal pressure and composition, denser than its surroundings. It therefore sinks back; a downward displacement likewise gives a restoring buoyancy force. If the inequality is reversed, the displaced parcel continues to rise. The stability condition is thus
Equality is the neutral boundary. There is no composition-gradient term because the ideal-gas parcel and the background have the same mean molecular weight.
In the thin upper stellar atmosphere, take a plane-parallel approximation with constant , optical depth increasing inward, and zero overlying pressure . Hydrostatic equilibrium in optical depth gives . The ideal gas equation converts the stipulated opacity into
Define , so the given grey atmosphere relation is . Integrating from the surface gives
This proves the required pressure law. In particular times a function only of and .
For the grey-atmosphere convection onset with density-linear opacity, differentiate the two profiles to obtain the radiative temperature gradient:
It increases monotonically from zero, since , and first reaches the adiabatic temperature gradient at . Therefore
Below this point the radiative stratification is unstable and convection must transport some of the outward energy flux.
At the outer convection boundary, depends only on , so and for fixed . The matching constant in is . Consequently
This temperature form of the polytropic relation follows from an ideal monatomic adiabat; its is not the density-polytrope coefficient in Question 3.
For the deep radiative boundary beneath a power-law-opacity convective envelope, negligible envelope mass and absent local energy generation give and . Divide the radiative temperature equation by the pressure equation to get the stellar radiative temperature gradient
This is the gradient that would be needed if radiation alone carried , and it is the quantity to compare with even while the actual region is convective. The deeper Kramers' opacity law gives . Substituting yields
As rises inward this decreases, so the transition back to a stable radiative interior is at
Thus . Using the outer matching and the Stefan–Boltzmann law,
The mass and radius cancel. Hence, at fixed opacity coefficients and composition,
This uses the same adiabat through the convection zone and distinguishes its inner transition temperature from its outer onset temperature.
For a regular, spherical star, hydrostatic equilibrium gives . Multiply by and integrate to the surface. Integration by parts gives
The central boundary term vanishes for finite central pressure. Since , the right-hand side of the integrated hydrostatic equation is
The last expression is the gravitational energy: adding a thin mass shell to the interior mass contributes , counting each interacting pair once. Hence
With the specified zero surface pressure this proves the stellar virial theorem in the required form,
To derive the nonrelativistic electron-degeneracy pressure in a hydrogen-helium mixture, take fully ionized hydrogen and helium with hydrogen mass fraction , the electron number density is
Here helium nuclei have mass in the approximation used in the question. Filling the two spin states of each momentum mode up to the Fermi momentum gives
The isotropic momentum flux gives pressure . For nonrelativistic electrons , so the electron degeneracy pressure is
Substitute the composition-dependent and simplify powers of two:
Comparing gives the polytropic index . The result assumes complete electron degeneracy, , and negligible nuclear mass corrections.
Now adopt the stated approximate partially degenerate equation of state, keeping composition and hence fixed during contraction. Set and define the mass-averaged temperature . The stellar virial theorem from part (a) gives
For an index- stellar polytrope, the permitted gravitational-energy formula is . Also, using ,
Thus, with ,
The normalized density is a fixed function of for this homologous polytropic structure. Therefore its integral and fix , which is independent of and for a common composition. The mean here is mass-weighted, rather than a volume average.
For the mass-averaged temperature maximum in a partially degenerate polytrope, put and , both positive. As a function of the temperature is the concave quadratic
It reaches its maximum when , with
At large radius, contraction increases the temperature as in a thermally supported star. Once electron degeneracy provides a sufficiently large fraction of the support, further contraction can lower the thermal temperature. The formal zero-temperature radius is , with . The continuation to negative temperatures would be unphysical; the adopted sequence only describes its positive-temperature part.
The contraction-limited central temperature has the same mass scaling. Indeed, if the total polytropic pressure is , the approximate equation of state gives . The fixed density profile therefore fixes . A sufficiently small mass never reaches the central temperature required for sustained hydrogen burning, producing a hydrogen-burning minimum mass. Such an object becomes a cooling, degeneracy-supported brown dwarf instead of settling onto the main sequence. More massive contracting objects can ignite hydrogen before reaching this degeneracy-limited maximum. This conclusion uses the nonrelativistic, fixed-composition model; it does not require extrapolation into the relativistic regime.
Treat the stars as point masses for the orbital motion, neglecting their spin angular momenta. Let . Their distances from the centre of mass are and . The sum of their orbital angular momenta is
Thus the circular-binary orbital angular momentum is
For conservative binary mass transfer, both and are fixed. In terms of the binary mass ratio , and , hence
Mass conservation also gives , so
Differentiating at fixed gives . Therefore the stipulated Roche lobe law has the Roche-lobe radius response exponent
Define the overfill by . On a time scale short compared with the slow change of , the donor's response is . Thus
During mass loss . If , mass loss increases the overfill and therefore drives more mass loss: this is positive feedback, giving more rapid Roche-lobe overflow. If , mass loss decreases the overfill and can balance the slow expansion due to stellar evolution. More explicitly, with , contact requires
A slow negative transfer rate can satisfy this only on the stable side . The rapid-transfer condition is therefore
Its time scale depends on the physical donor response represented by ; it need not be a dynamical time scale if is a thermal-equilibrium mass-radius exponent.
For Roche-lobe overflow feedback with a constant wind torque factor, retain mass conservation to the stated approximation but allow the prescribed angular-momentum loss. Since , differentiation gives
The new Roche-lobe radius response exponent is
The same overfill calculation gives
A loss of angular momentum has because , and it lowers the critical ratio by making the orbit contract more strongly.
For the specified solar-mass donor exponent , the answer is
Without the wind torque this becomes , so donors slightly less massive than their companions can already lie on the rapid-transfer side of this model. For the lower threshold is , and for every positive satisfies the strict inequality. There is no finite upper bound on from this criterion. Since the paper does not specify , the wind case has a conditional range rather than a unique numerical cutoff.
Let be the uniform mass density. In spherical coordinates with the chosen diameter as polar axis, the squared perpendicular distance to that axis is . The moment of inertia is therefore
This derives the moment of inertia of a uniform solid sphere about any central axis by spherical symmetry.
Put . The center of mass follows the circle of radius , with position . Let and . Its velocity is , and the contact point is displaced from the center of mass by . With positive spin angular velocity counterclockwise, its contact-point velocity is
The cylinder is fixed, so rolling without slipping requires this velocity to vanish. Consequently . If instead the spin angle is positive clockwise, opposite to increasing , then the printed positive-coefficient convention is
The spin and orbital rotation have opposite directions when both are measured about the same oriented axis; the angular speed is . This specifies the sign convention in the printed expression.
The instantaneous contact point is at rest, so the normal reaction and static friction do no work. Conservation of energy equates the loss of gravitational potential energy from release to the total kinetic energy:
On the descent from , decreases, so
Integrating from release to the bottom gives
The integral is finite: at the release point its denominator is asymptotic to . Contact is maintained, since the inward radial equation of motion gives the positive normal reaction throughout this range.
For frictionless sliding, the normal reaction passes through the sphere's centre and produces no torque; Newtonian gravity also exerts no central torque. With initial zero spin, the sphere remains nonrotating. All the lost gravitational potential energy becomes translational kinetic energy, so
A hollow spherical shell of the same mass and radius places its material farther from the centre, giving a larger moment of inertia. Under rolling without slipping, a given centre speed therefore requires more rotational kinetic energy; the same gravitational drop yields a smaller speed. Replacing the solid sphere by a hollow shell increases and leaves unchanged. The sliding case remains nonrotating and has the same center of mass path and gravitational drop. This is the general comparison for rolling inside a fixed circular cylinder.
During ejection, the remaining mass is , and the water is exhausted at . Consider the momentum of the octopus and the small portion of water expelled during . Initially it is . At the end of this interval, the remaining body has mass and velocity , while the expelled mass has laboratory velocity , to first order. Hence the total final momentum is
The external impulse from quadratic drag is . Equating this to the change in total momentum, and dividing by , proves the equation of motion
The outgoing water's momentum flux is essential: applying to the remaining body alone would omit it.
For dimensional analysis, the independent physical inputs are , with dimensions
The scales , and fix units of mass, time and length. In these units, the only remaining independent dimensionless variables are
For example, with and , the equation of motion becomes , with and ejection ending at . Thus the dimensionless terminal value depends only on , proving .
For its explicit value, put and . The initially stationary solution stays below , since the acceleration vanishes at and is positive below it. Separation of variables gives
The two integrals yield
This is constant mass-loss propulsion with quadratic drag. At , the remaining mass is , giving
Thus , verifying the required dimensional form. The limit as drag tends to zero is , the ordinary rocket equation, while the positive-drag answer is below the terminal velocity .
A renormalizable quantum field theory requires only finitely many independent field and parameter redefinitions to remove regulator dependence from its perturbative correlation functions, order by order. The allowed local counterterms must close within that finite set; it is not necessary that the unrenormalized integrals be finite.
In four spacetime dimensions the kinetic term fixes the mass dimension of to one. A coefficient multiplying therefore has dimension . Power counting excludes couplings of negative mass dimension for perturbative renormalizability. For a single scalar with a conventional nonderivative polynomial potential, the allowed form is
The constant controls vacuum energy and the linear term a possible tadpole. A field-reflection symmetry can remove the odd terms but is not necessary for renormalizability. Terms of degree above four generally require an infinite tower of counterterms when treated as fundamental interactions; they can instead be used in an effective field theory. Stability of the potential is a separate requirement, for example a positive quartic leading term, and should not be confused with the power-counting criterion.
For a graph with external scalar lines, internal lines, loops and vertices of degree , the superficial ultraviolet divergence degree is . The identities and give
When every , only the finitely many low-point local divergences, including the kinetic term, can require subtraction. Subdivergences are removed recursively by the same local counterterms. A vertex with instead permits arbitrarily high divergence degrees as more such vertices are added. This explains the restriction rather than merely asserting it.
The bare Lagrangian density is written in terms of regulator-dependent bare quantities, before taking the regulator away:
The coefficients in are related to the renormalized parameters by regulator-dependent redefinitions. The wave-function renormalization supplies a kinetic counterterm. In dimensional regularization , a quartic coupling has a bare relation of the form , with field factors incorporated according to the chosen definition. The bare parameters and bare correlation functions are held fixed when changing the renormalization scale.
Even without a physical mass, interacting loop subtractions require a reference scale to define the renormalized dimensionless coupling and field normalization. Finite logarithms then involve dimensionless combinations such as or . This scale is a subtraction convention, not a new physical mass parameter. Its explicit dependence is compensated by the running coupling and field normalization.
Let at separated points. Its bare counterpart satisfies . Define
Differentiating at fixed bare theory gives
This Callan-Symanzik equation expresses independence of the arbitrary subtraction scale. The renormalization-group beta function describes the compensating coupling change; the field-renormalization anomalous dimension describes the compensating rescaling of each insertion. Coincident composite insertions would require their own additional renormalization, rather than automatically obeying this separated-point equation.
For the two-point function, set and . We consider large positive spacelike , with continuation to other momenta as required. The dimensionless propagator factor obeys
Define and . The renormalization-group characteristic solution for a two-point function is
To verify it, the running-coupling flow has . Also, . Applying to the displayed expression therefore gives , exactly the required differential equation. The boundary value records the renormalization condition. This shows how the ultraviolet behaviour is determined by the coupling trajectory and its accumulated anomalous dimension: it can approach zero coupling, a nonzero fixed point, or a finite-scale singularity, depending on .
For with , direct integration gives
Consequently
For a regular free-theory normalization , the leading large- form is , times corrections from the small running coupling. This is asymptotic freedom with a logarithmic, rather than a fixed-power, propagator correction. The hypothetical beta and gamma functions are being used as given; this calculation does not assert that they are the actual beta and gamma functions of an arbitrary quartic scalar theory.
For the second beta function put , so . Its positive nonzero fixed point and linearized slope are
For the coupling increases with momentum, while for it decreases; each positive trajectory approaches in the ultraviolet. Negative initial couplings approach by the odd symmetry of the beta function, with the same limiting squared coupling. The exactly zero trajectory remains zero. Thus the nonzero ultraviolet fixed point of a cubic-quintic beta function replaces asymptotic freedom by a finite-coupling ultraviolet limit, with . If the anomalous dimension is regular there, the propagator factor has fixed-point behaviour when its fixed-point boundary value is nonzero. If is retained, the exponent is . Trusting this perturbative truncation as a statement about the underlying theory requires to be sufficiently small; the flow result itself follows from the specified beta function.
For a single chiral superfield, the supergravity F-term potential, with restored, is
The canonical Kähler potential has , so the Kähler covariant derivative of a superpotential is
Consequently the Polonyi model has the scalar potential
The negative term is essential; retaining only the square of the supergravity auxiliary field would give the global F-term scalar potential rather than the supergravity F-term potential.
Set and . In this notation , where
At , , one has and , so . To prove that this is a local minimum in both real directions, rather than checking only the real axis, put . Direct expansion gives the stable zero-energy Polonyi vacuum identity
Every coefficient outside the square is strictly positive, and the exponential factor is positive. Thus everywhere and equality requires , provided . This proves the stronger conclusion
a unique global minimum of this Polonyi model. If , the scalar potential vanishes identically and this uniqueness and supersymmetry breaking conclusion do not hold.
For an independent local stability check, the Hessian matrix of at this global minimum is diagonal:
Since and its first derivatives vanish there, differentiating the exponential introduces no additional Hessian matrix terms at the global minimum. For canonically normalized real fluctuations , the squared scalar masses are therefore
which are both positive.
Although the vacuum energy is zero, the supergravity auxiliary field is not. In the convention , its vacuum expectation value obeys
The matter fermion transformation contains a term proportional to , so no nonzero constant supersymmetry parameter leaves this vacuum invariant. Hence supersymmetry is spontaneously broken even though the vacuum energy vanishes: the positive auxiliary field contribution cancels the negative gravitational term. The matter fermion supplies the goldstino, which is absorbed by the gravitino in the super-Higgs mechanism.
Using the gravitino mass from a superpotential, and , gives
An auxiliary field has mass dimension two. Defining the supersymmetry breaking scale by , one obtains
Thus the gravitino mass is small compared with the supersymmetry breaking scale when that scale is sub-Planckian, as in the usual hierarchy . This smallness is gravitational suppression; it is not a conclusion valid for arbitrary choices of the mass parameter.
Use the scalar components of the three chiral superfields, and set , , and . We assume the physical sign and a nonzero coupling . The coefficient below is exactly the coefficient in the supplied D-term scalar potential; it should not be silently replaced by a differently normalized squared gauge coupling. The canonical Kähler potential has identity Kähler metric, and the three Kähler covariant derivatives of a superpotential are
Substitution in the supergravity F-term potential gives
In particular, the subtraction cancels only half of the six terms contributed by the three squares. The full scalar potential is
Its nonnegativity follows term by term from and . This is a special property of this cubic superpotential, not of a general supergravity F-term potential.
The D-term vacuum branches of a cubic superpotential must be described as branches: there are not two isolated field-space minima. First suppose . Vanishing of requires , and vanishing of the D-term scalar potential requires . Together these conditions give the vacuum expectation values
The phase of lies on a gauge orbit, and one may choose it positive and real. At this global minimum, all and the real auxiliary field of the vector multiplet vanish. It is a supersymmetric vacuum, even though the internal gauge symmetry is spontaneously broken. For , interchange and : the nonzero vacuum expectation value has magnitude . For , the supersymmetric vacuum is the entire family , with arbitrary.
The other branch has and arbitrary . Here but the real auxiliary field is nonzero for , and
Thus this branch has pure D-term supersymmetry breaking and a flat direction of a scalar potential in the neutral scalar. It is not stable at every neutral-field value. At a fixed , expansion in the charged scalars yields
Hence the two charged squared masses and the stability condition are
Above this threshold, both charged directions have positive mass squared. Continuity makes them uniformly positive in a sufficiently small neighborhood of the chosen neutral-field value, so the stationary branch is a genuine, non-strict local minimum valley. Its neutral direction stays flat. Below the threshold, a charged mass squared is negative and the branch is a saddle point. At equality it is a boundary of the stable valley, not a local minimum in the full field space: changing slightly downward and switching on the unstable charged field lowers . For example, for , put , , . Since has , one finds .
There are no further local-minimum branches under these generic assumptions. If , the expression for is strictly increasing in , so a stationary configuration with a charged condensate must have . At , an interior critical point with is impossible, because
On the boundary , the relevant scalar potential is ; on it is . Their nonnegative-domain minima give exactly the supersymmetric vacuum described above, or the origin of the charged-field branch. Thus, for nonzero , the two types are the zero-energy supersymmetric gauge orbit and the positive-energy D-breaking minimum valley above the threshold. If , the two types coalesce into the flat supersymmetric vacuum; the literal two-minimum assertion needs this parameter qualification.
In the global supersymmetry limit , the F-term scalar potential becomes and the D-term scalar potential retains its form. The same two types of branches occur, with the threshold . Supergravity instead introduces the exponential and additional positive terms. The goldstino of global supersymmetry breaking is a physical fermion; in local supersymmetry it participates in the super-Higgs mechanism. A positive vacuum energy also gravitates in supergravity. In particular, the positive-energy branch here has and hence vanishing gravitino mass from a superpotential; it is not the tuned zero-energy broken vacuum of Question 3.
There is a further gauge-consistency qualification. In conventional two-derivative matter-coupled supergravity, a constant Fayet–Iliopoulos term gauges an R-symmetry and requires a corresponding nontrivial transformation of the superpotential. The apparent charges make the supplied cubic superpotential neutral, so a nonzero independent constant does not by itself specify a consistent conventional supergravity action. The computed branches are those of the supplied scalar potential; a full embedding needs compatible charges or additional fields. This restriction is explained in Van Proeyen's derivation of the superpotential gauge-covariance condition.
More generally, the supergravity moment-map constraint on D-term breaking follows directly from gauge covariance. If and the moment map is , then at a point with ,
Consequently all chiral auxiliary fields vanishing implies the D-term also vanishes there. This relation has no analogue forcing that implication in an arbitrary global supersymmetry model. Its hypothesis matters: the formal D-breaking branch above has , so it does not contradict the identity. These conventional gauge-completion restrictions are additional to the elementary minimization calculation.
Take upward as positive and let be the change in rocket mass during . The exhaust emitted has mass and, to first order, upward velocity . The linear momentum of the rocket and just-emitted gas changes by
The external gravitational impulse is to first order. Applying the momentum balance to this material system, rather than applying a fixed-mass formula to the rocket alone, gives the rocket equation
The term is the upward thrust.
For the specified mass loss and exhaust law, and
Thus the condition for strictly positive upward acceleration at launch is
At equality the initial thrust just balances the weight; under this particular exhaust law the acceleration becomes positive immediately afterwards. If the initial thrust is smaller than the weight and the rocket initially rests on a launch surface, a nonzero supporting force persists until a later lift-off, so the free-flight equation does not describe its earlier supported motion.
For immediate free flight, integrate the rocket equation with :
This applies while and the nonrelativistic approximation remains valid. The formal divergence as the model's remaining mass approaches zero is not a prediction within that approximation.
For dimensional analysis, write the mass, length and time dimensions as . Then , , , , and . Both and have mass dimension, so is dimensionless; both terms in the boxed velocity have dimension . The launch inequality compares two forces, each of dimension .
Build the spherical system by bringing successive thin spherical shells in from infinity. The interior mass is , and the next shell has . By the spherical shell theorem, the interaction gravitational potential energy of that shell with the assembled interior is . Each pair of mass elements is counted exactly once in this construction. Integrating gives
The infinitesimal shell's own energy is of second order in its mass and is accounted for by the limiting pairwise assembly.
For spherical density projection, let run along the line of sight at projected radius . Spherical symmetry gives
On , use , so . The two halves of the line of sight contribute equally, yielding
To derive the inverse Abel transform, set , , and . The projection becomes . Compose this with the same square-root kernel and exchange the convergent integrals:
The inner integral equals , as follows from with . Write the left side as before differentiating; this avoids differentiating a singular moving endpoint. Under sufficient regularity and decay,
Finally put . Since , this gives the spherical Abel deprojection
Next, the strip density of a spherical system is the integral across the projected coordinate perpendicular to the strip:
For , the substitution on yields
There is a particularly simple way to compose these two projections: integrates the original mass density over the whole plane of coordinates . Polar coordinates in that plane, with , give
Differentiation now proves
The central density, if finite, is obtained by taking the limit as .
To obtain gravitational binding energy from strip density, use in the shell-assembly result, and note also that . Two integrations by parts then give
Here ; the calculation assumes convergent binding energy and sufficient regularity and decay for these boundary terms to vanish. Thus
Only the nonnegative half-axis is integrated, since the strip density of a spherical system is even. Dimensionally, is mass per length, so the last expression has the units of energy.
Use stellar homology: dimensionless density, temperature and luminosity profiles are taken to have the same shape, so the dimensional scale relations follow from the structure equations. Mass conservation and hydrostatic support give
The ideal gas equation then gives , with the dimensional constants held fixed. The specific stellar energy-generation rate has the prescribed hydrogen factor and temperature exponent, so
Independently, radiative diffusion in a star gives , hence
Equating the nuclear supply and radiative transport establishes composition-dependent CNO homology:
At one fixed composition, these reduce to and as required.
The Stefan–Boltzmann law gives , so . Thus the logarithmic Hertzsprung-Russell diagram slope along this zero-age main sequence is
Increasing mass moves upward and leftward on the conventional diagram, whose temperature axis is reversed.
For composition changes, neglect the small heavy-element contribution to the mean molecular weight and set . Then . The fixed trace catalyst abundance is absorbed into ; this is not an assertion that an exactly metal-free star can undergo CNO burning. Substitution into the homology laws gives
For a retained, nonnegligible fixed metallicity , the corresponding mean-molecular-weight factor is instead. The displayed requested formula uses the trace-metal approximation.
At fixed mass, differentiating with respect to at the formal reference gives
Consequently
and the initial fully mixed hydrogen-burning evolutionary track has slope
Hydrogen depletion has , so both luminosity and effective temperature initially increase. The star therefore moves upward and leftward, on a shallower logarithmic track than the mass sequence. The figure uses dimensionless luminosity and temperature relative to its initial reference star; the track is computed for declining at constant mass.
Figure 1.
Zero-age radiative homology sequence and the fixed-mass fully mixed hydrogen-burning track
. Details.
Ordinary cluster main-sequence turn-offs instead evolve toward brighter, cooler, more extended configurations and eventually the giant region. Normal stars are therefore not mixed throughout their entire mass during hydrogen burning: the core composition changes while the envelope retains hydrogen, and composition and entropy gradients develop. Core exhaustion and subsequent shell burning then produce very different structural evolution. Complete mixing is a special evolutionary scenario, not the standard explanation of ordinary giant formation.
These hydrogen-burning scalings must not be extrapolated literally to . Their vanishing nuclear coefficient would make the inferred radius tend to zero and the effective temperature singular; in reality hydrogen exhaustion requires a different structural and energy-generation model. A homogeneous remnant may be helium-rich, but its later structure is not described by the same CNO law. Also, the cover-sheet radiation-pressure expression has a temperature-power typo: the physical term is . Radiation pressure was neglected in this derivation, so that typo does not enter the calculation.
In the radiative layer, take the enclosed mass to be , the transported luminosity to be , and the mean molecular weight to be constant. The envelope's own gravity is neglected. Divide the hydrostatic and radiative temperature equations, and then eliminate density using the gas equation of state:
Put and . For , integration gives
Defining and yields
In the relevant positive-exponent case, writing produces the stated form whenever this integration constant is nonnegative. The exceptional values excluded above lead to logarithms: gives , gives , and both zero give .
For radiative-envelope convection matching, compute the logarithmic gradient of the pressure relation:
At the marginal boundary, the Schwarzschild criterion equates this with the ideal-gas adiabatic gradient . Hence
The positive real notation requires a nonnegative bracket; a general signed is the safer integration constant otherwise. For the later case and a monatomic gas, this bracket is , so is of the same order as .
Deep below the convective boundary, with the temperature-power term dominating , the pressure law becomes . Differentiating this approximation and using hydrostatic balance gives
The general temperature is a term proportional to plus a boundary-dependent constant. Well inside the boundary radius, where that term dominates, it reduces to
The inverse-radius profile is an inner approximation, not an assertion that the boundary term vanishes everywhere.
Now take , so , . Then , , and define so that . The luminosity supplied by the thin burning layer is
Because the integral is strongly concentrated near and , the outer term is negligible. Thus
The factor , rather than merely , is essential since the specified energy-generation rate is per unit mass. This integration uses the deep radiative profile as a leading thin-shell approximation.
Crucially, is not independent of luminosity. Its defining relation gives
Substitution into the shell integral therefore yields . The degenerate-core relation gives , so
This degenerate-core shell-burning mass-luminosity relation is specific to the adopted opacity, burning exponent and structural approximations.
Removing some envelope mass while retaining the hydrogen-burning shell and the deep-envelope regime leaves the leading luminosity approximately unchanged at fixed core mass. The envelope can instead change the radius and effective temperature. This independence is not valid for arbitrarily complete stripping: removing the fuel or the pressure-confined shell structure quenches or modifies shell burning, and the assumed giant model then fails.
Write and for the number densities of nickel-56 and iron-54. Both nuclei contain 28 neutrons, while their proton numbers differ by two. Dividing the positive-exponential nuclear statistical equilibrium formulas printed in the PDF therefore gives the nickel–iron nuclear equilibrium ratio
The neutron-density factors cancel. contains the isotope-dependent statistical and translational mass factors suppressed in the supplied proportionality; it is constant in this simplified calculation. The TeX conversion incorrectly reverses the exponential sign, so the PDF sign is essential here.
In the specified restricted composition, equality of total proton and neutron numbers means
Consequently
Differentiation, with the suppressed prefactors held constant, gives
Using the paper's supplied signed energy difference produces
With that negative value, tends to zero at both temperature endpoints, increases below this temperature and decreases above it. At the stationary point , so it is indeed a maximum for the printed toy data.
There is a genuine physical-data qualification: the stated negative energy difference is inconsistent with interpreting as the conventional positive total nuclear binding energy of these actual isotopes. Using atomic masses, whose electron masses cancel in this reaction,
The masses used are , and , from NIST hydrogen masses, NIST iron masses and NIST nickel masses. Small electronic binding corrections do not change the sign. With a positive , the simplified above is strictly decreasing for and has no finite maximum. Thus the numerical maximum follows from the supplied negative parameter, not from consistent physical binding energies; changing the TeX sign to recover a maximum would not repair the PDF's data.
Nuclear statistical equilibrium is appropriate when the relevant forward and inverse strong or electromagnetic reactions are fast compared with changes of temperature, density and composition. It is useful in sufficiently hot late silicon-burning material and in the hottest explosive burning regions of massive stars. Some burning zones reach only quasi-equilibrium, and cooling or expansion eventually freezes the abundances out; an equilibrium expression cannot by itself predict the full ashes when these timescale conditions fail. Temperature, density, conserved baryon number and the electron fraction in stellar matter must all be specified. Negligible free-particle contributions in bulk bookkeeping are not the same as exactly zero free densities in the equilibrium formula.
Strong and electromagnetic reactions preserve total charge and baryon number, so they conserve the proton-to-neutron ratio. Weak reactions need not: electron capture and beta-plus decay convert protons into neutrons and lower it. Such conversions occur in hydrogen-burning networks and in the nitrogen-to-neon chain during helium burning. At high density in advanced stellar cores and core collapse, electron capture becomes particularly important; neutrino emission also affects energy and lepton loss. Beta-minus decay changes the ratio in the opposite direction and can compete during later relaxation. Nuclear equilibrium for strong reactions does not automatically imply weak equilibrium or a fixed proton-to-neutron ratio throughout all stellar evolution.
Use for the neon-22 mass fraction, temporarily keeping it distinct from the earlier metallicity. Let be the baryon density, ignoring small nuclear mass defects. Equal masses of carbon and oxygen give
Carbon-12 and oxygen-16 each have equal proton and neutron numbers, while neon-22 has ten protons and twelve neutrons. The total proton and neutron numbers per baryon are consequently
Hence
The equal carbon/oxygen mass split cancels from this result: any mixture of these two symmetric nuclei has the same electron fraction. Taking the paper's metallicity symbol to mean gives its requested expression.
The literal pre-helium-burning nitrogen interpretation, however, needs a correction. Each initial nitrogen-14 nucleus receives two alpha particles, including eight additional baryons, before becoming neon-22; the intervening beta-plus decay also emits a neutrino. Thus, if the entire original metallicity is nitrogen-14 with mass fraction , conservation of the number of seed nuclei gives
This is neon-22 neutron excess from nitrogen-14. The two expressions are compatible only when their arguments are distinguished: the displayed law uses the final neon mass fraction, not literally the earlier nitrogen mass fraction. Treating those fractions as equal is an extra simplifying approximation, not a consequence of the stated capture chain.
A radioactively powered Type Ia supernova light curve rises to a broad maximum over roughly weeks as energy deposited in expanding ejecta diffuses out, then declines over the following weeks and months. Its power comes predominantly from the nickel-56 decay chain, first nickel to cobalt and then cobalt to iron. Gamma rays and charged decay products deposit heat, which is reprocessed into optical and infrared radiation. At later times, declining cobalt decay power and increasing gamma-ray escape shape the tail. It is not sustained nuclear fusion in the expanding ejecta. The amount of radioactive nickel, photon diffusion time and deposition efficiency all matter, so an uncorrected fixed luminosity is only an idealization.
For the requested two-isotope yield calculation, assume the same total processed iron-group mass, namely as inferred from the zero-neon reference. Treat the final freeze-out ashes as nickel-56 and iron-54, neglecting residual light particles and additional weak neutronization. If is the nickel mass fraction, conservation of electron fraction gives
Since , this reduces to the metallicity dependence of nickel-56 yield
Under the paper's intended identification ,
If instead denotes the earlier nitrogen mass fraction, retaining the actual capture stoichiometry gives and . The same temperature and density do not alone fix a total burned mass; equality of that mass is also implicit in the comparison. The two-isotope assumption is a bulk-ash approximation, not a complete equilibrium network with all free particles identically absent.
For otherwise comparable explosions, approximate the peak luminosity by a quantity proportional to the radioactive nickel mass. Let be the true distance to an extremely low-neon event, and use the solar calibration . Its measured flux is , but the inferred standard-candle distance is
The source is intrinsically brighter than the calibration, so the distance is underestimated by about , conventionally described here as about . Equivalently, , consistent with the supplied rounded factor . Using the literal nitrogen metallicity instead gives about a underestimate. These are conditional composition effects in this simplified model; real Type Ia distance estimation also standardizes light-curve shape and color.
Since the donor fills its Roche lobe, use and substitute its equilibrium mass-radius relation into the lobe approximation:
Kepler's third law then gives
so the Roche-lobe period-mass relation for a linear donor radius law is
The total binary mass cancels, rather than being assumed equal to the white-dwarf mass.
The orbital radii about the center of mass are and . Summing the two orbital angular momenta gives
With transferred matter retained by the accretor, total mass is fixed. Logarithmic differentiation gives
For , the separation response is . Since at fixed total mass,
The second relation uses the supplied equilibrium donor sequence. As the donor loses mass, . If , its Roche lobe shrinks fractionally faster than this donor radius, so the excess filling increases and drives more transfer: it is positive feedback. If , infinitesimal loss instead reduces overfill on that assumed sequence, so some continuing driver is needed to retain contact.
There are two limits to the literal dynamical interpretation. The given initial range is , so is an extrapolated instability test, not a state inside that range. Also, a true dynamical stability of binary mass transfer test needs the donor's adiabatic stellar radius response exponent , rather than its thermal-equilibrium exponent. In general,
Thus is the result for the assumed response. A fully convective idealized donor has , giving a dynamical threshold , already below . The equilibrium radius law alone should not be used to assert a universal dynamical threshold.
One continuing driver is magnetic braking of a binary star. A magnetized wind carries away the donor's spin angular momentum; tides that maintain near-synchronous rotation replenish that spin from the orbit, removing orbital angular momentum. With an external , the lobe response becomes
Equating this to the donor response proves the binary mass-transfer contact equation in this model:
A negative torque therefore drives mass loss while maintaining contact for , subject to the chosen donor response and braking assumptions.
To obtain the specified expansion, use the usual nova-ejection prescription: the escaping material is expelled isotropically from the white dwarf and carries its specific orbital angular momentum, with no additional torque. This is isotropic re-emission from a binary star. Circularity alone does not determine the angular momentum carried by the ejecta; the prescription is an additional physical assumption.
The accretor's orbital radius is , so its specific angular momentum satisfies
During ejection, take , and , where . Consequently
Differentiate to first order:
Using yields
The same lobe law now has changing total mass but unchanged donor mass, so
The donor radius is unchanged on the short eruption timescale, whereas its lobe expands. It therefore falls inside the lobe and Roche-lobe overflow stops in this idealized model. This is nova-induced binary detachment.
The need for the ejection prescription can also be seen algebraically. If the ejecta instead have mean specific angular momentum , the same circular-orbit differentiation gives
The requested coefficient is obtained only when . Other prescriptions can produce a different expansion or a contraction; the stated circularity condition does not rule them out. The ejection may be short compared with the transfer time while remaining slow compared with an orbit, consistent with circularity. If truly impulsive, an eccentricity is generally excited unless an additional circularization assumption is made.
Now put for the constant external torque. Keep the masses, mass ratio and torque coefficients fixed to leading order over one small- cycle; corrections in those coefficients affect the calculated times only at the next order. During detachment there is no transfer, so fixed component masses give
The fractional lobe gap closes in time
During the semi-detached accumulation phase, use the same external torque and the equilibrium donor contact equation derived in (a):
Accumulating the next layer therefore takes
Division cancels the layer mass, torque and angular momentum:
The result uses for a positive steady accumulation rate, unchanged donor structure during the eruption and the same external torque in the two phases. Here the given ensures a positive denominator. It describes an idealized nova cycle; extra frictional loss, irradiation-driven donor changes or a different ejection angular momentum would change the ratio.
All gas being used or expelled means , so in the formal leaky-box model of galactic chemical evolution and the final stellar mass is . Weight the stellar metallicity by formed stellar mass, using :
Here , by integration by parts. The final mass-weighted mean is
As a metal-budget check, newly produced metal mass is ; stars retain , while the galactic outflow carries . Their sum is exactly . An observed luminosity-weighted stellar metallicity need not equal this mass-weighted mean.