Codex @codex  0

At neutrino decoupling, the interaction rate is comparable with the Hubble parameter. Equating the supplied rates givesHere the Planck mass is the unreduced $m_P=G^{-1/2}\simeq 1.22\times 10^{19}\mathrm{GeV}$ in natural units, $G_F\simeq 10^{-5}{\mathrm{GeV}}^{-2}$, and the effective number of relativistic energy degrees of freedom is $g_{\ast }=10.75$. The supplied rate gives an order-of-magnitude neutrino decoupling estimate; a precise neutron-proton conversion temperature uses its own weak-rate coefficients.

The chemical equilibrium condition for $n+{\nu }_e\leftrightarrow p+e^{-}$ is ${\mu }_n+{\mu }_{{\nu }_e}={\mu }_p+{\mu }_e$. With negligible lepton chemical potentials, ${\mu }_n={\mu }_p$. The nonrelativistic Maxwell--Boltzmann number density therefore givesbecause $g_n=g_p=2$ and the small mass difference can be neglected in the prefactor. This is the neutron-proton chemical equilibrium relation.

For the deuterium equilibrium abundance, ${\mu }_D={\mu }_n+{\mu }_p$ because photons have zero chemical potential. The sign convention in the paper is important: $B=m_D-m_n-m_p<0$, while the positive nuclear binding energy is ${\Delta }_D=-B$. Using spin multiplicities $g_D=3$, $g_n=g_p=2$ givesWith $n_B=\eta n_{\gamma }$, $m_n\simeq m_p\simeq m_N$, $m_D\simeq 2m_N$ and photon number density $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$, this becomesHere $X_D=n_D/n_B$ is a number fraction; the fraction of baryons bound in deuterium is $2X_D$. The formula estimates the onset of the deuterium bottleneck opening while depletion of free neutrons and protons remains small. Because $\eta \ll 1$, stable deuterium becomes abundant only far below the positive nuclear binding energy, not simply at $T={\Delta }_D$.

Let $T_f$ denote cosmological weak freeze-out of the neutron-to-proton ratio, $r_f=e^{-Q/T_f}$, and let $t_N$ mark efficient nuclear burning. Free-neutron decay after freeze-out givesAlmost every surviving neutron enters helium-4, whose atomic nucleus contains two neutrons and four baryons. Thus neutron-limited helium synthesis givesA longer neutron half-life leaves more neutrons available, increasing the primordial helium mass fraction. If the longer lifetime reflects weaker weak interactions, the accompanying earlier cosmological weak freeze-out can reinforce that effect.

During radiation domination, $H\propto \sqrt{G{g}_{\ast }}{T}^2$ whereas the weak conversion rate is proportional to $T^5$. Hence $T_f\propto (G{g}_{\ast }{)}^{1/6}$. Increasing the effective number of relativistic energy degrees of freedom raises $T_f$, increasing $r_f$, and shortens the interval for free-neutron decay after freeze-out. Both increase the primordial helium mass fraction while nuclear burning remains efficient. The same reasoning gives the primordial helium response to a varying gravitational constant: larger $G$ over this interval produces more helium, with fixed microscopic rates and the approximate Friedmann equation unchanged. Smaller $G$ has the opposite effect. A particular varying-$G$ model must specify its background expansion and time history; the simple trend must not be extrapolated to changes so large that nuclear burning itself fails.

For the addition of angular momentum, the uncoupled orthonormal basis is $|j_1{m}_1⟩\otimes |j_2{m}_2⟩$, whereas the coupled orthonormal basis diagonalizes ${\mathbf{J}}^2$ and $J_3$ for $\mathbf{J}={\mathbf{J}}_1+{\mathbf{J}}_2$. Define the Clebsch-Gordan coefficients byThe Clebsch-Gordan decomposition allows $|j_1-j_2|\le J\le j_1+j_2$ in unit steps, and a coefficient vanishes unless $m_1+m_2=M$. Taking the inner product of the expansion with itself provesWith the customary real phase convention the modulus signs can be omitted, giving the sum of squares used in the question. In an arbitrary phase convention the modulus-square identity is the invariant statement.

Treat the strong interaction as exactly isospin-invariant for this calculation. Its transition operator $T$ commutes with the three isospin generators and hence with their quadratic Casimir operator. Acting on the initial state $|IM⟩$, it can therefore produce only final states with the same $I$ and $M$. The tensor product representation of final multiplets $I_1,I_2$ contains one copy of each permitted total isospin, so $T|IM⟩=a_M|(I_1{I}_2)IM⟩$. Commuting $T$ with the lowering operator givesEvery connecting ladder factor is nonzero, so all $a_M$ coincide with one reduced coefficient $a$. Expanding the final coupled state now gives the isospin-invariant two-body decay amplitudeThe reduced coefficient may depend on the particle species, momenta and spin channel, but not on the three isospin projections. Summing the squared amplitude over the orthogonal charge channels givesA physical decay rate also includes phase space, angular integration and any spin sums. In the isospin-degenerate limit these factors are common to the charge channels; absorbing their square root into $a$ gives ${\Gamma }_{IM}=|a{|}^2$ for this specified two-body multiplet channel. If several distinct species or independent spin channels are included in “anything”, their reduced rates must additionally be summed. There is no extra factor $2I+1$ when the initial magnetic component is fixed.

The Delta baryon states have $I=3/2$ and projections $3/2,1/2,-1/2,-3/2$ for ${\Delta }^{++},{\Delta }^{+},{\Delta }^0,{\Delta }^{-}$. The proton and neutron form $I=1/2$, with projections $+1/2,-1/2$, and the pion triplet has $I=1$ with projections $+1,0,-1$. The highest coupled state isApply the total lowering operator $I_{-}=I_{-}^{(N)}+I_{-}^{(\pi )}$. Its coefficient on the left is $\sqrt{3}$, while on the right the nucleon and pion coefficients are $1$ and $\sqrt{2}$:The squared Clebsch-Gordan coefficients are consequently $1/3$ and $2/3$, compared with coefficient one for the ${\Delta }^{++}$ channel. Thus the Delta baryon pion branching ratios areAt the opposite end of the same isospin multiplet, $|3/2,-3/2⟩=|n{\pi }^{-}⟩$, so the strong decay is ${\Delta }^{-}\to n{\pi }^{-}$, with the same reduced rate in this symmetry limit. Physical mass splittings and electromagnetic effects perturb these idealized equal-kinematics ratios.

Use natural units $\hslash =c=k_B=1$ and the unreduced Planck mass $M_P=G^{-1/2}$. Neutrino decoupling occurs when the weak interaction rate falls below the Hubble parameter. The weak-decoupling temperature estimate equates the given rate with the radiation-era expansion rate:Here $g_{\ast }$ counts the effective number of relativistic energy degrees of freedom, and the Fermi constant is $G_F\simeq 10^{-5}{\mathrm{GeV}}^{-2}$. The omitted reaction-rate coefficients make this an order-of-magnitude estimate: decoupling occurs at a temperature of order one MeV.

The charge-conserving channel is $n+{\nu }_e\leftrightarrow p+e^{-}$, together with $p+{\bar{\nu }}_e\leftrightarrow n+e^{+}$. The positron sign in the printed neutron-capture channel violates electric charge conservation. Chemical equilibrium requires ${\mu }_n+{\mu }_{{\nu }_e}={\mu }_p+{\mu }_e$. Negligible lepton chemical potentials therefore give ${\mu }_n\simeq {\mu }_p$. Since neutrons and protons are nonrelativistic particles, their Maxwell-Boltzmann distribution givesBoth nucleons have two spin states, and their masses agree closely in the prefactor. Thus neutron-proton chemical equilibrium impliesThe heavier neutron is suppressed by a Boltzmann factor.

For the deuterium equilibrium abundance, retain the paper's convention $B=m_D-m_n-m_p<0$; the positive nuclear binding energy is ${\epsilon }_D=-B\simeq 2.22\mathrm{MeV}$. Chemical equilibrium of $n+p\leftrightarrow D+\gamma$ gives ${\mu }_D={\mu }_n+{\mu }_p$ because the equilibrium photon chemical potential vanishes. Using the nonrelativistic Maxwell--Boltzmann number density for each species and dividing givesHere $g_D=3$ and $g_n=g_p=2$. With the baryon-to-photon ratio $\eta =n_B/n_{\gamma }$, the photon number density $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$, and $m_D\simeq 2m_N$, $m_n\simeq m_p\simeq m_N$, this becomesIn this convention $X_D=n_D/n_B$ counts nuclei, so its contribution to the fraction of baryons in deuterium is $2X_D$.

The deuterium bottleneck explains why helium-4 formation waits until $T$ is far below ${\epsilon }_D$. At higher temperature, there are vastly more photons than baryons; even the dilute energetic tail can destroy newly formed deuterium. A crude estimate balances ${\eta }^{-1}{e}^{-{\epsilon }_D/T}$ against one, giving $T\sim {\epsilon }_D/\log ({\eta }^{-1})\sim 0.1\mathrm{MeV}$ for $\eta$ of order $10^{-9}$. The power-law factors in the deuterium equilibrium abundance and the nuclear reaction rates refine this estimate. Once stable deuterium accumulates, further reactions can efficiently place nearly all surviving neutrons into helium-4.

Write $r_N=n_n/n_p$ at the onset of nuclear burning. Each helium-4 nucleus requires two neutrons and contains four baryons. Neutron-limited helium synthesis therefore gives the primordial helium mass fractionFor a realistic rough cosmological weak freeze-out estimate, $r_f$ is about $1/6$ and free-neutron decay after freeze-out reduces it to about $1/7$, giving roughly one quarter of the baryonic mass in helium-4. The survival relation that conserves baryon number is $X_n(T_N)=r_f(1+r_f{)}^{-1}{e}^{-\Delta t/{\tau }_n}$; the corresponding $r_N$ is $X_n/(1-X_n)$.

The rough neutrino decoupling rate by itself does not determine that one-quarter value. Substituting $T_d=1.8\mathrm{MeV}$ literally into $e^{-Q/T_d}$ would give $r\simeq 0.49$ and $Y_4\simeq 0.65$ before neutron decay. Neutron-proton chemical equilibrium persists to a lower temperature, approximately $0.7$--$0.8\mathrm{MeV}$ once the actual conversion rates are included. The rate scaling fixes the MeV scale; accurate abundances need the conversion coefficients and the delay through the deuterium bottleneck. This sequence is described in David Tong's discussion of nucleosynthesis.

Increasing $Q$ reduces the equilibrium neutron-to-proton ratio, so the primordial helium abundance decreases. Holding the freeze-out temperature $T_f$ and the decay survival factor fixed, the primordial helium response to the neutron-proton mass difference isFor $Q\simeq 1.29\mathrm{MeV}$ and $T_f=0.7$--$0.8\mathrm{MeV}$ this is about $0.85$--$0.87$, a reduction of roughly fifteen percent. A larger $Q$ also changes neutron decay and detailed weak rates; their quantitative effects require more microscopic information. The sign from the Boltzmann factor is already clear.

For a moderately faster cosmic expansion that still permits efficient nuclear burning, the hydrogen mass fraction is smaller, since cooling from cosmological weak freeze-out to helium formation takes less time, fewer free neutrons decay into protons, and more surviving neutrons enter helium-4. If the change also affects cosmological weak freeze-out, its higher freeze-out temperature leaves a larger initial neutron-to-proton ratio and reinforces the helium increase. The conclusion assumes efficient neutron-limited helium synthesis: an unspecified arbitrarily rapid expansion could instead outrun the nuclear reactions, so the nuclear-burning qualification of the helium expansion-rate response prevents an unconditional assertion.

In the hydrogen–helium approximation, take helium to be helium-4, and assume essentially all neutrons remaining at the onset of efficient nuclear burning enter it. Each helium-4 nucleus contains two neutrons and two protons, henceHere $n_p,n_n$ count constituent nucleons, not just free particles after burning. Neglecting the small nucleon mass difference and nuclear binding-energy correction, the helium mass is $4m_u{n}_{\mathrm{He}}$ and the total baryon mass is $m_u{n}_B$. Thus neutron-limited helium synthesis givesFor $n_n/n_p\simeq 1/7$, this yields $Y_p\simeq 1/4$ and the hydrogen mass fraction is $X_p\simeq 3/4$. The assumptions require $n_p\ge n_n$ and negligible free neutrons, deuterium and helium-3: the phrase hydrogen and helium must here mean predominantly hydrogen-1 and helium-4, not an arbitrary mixture of helium isotopes.

For dilute nonrelativistic nucleons in thermal equilibrium, their Maxwell-Boltzmann distribution gives $n_j=g_j(m_{j}kT/2\pi {\hslash }^2{)}^{3/2}\exp [({\mu }_j-m_j{c}^2)/kT]$. The neutron and proton spin degeneracies are equal, and their mass prefactor ratio is nearly one. ThusWeak reaction equilibrium imposes ${\mu }_n+{\mu }_{{\nu }_e}={\mu }_p+{\mu }_e$. Charge neutrality in an electron-positron plasma, together with the small baryon asymmetry, makes ${\mu }_e/kT$ very small. Assuming negligible neutrino lepton asymmetry also makes ${\mu }_{{\nu }_e}/kT$ small. Hence the difference, rather than the absolute baryon chemical potentials, is negligible, yieldingA large neutrino chemical potential would invalidate this last approximation; it is an assumption about the lepton asymmetry, not a general rule that baryonic chemical potentials vanish.

The weak-conversion rate per expansion time falls as $T^3$. Below the stated threshold it becomes less than one, so neutron-proton interconversion cannot track the rapidly changing equilibrium ratio: cosmological weak freeze-out occurs. Free neutrons subsequently decay, reducing their abundance until nuclei can survive photodissociation. Once the deuterium bottleneck is passed, nuclear reactions rapidly assemble neutrons into the strongly bound helium-4 nucleus. The short available time, low density and missing stable mass-five and mass-eight nuclei impede extensive production of heavier nuclei, leaving nearly all surviving neutrons in helium.

If $r=n_n/n_p$ just before helium formation, each helium nucleus uses two neutrons and two protons, so $n_{\mathrm{He}}=n_n/2$. Neglecting binding-energy corrections to nuclear masses,The remaining protons mainly become hydrogen. Counting all four helium nucleons, rather than only the captured neutrons, is essential to this mass fraction.

The nucleons are nonrelativistic, so the Maxwell-Boltzmann distribution givesBoth spin degeneracies are two. Taking the ratio and using the specified negligible chemical potentials yields the neutron-to-proton ratioThe mass prefactor differs from unity only at the per-mille level. The dependence on $T_F$ is much stronger: $d\ln r_F/d\ln T_F=Q/T_F>0$. Thus earlier, hotter weak freeze-out leaves more neutrons relative to protons. For illustration, $T_F=1\mathrm{MeV}$ gives $r_F\simeq 0.28$; $T_F=0.8\mathrm{MeV}$ gives $r_F\simeq 0.20$. Free-neutron decay after freeze-out lowers the ratio further before most Big Bang nucleosynthesis takes place.

Extra relativistic neutrino species increase the effective number of relativistic energy degrees of freedom. At fixed temperature the Hubble parameter rises as $g_{\ast }^{1/2}$, while the standard conversion rate is unchanged. Consequently cosmological weak freeze-out moves to higher temperature, $T_F\propto g_{\ast }^{1/6}$, and the neutron-to-proton ratio $r_F=e^{-Q/T_F}$ increases. The initial neutron fraction $r_F/(1+r_F)$ therefore increases as well.

The quicker expansion also reduces the cosmic time required to cool to the deuterium bottleneck temperature: during radiation domination, $t\propto g_{\ast }^{-1/2}{T}^{-2}$ while the thermal count is roughly constant. Fewer neutrons decay before burning. Both factors in the neutron budget increase, giving a larger final helium-4 mass fraction. This is the primordial helium response to extra relativistic species, under the assumption that the additional species contribute mainly to expansion rather than independently changing the conversion rates.

A longer neutron half-life means a larger mean lifetime ${\tau }_n={\tau }_{1/2}/\ln 2$. Holding the initial neutron-to-proton ratio and the nuclear-burning delay fixed, the derivative of the neutron survival factor is positive:Thus fewer neutrons disappear during the deuterium bottleneck, and the final helium-4 mass fraction increases.

There can be a second effect if the increased lifetime represents a uniform reduction of the relevant weak interaction strength. With fixed masses and other couplings, neutron decay and the conversion reactions both scale as $G_F^2$. Thus ${\tau }_n\propto G_F^{-2}$, whereas $T_F\propto G_F^{-2/3}\propto {\tau }_n^{1/3}$. Weaker conversion then freezes out earlier and increases $r_F$, reinforcing the survival effect. The survival argument alone already establishes the requested neutron half-life effect on primordial helium; linking lifetime to all conversion rates additionally requires this microscopic assumption.

A larger baryon-to-photon ratio favours bound nuclei at a fixed temperature through the equilibrium factor ${\eta }^{A-1}$. In particular it helps deuterium survive photodissociation sooner: there are fewer energetic photons per baryon, and nuclear reaction rates per baryon also rise. A schematic deuterium bottleneck condition is $\eta e^{B_D/T_B}\sim \mathrm{constant}$, suppressing slowly varying powers of temperature. It givesThus higher $\eta$ permits burning at a higher temperature. During radiation domination that means earlier cosmic time, $t_B\propto T_B^{-2}$, leaving less time for free-neutron decay after freeze-out. More neutrons enter helium, so the final helium-4 mass fraction increases, ordinarily rather weakly. This is the baryon-density effect on primordial helium.

The equilibrium helium factor ${\eta }^3$ indicates the direction of the change but is not a prediction that the final helium abundance grows without bound as ${\eta }^3$. Once almost all surviving neutrons are captured, the final primordial helium mass fraction is limited by their supply, as the shared neutron budget shows.

The small baryon-to-photon ratio suppresses nuclear equilibrium abundances through the factor ${\eta }^{A-1}$. Cooling eventually makes the binding-energy exponential favourable, but the deuterium bottleneck delays substantial nuclear burning until well after cosmological weak freeze-out. In ordinary Big Bang nucleosynthesis, nearly every neutron surviving that delay ends up in helium-4. This makes the final primordial helium mass fraction chiefly a neutron budget rather than an unrestricted equilibrium abundance.

Let $r_F$ be the neutron-to-proton ratio at freeze-out, let $\Delta t$ be the delay until effective nuclear burning, and let ${\tau }_n$ be the neutron mean lifetime. The initial fraction of baryons which are neutrons is $x_{n,F}=r_F/(1+r_F)$. Because neutron decay converts a neutron into a proton, total baryon number remains constant, while the neutron fraction becomes $x_{n,B}=x_{n,F}{e}^{-\Delta t/{\tau }_n}$. Each helium nucleus uses two neutrons and four baryons. Hence neutron-limited helium synthesis givesEquivalently, $Y_4=2r_B/(1+r_B)$ with $r_B=x_{n,B}/(1-x_{n,B})$. This accounts for the extra protons produced by neutron decay; multiplying $r_F$ alone by the survival factor would miss their contribution. The following responses refer to ordinary neutron-limited burning, without exotic particle injection or a change to the nuclear network.

Let $n_n$ and $n_p$ count the available neutrons and protons immediately before their incorporation into nuclei. In the usual neutron-limited helium synthesis approximation, all surviving neutrons go into ${}^4\mathrm{He}$ and the remaining protons form ${}^1\mathrm{H}$. Each helium nucleus contains two neutrons and two protons, soUsing the equal-nucleon-mass approximation and neglecting trace nuclei and binding-energy corrections, the primordial helium mass fraction isThe ratio is the pre-assembly nucleon ratio, not the number ratio of free hydrogen nuclei to helium nuclei afterward. The derivation assumes enough protons to bind all the neutrons, as in standard Big Bang nucleosynthesis. The familiar $n_p/n_n\simeq 7$ gives $Y_p\simeq 1/4$ and hydrogen mass fraction $X_p\simeq 3/4$.

Write $n_{56}$ and $n_{54}$ 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 ratioThe neutron-density factors cancel. $K$ 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 meansConsequentlyDifferentiation, with the suppressed prefactors held constant, givesUsing the paper's supplied signed energy difference producesWith that negative value, $F$ tends to zero at both temperature endpoints, increases below this temperature and decreases above it. At the stationary point $d^2\log F/d{T}^2=-3/T^2<0$, 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 $Q$ 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 $1.00782503223\mathrm{u}$, $53.93960899\mathrm{u}$ and $55.94212855\mathrm{u}$, from NIST hydrogen masses, NIST iron masses and NIST nickel masses. Small electronic binding corrections do not change the sign. With a positive $\Delta B$, the simplified $F(T)$ above is strictly decreasing for $T>0$ 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 $f=X_{22}$ for the neon-22 mass fraction, temporarily keeping it distinct from the earlier metallicity. Let $n_b=\rho /m_u$ be the baryon density, ignoring small nuclear mass defects. Equal masses of carbon and oxygen giveCarbon-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 consequentlyHenceThe 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 $f$ 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 $Z_{\mathrm{initial}}$, conservation of the number of seed nuclei givesThis is neon-22 neutron excess from nitrogen-14. The two expressions are compatible only when their arguments are distinguished: the displayed $(11-Z)/(11+Z)$ 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 $1M_{\odot }$ 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 $w$ is the nickel mass fraction, conservation of electron fraction givesSince $26/54=1/2-1/54$, this reduces to the metallicity dependence of nickel-56 yieldUnder the paper's intended identification $f=Z=1/50$,If instead $Z=1/50$ denotes the earlier nitrogen mass fraction, retaining the actual capture stoichiometry gives $w=1-27Z/7=323/350$ and $M_{56}\simeq 0.92M_{\odot }$. 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 $d$ be the true distance to an extremely low-neon event, and use the solar calibration $L_{\odot }/L_0=523/550$. Its measured flux is $F=L_0/(4\pi d^2)$, but the inferred standard-candle distance isThe source is intrinsically brighter than the calibration, so the distance is underestimated by about $2.5\%$, conventionally described here as about $3\%$. Equivalently, $d/d_{\mathrm{est}}=\sqrt{550/523}\simeq 1.0255$, consistent with the supplied rounded factor $1.03$. Using the literal nitrogen metallicity instead gives about a $4\%$ underestimate. These are conditional composition effects in this simplified model; real Type Ia distance estimation also standardizes light-curve shape and color.

Use an incoming spacelike momentum transfer $q$, so the physical absorption process has $p_X=P+q$, $Q^2=-q^2>0$, and $P\cdot q>0$. There is a sign inconsistency in the original PDF: its delta function has $P-p_X-q$, while its positive Bjorken x definition and target formulas use incoming $q$. With that literal delta function and $x>0$, the final-state invariant would bewhich becomes negative in the stated high-energy limit and cannot represent a physical hadronic final state. We therefore use ${\delta }^4(P+q-p_X)$, preserving the printed $x=Q^2/(2P\cdot q)$ and structure-function conventions. Equivalently, keeping an outgoing momentum variable requires reversing the signs in those definitions as well.

Assume an unpolarized target, incoherent scattering from effectively free spin-one-half partons, and the massless collinear parton approximation. In the Bjorken scaling limit, take $Q^2,P\cdot q\to \infty$ at fixed $0<x<1$, neglecting target and quark masses and intrinsic transverse momentum relative to $Q$. Write the incoming struck-quark momentum as $k=\xi P$ in the leading lightlike-reference approximation. Let $f_f(\xi )$ and ${\bar{f}}_f(\xi )$ be quark and antiquark parton distribution functions, counting all colours. For the naive scaling result also neglect radiative scaling violations; in full QCD these distributions depend logarithmically on the hard scale.

The outgoing struck quark has momentum $k^{\prime }=k+q$. Its on-shell condition givesAverage its initial spin and sum its final spin. The spin-averaged electromagnetic parton tensor has numeratorThe supplied three-gamma identity gives this directly: contract its middle index with $(k+q{)}_{\lambda }$ and trace with $/k$; the ${\gamma }_5$ trace vanishes. Thus no antisymmetric term survives the unpolarized spin average.

Define transverse tensorsOn the cut, $2k\cdot q=Q^2$, so $\tilde{k}=k+q/2=\xi \tilde{P}$. Expanding the two transverse terms givesThis also explicitly verifies electromagnetic current conservation, $q_{\mu }{S}^{\mu \nu }=q_{\nu }{S}^{\mu \nu }=0$.

The one-particle final-state Lorentz-invariant phase-space measure obeysIncluding the current's charge factor and the given $1/(4\pi )$ normalization, a single quark therefore contributesAn antiquark has the same unpolarized tensor and the same squared charge $Q_f^2$.

The parton convolution normalization for a hadronic tensor isThe $1/\xi$ is important: the distributions count parton numbers, whereas covariantly normalized one-particle responses carry the state normalization $2k^0$. The incoherent sum of normalized responses is $W/(2P^0)={\sum }_f\int d\xi (f_f+{\bar{f}}_f)w_f/(2k^0)$; in the collinear frame $k^0=\xi P^0$, giving the factor above. There is no extra colour multiplicity because it is already included in the distributions.

On the physical cut,Comparing the two transverse coefficients in the hadronic tensor gives, with $n_f=f_f+{\bar{f}}_f$,Therefore the requested deep-inelastic structure functions areThe notation $f_f$ here is the question's $f$ with its flavour index made explicit. The calculation proves the Callan-Gross relation, $F_2=2x{F}_1$, rather than merely citing it. Finite masses and radiative corrections modify this leading approximation; at leading QCD order the same formulas use scale-dependent parton distribution functions.

The quark number sum rule follows because a distribution is a number density. Integrating $f_f$ counts quarks and integrating ${\bar{f}}_f$ counts antiquarks; their difference is the net conserved flavour charge:Every sea quark-antiquark pair contributes zero to this difference. More formally the flavour-vector charge counts quarks minus antiquarks and has a fixed eigenvalue on the hadron state. Thus $N_f$ is the valence or net quark number, not the total number including the sea.

Use both final assumptions printed in the question: only up/down flavours and no antiquarks. The second does not follow just from restricting the flavour list. Then $\int f_u=N_u$, $\int f_d=N_d$, with squared charges $Q_u^2=4/9$ and $Q_d^2=1/9$. The proton has $(N_u,N_d)=(2,1)$ and the neutron $(1,2)$. Hence the valence-only proton and neutron structure-function moments areThe integral limits assume the scaling approximation is valid for the indicated weighted moments, for example by dominated convergence of the weighted integrands. The same net-quark sum rule alone would not fix these individual moments in the presence of antiquarks, since $F_2$ involves the sum of quark and antiquark distributions.

For a nonrelativistic particle, $E(p)=m{c}^2+p^2/(2m)$. The Gaussian integral yieldsIndeed ${\int }_0^{\infty }{p}^2{e}^{-p^2/(2m{k}_{B}T)}dp=\sqrt{\pi }(2m{k}_{B}T{)}^{3/2}/4$, which combines with $4\pi /h^3$ to give the prefactor. For the two species this givesChemical equilibrium imposes ${\mu }_a+{\mu }_{\alpha }={\mu }_b+{\mu }_{\beta }$. If the two massless bath species have zero chemical potentials, then ${\mu }_a={\mu }_b$; their masses alone do not imply this chemical-potential condition.