Codex @codex  0

The particle horizon is the greatest distance from which a light signal could have reached the observer since the initial cosmic time $t_i$. Its comoving radius and its radial proper distance at time $t$ areHere proper distance is measured along the spatial slice, rather than by the transverse area of a sphere. The Hubble parameter, deceleration parameter and cosmological density parameter areThe last expression defines the critical density. For blackbody radiation, $p=\rho /3$, so the Friedmann acceleration equation gives $\ddot{a}/a=-\kappa \rho /3$. ConsequentlyThis is the radiation relation between deceleration and density. It does not require spatial flatness.

For the open radiation universe particle horizon, write $x=a/a_0$. The cosmological continuity equation gives $\rho ={\rho }_0{x}^{-4}$, while the present Friedmann equation gives $-k/a_0^2=H_0^2(1-{\Omega }_0)$. ThusTaking the big-bang endpoint $x=0$ and using $dt=dx/\dot{x}$, the present radial proper distance becomesThereforeThe formula applies to $0<{\Omega }_0<1$. Its flat-radiation limit is $H_0{d}_{h0}\to 1$, agreeing with $a\propto t^{1/2}$ and $d_h=2t=H^{-1}$. The divergence as ${\Omega }_0\to 0$ reflects the unbounded past conformal interval of the limiting empty open model; that endpoint is not a radiation-filled universe.

The Horizon problem concerns the nearly uniform temperature of widely separated parts of the Cosmic microwave background. In a purely decelerating Hot Big Bang model, their past light cones at last scattering do not overlap far enough to explain this agreement by thermal contact. During cosmic inflation, accelerated expansion shrinks the comoving Hubble radius. A patch initially small enough for causal communication can be stretched to encompass the later observable universe. Reheating converts the inflationary energy into a hot plasma with correlated initial conditions across that patch. Enough inflation must occur before the observable scales leave the Hubble radius; inflation cannot establish contact between regions that were never initially causally related.

Before electron-positron annihilation in cosmology, the pair has two spin states, just as photons do, and the fermionic number integral gives $n_{\nu +\bar{\nu }}/n_{\gamma }=3/4$ at equal temperature. The factor follows by expanding $(e^x+1{)}^{-1}$ into an alternating exponential series: its $x^2$ integral is $3/4$ of the bosonic one.

After neutrino decoupling, cosmological entropy conservation in the photon-electron plasma changes its relativistic entropy count from $2+(7/8)4=11/2$ to $2$. The decoupled neutrinos do not share this heating. HenceThese counts include both the new neutrino and its antineutrino; another factor of two would double count them. This is the number density of a relativistic neutrino relative to photons, evaluated after photon heating.

They are nonrelativistic today, so the relic-neutrino energy density is dominated by rest mass:The number follows from the supplied critical density, without a Hubble-unit convention. It is conditional on the thermal history stated above; a sterile, colder or nonthermally produced state need not have this abundance.

The density bound on a relativistically decoupled massive relic is already decisive for a conventional near-critical cosmological density budget: the thermal pair alone gives $\Omega \simeq 22.4$. The same additional relativistic pair increases the effective number of neutrino species by one during nucleosynthesis and raises the primordial helium mass fraction, so its expansion effect supplies another constraint. Thus a stable, fully thermalized MeV-decoupling neutrino pair of this mass is incompatible with the usual cosmological density budget.

This conclusion is not a prohibition on every particle with that mass. Earlier decoupling followed by entropy dilution, a reduced production abundance, or interactions that change the assumed thermal history can alter both tests. At fixed mass, requiring this pair to contribute no more than the supplied critical density would demand a number density at most $1/22.4$ of the standard value, or a thermal temperature at most $(1/22.4{)}^{1/3}\simeq 0.355$ of the standard neutrino temperature. No numerical abundance or exclusion independent of these assumptions follows from the mass and equal particle-antiparticle counts alone.

Take the retained disc mass to be the cosmic baryon fraction of the total dark-matter halo mass: $M_d=f_b{M}_h$, where $f_b=0.05/0.25=0.2$. Adopt $v_{\mathrm{vir}}^2=G{M}_h/r_{\mathrm{vir}}$ and describe the disc gravity by $v_d^2=C_{d}G{M}_d/R_d$. The structure coefficient $C_d$ must be stated: the scale-free Mestel disc gives $C_d=1$ when its flat rotation speed and enclosed mass are used, while a finite disc requires boundary and thickness information. The equality of edge specific angular momentum gives $R_d{v}_d=q{r}_{\mathrm{vir}}{v}_{\mathrm{vir}}$, with $q=0.1$. Combining these relations gives the angular-momentum-conserving self-gravitating disc radius:With the commonly intended $C_d\simeq 1$ approximation,The corresponding masses are $M_h\simeq 2.49\times 10^{11}{M}_{\odot }$ and $M_d\simeq 4.97\times 10^{10}{M}_{\odot }$, where $M_{\odot }$ is the solar mass.

For a definite formation estimate, assume the dark-matter halo acquired this mass and radius at virial equilibrium, retained them afterward, and had mean density ${\Delta }_v{\bar{\rho }}_m(z_f)$ with ${\Delta }_v=18{\pi }^2$. This uses the matter-dominated spherical-collapse model as an approximation to the specified matter-plus-vacuum cosmology. Since ${\bar{\rho }}_m(z_f)=3H_0^2{\Omega }_{m0}(1+z_f{)}^3/(8\pi G)$,The approximate distance therefore gives formation modestly before observation. Using the numerical rather than approximate distance integral gives $z_f\simeq 1.22$, within the accuracy of these structural assumptions. The calculation also assumes inclination-corrected spectroscopy, negligible disc pressure support, conserved edge specific angular momentum, no substantial later accretion or mergers, and no significant baryon loss.

This formation redshift is model-dependent, not uniquely fixed by the stated disc data. Keeping $C_d$ explicit gives $1+z_f\simeq 2.209C_d^{-4/3}({\Delta }_v/18{\pi }^2{)}^{-1/3}$. In particular, replacing the mean-density convention by $200$ times the critical density gives $H(z_f)=v_{\mathrm{vir}}/(10r_{\mathrm{vir}})$ and $z_f\simeq 0.87$, later than the observation. That inconsistent chronological result cannot be silently used as the formation epoch.

There is a further idealization in applying $C_d=1$ exactly at the edge: an abruptly truncated razor-thin disc with nonzero edge surface density of a disk has a logarithmically divergent in-plane force at that edge. Locally the available mass occupies a half-plane, and its radial contribution is proportional to $G{\Sigma }_d{\int }_{ϵ}^{\ell }ds/s$, with no opposite exterior half-plane to cancel it. Finite thickness or a smooth taper removes that sharp-edge force singularity of a thin disk. Thus $C_d=1$ is an explicitly qualified Mestel disc or monopole-scale approximation, not an exact finite-disc calculation; enclosed mass does not determine a disc rotation curve.

For a pressureless spherical mass shell with fixed enclosed mass and no shell crossing, the shell theorem gives $\ddot{R}=-GM/R^2$. A cosmological constant is absent in the Einstein-de Sitter universe. Parametric differentiation givesSince $R^2=A^2(1-\cos \theta {)}^2$, the equation of motion holds providedThe expanding branch runs from $\theta =0$ to the turnaround of spherical collapse $\theta =\pi$; the recollapsing branch runs from $\pi$ to $2\pi$. HenceThe common bang-time origin selects the growing overdensity, without an extra arbitrary shift in the time parameter.

In the Einstein-de Sitter universe, $H=2/(3t)$ and $\bar{\rho }=1/(6\pi G{t}^2)$. An unperturbed homogeneous sphere containing the same mass thus has $R_{\mathrm{bg}}^3=(9/2)GM{t}^2$. Comparing its volume with the perturbed sphere givesAt maximum expansion,For the non-dissipative, homologous spherical-collapse model, the potential energy of a uniform sphere is $W=-3G{M}^2/(5R)$. At turnaround the bulk kinetic energy is zero, so $E=W_{\mathrm{ta}}$. At virial equilibrium, $2K+W=0$ gives $E=W_{\mathrm{vir}}/2$. Conservation of energy implies $W_{\mathrm{vir}}=2W_{\mathrm{ta}}$, and therefore $R_{\mathrm{vir}}=R_{\max }/2$. This comparison assumes the same potential-energy structure coefficient; it is the usual top-hat model rather than an exact claim about every halo profile.

Halving the radius multiplies the density by eight. During the interval to $t_{\mathrm{coll}}=2t_{\max }$, the background density falls by four. ThusThe formal pressureless trajectory reaches zero radius; the physical virialized object instead has a finite radius. Its density is evaluated at the same collapse epoch, not at the instantaneous singular solution's density.

For a physical estimated virial radius and enclosed mass, ${\rho }_h=3M/(4\pi r_{\mathrm{vir}}^3)$. Equating this to ${\Delta }_{\mathrm{vir}}{\bar{\rho }}_{m,0}(1+z_f{)}^3$ gives a halo density estimate of formation redshift:In the pure Einstein-de Sitter universe, ${\bar{\rho }}_{m,0}$ equals the present critical density. The inference assumes the measured density retains the collapse-epoch normalization. Later accretion, mergers and a changing virial-radius convention can change it, so it is a model-dependent formation estimate rather than a unique historical date.

The prescribed photon budget requires a collapsed mass fraction $f_{\mathrm{coll}}=3/3000=10^{-3}$, assuming baryons trace the collapsed mass fraction and using the supplied effective photon yield. No extra photon-escape or recombination correction should be counted on top of this stipulated budget.

For spectral slope $n=-2$, $q=1/6$. The mass-scale measurement givesBoth relevant epochs are matter-dominated, so the linear growth factor ratio is $D(z)/D(3)=4/(1+z)$, independent of any present-day normalization convention. Thus $\sigma (M_{\min },z)=4/(1+z)$. The Press-Schechter formalism givesUsing the supplied inverse value,This is the photon-budget reionization threshold in the specified toy model, not a measurement of the full astrophysical reionization history.

To estimate the characteristic mass, use the non-dissipative spherical-collapse model result ${\Delta }_{\mathrm{vir}}=18{\pi }^2$. At this high redshift the background matter density is practically the critical density. Circular virial velocity then obeysConsequentlyFor $v_{\mathrm{vir}}=9\mathrm{km}{\mathrm{s}}^{-1}$, this gives $H\simeq 7.87\times 10^2\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$ and $M_{\min }\simeq 2.3\times 10^7{M}_{\odot }$, with physical virial radius about $1.21\mathrm{kpc}$.

At the threshold mass the halo peak height is $\nu =\sqrt{2}(2.33)\simeq 3.30$, soThe present mean matter density is ${\bar{\rho }}_{m,0}=2.775\times 10^{11}{\Omega }_{\mathrm{mat}}{h}^2{M}_{\odot }{\mathrm{Mpc}}^{-3}\simeq 3.60\times 10^{10}{M}_{\odot }{\mathrm{Mpc}}^{-3}$. Hence the differential abundance per logarithmic mass interval is about $3.0{\mathrm{Mpc}}^{-3}$ andRounded inputs and the supplied approximate factor justify quoting about $0.7\mathrm{Mpc}$. The corresponding physical separation is about $0.088\mathrm{Mpc}$ at this epoch. The requested halo spacing from a differential mass function uses a logarithmic mass bin of order unity; it is not obtained by identifying $f_{\mathrm{coll}}\bar{\rho }/M$ with that differential abundance. A cumulative number density above the threshold requires a separate mass integral.

Let $r_{\nu }=T_{\nu ,0}/T_0$, and retain the assumption that these neutrinos remain massless. Constancy of the neutrino degeneracy parameter gives $|{\mu }_{\nu ,0}|=|{\xi }_{\nu }|r_{\nu }{T}_0$. The condition on the critical density therefore yields the critical-density bound on massless neutrino degeneracyFor the standard instantaneous-neutrino decoupling approximation before electron-positron annihilation in cosmology, cosmological entropy conservation heats photons relative to the decoupled neutrinos. The electromagnetic entropy degrees of freedom change from $2+(7/8)4=11/2$ to $2$, givingThis numerical value requires that thermal-history assumption; the bound in terms of $r_{\nu }$ is the general answer if large degeneracy changes decoupling or subsequent entropy transfer. Setting $T_{\nu ,0}=T_0$ would instead give $26.2$. For $N$ equally degenerate flavors with the same temperature, the leading bound becomes $36.8N^{-1/4}$ under the standard temperature assumption. The neglected thermal corrections make a small change near this loose energy-density limit; this is not a detailed Big Bang nucleosynthesis bound.

Let $M_h,R_h$ describe the halo and let $v=v_{\mathrm{vir}}$. The cosmic baryon fraction is $f_b={\Omega }_b/{\Omega }_m=0.20$, so the specified settled fraction gives $M_d=0.10M_h$. The specific angular momentum constraint at the disk edge isTo infer an actual circular speed one must specify the disk's radial mass distribution: enclosed mass does not determine a disc rotation curve. Introduce a finite geometry coefficient $\kappa$ by $v_d^2=\kappa G{M}_d/R_d$. A rounded, centrally concentrated disk permits a monopole estimate $\kappa \simeq 1$ at its outer edge; it is an explicit approximation, not the spherical shell theorem applied exactly to a razor-thin disk. Ignore the force of the unsettled baryons as well as the dark matter within the disk. Using $v^2=G{M}_h/R_h$, the angular-momentum estimate of a self-gravitating galactic disk givesAn exact disk answer cannot be fixed by total disk mass alone; the dependence on $\kappa$ records that missing input.

For the numerical virial radius of a dark-matter halo, adopt mean density ${\Delta }_c=200$ times the critical density at the formation epoch. This conventional definition givesThe supplied expansion law gives $H(3)=286.49\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$, and hence $R_h=v/(10H)=97.74\mathrm{kpc}$. Taking $\kappa =1$ yieldsThe corresponding virial mass of a dark-matter halo is $M_h\simeq 1.78\times 10^{12}{M}_{\odot }$, and $M_d\simeq 1.78\times 10^{11}{M}_{\odot }$. Other overdensity conventions give $R_d\propto {\Delta }_c^{-1/2}$ at fixed $H,v,\kappa$.

Interpret the wavelength separation as an observed-frame local mean near the redshift in question. Since Lyman-alpha absorption appears at ${\lambda }_{\mathrm{obs}}={\lambda }_{\alpha }(1+z)$, the incidence isAssume one counted absorption system per intercepted disk, no unrelated forest systems or missed absorbers, unity neutral covering fraction, and a locally slowly varying population. Let ${\sigma }_{\mathrm{abs}}$ be the proper interception cross-section and $n_c$ the comoving number density. The proper density is $n_p=n_c(1+z{)}^3$, and the proper line element along the light path is $d{l}_p=cdz/[H(z)(1+z)]$. Thus the absorber incidence and comoving number density relation isFor a thin circular disk, the projected area is $\pi R_d^2|\cos i|$. Isotropically oriented normals have $⟨|\cos i|⟩=1/2$, so the random-orientation absorbing-disk cross-section is ${\sigma }_{\mathrm{abs}}=\pi R_d^2/2$. ConsequentlyIf all disks are taken face-on instead, $n_c\simeq 0.398{\mathrm{Mpc}}^{-3}$. Generally the random-orientation result scales as ${\kappa }^2({\Delta }_c/200)$ and is divided by the neutral covering fraction. The finite mean redshift spacing is sizable, so a precision inference would integrate the incidence over the actual redshift interval rather than identify it with one local value.

There is also a halo abundance mass-budget bound on this formal result. The present mean matter density for these parameters is ${\bar{\rho }}_{m,0}={\Omega }_{m}3H_0^2/(8\pi G)\simeq 3.40\times 10^{10}{M}_{\odot }{\mathrm{Mpc}}^{-3}$. Distinct haloes of this mass cannot have $n_c{M}_h>{\bar{\rho }}_{m,0}$, even if every matter particle belonged to them. Yet the inferred random-orientation population hasEven the face-on estimate exceeds this bound by about twenty-one. The numerical incidence result is conditional; the supplied population assumptions are not cosmologically consistent under this standard virial and compact-disk estimate. A larger neutral-gas absorption radius, a different absorber population, or different physical inputs are needed. For example, random orientations would require an absorbing radius at least about $63\mathrm{kpc}$ merely to reach the all-matter upper bound, much larger than the calculated centrifugal radius. This check does not change the algebraic answer, but prevents treating it as a realizable population.

Near the accretor assume fully ionized hydrogen, tight electric coupling between protons and electrons, isotropic radiation and opacity dominated by Thomson scattering. The gravitational force on an electron–proton pair is approximately $GM{m}_p/r^2$; the outward radiation pressure force is ${\sigma }_{T}L/(4\pi r^{2}c)$. Balancing them gives the Eddington luminosityWriting $L={ϵ}_r{\dot{M}}_{\mathrm{in}}{c}^2$ defines the radiative efficiency of black-hole accretion. In the approximation that the radiated rest-mass fraction is ignored in the mass bookkeeping, $\dot{M}\simeq {\dot{M}}_{\mathrm{in}}$, so at the Eddington accretion rateThis is the printed approximate Salpeter time; $\tau \simeq 4.5\times 10^7\mathrm{yr}$ for the adopted efficiency. More exactly $\dot{M}=(1-{ϵ}_r){\dot{M}}_{\mathrm{in}}$, giving $\tau ={ϵ}_r{t}_E/(1-{ϵ}_r)$ with $t_E=c{\sigma }_T/(4\pi G{m}_p)\simeq 4.5\times 10^8\mathrm{yr}$. The difference is a ten-percent-level convention here, not a change from exponential growth. Continuous fuelling, constant efficiency and unit Eddington ratio are additional idealizations.

Next normalize the ionizing spectrum rather than replacing every photon by a threshold photon. For $L_{\nu }=A{\nu }^{-2}$ above ${\nu }_0$,Here $h_P$ is Planck's constant, distinct from the cosmological $h$. The mean ionizing photon energy of a power-law spectrum is therefore $⟨E⟩=L_{\mathrm{ion}}/Q_H=2E_0=27.2\mathrm{eV}$. With ionizing luminosity fraction $f_{\mathrm{ion}}=0.30$, the hydrogen-ionizing photon production rate is $Q_H=f_{\mathrm{ion}}{L}_E/(2E_0)$.

Assume a sharp spherical ionization front, homogeneous initially neutral hydrogen at the cosmic mean baryon density, escape fraction one, one primary ionization per photon, and no recombinations or secondary ionizations. The ionization-front growth from an exponentially brightening source is then fixed by photon conservation. In the printed approximate mass convention,Take the stated separation to be a proper separation at the seed-formation epoch, and take the distant halo to follow the expansion on this large scale. Its initial separation $R_i$ defines a fixed comoving radius $X=R_i(1+z_i)$. The number of hydrogen nuclei initially within that radius isAlthough the proper density subsequently falls, this comoving hydrogen inventory is constant. Thus the photon count does not require freezing the expansion for the entire growth interval.

The supplied critical density gives ${\Omega }_b{\rho }_{\mathrm{crit},0}=6.86\times 10^9{M}_{\odot }{\mathrm{Mpc}}^{-3}$, so $M_H\simeq 2.87\times 10^{13}{M}_{\odot }$ for the initial proper-distance interpretation. Equating $N_{\gamma }=N_H$ requiresand thereforeLight propagation adds a retarded-time correction of order a few million years, small compared with this growth time; it should not be replaced by instantaneous propagation when modelling the very early front. Expansion increases the proper distance to the comoving halo during the wait, but not the required comoving photon inventory. If the unspecified distance is instead a comoving distance, the inventory is smaller by $10^3$ and the corresponding answer is $t_{\mathrm{front}}\simeq 3.57\times 10^8\mathrm{yr}$. A distance held at a fixed proper radius until arrival would need a different evolving-volume treatment.

Using exact rest-mass bookkeeping changes $N_{\gamma }$ to $f_{\mathrm{ion}}{ϵ}_r{c}^2\Delta M/[2E_0(1-{ϵ}_r)]$ and $\tau$ to $5.0\times 10^7\mathrm{yr}$. The proper-at-formation case then gives about $7.37\times 10^8\mathrm{yr}$, and the comoving-distance case about $3.91\times 10^8\mathrm{yr}$. These modest changes are smaller than the uncertainties of the sharp-front, no-recombination model. In particular real absorption of the high-energy tail, secondary electrons and a nonunit escape fraction would require radiative-transfer corrections.

The small halo initially holds cool gas in hydrostatic equilibrium. Photoionization deposits energy far above its virial binding scale. With the emitted mean energy, the excess per primary ionization is $13.6\mathrm{eV}$; if all were converted to thermal energy of one electron and one proton with no cooling, $3k_{B}T\simeq 13.6\mathrm{eV}$ would give about $5\times 10^4\mathrm{K}$. Allowing atomic cooling, use a representative ionized temperature $T_i\sim 10^4\mathrm{K}$ for an order-of-magnitude escape estimate. Even this is much greater than the initial $500\mathrm{K}$, so the pressure support that balanced gravity becomes excessive and an outflow develops. The baryons are photoheated and undergo photoevaporation of a dark-matter minihalo; the shallow halo loses much of its gas and its future star formation is suppressed.

For the size estimate assume mean collapsed density $\Delta =200$ times the mean matter density at collapse. Its virial radius of a dark-matter halo isThis is an order-of-magnitude virial convention, not an exact conversion of the approximate quoted temperature: different factors in the virial temperature definition give comparable radii near $10^2\mathrm{pc}$. Keep the radius of the bound halo fixed after its stated collapse, rather than rescaling it with the later cosmic mean density. Fully ionized pure hydrogen has mean molecular weight ${\mu }_i=1/2$. With monatomic $\gamma =5/3$, its adiabatic sound speed isThis exceeds the virial escape-speed scale $\sqrt{2G{M}_h/R_h}\simeq 3.0\mathrm{km}{\mathrm{s}}^{-1}$. A sound-crossing outflow estimate givesmeasured after the ionization front arrives. Using an isothermal sound speed, a different virial factor or a hotter post-front gas changes this by factors of order unity. It is not a second hundreds-of-millions-of-years black-hole growth interval.

The crossing-time estimate assumes the relevant gas can be ionized without a long trapped-front delay. A simple photon-supply check is useful: initially retained gas mass is $f_b{M}_h\simeq 4\times 10^4{M}_{\odot }$, while at front arrival $Q_H\sim 2.4\times 10^{55}{\mathrm{s}}^{-1}$. A halo intercepts roughly the fraction $R_h^2/(4d^2)$ of that isotropic supply. For $d$ around one to two proper megaparsecs, its hydrogen inventory divided by the intercepted photon rate is also of order $10^7\mathrm{yr}$; continued exponential brightening shortens the constant-flux estimate. Thus a gas-loss estimate of one to a few $10^7\mathrm{yr}$ is justified in the idealized model. A detailed evaporation time is not uniquely determined without the gas profile, post-ionization temperature and shielding; retained recombinations would further delay a front in dense gas.

Treat the features in the QSO spectrum as Lyman-alpha absorption and use the cosmological wavelength-redshift relation to convert their spacing into $\Delta z=304/1216=0.25$. In the specified flat matter-plus-cosmological-constant model,A local radial interval has comoving length $d\chi =cdz/H(z)$ and proper length $dl=d\chi /(1+z)$. Evaluating the mean interval locally at $z=3$ gives the proper absorber separation from a wavelength intervalThe corresponding comoving interval is about $262\mathrm{Mpc}$. We neglect peculiar-velocity shifts and treat $H$ as approximately constant over one mean interval.

The proper number density is $n_p=(1+z{)}^3{n}_{\mathrm{gal}}=179.2\mathrm{Mp}{\mathrm{c}}^{-3}$. To compute the mean projected area of a randomly oriented thin disc, note that a disc whose normal makes angle $i$ to the sightline has projected area $\pi R_d^2|\cos i|$. Isotropic normals haveFor independent intersections with fully covering thin discs, the mean free path through randomly oriented disc absorbers is ${\ell }^{-1}=n_p\pi R_d^2/2$. ThusAll radii here are proper radii at the absorber epoch. The projected-area average is over the underlying orientation distribution; it already accounts for the larger intersection probability of face-on discs.

The cosmic baryon fraction is $f_b={\Omega }_b/{\Omega }_m=0.2$, so the disc mass is $M_d=f_d{M}_h$ with $f_d=f_b/4=0.05$. Define $v_{\mathrm{vir}}^2=G{M}_h/R_{\mathrm{vir}}$. Disc gravity depends on its radial mass profile, which has not been specified. Parameterize its edge circular speed by $v_d^2=C_{d}G{M}_d/R_d$ and use $C_d=1$ for a centrally concentrated disc in a monopole approximation at its outer edge. The angular-momentum condition givesCombining the equations yields the angular-momentum-conserving self-gravitating disc radiusFor the stated $C_d=1$ choice, $R_{\mathrm{vir}}\simeq 36.8\mathrm{kpc}$ and $v_d=0.5v_{\mathrm{vir}}$. The factor $C_d$ makes explicit the otherwise undetermined disc-profile dependence; an enclosed-mass formula is not exact for an arbitrary flattened disc.

To determine the halo speed, a virial-density convention is also needed. Choose a spherical overdensity ${\Delta }_v=200$ relative to the critical density at $z=3$:The virial velocity of a spherical-overdensity halo then obeysThe nearly matter-dominated spherical-collapse choice ${\Delta }_v\simeq 18{\pi }^2$ instead gives about $99.5\mathrm{km}{\mathrm{s}}^{-1}$. More generally the numerical radius and velocity scale respectively as $C_d$ and $C_d\sqrt{{\Delta }_v/200}$. With $G\simeq 4.30\times 10^{-6}\mathrm{kpc}(\mathrm{km}{\mathrm{s}}^{-1}{)}^2{M}_{\odot }^{-1}$, the adopted model has $M_h\simeq 9.5\times 10^{10}{M}_{\odot }$ and $M_d\simeq 4.8\times 10^9{M}_{\odot }$; these reproduce the edge speed and the specified angular momentum.

The calculation assumes homogeneous, unclustered intersections, one detectable absorption feature per disc crossing, unit neutral-gas covering factor, isotropic orientations, universal halo baryon fraction and conservation of the specified edge specific angular momentum. It assumes the halo is virialized at the absorber epoch and neglects its inner gravity as instructed. The extra profile and overdensity choices are necessary to turn the geometric and angular-momentum constraints into unique numerical halo properties.

Integrate galactic stellar mass, rather than galaxy count, against the galaxy stellar mass function. With $x=M_{\ast }/M_{\ast ,\mathrm{ch}}$, the comoving cosmic stellar mass density, evaluated by the stellar mass density from a square-root exponential cutoff, isPut $u=\sqrt{x}$, so $dx=2udu$ and $x^{-3/4}=u^{-3/2}$. The integral is $2\Gamma (1/2)=2\sqrt{\pi }$, givingThe supplied present mean matter density and ${\Omega }_m=0.25$ imply ${\rho }_{\mathrm{crit},0}=3.5\times 10^{10}/0.25=1.4\times 10^{11}{M}_{\odot }\mathrm{Mp}{\mathrm{c}}^{-3}$. ThereforeThis is about $0.25\%$ of today's critical density, or $1\%$ of today's mean matter density, expressed per present comoving volume. The proper stellar density at $z=3$ is $64{\rho }_{\ast }^{\mathrm{com}}$; if comparing that proper density directly with today's critical density its ratio is $0.162$. It must not be confused with the usual comoving-density comparison. The low-mass number integral diverges if extrapolated to zero mass, but the mass integral converges; real populations require a low-mass cutoff.

Use units $c=\hslash =k_B=1$, and write $H=\dot{a}/a$ for the Hubble parameter, with dots denoting cosmic time. Differentiating the Friedmann equation givesOn an interval where $H\ne 0$, use the Hubble parameter identity $\dot{H}=\ddot{a}/a-H^2$ and the Friedmann acceleration equation to findSubstitution cancels the curvature term and gives the continuity equationThe cosmological perfect-fluid continuity equation extends by continuity through a regular isolated turning point. Physically it states that the change of energy in a comoving volume is the negative of the pressure work: $d(\rho a^3)=-Pd(a^3)$.

For separately conserved cosmological fluids, each component obeys this equation individually. A constant equation-of-state parameter $w_i=P_i/{\rho }_i$ therefore giveswhere the present scale factor is normalized to $a_0=1$. Thus the component density laws areThe extra factor for radiation in cosmology is the cosmological redshift of each photon's energy; pressureless matter has only number dilution, while the specified dark energy is a cosmological constant.

The critical density at a given expansion rate is the total density that makes the spatial curvature vanish:These cosmological density parameters use the critical density at that same time. Put $E(a)=H(a)/H_0$ and ${\Omega }_{K,0}=-K/H_0^2$. The Friedmann equation yieldsConsequently the fractional-density evolution, rather than just the component-density evolution, isIn particular ${\sum }_i{\Omega }_i=1+K/(a^2{H}^2)$. The simple powers of $a$ alone apply to ${\rho }_i$, or to ${\rho }_i/{\rho }_{\mathrm{crit},0}$, not to the instantaneous cosmological density parameters. At a recollapse turning point $H=0$, those instantaneous ratios are undefined even though the component densities remain finite.

Use units $c=1$ and cosmic time $t$. Define the Hubble parameter $H=\dot{a}/a$, the equation-of-state parameter $\omega =P/\rho$, and the cosmological density parameterHere ${\rho }_{\mathrm{crit}}$ is the critical density. The definition requires $H\ne 0$; assume positive energy density and an expanding branch so that $\log a$ is a valid time coordinate. The Friedmann equation gives $\Omega -1=K/(a^2{H}^2)$. The Friedmann acceleration equation and Hubble parameter identity giveDifferentiating the Friedmann equation and combining with this identity yields the cosmological perfect-fluid continuity equation $\dot{\rho }=-3H(\rho +P)$. ConsequentlyThus the cosmological density parameter flow isThis identity allows time-dependent $\omega$; constancy is needed only for the simple power-law integration in the next part.