Gravitational Hydrodynamics of Large Scale Structure Formation
Total Page:16
File Type:pdf, Size:1020Kb
epl draft Gravitational hydrodynamics of large scale structure formation Th. M. Nieuwenhuizen1(a), C. H. Gibson2(b), and R. E. Schild3(c) 1 Institute for Theoretical Physics, Valckenierstraat 65, 1018 XE Amsterdam, The Netherlands 2 Mech. and Aerosp. Eng. & Scripps Institution of Oceanogr. Depts., UCSD, La Jolla, CA 92093, USA 3 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA PACS 98.80.Bp { Origin and formation of the universe PACS 95.35.+d { Dark matter PACS 98.20.Jp { Globular clusters in external galaxies Abstract. - The gravitational hydrodynamics of the primordial plasma with neutrino hot dark matter is considered as a challenge to the bottom-up cold dark matter paradigm. Viscosity and turbulence induce a top-down fragmentation scenario before and at decoupling. The first step is the creation of voids in the plasma, which expand to 37 Mpc on the average now. The remaining matter clumps turn into galaxy clusters. At decoupling galaxies and proto-globular star clusters arise; the latter constitute the galactic dark matter halos and consist themselves of earth-mass MBD. Frozen MBDare observed in microlensing and white-dwarf-heated ones in planetary nebulae. The approach explains the Tully-Fisher and Faber-Jackson relations, and cosmic microwave temperature fluctuations of micro-Kelvins. (a) [email protected] the first stars should be 100-1000 M population III su- (b)[email protected] perstars. OGCs do not spin rapidly so they cannot be (c)[email protected] condensations, and their small stars imply gentle flows in- consistent with superstars. Other difficulties are posed by Introduction. { Structure formation in the Universe empty supervoids with size up to 300 Mpc reported from starts in the plasma of protons, electrons, He atoms and radio telescope measurements [1], dwarf galaxies with a lot neutrinos, that exists up to some 400,000 yr after the Big of dark matter [2] and a preferred axis of evil spin direc- Bang, the time of decoupling (dc) of photons from matter tion (AE) that appears at scales extending to 1.5 Gpc, (last scattering (L) or recombination). Then the plasma a tenth of the horizon scale [3]. Nearly every month transforms to a neutral gas of H and 24% by weight 4He, new observations arise that pose further challenges to the with the neutrinos remaining free streaming. As this oc- ΛCDM paradigm: Correlations in galaxy structures [4]; curs at about four thousand degrees Kelvin, a moderate absence of baryon acoustic oscillations in galaxy-galaxy plasma temperature, we shall seek an explanation in terms correlations [5]; galaxies formed already when the universe of plasma physics and gravitational hydrodynamics alone. was 4 { 5 billion years old [6]; dwarf satellites that swarm This embodies a return to the top-down scenario of large our own galaxy just like its stars [7]. scale structure formation. The recent conclusion by one of us that dark matter Currently it is assumed that cold dark matter (CDM) particles must have mass of a few eV and probably are 1.5 also exists and, clustered before decoupling, has set seeds eV neutrinos [8], means that dark matter is hot (HDM), for baryon condensation. The so-called the concordance urging once more for an explanation of structure formation or ΛCDM model involves also a cosmological constant or from baryons alone, without a cold dark matter trigger. dark energy. It describes a hierarchical bottom-up ap- We shall discuss such baryonic clustering due to a vis- proach to structure formation, stars first, then galaxies, cous instability in the plasma, overlooked by the currently clusters, and, finally, voids. popular linear models of structure formation. CDM is But observations of dense clumps of ancient small stars assumed not to exist, while HDM, though initially impor- in old globular clusters (OGCs) in all galaxies contradict tant to maintain the homogeneity of the plasma, has no the ΛCDM predictions that star formation should begin role in the structure formation. Central in our discussion only after about 300 million years of dark ages and that will be the huge plasma viscosity ν ∼ 6 1027m2s−1 arising p-1 Nieuwenhuizen, Gibson & Schild from photons that scatter from free electrons. This makes small scales to large, contrary to standard turbulence the- the plasma increasingly viscous, while it is also expanding ories that include in turbulence also irrotational flows and with space. At some age before the decoupling an insta- motions dominated by other forces. bility creates the first structures, proto-voids and proto- Gravitohydrodynamics (GHD). { In his approach galaxy-clusters. At decoupling the viscosity drops to hot with hydrodynamic and diffusive modelling, Gibson 1996 gas values ∼ 1013m2s−1, which creates further structures [10] derives several gravitational Schwarz length scales for at the Jeans scale and at the new, small viscous scale. structure formation by Kolmogorian dimensional analy- The plan of this Letter is to review modern theory of p sis. With τ = 1= ρG the gravitational free fall time gravitational hydrodynamical structure formation, to eval- g in the Jeans length L = V τ , there occur first the vis- uate the estimates for various fragmentation scales within pJ S g cous length L = τ γν, where ν the kinematic viscos- Friedman cosmology, and to compare with observations. SV g ity and γ the rate of the strain, i. e., the magnitude of 1 Hydrodynamics. { The description of gravitational eij = 2 (@vi=@xj + @vj=@xi). Second, there is the turbu- 3 1=2 structure formation starts with Jeans 1902. He proposes lent length LST = (ετg ) , where " is the rate of energy that scales for gravitational condensation of a uniform dissipation per unit mass, and third, the diffusive length p fluid of mass density ρ must be larger than thep Jeans LSD = Dτg, where D is the diffusion coefficient. Within 1=2 ac acoustic scale LJ = VS=(ρG) , where VS ∼ c= 3 is the the acoustic horizon scale, dH ∼ ct structures can form at −11 3 2 ac sound speed of the plasma and G = 6:67 10 m /kg s is scale L if dH ≥ L ≥ max(LSV ;LST ;LSD). Newton's constant. As the plasma Jeans scale is always We shall evaluate these scales within the flat Friedman larger than the horizon scale of causal connection, this for- metric ds2 = c2dt2 −a2(t)dr2. The Friedman equation for bids gravitational structure formation before decoupling. baryonic and neutrino matter reads The Jeans criterion reflects a linear gravitational instabil- ity from acoustics, but it neglects the fact that self grav- 2 itational instability of a gas is absolute [9]. All density a_ ΩB + Ων Ωγ 2 2 ≡ Ω(a) = ΩΛ + 3 + 4 : (2) variations will grow or decrease unless prevented by the H0 a a a viscous forces, turbulent forces and diffusion effects. p d R 1 p −τ 1+x Conservation of the specific momentum in a fluid is ex- where Ων = −Ων;0 dτ log 0 dx xe with τ = 2 amν c =kBT0 becomes a-dependent above the Compton pressed by the Navier-Stokes equation, p temperature ∼ 17; 000 K. With dt = da=H0a Ω(a) the @~v age is R t dt, while the `angular' distance to an object at = rB + ~v× ~! + F~viscous + F~other; (1) 0 @t R t0 redshift z = 1=a − 1 reads dA(z) = [c=(1 + z)] t dt=a. averaged over system control volumes exceeding the mo- We adopt the Hubble parameter h = H0=[100km=s Mpc] 3=4 1=2 2 mentum collision length scale. B is the Bernoulli group =0.744 favored in Ref. [8], so that mν = 2 (GF ) me = 1 2 3=2 of mechanical energy terms B = p/ρ + 2 v + lw and the 1:4998 eV and Ων = 0:111=h = 0:173. For baryons we 2 1 2 viscous force is F~viscous = νsr ~v + ( νs + νb)r · (r · ~v), take ΩB = 0:02265=h = 0:0409 from WMAP5 [13], while 3 −5 2 −5 with kinematic shear viscosity νs = η/ρ and bulk viscosity for photons Ωγ = 2:47 · 10 =h = 4:46 · 10 . Finally, νb = ζ/ρ, while other fluid forces may arise. The inertial- ΩΛ = 1 − Ων − ΩB − Ωγ = 0:786 assures a flat space. vortex force per unit mass ~v ×~!, with ~! = r×~v, produces Viscous instability in the plasma. { GHD starts turbulence if it dominates the other forces; for example, with acknowledging the importance of the photon viscos- ~ Re ≡ j~v × ~!j=jFviscousj is the Reynolds number. A large ity. Because it strongly increases in time, already before viscosity corresponds to a small Reynolds number, with c decoupling the plasma becomes too viscous to follow the universal critical value Re ∼ 25 − 100. For adiabatic flows expansion of space [10]. Thus a gravitational instability the \lost work" term lw due to frictional losses is negligible occurs, that tears the plasma apart at density minima, so the enthalpy p/ρ decreases or increases to compensate thus creating voids. Cosmic (super)voids surround us at for changes in the kinetic energy per unit mass 1 v2. 2 any distance and the furthest observable ones are located The turbulence problem is put in a new perspective by at the decoupling redshift. Presently, voids have a 20 times Gibson [10]. Universal similarity laws are explained in under-density with respect to the critical density. In be- terms of the inertial vortex forces ~v × ~!. From Eq. (1), tween voids the galaxy clusters are located on \pancakes" turbulence is defined as an eddy-like state of fluid motion that join in superclusters. where the inertial-vortex forces of the flow are larger than 2 The shear viscosity ν ≡ νs reads for kBT mec [14] any other forces that tend to damp the eddies out [11,12].