Peculiar velocity fields from analytic solutions of General Relativity

Peculiar velocities are analyzed through cosmological perturbations in the Newtonian longitudinal gauge characterized by irrotational shear-free congruences in an Eulerian frame. We show that non-trivial peculiar velocity fields can be generated through Lorentzian boosts in the non-relativistic limit, where the Eulerian frame is obtained from analytic solutions of Einstein’s equations sourced by an irrotational shear-free fluid with nonzero energy flux. This approach provides a physically viable interpretation of these analytic solutions, which (in general) admit no isometries, thus allowing, in principle, for modeling time and space varying 3-dimensional fields of peculiar velocities that can be contrasted with observational data on our local cosmography. As a ``proof of concept’’ we examine the peculiar velocities of varying dark matter and dark energy perfect fluids with respect to the CMB frame using a simple, spherically symmetric particular solution. The resulting peculiar velocities are qualitatively compatible with observational data on the CMB dipole.

💡 Research Summary

The paper presents a novel relativistic framework for modelling cosmic peculiar velocities by exploiting exact solutions of Einstein’s equations sourced by an irrotational, shear‑free fluid with a non‑zero energy‑flux vector. Traditional approaches treat peculiar velocities either as low‑redshift Doppler shifts or as linear perturbations in the Newtonian longitudinal gauge, but these become inadequate on scales approaching the Hubble horizon where relativistic effects matter. The authors therefore turn to the most general analytic solution known since 1993 (Sussman et al.), characterized by a metric ds² = −N²dt² + L²δᵢⱼdxⁱdxʲ and a 4‑velocity uᵃ that is irrotational and shear‑free. The solution contains up to ten free functions of time, encoded in a generating scalar Φ(t, xᵢ) and an arbitrary function J(Φ).

A key insight is the reinterpretation of the off‑diagonal energy‑momentum component qᵃ. In earlier literature qᵃ was identified with heat conduction, an interpretation unsuitable for the long‑range, collisionless gravitational dynamics of the universe. The authors instead consider two perfect fluids moving with distinct 4‑velocities ũᵃ and uᵃ, related by a Lorentz boost ũᵃ = γ(uᵃ + vᵃ), where vᵃ is orthogonal to uᵃ and represents the peculiar velocity field. Transforming the perfect‑fluid energy‑momentum tensor under this boost yields Tᵃ_b = (ρ + p)uᵃu_b + p gᵃ_b + 2(ρ + p)v_(a u_b). In the non‑relativistic limit (v²≪c²) the Lorentz factor γ≈1, and the energy‑flux vector reduces to qᵃ≈(ρ + p)vᵃ. Consequently, the exact solution’s qᵃ can be directly identified with a physically observable peculiar velocity field, and the free functions of the solution become the degrees of freedom that shape a three‑dimensional velocity map.

To demonstrate feasibility, the authors select the simplest subclass with J(Φ)=0, which forces the Weyl tensor to vanish (Petrov type O) and yields N = L, eliminating 4‑acceleration. Imposing spherical symmetry by setting the linear terms in Φ to zero gives Φ = α(t)+β(t)r². The resulting metric (eq. 34) resembles a closed FLRW line element multiplied by a factor L = 1 + b(t)F(r), where F(r)=1−cos r. The Einstein equations then relate the scale factor a(t), the function b(t), and the energy‑flux component Q = −2 ḃ f a⁻² (with f(r)=sin r). Because Q∝ḃ, any time variation of b(t) generates a radial peculiar velocity v_r ∝ ḃ r.

The authors expand the dynamics around a closed FLRW background, showing that for small b(t) the density and pressure acquire corrections that match those obtained in standard linear perturbation theory. They then construct two physically motivated scenarios: (i) a mixture of matter (baryons + CDM) and a varying dark‑energy fluid, each with its own 4‑velocity, and (ii) separate dark‑matter and dark‑energy perfect fluids that are not comoving. By choosing ḃ(t) of order 10⁻³ H₀, the model predicts radial velocities of 300–400 km s⁻¹ at distances of ~100 Mpc, in good qualitative agreement with the observed CMB dipole (≈370 km s⁻¹) and bulk‑flow measurements.

The paper concludes that irrotational shear‑free exact solutions, when reinterpreted as tilted fluid mixtures, provide a powerful non‑perturbative tool for generating realistic, time‑ and space‑dependent peculiar velocity fields. The framework accommodates arbitrary 3‑D velocity configurations through its free functions, opening the door to direct fitting of velocity‑survey data. Future work is outlined: exploring more general choices of J(Φ), extending beyond spherical symmetry, and performing quantitative comparisons with modern peculiar‑velocity catalogs. Appendices address regularity conditions, coordinate transformations, and the accuracy of the series expansion versus the full exact expressions.