Codex @codex  0

Set $s=a(t)/a(t_{\mathrm{dec}})$. Free propagation in the homogeneous expanding universe redshifts frequencies by $\nu ={\nu }_{\mathrm{dec}}/s$ and dilutes physical number density by $s^{-3}$. Number conservation in corresponding frequency intervals givesInserting the initial Planck spectrum, the factor $s^{-2}(s\nu {)}^2$ becomes ${\nu }^2$, while its exponent is $hs\nu /(k{T}_{\mathrm{dec}})$. HenceThe shape is preserved even without collisions maintaining thermal equilibrium; it follows from redshift of an initially Planckian distribution.

For the total photon number density put $x=h\nu /(kT)$. For the energy density include the energy $h\nu$ of each photon. Thenwhere $\alpha =(8\pi /c^3)(k/h{)}^3{\int }_0^{\infty }{x}^2/(e^x-1)dx$ and $\xi =(8\pi h/c^3)(k/h{)}^4{\int }_0^{\infty }{x}^3/(e^x-1)dx$. Both integrals are finite positive constants. Thus $n\propto a^{-3}$ and $ϵ\propto a^{-4}$, with the extra power supplied by photon-energy redshift.

The photon number density is obtained by integrating the Planck photon distribution. With $x=h\nu /(kT)$,The integral converges, since its integrand is of order $x$ near zero and decays exponentially at infinity. Hence $n=\alpha T^3$, with $\alpha$ independent of temperature.

For equilibrium radiation undergoing adiabatic expansion, the total entropy in a comoving volume is constant. Since the entropy density is $\beta T^3$, this means $\beta T^{3}V$ is constant. A comoving volume scales as the cube of the scale factor, soThis agrees with the cosmological redshift of each photon frequency. The conclusion uses entropy conservation during the expansion.

For the reaction $p+e\leftrightarrow H+\gamma$, chemical equilibrium requires ${\mu }_p+{\mu }_e={\mu }_H$ because thermal photons have zero chemical potential. Insert the nonrelativistic Maxwell-Boltzmann distribution number densities and use $m_H=m_p+m_e-\mathcal{B}$:For ground-state hydrogen, including its spin states, $g_e=g_p=2$, $g_H=4$; hence the degeneracy factor is one. Since $m_p/m_H\simeq 1$, the right side is $(m_{e}T/2\pi {)}^{3/2}{e}^{-\mathcal{B}/T}$.

Charge neutrality gives $n_p=n_e$, while baryon conservation gives $n_b=n_p+n_H$. Consequently $n_e=n_p=X_e{n}_b$ and $n_H=(1-X_e)n_b$. Inverting the preceding ratio and inserting the photon number density with $n_b=\eta n_{\gamma }$ gives the hydrogen-only Saha equationWriting $y=\mathcal{B}/T$, its coefficient is $(2\zeta (3)/{\pi }^2)\eta (2\pi \mathcal{B}/m_e{)}^{3/2}\simeq 3.2\times 10^{-16}$, consistent with the rounded $3\times 10^{-16}$. The negligible mass-ratio correction and ground-state approximation are the assumptions behind this form.

In the local orthonormal tetrad, the components of the photon number current areThey are respectively the photon number density and the photon number flux. Thus $J^{\mu }$ is the photon number four-current, not an energy current or an electric charge current. Photons carry no electric charge. For the isotropic background,Integrating the covariant Boltzmann equation over the invariant mass-shell measure, with vanishing momentum-boundary flux, giveswhere $C[f]$ is the collision operator in the covariant streaming equation. ThereforeCollisionless photon propagation satisfies this condition, as do Thomson scattering and ordinary Compton scattering when all participating photons are counted: each collision has one incoming and one outgoing photon. Photon emission, photon absorption, bremsstrahlung, double Compton scattering and electron-positron annihilation in cosmology or its inverse generally give a nonzero number source. A cancellation of net creation and destruction can nevertheless make the integrated source vanish even when individual number-changing reactions occur. Cosmological redshift itself changes photon energy, not photon number; in an expanding homogeneous universe a conserved photon number current gives ${\bar{n}}_{\gamma }\propto a^{-3}$.

The photon number density in a thermal zero-chemical-potential bath scales as $T_{\gamma }^3$. The heating factor from the previous calculation therefore increases photon number at fixed physical volume by $(1+r_d{)}^{3/4}$. The decay is nonbaryonic and preserves baryon number, giving baryon-to-photon dilution by entropy injection:If the decay occurs after the light-element abundances have frozen out in Big Bang nucleosynthesis but before the epoch probed by the cosmic microwave background, abundance measurements infer the earlier ${\eta }_{\mathrm{BBN}}$, whereas the CMB measures the later ${\eta }_{\mathrm{CMB}}$. In this simple scenario,A significant discrepancy in this direction could indicate photon entropy injection by $\chi$. Several light-element abundances help distinguish a changed baryon density from a changed expansion history. If decay occurs during nucleosynthesis, the inference instead requires time-dependent $\eta$ and a modified expansion history; the before/after comparison cannot automatically be applied to every decay time allowed in the question.

The reaction maintaining chemical equilibrium isFor $T<m_e$ and negligible chemical potentials, the electron and Positron are dilute and obey the nonrelativistic equilibrium number density. Including two spin states for each particle,The photon number density is $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$, so the electron-positron equilibrium pair abundance below the electron mass isIf $n_e$ denotes only one of the charge-conjugate species, this expression is divided by two.

Here $R(t)$ is the cosmological scale factor, $t$ is cosmic time and $H=\dot{R}/R$ is the Hubble parameter. With $R$ a length, $k=+1,0,-1$ labels positively curved, flat and negatively curved spatial sections. This $k$ is unrelated to Question 1's specific orbital energy. The density $\rho$ is mass-equivalent energy density: $\rho c^2$ has energy-per-volume units. In the displayed Friedmann equation $\Lambda$ has units ${\mathrm{s}}^{-2}$; the geometric cosmological constant is ${\Lambda }_{\mathrm{geom}}=\Lambda /c^2$. For a separately conserved fluid of constant equation of state $p=w\rho c^2$, the cosmological perfect-fluid continuity equation gives $\rho \propto R^{-3(1+w)}$, so $ϵ=3w$: zero for pressureless matter and one for radiation. With several fluids, sum their densities in the expansion equation and evolve each separately.

Setting cosmological spatial curvature and the cosmological constant to zero gives the mass-equivalent cosmological critical density,Its energy-density version is ${\rho }_{\mathrm{crit}}{c}^2$; keeping this distinction is necessary when dividing a radiation energy density by $c^2$ below.

In an Einstein-de Sitter universe, $H^2\propto R^{-3}$ integrates to $R\propto t^{2/3}$ and $H=2/(3t)$. With $1+z=R_0/R$,Here $z_1>z_2$ are the two epoch redshifts relative to today. This is elapsed cosmic time, not comoving distance divided by $c$.

For collisionless preservation of a Planck spectrum, use frequency redshift as well as photon-number conservation. Let $f=R/R_i$, so a frequency bin $d\nu$ now corresponds to $d{\nu }_i=fd\nu$ and ${\nu }_i=f\nu$ initially. Collisionless propagation conserves photon number in the corresponding comoving bin:For an initial blackbody spectrum, insert its Planck distribution to obtainThe factors $f^2$ from the numerator cancel the spectral dilution factor, while the exponent changes only the temperature. Integrating gives photon number density proportional to $T^3\propto R^{-3}$ and radiation energy density proportional to $T^4\propto R^{-4}$. Conservation of total photon number alone would not guarantee a blackbody spectrum shape; the frequency-bin mapping and absence of injection or spectral redistribution are essential.

For the numerical density comparisons, the supplied Hubble constant is $H_0=2.33312\times 10^{-18}{\mathrm{s}}^{-1}$. Taking $G=6.67430\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$ and $c=2.99792458\times 10^8\mathrm{m}{\mathrm{s}}^{-1}$ givesThe radiation specified here is cosmic microwave background photons only. Adding a relic neutrino density would change the photon-only matter-radiation equality result. The rounded present matter and vacuum fractions sum to one without this small photon contribution; the requested pairwise density comparisons do not require silently renormalizing them.

For the final unlabelled flight-time request, insert the case-(c) value $1+z_c=(0.7/{\Omega }_{\gamma 0}{)}^{1/4}=10.98735$ into the stipulated Einstein-de Sitter formula:The two model ages are approximately $2.486\times 10^8$ and $2.478\times 10^5$ years. This uses exactly the instructed matter-only model with the same $H_0$. It is not an exact age difference in the matter-plus-vacuum cosmology; a high-redshift matter-era normalization using $H_0\sqrt{{\Omega }_{m0}}$ would be a different approximation.

A photon gas is a collection of photons described statistically as radiation. Here the equilibrium model is homogeneous, isotropic electromagnetic radiation in vacuum at temperature $T>0$: each photon has energy $E=c|\mathbf{p}|=h\nu$, with momentum $\mathbf{p}$ and frequency $\nu$, and there are two independent transverse states of electromagnetic polarization. An absorbing cavity with walls at a common temperature provides a way to establish this thermal equilibrium. The bulk formulas assume cavity dimensions much larger than the typical wavelength $hc/(k_{B}T)$, so the electromagnetic mode density in a box can replace individual cavity modes by a continuum. ^[ref]

The crucial constraint is whether photon emission and photon absorption, or other number-changing processes, can equilibrate the photon population. When they can, photon number is not an independently conserved quantity and the equilibrium photon chemical potential is ${\mu }_{\gamma }=0$. For each polarization and frequency, the mean occupation number then follows the Bose-Einstein distribution:If number-changing processes are ineffective but number-conserving interactions redistribute energy, a constrained Bose-Einstein distribution can instead have $f=[\exp ((h\nu -{\mu }_{\gamma })/(k_{B}T))-1{]}^{-1}$ with chemical potential ${\mu }_{\gamma }\le 0$ fixed by the number constraint. Thus thermal equilibrium with respect to energy exchange alone does not establish ${\mu }_{\gamma }=0$. The zero value can also persist in an already prepared Planck photon distribution during suitable subsequent evolution; continual photon emission is not required to preserve it. ^[ref]

Combining $f$ with the two-polarization density of states, $g_{\nu }/V=8\pi {\nu }^2/c^3$, gives the Planck photon distribution and the spectral photon energy density:Here $B_{\nu }$ is the spectral radiance in Planck's law; $u_{\nu }$ is stored energy density per unit frequency, not outgoing power per area. Their different units are connected by the full-sphere solid angle and the speed of light $c$. The Boltzmann constant is $k_B$, and $h$ is the Planck constant. ^[ref]

The equilibrium equation of state is particularly simple:The radiation constant $a_R$ is the same constant called $a_{\mathrm{rad}}$ in black-body radiation; $\sigma$ is the Stefan-Boltzmann constant. An outline of the energy density calculation substitutes $x=h\nu /(k_{B}T)$ and uses ${\int }_0^{\infty }{x}^3/(e^x-1)dx={\pi }^4/15$. The radiation pressure follows because the angular average of the squared directional component of photon momentum is $1/3$ for isotropy; the thermodynamic derivation of blackbody radiation pressure supplies an independent route. These formulas describe the radiation contribution, not the cavity walls. The isotropic massless relation $p=u/3$ survives a different isotropic spectrum, whereas the coefficient in $u=a_R{T}^4$ requires the zero-chemical potential Planck photon distribution.

For a computed example with stipulated $T=2.725\mathrm{K}$, the photon energy density is $u\simeq 4.17174\times 10^{-14}\mathrm{J}{\mathrm{m}}^{-3}$ and the radiation pressure is $p\simeq 1.39058\times 10^{-14}\mathrm{Pa}$. The photon number density in the International System of Units isHere $\zeta$ is the Riemann zeta function. The mean photon energy differs from the energy at the peak of a spectral curve: counting photons and weighting them by energy emphasize different frequencies. The calculations use the exact International System of Units values of the Planck constant, Boltzmann constant and speed of light. ^[ref]

A motivating observation is the cosmic microwave background spectrum measured by the Far Infrared Absolute Spectrophotometer (FIRAS) aboard COBE. Its apparatus used a polarizing Michelson interferometer, with one input viewing a roughly seven-degree sky patch and the other an internal blackbody. Scanning mirrors changed the optical path length difference; heat-sensitive detectors recorded the varying power, and Fourier transform spectroscopy recovered the spectral difference. Operation outside the terrestrial atmosphere removed atmospheric emission and absorption, while shielding and cooling the optics in a helium cryostat suppressed the instrument's own radiation. ^[ref]

Absolute calibration came from moving an external blackbody into the sky aperture, filling the same beam used for the sky. This was an iron-loaded epoxy absorber shaped as a re-entrant cone. Repeated reflections made it a close approximation to a blackbody, while thermometers supplied its temperature. Mather and colleagues explicitly evaluated reflections, temperature gradients, radiation leaking through the material and leakage around its edge. These are limits on how faithfully the calibrator realizes thermal equilibrium radiation, rather than details that can be eliminated merely by calling it black. ^[ref]

During sky observations the internal blackbody temperature was adjusted to nearly null the spectral difference, reducing gain and drift errors. Fixsen and colleagues' 1996 analysis used calibrated spectra over wavenumbers $2$ to $21\mathrm{c}{\mathrm{m}}^{-1}$, about $60$ to $630\mathrm{GHz}$. They separated the sky's spatially uniform component, its dipole and Galactic emission, and varied sky cuts and foreground templates as controls. The resulting cosmic microwave background spectrum agreed with a single Planck function, with reported RMS deviations below 50 parts per million of its peak brightness. This is a residual statistic over the analyzed spectrum, not a pointwise relative error bound at every frequency. ^[ref]

The observation motivates using a photon gas as a physical model rather than merely a formal counting exercise. Its interpretation has a clear boundary: agreement with Planck's law constrains departures from an equilibrium spectral shape within the calibrated band and foreground assumptions; it does not directly measure radiation pressure or entropy density of a photon gas. Nor does it establish that today's cosmic microwave background is continually re-equilibrating with matter. A previously thermal spectrum can retain its Planck function shape as expansion of the universe lowers every photon frequency in the same proportion. The usefulness of the model is that one measured spectral temperature, together with its equilibrium assumptions, predicts energy density, radiation pressure, photon number density and entropy density of a photon gas consistently.

In an orthonormal tetrad, $J^{\hat{0}}=\int f{d}^{3}p$ is the photon number density and $J^{\hat{i}}=\int f{e}^i{d}^{3}p$ is its number flux. Integrating the Boltzmann equation gives ${\nabla }_{\mu }{J}^{\mu }=\int C[f]d^{3}p/E$ when momentum-boundary terms vanish. Collisionless propagation and photon-number-conserving scattering give zero divergence. Photon emission, photon absorption, electron-positron annihilation in cosmology, bremsstrahlung and double Compton scattering generally do not. This is a number current, not an electric charge current: photons are neutral.

For hydrogen ionization fraction $X_e$ of order one, the hydrogen-only Saha equation gives $1\sim n_b(2\pi /(m_{e}T){)}^{3/2}{e}^{B/T}$. The baryon-to-photon ratio and photon number density imply $n_b\sim \eta T^3$. Therefore $1\sim \eta (T/m_e{)}^{3/2}{e}^{B/T}$ and the displayed logarithmic balance follows. A small $\eta$ supplies many ionizing photons per baryon, so neutral hydrogen forms at a temperature much smaller than its Hydrogen binding energy. The numerical temperature depends on the ionization-fraction convention and on whether equilibrium remains valid.