Cosmological continuity equation ($\dot{\rho }+3H(\rho +P)=0$)

In natural units, conservation of the stress-energy tensor for a homogeneous perfect fluid in an FLRW metric gives the displayed equation for its energy density and pressure. It is equivalent to $d(\rho a^3)=-Pd(a^3)$. For a constant equation-of-state parameter $w=P/\rho$, separation gives $\rho \propto a^{-3(1+w)}$. If $\rho$ instead denotes mass density, the pressure term is $P/c^2$.

This is local stress-energy conservation; it does not in general define a globally conserved total energy for an expanding universe. The adiabatic work in a comoving volume makes this distinction explicit.

For a homogeneous cosmic fluid, the physical volume bounded by fixed comoving coordinates obeys $V\propto a^3$. The cosmological continuity equation therefore gives $d(\rho c^{2}V)=-PdV$. This is the first law of thermodynamics with conserved entropy and no additional chemical work. The internal energy of that volume need not be constant: a component with negative pressure gains internal energy as the universe expands. Local stress-energy conservation does not imply a global conserved energy for an arbitrary expanding spacetime.

A cosmological fluid with constant equation-of-state parameter $w$ has $P=w\rho$ and $\rho \propto a^{-3(1+w)}$.

Quintessence is dynamical dark energy with negative pressure, conventionally $-1<w<-1/3$.

A homogeneous scalar with $V=m^2{\varphi }^2/2$ obeys $\ddot{\varphi }+3H\dot{\varphi }+m^2\varphi =0$. If it is to provide potential-dominated dark energy today rather than rapidly oscillating matter, $m\lesssim H_0\sim 10^{-33}\mathrm{eV}$ is required. Its present density gives ${\varphi }_0\simeq \sqrt{6{\Omega }_X}{M}_{\mathrm{Pl}}{H}_0/m$. Taking $m\sim H_0$ gives a Planck-scale field; a much smaller mass requires a larger field value, so density alone does not fix a unique mass.

// Target: cosmology.bigb

A flat universe dominated by a canonical homogeneous scalar field with constant equation-of-state parameter $w>-1$ has $\rho ={\rho }_0(a/a_0{)}^{-3(1+w)}$ and ${\dot{\varphi }}^2=(1+w)\rho$. The positive rolling branch gives the displayed relation and $V(\varphi )=(1-w){\rho }_0\exp [-\sqrt{24\pi G(1+w)}(\varphi -{\varphi }_0)]/2$. The opposite rolling branch reverses the exponential slope. For $w=-1$ the field is constant; $w<-1$ cannot be realized by a real scalar with this positive kinetic term. An exponential slope $\lambda \sqrt{8\pi G}$ yields $w={\lambda }^2/3-1$, with accelerating nonphantom expansion when $0<{\lambda }^2<2$.

For $w=-2/3$, zero cosmological constant, and negative curvature, $a(t)=\sqrt{\beta }t+\gamma t^2/4$ with density and curvature parameters $\gamma ,\beta >0$.

For $w=-2/3$, a quintessence-dominated scale factor grows quadratically and $tH$ approaches two.