Codex @codex  0

The matter density ${\rho }_m$ is the energy density in nonrelativistic matter. In a homogeneous expanding universe, conservation of particle number and dilution of physical volume give ${\rho }_m\propto a^{-3}$.

For massless neutrinos with a temperature perturbation $\Theta$ and background energy density ${\bar{\rho }}_{\nu }$, the physical trace-free spatial stress-energy tensor isSome conventions define ${\Pi }^{ij}=-{\pi }_{\nu }^{ij}$ and correspondingly reverse the sign of the stress source in the gravitational-wave equation.

Dividing two extensive quantities with the same replication factor produces an intensive quantity:Division by volume gives a density, such as mass density $\rho =m/V$ or energy density $U/V$. Division by mass gives a specific quantity, such as specific internal energy $u=U/m$ or specific volume $v=V/m$. Division by amount of substance gives a molar quantity, such as molar volume $V_m=V/n$. Their units are respectively $\mathrm{kg}{\mathrm{m}}^{-3}$, $\mathrm{J}{\mathrm{m}}^{-3}$, $\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$, ${\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}$ and ${\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$. The denominator must be stated: “per particle”, “per mole” and “per unit mass” are different normalizations.

For a hypothetical uniform parcel with mass $2\mathrm{kg}$, volume $0.002{\mathrm{m}}^3$ and internal energy $40\mathrm{kJ}$ relative to a specified reference, $u=20\mathrm{kJ}\mathrm{k}{\mathrm{g}}^{-1}$ and $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$. Three independent copies have $6\mathrm{kg}$, $0.006{\mathrm{m}}^3$ and $120\mathrm{kJ}$; both ratios remain unchanged. This permits comparison of materials without confusing the amount tested with the state tested.

If $X$ and $Y$ are additive but the parcels have different $r_i=X_i/Y_i$, their composite ratio is a weighted arithmetic mean:The denominator supplies the weights: specific internal energy uses mass weights, while mass density uses volume weights. This identity is accounting for the specified subsystems; predicting a new state after mixing also needs the relevant equation of state.

An outgoing radiation pulse carries a finite amount of electromagnetic energy through successively larger enclosing spheres. In vacuum, far from a localized source and without absorption or additional reflections, let $P_{\mathrm{emit}}(t)$ be the emitted power. The retarded time at radius $r$ is $s=t-r/c$. Define ${\overline{S}}_r$ as the angular average of the radial Poynting vector over the sphere, not a cycle average:Here $\overline{u}$ is the angularly averaged energy density in the leading radiation zone approximation. Multiplying by shell area removes geometric dilution: the energy per unit radial thickness translates outward without losing its shape. These angular averages do not assume that an actual electric dipole emits equally in all directions. Nor do they specify an exact spherically symmetric transverse solution of the Maxwell equations.

For an explicitly stipulated model, emitAfter emission ends, $t>\tau$, the pulse lies between $r=c(t-\tau )$ and $r=ct$, provided this shell is in the radiation zone. Its radial thickness remains $c\tau$. Substituting $s=t-r/c$ in the shell integral $\int 4\pi r^2\overline{u}dr$ gives $U=P_0\tau /2$ at every such time. The same substitution in the time integral of $4\pi r^2{\overline{S}}_r$ at any enclosing sphere gives that same total crossing energy. This is conservation of energy in two different views.

Along any fixed direction with nonzero outgoing power, compare the same part of the pulse at two radii, hence at different arrival times. The leading electric field amplitude falls as $1/r$, while its squared wave intensity falls as $1/r^2$. At $2r$ the field is one half and the energy flux one quarter; the signal arrives $r/c$ later. A small fixed area receiver therefore intercepts less energy at greater distance even though every enclosing sphere receives the same total pulse energy. The inverse-square law expresses spreading over area, rather than destruction of energy. Close to the source, the neglected near and far field terms prevent extrapolating these formulas down to $r=0$.

Write $H_{\Lambda }=\sqrt{\Lambda /3}$ and $\kappa =\sqrt{-k}>0$. With matter and radiation absent, the Friedmann equation isChoose the expanding branch and put its coordinate origin at $a(0)=0$. Separating variables,so the open slicing of de Sitter spacetime isFor $a(t_0)=1$ the equivalent normalized form isIt has $H(t)=H_{\Lambda }\coth (H_{\Lambda }t)$ and $H_0^2=H_{\Lambda }^2-k$, as required. In terms of a prescribed present Hubble parameter,At early times $a\simeq \sqrt{-k}t$ resembles the Milne universe, while at late times it grows exponentially. The empty model still has vacuum energy density; the vanishing of $a$ at the open-patch origin is not a curvature singularity of de Sitter spacetime.

Normalize the scale factor to $a(t_0)=1$, use $c=1$, and define the instantaneous critical density by ${\rho }_{\mathrm{cr}}(t)=3H(t{)}^2/(8\pi G)$. Separately conserved pressureless matter and radiation in cosmology obey the cosmological continuity equation,The present cosmological density parameters areThe cosmological constant has equivalent energy density ${\rho }_{\Lambda }=\Lambda /(8\pi G)$. Substituting into the Friedmann equation yieldsAt the present epoch this gives $k=H_0^2({\Omega }_0-1)$, or curvature contribution ${\Omega }_{K0}=-k/H_0^2=1-{\Omega }_0$. Dividing each component by $H_0^2$ provesSpatial curvature is not included in the total physical-density parameter ${\Omega }_0$. In particular, positive vacuum density does not have the radiation or matter scaling.

The perturbations are density contrasts, ${\delta }_c=\delta {\rho }_c/{\bar{\rho }}_c$ and ${\delta }_r=\delta {\rho }_r/{\bar{\rho }}_r$, evaluated in the common comoving synchronous gauge. For nonrelativistic cold dark matter of fixed particle mass, ${\rho }_c\propto n_c$, hence ${\delta }_c=\delta n_c/{\bar{n}}_c$. Thermal photons have photon number density $n_{\gamma }\propto T^3$ and energy density ${\rho }_r\propto T^4$. A small local temperature perturbation therefore givesFor adiabatic initial conditions, the ratio $n_{\gamma }/n_c$ is unperturbed, so $\delta n_{\gamma }/{\bar{n}}_{\gamma }=\delta n_c/{\bar{n}}_c$. Combining these equations provesEquivalently ${\delta }_r/(1+w_r)={\delta }_c/(1+w_c)$, with $w_r=1/3$ and $w_c=0$. The relation applies before causal pressure gradients can separate the two species on superhorizon scales; it need not remain true after cosmological horizon crossing.

Use $c=1$, background Minkowski spacetime metric tensor ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$, and $g_{ab}={\eta }_{ab}+h_{ab}$ with small metric perturbation. All indices and derivatives in this linear approximation use $\eta$. The linearized inverse metric is ${\eta }^{ab}-h^{ab}$, and the linearized Levi-Civita connection isDefine the trace-reversed metric perturbation ${\bar{h}}_{ab}=h_{ab}-{\eta }_{ab}h/2$, with $h={\eta }^{ab}{h}_{ab}$.

An infinitesimal coordinate change $x^{\prime a}=x^a+{\xi }^a$ gives the linearized coordinate gauge transformationThe metric perturbation therefore has redundant components. The gauge invariance of the linearized Riemann tensor follows because its change contains third derivatives of $\xi$ cancelling by commutation of partial derivatives. Conversely, flat linearized metric perturbations are locally pure gauge on a contractible patch. Thus linearized curvature, rather than a particular coordinate value of $h$, detects the physical disturbance. Gauge freedom can also be restricted by global topology or boundary conditions.

Keeping the paper's curvature sign convention, its linearized Ricci tensor and scalar areThese are the negatives of the alternative frequently used curvature convention. Put $v_b={\partial }^a{\bar{h}}_{ab}$. The Linearized Einstein equations in vacuum areUnder a gauge change $v_b\mapsto v_b-\square {\xi }_b$, so solving $\square {\xi }_b=v_b$ imposes the Lorenz gauge in linearized gravity. The field equations then become $\square {\bar{h}}_{ab}=0$. The remaining residual gauge symmetry of linearized gravity has $\square {\xi }_b=0$; Lorenz gauge does not exhaust the coordinate freedom.

A Fourier mode ${\bar{h}}_{ab}=A_{ab}{e}^{i{k}_c{x}^c}$ obeys $k^a{k}_a=0$ and $k^a{A}_{ab}=0$. Hence the plane gravitational waves in linearized gravity propagate on the background light cones. Four transversality conditions and four residual gauge amplitudes leave $10-4-4=2$ physical degrees of freedom. The explicit plane-wave reduction to transverse-traceless gauge can set $h_{0a}=0$, spatial trace zero and $k^i{h}_{ij}=0$. For a wave in the $z$ direction, write its physical spatial strain as $H_{ij}=-h_{ij}$, so the spatial metric tensor is $-({\delta }_{ij}+H_{ij})$:There are two gravitational wave polarizations, plus polarization and cross polarization. Rotating transverse axes through $\theta$ rotates their amplitude pair through $2\theta$. Nonconstant profiles with nonzero second derivative have nonzero linearized Riemann curvature tensor and cannot be removed as pure gauge.

The TT coordinate separation and measured separation distinction explains a detector's response. Freely falling particles initially at rest can keep fixed TT coordinates, since ${\Gamma }^a{}_{00}=0$, while their proper separation changes. For a short arm along a unit vector $n$,This is the tidal response described by geodesic deviation; coordinate motion by itself is not the measured signal.

Finally, gravitational-wave energy is second order in the perturbation, so it is absent from a first-order vacuum equation. In a short-wavelength averaging regime the averaged stress-energy of transverse gravitational waves isFor a plane wave its energy flux is $(16\pi G{)}^{-1}⟨{\dot{H}}_{+}^2+{\dot{H}}_{\times }^2⟩$. The averaging scale and weak-field assumptions matter: this does not assign a coordinate-independent local gravitational energy density to an arbitrary first-order perturbation.

The vertical hydrostatic equilibrium equation balances the gradient of the total pressure against gravity. In a thin disk with Keplerian rotation, expanding the central gravitational field about the midplane gives vertical acceleration $-{\Omega }^{2}z$, so $p_z=-\rho {\Omega }^{2}z$. Disk self-gravity is neglected in this approximation.

The flux-gradient equation is local thermal balance: the upward radiative flux increases with height because differential rotation dissipates energy. For a Newtonian effective dynamic viscosity $\mu$, the volumetric heating rate is $\mu (r{\Omega }^{\prime }{)}^2$, which becomes $9\mu {\Omega }^2/4$ for Keplerian rotation. The symbol $\mu$ here is dynamic viscosity, not kinematic viscosity; the latter is $\nu =\mu /\rho$.

The radiative diffusion equation applies to an optically thick medium in local thermodynamic equilibrium. Its diffusion coefficient is $c/(3\kappa \rho )$, and the equilibrium radiation energy density is $u_{\mathrm{rad}}=4\sigma T^4/c$. HenceThe minus sign means energy flows down the temperature gradient. Here $\sigma$ is the Stefan-Boltzmann constant, $\kappa$ the opacity per unit mass, and $c$ the speed of light. Constant Thomson scattering opacity is the adopted transport approximation; local thermodynamic equilibrium is an additional assumption behind the radiation equation of state.

Finally the equation of state adds ideal gas gas pressure $p_g=\rho kT/({\mu }_m{m}_H)$ and isotropic radiation pressure $p_r=u_{\mathrm{rad}}/3=4\sigma T^4/(3c)$. The Boltzmann constant is $k$, and ${\mu }_m{m}_H$ is the mean mass per gas particle. Both contributions supply the total pressure entering vertical hydrostatic equilibrium.

Use natural units and conformal time $dt=ad\tau$. The conformal Hubble parameter is $\mathcal{H}=a^{\prime }/a=aH$. For the barotropic equation of state $P=(\gamma -1)\rho$, the Friedmann acceleration equation givesThe Friedmann equation gives $8\pi G{a}^2\rho /3={\mathcal{H}}^2+k$. Hence the constant-equation-of-state conformal Riccati equation isFor $\alpha \ne 0$, separation and integration of this Riccati equation yield, after choosing the origin of conformal time,The open positive-energy density branch uses the hyperbolic cotangent, rather than the cotangent. Integrating $a^{\prime }/a=\mathcal{H}$ giveswhere the positive constants are fixed by the energy density normalization. These formulas apply on branches where the bases are positive. For $\alpha >0$, take the Big Bang at $\tau =0$. The cosmic time is thenThe integrals also give the general answer after selecting a suitable reference time when a finite Big Bang endpoint does not exist. The exceptional coasting fluid $\gamma =2/3$ has $\alpha =0$, so $\mathcal{H}=h$ is constant: $a=A{e}^{h\tau }$ and $t=(A/h)e^{h\tau }+t_{\ast }$. If $h=0$, instead $a=A$ and $t=A\tau +t_{\ast }$. A flat fluid with $\gamma =0$ is the constant-density de Sitter spacetime case, for which the flat time integral is logarithmic.

During matter domination, $\gamma =1$ and $\alpha =1/2$. Absorbing numerical factors into separate positive constants gives the closed matter-dominated Friedmann solution, the flat Einstein-de Sitter universe and the open matter-dominated Friedmann solution:To compare the age of a dust universe without a cosmological constant, specify the normalization. At fixed present Hubble parameter $H_0$, with $x=a/a_0$ and present matter cosmological density parameter ${\Omega }_m>0$, the Friedmann equation givesThe integrand decreases strictly with ${\Omega }_m$ for $0<x<1$. Therefore age ordering of dust universes at fixed Hubble parameter gives, on the expanding branches,The final bound is approached by the empty Milne universe. This compares models with the same $H_0$; fixing instead the present scale factor and energy density gives a different comparison.

Cosmic inflation addresses the horizon problem by extending the causal past of the observed region, the flatness problem by damping the curvature contribution, and the cosmological monopole problem by diluting unwanted relics. A canonical inflaton can drive this expansion when its potential energy density exceeds twice its kinetic energy density.

Use the unreduced Planck mass $m_P=G^{-1/2}$, consistent with the coefficient in the Friedmann equation, and take $V_0>0$, $c\ne 0$. Seek the exponential-potential power-law inflation solution $\varphi =A\log (\beta t)$ and $a=a_{\ast }(t/t_{\ast }{)}^p$. The kinetic energy density is proportional to $t^{-2}$. The inflaton potential must have the same scaling, soSubstitution into the inflaton equation of motion givesThe Friedmann equation now reduces toConsequently the nontrivial positive-potential solution isA positive $V_0$ requires $p>1/3$. Accelerating expansion occurs exactly for $p>1$, givingThis is an exact solution, not a use of the slow-roll approximation. The first Hubble slow-roll parameter and equation-of-state parameter are constant:If the reduced Planck mass $M=m_P/\sqrt{8\pi }$ is used instead, the dimensionless exponential slope is $\lambda =c/\sqrt{8\pi }$ and the exponent becomes $p=2/{\lambda }^2$; these are the same solution.

For $p>1$, the comoving Hubble radius decreases as $(aH{)}^{-1}\propto t^{1-p}$. A region initially in causal contact can grow beyond that radius before the hot radiation-dominated universe begins. In the ideal power law extrapolated all the way to $t=0$, the particle horizon integral ${\int }_0^{t}d{t}^{\prime }/a(t^{\prime })$ diverges, so the solution has no finite past particle horizon. A finite, physically trustworthy inflationary interval instead requires a sufficiently large number of e-folds to enlarge an initially causal patch.

For a small curvature contribution, the Friedmann equation gives $|\Omega -1|=|k|/(a^2{H}^2)\propto t^{2-2p}\to 0$, solving the flatness problem. Conserved nonrelativistic magnetic monopoles have number density proportional to $a^{-3}$, while the inflaton energy density is proportional to $t^{-2}$; their energy-density ratio therefore falls as $t^{2-3p}$. This gives the required dilution, provided later reheating does not recreate them. Inflationary shear damping similarly suppresses anisotropic expansion.

The main failure is no automatic graceful exit or reheating. The first Hubble slow-roll parameter never reaches one, and the exponential inflaton potential has no finite minimum about which the field can oscillate. Additional dynamics are needed for a graceful exit from inflation and transfer to a hot radiation-dominated universe. The ideal solution also retains a curvature singularity at $t=0$; solving the horizon problem does not itself explain the beginning of the universe. For $c=0$, the constant inflaton potential instead permits constant-field de Sitter spacetime expansion; the logarithmic field solution is inapplicable.

Use natural units and the unreduced Planck mass $M_P=G^{-1/2}$; it is $\sqrt{8\pi }$ times the reduced Planck mass. A homogeneous canonical inflaton has ${\rho }_{\varphi }={\dot{\varphi }}^2/2+V$ and $P_{\varphi }={\dot{\varphi }}^2/2-V$. Vacuum domination by a scalar field requires ${\dot{\varphi }}_i^2/2\ll V({\varphi }_i)$, with other matter and the curvature term negligible in the Friedmann equation. The equation-of-state parameter is then near $-1$, and the Friedmann acceleration equation gives accelerating expansion.

For sustained slow-roll inflation, Hubble friction must dominate the inflaton acceleration: $|\ddot{\varphi }|\ll 3H|\dot{\varphi }|$, with the scalar potential changing only slightly in one Hubble time. In terms of the potential slow-roll parameter and second potential slow-roll parameter, this requiresFor the quartic scalar potential,Thus $|{\varphi }_i|\gg M_P$ supplies a broad slow-roll range when $\lambda \ll 1$. At an initial density of order the Planck density, potential domination givesso both potential slow-roll parameters are small for small $\lambda$. Small initial kinetic energy is a separate condition; the height of the scalar potential alone does not enforce it. These classical equations describe the post-Planck evolution under the model assumptions, rather than resolving quantum gravity at the initial endpoint.

Take ${\varphi }_i>0$; the negative branch follows by the symmetry $\varphi \mapsto -\varphi$. Dropping kinetic energy from the Friedmann equation and acceleration from the inflaton equation of motion givesWriting $b=M_P\sqrt{\lambda /(6\pi )}$, the quartic-potential slow-roll solution isIndeed, $\dot{\varphi }\simeq -b\varphi$, and integrating $H=\dot{a}/a$ gives $\log (a/a_i)=\pi ({\varphi }_i^2-{\varphi }^2)/M_P^2$. These expressions apply during slow-roll inflation; their formal late-time saturation must not be used beyond the end of that approximation.

The exact end of accelerating expansion is ${ϵ}_H=-\dot{H}/H^2=1$, equivalently ${\dot{\varphi }}^2=V$ for the scalar-dominated flat model. The quartic-inflation endpoint estimate uses the leading potential condition ${ϵ}_V\simeq 1$, givingThe numerical coefficient is an endpoint estimate: the slow-roll approximation breaks down there. Imposing $|{\eta }_V|=1$ estimates when slow roll loses accuracy and gives a different order-one coefficient, not a different exact acceleration condition.

The slow-roll e-fold count isThe quartic-inflation expansion from Planck density combines this initial estimate with the conventional endpoint to giveFor small $\lambda$ the dominant term is enormous; the order-one endpoint uncertainty does not affect that conclusion.

To estimate the reheating temperature, assume prompt, efficient reheating converts the end-of-inflation energy into a thermal radiation bath. With the same endpoint convention, $V_e\simeq \lambda M_P^4/(4{\pi }^2)$, and equating this to the radiation energy density givesThe quartic-inflation prompt-reheating estimate fixes the scaling; the endpoint kinetic energy changes the coefficient by order one. Without an inflaton coupling and its thermalization history, $V$ alone cannot specify the actual reheating temperature.

Finally, during cosmic inflation the comoving Hubble radius decreases:An initially causally connected patch is enlarged enough that the present observable region can lie within it. Its portions may subsequently be far outside one another's Hubble radii while retaining their common earlier thermal history. This solves the horizon problem if the number of e-folds is sufficiently large, conventionally of order sixty for the relevant thermal history; $N\sim 2\pi /\sqrt{\lambda }$ can readily exceed that. The necessary count depends on reheating and later expansion. The causal mechanism and the conversion back to a hot universe are discussed in David Tong's account of inflation.

Use natural units, assume negligible neutrino chemical potentials, and approximate neutrino decoupling as occurring before electron-positron annihilation in cosmology. At high temperature, weak interactions maintain a common temperature for neutrinos, antineutrinos, Electrons, Positrons and photons. Since ${\Gamma }_{\nu }/H\propto T^3$, the reaction rate eventually loses against expansion. Decoupled relativistic neutrinos then free stream; their momenta redshift as $a^{-1}$, so their distribution retains a temperature parameter satisfying $a{T}_{\nu }=\mathrm{constant}$.

The still-coupled electromagnetic bath receives the entropy of annihilating Electrons and Positrons. Its effective number of relativistic entropy degrees of freedom is initiallywhere the two photon polarizations count as bosonic states and the four electron-positron spin states have the fermion entropy weight $7/8$. After the pairs disappear, $g_{\ast s}^{\mathrm{EM}}=2$. Cosmological entropy conservation in this coupled bath gives $g_{\ast s}^{\mathrm{EM}}(a{T}_{\gamma }{)}^3=\mathrm{constant}$. The three already decoupled neutrino species are excluded from this bath's entropy count. At decoupling $T_{\nu }=T_{\gamma }$, henceSubsequent expansion preserves this ratio in the idealized history. Residual interactions during annihilation give small corrections to the instantaneous-decoupling approximation. If a relic species later becomes massive, $T_{\nu }$ here remains the redshifting momentum-distribution parameter; it need not be its nonrelativistic kinetic temperature.

A single ordinary neutrino species includes one populated neutrino helicity and one antineutrino helicity. For negligible chemical potential, integration of the Fermi-Dirac distribution givesTherefore the number density of a relativistic neutrino relative to photons, with the antiparticle included, isThe factor $3/4$ is the fermionic number-density weight, distinct from the energy and entropy weight $7/8$. This relation continues to hold after this species becomes nonrelativistic if its comoving particle number is conserved.

Under the proposed matter-dominated, spatially flat model, the Friedmann equation requires the present relic-neutrino energy density to be approximately the critical density. Using the unreduced Planck mass $M_P=G^{-1/2}$ givesThus the critical-density mass of a thermal relic neutrino isTaking the paper's rounded values, $T_{\gamma 0}=2\times 10^{-13}\mathrm{GeV}$, $H_0=h/(4.5\times 10^{41}{\mathrm{GeV}}^{-1})$ and $M_P\sim 10^{19}\mathrm{GeV}$, givesFor $0.5<h<1$, this is of order $30$--$100\mathrm{eV}$ with these approximate constants. More precise constants change the coefficient, not the requested mass scale. This mass is much larger than the present momentum scale $T_{\nu 0}\simeq 1.4\times 10^{-4}\mathrm{eV}$, validating the nonrelativistic energy density approximation. The result assumes that one thermal relic species dominates the matter density; it is conditional on that model rather than a measurement of a neutrino mass.

In the Newtonian dust calculation, use a comoving spherical shell of physical radius $R(t)=a(t)\chi$. Homogeneity gives its enclosed mass $M=4\pi \rho R^3/3$, and the shell theorem removes the gravitational force of exterior shells. Assume pressureless matter, conserved enclosed mass, no cosmological constant, and no shell crossing. Its specific conserved mechanical energy isHomogeneity permits the energy per ${\chi }^2$ to be a constant, written $-k{c}^2/2$. Dividing by $a^2$ givesThis is a Newtonian derivation for dust. For radiation the enclosed gravitating mass is not conserved, so this same mechanical-energy argument cannot simply be reused; general relativity supplies the same Friedmann constraint with the total energy density divided by $c^2$, together with pressure-dependent fluid conservation.

The numerical fractions describe physically very different components. Use $c=1$ so that pressure and energy density have the same units. The stress-energy conservation equation for a homogeneous component, together with its equation of state $p_i=w_i{\rho }_i$, givesNonrelativistic baryons and cold dark matter have $w\simeq 0$ and dilute as $a^{-3}$; radiation has $w=1/3$ and dilutes as $a^{-4}$; vacuum energy has $w=-1$ and remains constant. The curvature term scales as $a^{-2}$ in the Friedmann equation but is a geometric contribution, not an additional material fluid.

There is consequently a coincidence in the present mixture. Normalize $a_0=1$. ThenAt early times vacuum is negligible, while at late times it dominates overwhelmingly. Observing order-unity matter and vacuum fractions during their comparatively recent crossover motivates the cosmic coincidence problem: why should this crossover occur near our observing epoch?

The absolute scale gives a second puzzle. Restoring SI units, $H_0\simeq 2.27\times 10^{-18}{\mathrm{s}}^{-1}$ yieldsIts energy density is ${\rho }_{\Lambda ,0}{c}^2\simeq 5.8\times 10^{-10}\mathrm{J}{\mathrm{m}}^{-3}$. Quantum vacuum contributions associated with ordinary particle-physics scales can be far larger than this observed value. Explaining the tiny net value, including cancellations among contributions, is the cosmological constant problem. The required negative pressure also differs qualitatively from ordinary dilute matter.

The matter budget supplies another explanatory gap: only $0.045$ is baryonic, while $0.255$ must be nonbaryonic matter in this benchmark. Together with $0.7$ in vacuum energy, most of the cosmic density is in components whose microscopic explanation is not supplied by fitting their fractions. Finally, near-zero curvature raises the flatness problem. For nonzero $k$,During radiation domination $H\propto a^{-2}$, so $|{\Omega }_k|\propto a^2$; during matter domination $H\propto a^{-3/2}$, so $|{\Omega }_k|\propto a$. Near-flatness today therefore requires extremely small early curvature in an unmodified decelerating history. Cosmic inflation offers a mechanism that reduces this curvature fraction through accelerated expansion.

The apparent contrivances are the recent matter-vacuum coincidence, the tiny vacuum-energy scale, the near-flat initial geometry and the unexplained microscopic dark components. Their observationally inferred values impose concrete requirements on a deeper physical account.

Let ${\rho }_0$ be a constant reference density and define total buoyancy acceleration by $\sigma =-g(\rho -{\rho }_0)/{\rho }_0$. Absorb the reference hydrostatic pressure into a kinematic pressure $p$. The ideal nonrotating Boussinesq equations areFor a resting density profile ${\rho }_b(z)$, the background satisfies $p_b^{\prime }(z)={\sigma }_b(z)$ and the squared buoyancy frequency isStable stratification has $N^2>0$. The Boussinesq approximation requires fractional density variations to be small over the domain and parcel excursions, so density can be replaced by ${\rho }_0$ in inertia while its small variation is retained in the gravitational force. For an approximately uniform density gradient over height $H$, this requires $N^{2}H/g\ll 1$. The stipulated incompressibility alone does not justify neglecting density variation in inertia.

Put $b=\sigma -{\sigma }_b$ and $p^{\prime }=p-p_b$. To linear order the two-dimensional internal gravity wave equations areThe exact impermeability condition on $z=\eta (x,t)$ is the kinematic boundary conditionTaylor expansion about $z=0$ permits $w(x,0,t)={\eta }_t$ at leading order if the boundary slope and its displacement relative to the vertical wave scale are small: $|kϵ|\ll 1$ and $|mϵ|\ll 1$. The second condition controls both the omitted evaluation correction $\eta w_z$ and the product $u{\eta }_x$, since continuity gives $u\sim (m/k)w$. Hence for the harmonic forcing,The propagating calculation assumes $k\ne 0$ and $0<|\omega |<N$; a time-independent boundary does no oscillatory work on fluid initially at rest. If $k=0$ but $\omega \ne 0$, incompressibility forces the vertical velocity to be independent of height, so the unbounded column has no finite-energy radiating internal-wave solution. The propagating wave-generation problem therefore needs the stated nonzero horizontal wavenumber.