Solution of the dark matter riddle within standard model physics: From black holes, galaxies and clusters to cosmology

It is postulated that the energy density of the (quantum) vacuum acts firstly as dark energy and secondly as a part of dark matter. Assisted by electric fields arising from a small charge mismatch in the cosmic plasma, it can condense on mass concentrations. No longer participating in the cosmic expansion, this constitutes ``electro-aether-energy'' (EAE), ``electro-zero-point-energy'' or ``electro-vacuum-energy'', which solves the dark matter riddle without new physics. A radial electric field of 1 kV/m is predicted in the Galaxy. For proper electric fields, EAE can cover the results deduced with MOND. An instability allows a speedy filling of dark matter cores. Hydrostatic equilibrium in galaxy clusters is obeyed. Flowing in of EAE explains why black holes become supermassive, do not have mass gaps and overcome the final parsec problem. Rupture of charged clouds reduces, e. g., the primordial baryon cloud to the cosmic web. The large coherence scale of the electric field acts as a scaffold for gentle galaxy formation and their vast polar structures. In galaxy merging and bars, there occurs no dynamical friction. At cosmological scales, EAE acts as pressureless dark matter. Its amount increases in time, which likely solves the Hubble tension by its late time physics. A big crunch can occur. Of the large cosmological constant injected at the Big Bang, a small part kept that form, without fine-tuning.


TABLE OF CONTENTS
1. Introduction 1.1 Short history of dark matter 1.2 Present status 1.3 A new approach to the problem of dark matter 1.4 Nomenclature for the proposed form of dark matter 1.5 Setup 2 THE STANDARD MODEL AND ITS ZERO POINT ENERGY 2.1 The essence: the Casimir effect 2.2 Zero point/aether energy in the black hole interior 2.3 Zero point energy as a physical entity 2.4 Zero point energy as the dark energy 2.5 Electric field as partner in zero point energy condensation 3 THEORETICAL FRAMEWORK 3.1 Full Einstein equations 3.2 Coulomb electrostatics in electrogravity 3.3 Linearized Einstein equations 3.4 Electro-zero-point energy (electro-aether energy) as the dark matter 3.5 A toy galaxy 3.6 The matter components 3.7 The X-ray gas 3.8 Hydrostatic equilibrium 4 PHYSICAL ESTIMATES 4.1 Charge mismatch in the plasma 4.2 Galactic electric field 4.3 Seeding of galactic magnetic fields 4.4 Implementation of hydrostatic stability 4.5 Metastability vs instability 5 ELECTRO-AETHER-ENERGY IN BLACK HOLES 5.1 Black hole metrics with a macroscopic core 5.2 Super-Eddington accretion by aether energy inflow 5.3 The final parsec problem under EAE accretion 6 ELECTRO-AETHER-ENERGY IN GALAXIES 6.1 Relation to MOND 6.2 Stability analysis 6.3An instability leading to explosive core growth 6.4 Evidence for constant-density, non-cusped cores 6.5 Dissolution of galactic cores 6.6 The electric field scaffold and the vast polar structure

Short history of dark matter
The matter in the world and skies we experience is called "normal matter" by specialists.It consists of protons and neutrons, that are bound by the strong force in atomic nuclei, and electrons that encircle the nucleus due to the Coulomb attraction between the positive protons and the negative electrons.
Nowadays it is understood, however, that normal matter makes up only some 5% of the total mass budget in the Universe.In fact, the stars that we observe in the night sky make up only a modest 4 of the total (Rich, 2009); most of the 5% lies in hydrogen clouds and hydrogen bridges between galaxies, which can be observed in the 21 cm radio line due to spin flipping in the hydrogen nucleus, as was discovered by Hendrik van de Hulst (Cook, 2001).
The remaining 95% of the total is matter that we do not perceive directly, even though its existence has been established rather firmly.After long suspicion of dark stars and the suggestion of dark matter based on stellar velocities by Jacobus Kapteyn in 1922(Kapteyn, 2013), the existence of dark matter was established by Fritz Zwicky in 1933 (Zwicky, 1933) and much support for it has emerged.An account of the history of dark matter is given in (De Swart et al., 2017).Actually, this matter is not dark but transparent.Its French name "matière obscure" translates as "obscure matter", hidden or unexplained matter, which does more justice to its nature.
Observations by Vera Rubin and Kent Ford in the seventies demonstrated that in the outer part of galaxies, circular orbits have nearly the same rotation speed (Rubin and Ford Jr, 1970), constituting "flat rotation curves".This has led to general acceptance of dark matter's existence.
In 1998, it was established that there also exists dark energy (Riess et al., 1998;Perlmutter et al., 1999), which constitutes some 70% of the mass budget in the Universe, while dark matter makes up some 25%.Of all the matter/energy in the Universe, 95% is unexplained so far; it will be the focus of this work.

Present status
Dark matter must be cold, that is: slowly moving, in order to account for the creation of galaxies.Many searches have been carried out.The historical candidate, MACHOs, massive astrophysical compact halo objects such as dark stars or planets, has been ruled out as the main contributor (Tisserand et al., 2007;Alcock et al., 2000).The next candidate is the WIMP (weakly interacting massive particle) proposed by Jim Peebles and followers (Peebles, 1982;Bond et al., 1982;Blumenthal et al., 1982Blumenthal et al., , 1984)).Being massive and moving slowly, it is termed "cold" and leads to the present paradigm of Lambda cold dark matter (ΛCDM).Intensive searching for possible candidates has not yielded a discovery (Arcadi et al., 2018), and its 10 years window since 2010, during which various new searches should find it if it exists, has all but closed (Bertone, 2010).Nowadays, the focus is shifting to axions and axion-like-particles (Wilczek, 1978;Weinberg, 1978;Sikivie, 1983), to warm DM (Boyarsky et al., 2014), and to dark photons, see, e.g., (Pignol et al., 2015).
Of the many other approaches we mention theories without a new particle, like Modified Newtonian Dynamics (MOND) (Milgrom, 1983) and entropic (Verlinde, 2011) and emergent gravity (Verlinde, 2017), which leads to related predictions (Nieuwenhuizen, 2017;Brouwer et al., 2017).Here the Newton law is modified for weak acceleration.A definite test of general relativity versus extended theories of gravity has been worked out for gravitational wave observation (Corda, 2009).
The dark matter (DM) problem is still an outstanding riddle, but various aspects challenge ΛCDM: absence of dynamical friction in galaxy merging (Kroupa, 2015;Oehm and Kroupa, 2024) and in galactic bars (Roshan et al., 2021); as per MOND, structures in the baryons correspond to structures in the rotation velocity.The Hubble tension expresses the difference between the Hubble constant H 0 ≈ 68 km/s Gpc from the cosmic microwave radiation (Aghanim et al., 2020) versus the 73 km/s Gpc from supernovae in the nearby Universe (Brout et al., 2022).This is not solvable (within ΛCDM) with new physics at low redshift (Keeley and Shafieloo, 2023).Likewise, the Lithium-7 problem expressing a factor ∼ 3 larger nucleosynthesis prediction than its observed Spite plateau (Spite and Spite, 1982;Coc and Vangioni, 2005;Pitrou et al., 2018), has remained open.
Observation by the James Webb telescope poses further challenges for ΛCDM.Black holes became supermassive early on.The black hole ULAS J1342+0928 has a mass of 8 10 8 M ⊙ at redshift z = 7.54 (Bañados et al., 2018).The one of GN-z11, an exceptionally luminous galaxy at z = 10.6, weighs several million solar masses and has accreted at about 5 times the Eddington rate (Maiolino et al., 2024).Early galaxies do not look chaotic but are more mature and heavier than predicted.A set of 12 of quiescent galaxies at redshifts 3 − 4 has masses in the 10 11 M ⊙ range, comparable to massive galaxies in the local Universe, and are quenched for more than a billion years (Nanayakkara et al., 2024).For redshifts z = 0.5 -8, many dwarf galaxies are prolate (van der Wel et al., 2014;Zhang et al., 2019), now termed "bananas" (Pandya et al., 2024).The current record holder is JADES-GS-z13-0, a galaxy with spectroscopically confirmed redshift of 13.2, which we observe as it existed 350 million years after the Big Bang (Curtis-Lake et al., 2023).As to clusters, an overdensity related to a galaxy protocluster is observed already at redshift z = 7.88 (Hashimoto et al., 2023).
These challenges put forward to abandon ΛCDM and start anew.

A new approach to the problem of dark matter
We propose an approach1 towards many such riddles without new physics; it will suffice to take a new look at the vacuum of the quantum field theory for the Standard Model of particle physics.Let us consider an analogy with the atmosphere.A "particle" in the atmosphere can be identified with a tornado.Without "particles", the atmosphere is in its "vacuum" state.But this is not a trivial state: the weather can be gentle or involve winds, rains, storms, . . .We shall take this notion over to the (quantum) vacuum, a state without matter, but not without energy.Because of Einstein's relation E = mc 2 , vacuum energy condensed on galaxies has similar (but not the same) gravitational effects as the purported particle dark matter.
Rather than invoking a new particle, we postulate that the zero point energy density of the quantum aether has specific properties: next to being uniform and acting as dark energy, it can be inhomogeneous and flow; it can condense on mass concentrations when assisted by electric fields; locally, it can be positive or negative.
Alternatively, we may view this vacuum as a classical aether (though not the historic fixed aether) with the said properties of energy.Either way, the cosmological "constant" can slowly vary in space, and in time, while its local value can have either sign.The gradient of the related negative or positive pressure acts as a genuine force density that can counteract, e.g., the electrostatic force density.In this picture, dark energy partly acts, in combination with electrostatics, as the dark matter.

Nomenclature for the proposed form of dark matter
For our theory of dark matter as a combination of electrostatic energy and zero-point or aether energy, (or "ether of space", an old term from Sir Oliver Lodge), we have initially employed the term "electro-zeropoint energy" (EZPE) (Nieuwenhuizen, 2023b).The drawback is that this refers to a quantum concept, while such a connection is, for now, merely an assumption.A more neutral term is "electro-vacuum energy" (EVE) or "electro-aether energy" EAE.Another one is (the energy density related to) a local cosmological constant (LCC), a term we employed in our earlier work on black holes with a smooth (singularity-free) interior (Nieuwenhuizen, 2023a).
We are thus hampered by traditions for the use the term "vacuum" as energetically either empty or nonempty.Wishing to distinguish the standard vacuum of quantum field theory, we propose to call that "the vacuum", and use "aether" for the (substance that sustains a) nontrivial vacuum, one with nonzero energy density and/or flow.Henceforward, we will use the acronym EAE.

Setup
The setup of this paper is as follows.Aspects of zero point energy are discussed in section 2. The theoretical framework is presented in section 3. The sizes of various effects are estimated in section 4. Applications to black holes are discussed in section 5.The working in galaxies and comparison with MOND is treated in section 6. Section 7 analyses the application to galaxy clusters.An implementation for cosmology in various epochs is worked out in section 8.The work closes with a conclusion, a summary and an outlook in sections 9, 10 and 11.

THE STANDARD MODEL AND ITS ZERO POINT ENERGY
The standard model of particle physics (SMPP) is a quantum field theory for the U(1) × SU(2)× SU(3) gauge group with three families of quarks and leptons.Since its conception in the 1960's and 1970's, all its particles have been established, the latest ones being the top quark in 1995 and the Higgs boson in 2012.Like in the decades before, SMPP has been capable to explain all experiments so far, the latest success being to rule out the breaking of lepton universality which was suggested in earlier experiments (Collaboration et al., 2022).
In quantum mechanics there is the notion of zero point energy (ZPE).A harmonic oscillator has a ZPE of 1 2 ω, where its oscillation frequency reads ω = k/m with k the spring constant and m the mass.With n quanta, it has energy levels (n + 1 2 ) ω.While the energy differences between the levels have clear experimental meaning, the meaning of ZPE is less obvious.Adiabatically changing either k or m, changes the ZPE; the difference has a physical meaning.
A quantum field decomposes in a large set of harmonic oscillators.Their total ZPE is formally a sum of ± 1 2 ω k (+ for bosons, − for fermions) over the 3d momenta k, a quartically divergent expression.If the momenta are cut off at the Planck scale, there remains a result which is about 123 orders of magnitude2 larger than the dark energy density in the present Universe.This enormous mismatch is a ground for unease to connection with ZPE.
The question we put forward is however: Isn't there a more prominent role for the ZPE? Haven't we, by focussing on the particles and taking the ZPE for granted, been picking out the raisins while overlooking the pudding?Is there room within the standard model to address its deficiencies like the description of the baryon asymmetry, dark matter and the dark energy related to the Universe's accelerated expansion?

The essence: the Casimir effect
In 1948 Hendrik Casimir discovered that two parallel conducting plates of area A at distance d have an energy −π 2 cA/720d 3 (Casimir, 1948;Milton, 2003;Balian and Duplantier, 2004).It is generally understood that this energy is gained from the quantum aether when bringing the plates from infinity to distance d.Interpretations of the effect differ, however, since energy can not be localized in electrodynamics and neither in gravitation.
Upon adiabatically moving the plates to a distance d ′ < d (or d ′ > d), the aether energy changes, so one can say that more aether energy flows out (in).When taking them back to infinity, this energy has to be resupplied as work to readjust the boundary conditions at the surface of the plates: the Casimir effect works as an ideal battery (Nieuwenhuizen, 2023a).
While the Casimir energy for parallel conducting plates is negative, it is positive for a conducting spherical shell, which has a tendency to expand for small and large radii (Boyer, 1968).While the Casimir energy is always positive, a bump in the energy as function of the radius exposes an intermediate regime with tendency to contract (Cetto and De La Peña, 1993).

Zero point/aether energy in the black hole interior
The assumption that aether energy can behave as a fluid is the corner stone for the present work.It arose recently as sine qua non element in our singularity-free solution for the black hole interior (Nieuwenhuizen, 2021(Nieuwenhuizen, , 2023a)).We propose that there is an extended core with a net (positive) charge.In the formation process, binding energy released by dissolving the protons and neutrons into free quarks is partly employed for electrostatic energy.Another part of the binding energy is proposed to be inserted in the quantum aether (without creating particles), which acts as a local cosmological constant (LCC).
The Einstein equations impose that the LCC and the electric field coexist.A change in the charge distribution changes the electric field and imposes a change in the LCC

Zero point energy as a physical entity
In the black hole problem, the zero point energy was treated as physical.Let us motivate this.As said, the naive (bare) expression for the energy density of the quantum vacuum is divergent.That a constant can be subtracted from the bare value is compatible with the Callan-Symanzik equation, the renormalization equation in quantum field theory (Peskin and Schroeder, 2018).This allows to define the renormalized aether energy density as the physical energy density.It can be positive, negative or zero.

Solution of the dark matter riddle within standard model physics
As for the single harmonic oscillator mentioned above, the aether modes can be deformed which relates to a finite aether energy density, positive or negative.For cosmological applications, no boundaries like Casimir plates or shells are present, but an electric field due to a mismatch of + and − charges may take their role.
In the application to the Casimir effect, to BHs and to the cosmos in the present work, the physical zero point energy does not vanish.In each of these cases, this is due to a physical effect.For the cosmos, we shall assume that ZPE was injected in the aether during the Big Bang, to turn, in the course of time, partly into particles, fields and dark matter.

Zero point energy as the dark energy
As a first step to give zero point energy a more prominent role, we identify the ZPE density of the quantum fields with the sought cosmic dark energy.In doing so, we follow our teacher Tini Veltman (Veltman, 1975) at the University of Utrecht, who cites Andrei Linde (Linde, 1974) and Joseph Dreitlein (Dreitlein, 1974), and our colleague Sander Bais with coauthor Robert Russell (Bais and Russell, 1975).The idea itself goes back to Yakov Zeldovich (Zel'dovich, 1967;Zel'dovich, 1968).The ZPE density is non-zero, having the "renormalized" (physical) value of 70% of the critical cosmic mass density ρ c ≃ 9 10 −30 gr/cm3 , not the "bare" (unphysical) one which is 123 orders of magnitude too large -but see section 8.8 for a fresh view.
In making this step, we explain -by default -the dark energy.Its merit is that the case is based on known physics and that we can postulate that the sought dark matter originates as a specific part of the very same ZPE.In section 8.8 we argue that the present small value of the cosmological constant is a dynamical effect, not due to fine-tuning.

Electric field as partner in zero point energy condensation
In the plasma of galaxies and clusters, the first actors are the free charges.The electrons, protons and ions occur at high density ∼ 0.01/cm 3 , and in principle compensate each other.Small local mismatches produce an electric field with its negative longitudinal pressure, which has to be compensated by the ZPE in the way prescribed by the Einstein equations.Hereto, we view the ZPE as a dynamical quantity: energy stored in the aether can flow and partly condense on mass concentrations such as black holes, galaxies and clusters.The free charges are not an accidental property; they provide a skeleton to which the ZPE is attached.
Together, the ZPE and the electric field form a scaffold, a large correlated structure, for normal matter.It constitutes our proposal for the dark matter, which gets supported by various arguments.

Full Einstein equations
We start in general relativity and express the invariant line element ds 2 = g µν dr µ dr ν in spherical spatial coordinates, r µ = (t, r, θ, φ) with µ = 0, 1, 2, 3, as 3 (Nieuwenhuizen, 2021) with functions N(r) and S(r) to be specified later.As generalization of the Schwarzschild metric, it is the most general spherically symmetric one, see e. g.Weinberg's book (Weinberg, 1972).It leads to a diagonal Einstein tensor G µ ν with elements where G 2 2 = G 3 3 due to the spherical symmetry.We express the stress energy tensor T µ ν in parameters ρ λ and ρ E to be specified later, and parameters ρ m , p m r and p m θ = p m φ = p m ⊥ connected to normal matter (m), possibly thermal, with, in principle, an anisotropic material pressure, with all coefficients functions of r.They are fixed by the Einstein equations G µ ν = 8πGT µ ν and may be presented in a mixed way as where (3.7) In case the matter has an isotropic pressure, these reduce to When material energy and pressure can be neglected, one simply has ρλ = ρ λ , ρE = ρ E and ρm = 0.It implies that N(r) = 1, which simplifies the analysis.
Combining the above expressions, or by considering the T 0 0 component, eq.(3.4) can be replaced by This can be integrated to solve S(r) with vanishing value of S(0) as (3.10) (The singular behavior ρ tot (r) ∼ 1/r 2 for r → 0 leads to a finite S(0) and a still regular metric.)Using (3.6), the function N can be eliminated from (3.5), to yield an equation for S alone, (3.11) The strategy to solve these equations is as follows.The energy densities and pressures of matter and electrostatics are assumed to be known, so that S can be solved from (3.11).It will determine M tot and ρ tot via (3.10), next ρ λ via (3.9); lastly, N can be solved from (3.6).After these steps, ρλ can just be read off from (3.4); being determined by the electric charge density and normal matter, it is "enslaved".
From (3.10) we can consider the rotation speed for circular orbits, hence any result for S can be expressed in terms of v rot or M tot .
The task simplifies when the energy and pressure (but not the charges) of normal matter can be neglected, so that N(r) = 1.For the application to singularity-free, cored black holes, this yields a class of exact solutions (Nieuwenhuizen, 2021).

Coulomb electrostatics in electrogravity
We consider a static potential A µ = δ 0 µ A 0 (r) and define the radial electric field as E(r) = −A ′ 0 (r)/N(r).It equals Here Q(r) is the total charge enclosed within radius r, given the charge density ρ q .Since the metric (3.1) is diagonal, the stress energy tensor keeps its special relativistic form, as employed in (3.3).T E is traceless, and involves a positive energy density and transversal pressures, but a negative longitudinal presssure; all equal in magnitude.

Linearized Einstein equations
For application to galaxies and clusters we can linearize around the vacuum N = 1 and S = 0. Indeed, the observed rotation speeds v rot are much smaller than the speed of light c = 1, so that S(r) = 2v 2 rot (r) ≪ 1. Eq. (3.6) yields In linearized form, Eq. (3.11) reads The homogenous solutions are S = r 2 and 1/r.The proper inhomogeneous solution, obeys M tot (0) = 0 and S(r) ∼ r 2 at small r.The involved functions are where Λ, the standard cosmological constant, enters as an integration constant.Surprisingly, the last term in (3.17) involves the outer region u ≥ r. (This would be avoided by subtracting its integral from 0 to ∞, which connects to a homogenous solution S ∼ r 2 , but, acting at large r as a huge negative cosmological constant, that is not acceptable.) Equating (3.17) to the integral of ρ tot in (3.9), viz.
is assumed to have flown in from infinity in the dynamics towards this phase.Notice that even without electric field, the aether already supplies the pressure term M p m , often called a "relativistic correction" to the matter term M ρ m .
The derivatives M tot and M λ yield The total pressure is always anisotropic.After eliminating ρ λ , the pressures (3.21) read The cosmological constant ρ Λ has been included for completeness; for applications to galaxies and clusters it is negligible.Outside a region with matter and a net charge Q, ρ E = ρ > E = Q 2 /8πr 4 leads to ρ λ = ρ Λ + p m .For galaxies and clusters, the net included charge will vanish at some radius, which can be identified with their outer border.

T.M. Nieuwenhuizen
Solution of the dark matter riddle within standard model physics

Electro-zero-point energy (electro-aether energy) as the dark matter
We have now arrived at the point to introduce EAE as the combination of the electric field energy and the ZPE in the LCC.In this paradigm we can identify the DM parts, For small r, the expression for ρ dm typically leads to a constant value, that is to say, EAE naturally involves constant-density dark matter cores.For empirical support of this, see sec.6.4.Nevertheless, powerlaw behavior S ∼ r 2n at small r with n ≥ 0 is possible, with The extreme case is n = 0, where S( 0) is a constant between 0 and 1.
The build up of a net positive charge inside a central region of a galaxy or galaxy cluster, implies that some electrons must move outwards.Let us introduce a crossover radius by R co ; beyond it, the expelled electrons make up a net negative charge density.The included net charge Q(r) grows up to R co but decays beyond it, making ρ E = Q 2 /8πr 4 decay quicker than 1/r 4 at moderately large r.The radial size of the galaxy (cluster) R g can be defined as the radius where Q(r) → 0; due to the large strength of electrostatics, this is a sharp boundary.
Approximating v 2 (r) here by GM dm (r)/r (dark matter only) and denoting v 2 (r) = v 2 (t) and Q 2 (r) = r 2 q 2 (t)/4G with t = log(r/kpc), it follows that The crossover radius where Q ′ (R co ) = 0, now set by the condition 2q 2 + q2 = 0, is related to the rotation curve.In integral form Eq. (3.25) reads These relations, extended to include normal matter, offer a test between the electric charge (and field) profile and the rotation curve.

A toy galaxy
Let us consider a toy galaxy, upon neglecting ρ Λ , ρ m and p m .It has normalized radius R g = 1 and normalized net charge density, which is positive for 0 < r < R co = 3 5 and negative up to R g .The included electric charge and energy density are The local cosmological constant ρ λ starts with positive value but has a negative tail, as shown here for a toy galaxy with charge distribution (3.27).While ρ dm is strictly positive and p r dm strictly negative, p ⊥ dm starts out negatively and has a positive tail.
with Q(R co ) = 432π/3125 = 0.4343.The total charge vanishes, Q(1) = 0. Eq.(3.23) yields While ρ dm is strictly positive and p r dm strictly negative, ρ λ is positive up to r 0 = 0.4849 and then negative beyond.Notice that r 0 < R co = 0.6.Conversely, p ⊥ dm is negative up to r = 0.3449 and positive beyond.The profiles (3.29) vanish at R g = 1.They are plotted in fig. 1.That ρ λ has generally a negative tail follows from eq. (3.20) near R g , where ρ E vanishes.
The accumulated zero point energy, has total M λ ≡ M λ (1) = 4π 2 /945 = 0.04177.Its positive density part between 0 and r 0 contains M λ (r 0 ) = 1.652M λ , while the tail between r 0 and 1 has value −0.652M λ .With respect to M λ , the ZPE imported from remote surroundings, this exhibits an additional 65.2% taken out of the vacuum in the outer region r 0 < r < 1, and transferred to the inner region r < r 0 .With vanishing net total charge, Q(r) = 0 for r ≥ R g , the total electrostatic energy M E = 4π 1 0 dr r 2 ρ E (r) equals 3M λ .These properties reflect our main assumption: zero point (vacuum) energy can be taken out or put into the (quantum) aether, at amounts governed by the Einstein equations.

The matter components
In a galaxy there is normal matter such as stars, planets, free floating planets, hydrogen gas clouds, and a plasma of protons, ions and electrons.
In a galaxy cluster, of all the galaxies, only the brightest cluster galaxy (bcg), located in the center and often the brightest X-ray source in the cluster, brings a relevant contribution to the enclosed mass M tot (r); it can be described by its mass density ρ bcg (r).This galaxy will also contain dark matter, which is a part of ρ dm .Our fit profile ρ bcg in section 7.1 combines normal and dark matter; as of now, there are no data for the separate parts.For simplicity of notation, we shall nevertheless write (3.31)where ρ g and p g describe the X-ray gas, see eq. (3.32) below.

The X-ray gas
Free electrons, accompanied by protons and ions, occur in galaxies and clusters.The mass density ρ g = n e mN involves the local free electron density n e and the average molecular weight mN ≈ (7/6)m N with m N the nuclear weight (Morandi et al., 2010).The pressure is p g = n e T , where the typical scale of T (r) is several keV.So the X-ray gas involves ρ g (r) = n e (r) mN , p g (r) = n e (r)T (r). (3.32) With T in the keV regime, the gas is hot but the protons and nuclei are nonrelativistic, with p g /ρ g ∼ T /m N 0.01.
In the EAE approach the X-ray plasma is indispensable, since by setting up a small charge mismatch (estimated in sec.4.1), it creates the electric field that combines with ZPE as "dark matter".Turned around, in this paradigm the presence of DM requires the presence of a net charge density and hence a (partial or full) plasma.

Hydrostatic equilibrium
For a fluid element in equilibrium, the balance of forces is expressed by hydrostatic equilibrium.Expressed by the energy conservation T µ ν;µ = 0, this is automatically satisfied for a bona fide solution of the Einstein equations, since G µ ν;µ = 0 by construction.For the stress energy tensor (3.3) in the metric (3.1) this leads to the exact force balance where we employed the cosmological pressure p λ = −ρ λ .In exact form, the respective force densities read, using Here F E = ρ q E is the Coulomb force on a unit volume with charge density ρ q in an electric field E. Eq (3.20) allows to replace the combination ρ In EAE, the leading terms are F λ and F E ; neglecting F G , one recognizes the derivative of ρ λ in eq.(3.20) as the solution for hydrostatic equilibrium.The exact approach adds the nonlinear, Newtonian F G term, which is small, see sec.4.4.Integrating (3.33) from r to ∞ yields the full nonlinear correction to ρ λ ; it becomes a self-consistent relation, since ρ λ enters ρ tot through eq.(3.20), and hence through M tot .
The fact that F G contributes to ρ ′ λ , exhibits the malleability (enslavement) of the ZPE/AE, doing just the right thing in the situation at hand.

Charge mismatch in the plasma
For the rotation speed v 2 = GM/r, eq.(3.17) estimates the electric field as E ∼ v/cℓ P r.In the galaxy clusters A1689 and A1835, the free electron density in the center is of order n e (0) = 0.05 − 0.09/cm 3 (Nieuwenhuizen et al., 2021).The density n p of plus charges (protons, ions) is very close to this.Equating E to e(n p − n e )r yields the relative mismatch δ q ≡ (n p − n e )/(n p + n e ) ∼ v/cen e ℓ P r 2 ∼ 6 10 −15 for v = 1000 km/s and typical scale r = 100 kpc.
In our Galaxy near the Sun n e ∼ 0.02/cm 3 leads for r = 10 kpc and v = 200 km/s to δ q ∼ 10 −13 .In 2013, the Voyager 1 spacecraft observed electron plasma oscillations corresponding to an electron density n e = 0.08/cm 3 , very close to the value expected in the interstellar medium and confirming the order of magnitude (Gurnett et al., 2013).The measured electric field of about 1 µV/m at kHz frequencies does not refer to a static field.
These tiny mismatches express that the plasma is not completely neutral, but slightly charged.The electric field is an average of spatially and temporally fluctuating E fields with their related B fields caused by the relatively high density of moving protons, ions and electrons.Since they are so light, it is easier to push out electrons than protons and ions; hence we assume that the inner galaxy is positively charged, and the outskirts negatively, separated by a crossover radius R co .Due to the strength of the Coulomb force, a galaxy has has a sharp boundary R b , where the included charges cancel each other.
In reverse, the locally available amount of ZPE sets the size of the static E field.Such a correlation can be searched for via the determination of the dark matter, also in the intragalactic medium.

Galactic electric field
The possibility of a galactic electric field has been considered in the literature.Ref. (Reucroft, 2014) estimates from rotation curves a field of 1 V/m at the solar location.To keep most electrons in the stellar interior, ref. (Chakraborty et al., 2014) estimates a center-to-surface potential difference of 1000 V; giant galaxies to have a similar potential difference, and rich clusters ∼10 kV.
In sec.7.2.2 and figure 4 we predict similar electric fields of ∼ 1.5 kV/m and 0.5 kV/m in the central ∼ 500 kpc of the galaxy clusters A1689 and A1835, respectively.

Seeding of galactic magnetic fields
Large scale, quasi static magnetic fields abound in the cosmos and play an important role in the evolution of galaxies, but their origin is still debated (Beck et al., 2013).In EAE theory, the formation of a dark matter core is a dynamical effect with a slowly varying electric field, which, according to Maxwell law ε 0 Ė = ∇ × H, generates cosmic magnetic fields.In the local Galaxy and in the fat clusters A1689 and A1835 we predict electric fields in the 1 kV/m range, which corresponds in strength to a magnetic field of 0.029979 G, while observed magnetic fields lie typically in the 1-10 µG regime.
As to orders of magnitude, Maxwell's law yields the estimate B ∼ LE/ct.For the Galaxy, E ∼ 1 kV/m, L ∼ 100 kpc and t ∼ 10 Gyr imply B ∼ 1µG, and for clusters with L ∼ 1 Mpc it yields B ∼ 10µG.As these rough estimates produce the right order of magnitude, this connection deserves a detailed analysis.

Implementation of hydrostatic stability
The net-charge ratio δ q ∼ 10 −13 − 10 −14 needed in the EAE approach is actually relatively large and unexpected.A standard argument is to consider a sphere with mass M = Nm N and charge Q = δ q Ne.In principle, the ratio of Coulomb and Newton forces at the surface, Q 2 /GM 2 = (Qm P /M) 2 , can not exceed unity, which occurs for δ q = m N /2em P = 0.5 10 −18 , quite smaller than the above estimates.Including dark matter by setting M ≈ 2Nm N only alleviates this to 10 −18 .While a huge charged cosmic object seems unrealistic, in the standard Einstein-Maxwell theory there cannot be a stable shell with negative charge density around a positively charged core, as it would add to rather than compensate the Newton attraction.Nevertheless, the typical value of only δ q ∼ 10 −18 looks problematic for EAE theory, so let us analyze the situation.
In the EAE setting, the galaxy or cluster has a large core with net positive charge density, surrounded by a halo of negative net charge density.As a whole, it is not charged.One may define the boundary of a galaxy or cluster R b as the radius where Q(R b ) = 0.
In the interior r < R b , our starting point is eq.(3.33), the hydrostatic equilibrium condition for the balance of forces acting on a fluid element.In terms of the local cosmological pressure p λ = −ρ λ , this takes to leading order the form In absence of ρ E and p ′ λ as in ΛCDM, the pressure must decay in such a way that F G , the Newton attraction per unit volume, is balanced.But notice that, according to (3.31), ρ m = ρ bcg + ρ g includes both the gas and the normal matter mass densities.As we discuss in section 7.2, this approach does not work out properly in clusters.
In EAE, the finite F E = Eρ q is the local Coulomb force density, involving a positive E = Q(r)/r 2 and a ρ q that is positive for r < R co and negative for r > R co .
We can estimate F G /F E as (ρ m /ρ E )(GM tot /r) (GM tot /r) ∼ 10 −7 for a galaxy with M tot = 10 11 M ⊙ and r = 100 kpc, and ∼ 5 10 −5 for a fat cluster with M tot = 10 15 M ⊙ and r = 1 Mpc.These values play the role of the above factor (m N /δ q em P ) 2 , showing that δ q is relatively large, and F G relatively small, in the presence of EAE.Clearly, in eq. ( 4.2) the combination p ′ m + p ′ λ should be balanced, not p ′ m alone.Because p ′ m ∼ F G can even be neglected, the balance is in essence p ′ λ = F E , obeyed grace to eq. (3.20).
In conclusion, the above Newtonian argument fails, since it overlooks the stress (pressure) related to a positive (negative) aether energy density.On a fluid element of EAE, the strong Coulomb force (repulsive for r < R co and attractive for r > R co ) is balanced by an equally strong (inward c.q. outward) force from the zero-point/aether pressure gradient.As a result, galaxies and clusters with charge ratios δ q well exceeding 10 −18 are allowed in EAE.

Metastability vs instability
Hydrostatic equilibrium for patches of matter describes stability at the macroscopic scale, but it does not automatically imply equilibrium at the microscopic scale.Indeed, an individual ion is affected by the strong outward Coulomb force but not know to involve a restoring force from the cosmological pressure, and neither an effect of the outer shell of negative charges.Likewise, an electron is strongly attracted inwards, the ones moving from r > R co to r < R co will lessen the charge mismatch.While for a reshuffling of the charge distribution, an accompanying reshuffling of the ZPE is required, an inherently unstable situation may remain.This observation suggests that EAE cores are unstable, which in principle demonstrates the failure of the theory.But the involved timescales need not coincide, and may be cosmological, so that metastable cores on Gyr timescales are compatible with EAE.
We estimate the effects by connecting to observations.While galaxies and clusters clearly exist with dark matter supposedly arising from EAE theory, there are also indications for their subsequent dissolution.For galaxies these are discussed in sec.6.5.

Black hole metrics with a macroscopic core
The present work was inspired by our solutions for black holes (BHs) with a smooth core enclosed by the inner horizon and an empty mantle up to the event horizon.Assuming that in the stellar collapse a bit more electrons than protons were ejected, the core-collapsed BH will be positively charged.The protons may dissolve into up and down quarks, thereby releasing their binding energy, 98% of their mass.Upon neglecting the quark and electron masses, analytical solutions were presented based on what we now call EAE theory.The stress energy tensor takes the form (5.1) Given the charge distribution and hence ρ E , a class of exact solutions of the Einstein equation (3.5), viz.2S − r 2 S ′′ = 32πGr 2 ρ E , was presented, after which the LCC (ZPE density) ρ λ follows from (3.4), ρ λ = (2S +4rS ′ +r 2 S ′′ )/32πGr 2 , under the condition that it matches the trivial aether ρ λ = 0 at the inner horizon.Special cases were worked out, and it was found that always M λ = 1 4 M, M E = 3M λ = 3 4 M. Next, accounting for quark and electron masses was carried out in a numerical analysis.It was found that the fluctuation spectrum has oscillating modes, but no growing (unstable) ones.
Notice that with ρ tot = −p r tot = ρ λ + ρ E the the first law of thermodynamics is locally satisfied in the form dU(r) ≡ ρ tot (r)dV = −p r tot (r)dV with dV = 4πr 2 dr, confirming that neither a temperature nor a chemical potential is connected to EAE.

Super-Eddington accretion by aether energy inflow
Black holes grow by mass accretion.It falls firstly on the accretion disk, and from there into the back hole or get spit out in the jets.Elsewhere, we show that black holes can also become more massive by inflow of aether energy.Here neither an accretion disk nor angular momentum plays a role.But there is a caveat.In the limit where the energy density (but not the charge density) of the normal matter can be neglected, we have presented a class of exact solutions.They involve a charge distribution with a net charge with charge-to-mass ratio When aether energy is absorbed, without additional charges, this limit will be reached.Further mass increase needs charge increase, which can come from the surroundings.This hints at a common growth of the supermassive black hole and the galaxy.Such has been detected and known as the M − σ relation.
The largest known black hole mass, 66 10 9 M ⊙ for Ton 618, is likely determined by super-Eddington accretion assisted by or even dominated by EAE absorption.
Recent observations support BH growth related to ZPE/AE .Ref. (Farrah et al., 2023b) reports statistics on hundreds of elliptical galaxies in 3 redshift bins in the domain 0 < z 2.5, showing that the supermassive BHs in massive, red-sequence elliptical galaxies must have grown in mass by a factor of 7 from z ∼ 1 to z = 0, and a factor of 20 from z ∼ 2 to z = 0.
GN-z11 is an exceptionally luminous galaxy at redshift z = 10.6.Its spectral features indicate that GN-z11 hosts an accreting black hole.The assumption of local virial relations leads to a black hole mass of log 10 (M BH /M ⊙ ) = 6.2 ± 0.3, which amounts to accretion at about 5 times the Eddington rate (Maiolino et al., 2024).
Both the existence of early black holes, even primordial ones (see later sections), with their super-Eddington accretion, is a central tenet of EAE theory.

The final parsec problem under EAE accretion
For two widely separated black holes to become bound and finally merge, potential energy must be lost.This can occur by dynamical friction, whereby kinetic energy is transferred to nearby matter.For example, a bypassing star can get a slingshot in which it gains kinetic energy and the BH pair becomes more bound.When the pair has a separation of a few parsecs, there is not enough matter to effectively continue this process, while gravitational radiation becomes relevant at distances of 0.001-0.01pc.This is called the final parsec problem (Milosavljević and Merritt, 2003).Various ways out have been proposed, including merging with help of further stars or a third BH.Also disk accretion during the merger of supermassive black hole binaries in galactic nuclei works for them (Armitage and Natarajan, 2002).
In the EAE paradigm, merging happens in a galaxy with a dark-matter core, which gets continuously filled up while being depleted by the BHs.For a BH pair this feeding increases both masses and also diminishes their distance, as we now discuss by a standard analysis.
For a two body problem, like a BH pair, the kinetic and potential energies are (5.2) When m 1,2 depend on time, this reads in terms of the mutual position r = r 1 − r 2 and the barycentre with, as usual, the total and reduced masses (5.4) The center of mass motion involves the conserved momentum so that the centre-of-mass energy P 2 i /2M decreases on mass inccrease.In the frame where P i = 0, the energy reads as for constant masses, which leads to a μ term in the equation of motion For the time-dependent masses, it results in the rate of change of energy with a minus sign in the kinetic term.As expected, only derivatives of the masses, not of the orbit, occur.The angular momentum in the rest frame is conserved due to (5.7).
If m 1,2 change negligibly during one period, circular orbits involve K = −E, V = 2E, leading to Ė/E = 3 μ/µ + 2 Ṁ/M.Integration from t i to t f and E = −GµM/2r yields (5.10) with i denoting initial and f final values, showing tighter binding upon mass accretion, r f ∼ 1/m 3 f .For Kepler orbits, the relations for E f /E i holds due to the virial theorem K = −E, V = 2E.The relation between r f and r i still holds for the epicentre and the apocentre.
The Lagrange-Runge-Lenz ellipticity vector ε and ellipticity parameter ε, with 0 (5.11) From (5.9) and (5.10) it is seen that the ellipticity is essentially constant for slow mass increase.
For m 1 = m 2 = m one has µ = m/2.Merging can be estimated to happen at the Schwartzschild value r f = 2Gm f .This leads to a final mass m f = m i (r i /2Gm i ) 1/4 .In terms of m i = mi M ⊙ and r i = ri kpc, this reads m f = 10 4 ( m3 i ri ) 1/4 M ⊙ , an appreciable mass gain from the aether.Relativistic corrections and gravitational radiation may bring corrections of order unity.

Relation to MOND
Observations by Vera Rubin and Kent Ford in the seventies demonstrated that in the outer part of galaxies, circular orbits have nearly the same rotation speed (Rubin and Ford Jr, 1970), constituting "flat rotation curves".To explain this, Modified Newtonian Dynamics (MOND) was introduced by Milgrom in 1983 (Milgrom, 1983).It assumes that the Newton force decays as 1/r at large r, which can be expressed as an effective enclosed mass behaving as M(r) ∼ r.For recent reviews, see (Lelli et al., 2017;Banik and Zhao, 2022).
For a proper amount of EAE, our approach can explain the MOND results.Let the ρ E profile take the form of an isothermal sphere, E 2 (r) = v 2 /Gr 2 .This leads to ρ λ ≈ ρ E and GM/r ≈ v 2 , so that v(r) = v exhibits the flat rotation curve.The acceleration g consists of the Newton term g N = GM b (r)/r 2 of the baryons (stars and hydrogen clouds) enclosed within r, and the DM term; they combine essentially as (Banik and Zhao, 2022) (The true MOND relation g = g N f (g N /a 0 ), with f (x) = √ x at small x and f = 1 at large x, employed in applications is more rich than this form (Milgrom, 1983;Lelli et al., 2017;Banik and Zhao, 2022).)It was supposed that a 0 ∼ 1.2 10 −10 m/s 2 is universal; empirical values fluctuate around this (Rodrigues et al., 2018), but confirm the baryonic Tully-Fisher relation In our EAE approach, the amount of DM in a galaxy depends on its history, and there is no universal a 0 .The next subsection discusses that constant-density cores are favored in EAE, and next we discuss various empirical evidence for that.
Renzo's rule states that "for any feature in the luminosity profile" (due to stars, X ray gas or hydrogen clouds) "there is a corresponding feature" (a bump, a dip) "in the rotation curve and vice versa" (Sancisi, 2004).The EAE approach offers a mechanism for this: extra structures in the electric field of extra mass structures.However, data for the electric field is needed to give EAE predictive power.Presently, this is not available.
While MOND suffers from its conflict with wide binaries (Banik et al., 2024), EAE is just Einsteinian/Newtonian in this regime.

Stability analysis
The stability of the theory needs to be considered.Following ref. (Nieuwenhuizen, 2023a), we consider a perturbation of the electrostatic potential and of the accumulated charge, δA 0 (r, t) = εa 0 (r)e −iωt , a ′ 0 (r) = −j(r)E(r), δQ(r, t) = εj(r)Q(r)e −iωt , (6.2) with 0 < ε ≪ 1 and yet unspecified profiles j(r), while we again neglect the matter density, but not the charge density, keeping N(r) = 1.This induces radial fluctuations of the metric, δg µµ (t, r) = 2εg µµ (r)h µ (r)e −iωt with h 3 = h 2 for spherical symmetry.The Coulomb energy density picks up metric fluctuations, 3) The remaining Einstein equations correspond to perturbations in the coefficients of (3.3), (6.4) where U µ = δ µ 0 / √ g 00 with U ν = g 0ν / √ g 00 is the time-like unit vector.A first effect is the appearance of elements G 0 1 ∼ G 1 0 ∼ ω, which must vanish since there is no energy current.(While the radial electric field is enhanced by a radial charge current, viz.Ė = −J q , this does not generate a magnetic field.)This imposes a relation best employed only near the end of the analysis.The remaining Einstein equations correspond to the first order perturbations in the coefficients of (3.3), Since the baryons make up a minor part of the energy, we first omit their mass and pressure, neglecting σ m and δσ m , but keeping their net charge distribution.This sets N(r) = 1.One can eliminate h 0 (r) in favor of a nucleus h(r), after which the h i functions follow as (6.8)This finally leads to an equation for h sourced by j, (6.9) The E ′ term enters through S ′′′ 1 from G µ ν , while the other ones come from T µ em ν .
In our nonrelativistic application we have S, S ′ ≈ 0, S ≈ −1, so the leading terms are (6.10) and Eq. ( 6.11) has a well behaved inhomogeneous solution, related to charge relocation coded by δQ ∼ j(r)Q(r) from eq. ( 6.2). h behaves at large r as h(r To determine it in general, let us denote the two real valued homogeneous solutions of (6.11) as J E (r) and Y E (r).Their Wronskian reads where E * > 0 follows by fixing the amplitudes and signs of J E and Y E .There is an inhomogenous solution of (6.11) that is well behaved at large r, reading (6.13)

An instability leading to explosive core growth
In our picture, dark matter stems from ZPE that flows in from infinity, and thereby creating a small mismatch between the charged particles so as to set up an E field.Does the inflowing DE fill up profiles ρ E ∼ 1/r 2 and hence M tot ∼ r and S(r) ∼ r 0 at large r, and if so, how does this happen?
To investigate this, we assume S(r) = 2v 2 n r 2n for small and moderate r, for which the rotation speed equals v(r) = v n r n .While n r 2n−2 /8πG must be non-negative, allowing − 1 2 ≤ n ≤ 1, with ρ E → 0 for n → 1.For small r, our interest is 0 ≤ n < 1 where S ≪ 1 vanishes for r → 0 (it remains bounded for n = 0).For relatively large r, the regime of interest is − 1 2 < n ≤ 0. On this background, we consider an electric field perturbation δE(r, t) = εj(r)E(r)e −iωt with ε ≪ 1 and some profile j(r).With the connection ρ E = E 2 /8π = Q 2 (r)/8πr 4 this is a special case of the treatment in subsection 6.2, with E ∼ r n−1 .Eq. (6.11) now reads For real ω, this has modes oscillating in time, involving modified Bessel functions I ν (ωr) and K ν (ωr).
More interesting are the unstable (growing) modes j(r)e ̟t , h(r)e ̟t with ̟ > 0 and ω = i̟, where eq.(6.2) involves δQ(r, t) = εj(r)Q(r)e ̟t .The case j(r) > 0 connects to increase of positive charge in the region r < R co , making the region r > R co more negatively charged; the case j(r) < 0 describes the reverse.Eq. (6.13) reduces to  6.16) for a toy galaxy with only ν = −1 (n = 0, ρ q = cst/r 2 ) dark matter and enclosed-charge growth function j(r) = 1/ √ r 2 + 1.The temporal growth rates are ̟ = 2 k with k from −3 (top) to 3 (bottom).which involves the ordinary Bessel functions J ν and Y ν .Compared to eqs.(6.12) and (6.11), it holds that J E (r) = r −ν J ν (̟r) and Y E (r) = r −ν Y ν (̟r), while W E (u) = 2/πu 1+2ν .The homogeneous solutions J E and Y E behave as r (1−n)/2 cos(̟r−φ) at large r and have to be omitted.Conversely, h i is well behaved and does not oscillate.With δρ λ /ρ λ = δρ E /ρ E , it follows that the increase of energy is exponential in time, With n ≤ 1 and, if j(r) > 0, h ′ i < 0, it is seen that G(r) > 0 describes mass accretion triggered by the increasing enclosed charge.The profile G(r) is plotted in Fig. 2 for the case j(r) = 1/ √ r 2 + 1 and n = 0, so that ρ q = cst/r 2 .It shows that the "explosive" filling of the DM profile is most effective in the center.
It came as a surprise to us that MOND, corresponding to n = 0 at large r, is not a limiting case, the allowed interval being − 1 2 ≤ n ≤ 1.Actually, the n = 0 case at small r, corresponding to a cusp ρ q ∼ 1/r 2 , E(r) ∼ 1/r, ρ E ∼ ρ λ ∼ 1/r 2 , is an extremal case (Nieuwenhuizen, 2024).
There is evidence for the n = 1, constant-density case, see next section.
When j(r) < 0 in the galaxy or cluster center, the above argument exhibits a case of EAE flowing out of the core.This may be connected with expanding and diluting cores of galaxies, see sec.6.5, and the cooling of cluster cores.

Evidence for constant-density, non-cusped cores
The NFW profile for ΛCDM (Navarro et al., 1997) has an 1/r singularity at the origin, expressing the mutual gravitational attraction of the CDM particles, playing out at their low ("cold") speeds.Whether the CDM consists of axions, WIMPs or MACHOs is not relevant on the scale of galactic cores.The 1/r divergence is called a cusp.However, dark matter cores are often observed to be flat, and the issue is called the "cusp-core problem"; for recent reviews, see (Del Popolo and Le Delliou, 2021;Boldrini, 2021).Dark matter halos seem to have stopped growing: those of nearby quasars are not heavies than those at z ∼ 6 (Arita et al., 2023).
In galaxies, the nearly flat rotation curves occur since normal matter (stars, gas clouds) at small r adds to the DM which dominates at large r.In NGC 3626 the rotation curve keeps growing beyond 8 kpc, for data up to 18 kpc.After accounting for stars and clouds, ref. (Shelest and Lelli, 2020) adds a DM component with v dm ∼ r n for n ≈ 1 or n = 1 in 3 galactic models.The same is done for NGC 2824 and NGC 6176, which do show flattening.The DM with index close to or equal to our limit n max = 1 supports the above analysis.
In the Triangulum galaxy (M33, Cartwheel galaxy) the rotation speed increases as constant+linear beyond 3 kpc up to 15 kpc.It was modelled by determining the contribution from stars and gas, while DM was modeled by an NFW profile (López Fune et al., 2017).As in previous cases, modelling by an n = 1 (constant density) EAE profile may work.
A variety of observations are at odds with the presence of a cusp (Palunas and Williams, 2000;Salucci and Burkert, 2000;De Blok et al., 2001;Karukes and Salucci, 2017;Di Paolo et al., 2019); they favor a constant-density core of a few kpc in size.Ref. (Di Paolo et al., 2019) mentions a mysterious entanglement between the properties of the luminous and the dark matter, which has similarity to Renzo's rule of section 6.1.In EAE theory this occurs, since its ZPE dark matter is regulated by an electric field, due to a net charge distribution that has to adjust itself.Our EAE explanation for the DM in the Galactic center is an indirect support for the binary milisecond pulsar interpretation of the 511 keV Fermi-LAT line (Bartels et al., 2018).

Dissolution of galactic cores
Another piece of evidence is the observational evidence of evolving constant-density dark matter profiles (Sharma et al., 2022).By subtracting the contributions from normal matter, these authors study the dark matter halos of a set of 256 star-forming disk-like galaxies at redshift z ∼ 1.They find constant-density DM cores, as expressed by eq.(3.23) and corresponding to our above case n = 1.But, statistically, the DM cores at z ∼ 1 are denser by 1.5 dex than current ones and smaller by a factor of 0.3 dex compared to present.
Within EAE theory, this can be explained by a diminishing of the net charge ratio, related to the instability of cores mentioned in sec.4.5.It leads to ZPE moving outwards after having moved inwards; the same effect may also have caused the adiabatic expansion of the cluster gas in cool cores in clusters.In principle, this outflow can be described by an approach like the one in sec.6.3, with negative source function j(r), underlining once more the ZPE's fluid character.

The electric field scaffold and the vast polar structure
It has been established that the Milky Way galaxy is surrounded by a vast polar structure of subsystems: satellite galaxies, globular clusters and streams of stars and gas, spreading from Galactocentric distances as small as 10 kpc out to 250 kpc (Metz and Kroupa, 2007;Pawlowski et al., 2012Pawlowski et al., , 2014)).A similar structure occurs in Andromeda (Kroupa, 2014).While an explanation was given as tidal tails of material expelled from interacting galaxies, the predicted EAE scaffold offers a fresh viewpoint.It presents a rather strong structure with which the matter inside it correlates.This picture offers an alternative for the fortunate temporal alignment put forward recently for Gaia data interpreted within ΛCDM (Sawala et al., 2023).(lower) with their fits to the truncated isothermal profile of Eq. ( 7.1) with index n co = 2.

Modified isothermal sphere as a fit for lensing
To apply the idea of isothermal spheres to the galaxy clusters, we consider the clusters A1689 and A1835.In ref. (Nieuwenhuizen et al., 2021), precise strong lensing and gas data and fits to them were presented.Here we consider a different modelling for the DM.A regularization of an isothermal sphere is with r 0 ≪ r 1 ≪ r 2 .At small r it yields E(r) ∼ r and a finite central charge density ρ q (0).In the middle region, r 0 ≪ r ≪ r 2 , it acts as a truncated isothermal sphere.At large r ≫ r 2 , it exhibits incomplete build up (underfill) with index n co .
The data for the cylindrical mass M 2d (r) are expressed in the cylindrical mass density Σ(r) = M 2d (r)/πr 2 , which derives from the 3d mass density as We consider the profile (7.1); the case n co = 2 works well; an analytic expression for its contribution to Σ can be derived.For the brightest cluster galaxy (bcg) with mass M bcg , we add a stretched exponential (se) profile ρ bcg = [M bcg /4πΓ(3n se )n se R 3 se ] × exp[−(r/R se ) 1/nse ], while the X-ray gas has also been modelled in (Nieuwenhuizen et al., 2021).In Fig. 1 we present the data for Σ(r) and fit this to Σbcg + Σdm + Σgas with r 0 → 0. The further parameters are for A1689: (7.3) and for A1835: These values work well, but are not optimized; the error bars will be comparable to the ones in related fits (Nieuwenhuizen et al., 2021).Since the bcg is poorly constrained, other shapes may function as well.
The "underfill" for r r 2 expresses that the surplus electrons pushed outwards are dominant there and diminish the net enclosed charge Q(r) and hence ρ E .
According to (3.20), the LCC ρ λ (r) equals ρ > E (r) − ρ E (r).For ρ E in eq. ( 7.1) with n co = 2, the ρ > E term can be solved analytically, As in the toy galaxy of sec.3.5, ρ λ is positive for small and moderate r, while it has a negative tail −r 4 2 v 2 /24πGr 6 .Its zero crossing lies at 1.2 Mpc for A1689 and at 1.7 Mpc for A1835.
The crossover radius which separates the inner region with positive net charge density ρ q from the outer region with a negative one, occurs for d(r 4 ρ E )/dr = 0, implying R co ≈ r 2 for n co = 2, so that R co = 1.5 and 2.1 Mpc for A1689 and A1835, respectively.

The hydrostatic equilibrium puzzle in clusters
In the Earth's atmosphere, hydrostatic equilibrium is broken by lightning, after which its restoration leads to thunder.In studies of clusters, hydrostatic equilibrium of the X-ray gas is investigated, but found to be dissatisfied (Morandi et al., 2010;Lemze et al., 2011;Morandi et al., 2012); it leads to a ∼ 40% "nonthermal pressure" component in the center of A1689 (Molnar et al., 2010), supposedly due to turbulence, spurious gas dynamics or the dynamical build up of the cluster.. Figure 5 of our ref.(Nieuwenhuizen and Morandi, 2013) shows that for hydrostatic equilibrium the gas temperature (and with it, the pressure) should be larger than observed by a factor ∼ 1.5.This riddle can be solved in EAE theory.The condition for hydrostatic equilibrium, Eq. (3.33), reads

Hydrostatic equilibrium in ΛCDM
In ΛCDM one has ρ E = 0 and p λ (r) = −Λ/8πG, a constant; employing ρ m = ρ bcg + ρ g and p m = p g from (3.31), eq. ( 7.6) results in This exhibits the ΛCDM hydrodynamic equilibrium condition for the gas pressure, but notice that ρ m = ρ bcg + ρ g also involves the bcg, a point generally overlooked.
A consistent approach in ΛCDM simulations should satisfy (7.7), but observed clusters need not, and, as mentioned above, the analyzed ones do not satisfy it.Away from the bcg, where ρ bcg → 0, the relation remains violated.

Hydrostatic equilibrium in EAE
In the EAE approach, the ρ E and p λ terms are present with |p λ | ∼ ρ E , and they are by far the dominant terms, making (7.6) a relation between them with the matter terms as spectators.As we realized in our fitting of the lensing profiles of the clusters A1689 and A1835 (Nieuwenhuizen and Morandi, 2013;Nieuwenhuizen, 2017Nieuwenhuizen, , 2020;;Nieuwenhuizen et al., 2021), the gas is just a spectator, with right nor need for "its own" hydrostatic equilibrium.
In the description of section 3.3 the nonlinear F G terms were left out, which led to the ρ λ in eq.(3.20) as solution for hydrostatic equilibrium.F G indeed acts as a small correction; compared to F E , it is of relative size 10 −6 − 10 −5 .Integrating it from r to ∞ yields a correction to ρ λ , which exhibits the malleability of the ZPE, doing just the right thing in the situation at hand.The respect for hydrostatic equilibrium in EAE implies that no big effects of turbulence or other (gas) dynamics are to be sought for.Rather, it underlines that its violation in ΛCDM is a real deficit.
While ρ λ + ρ E − p m can be identified with the empirical ρ dm profile, data for ρ E are needed to test the hydrostatic equilibrium condition.Unlike the ZPE density, that can only be inferred gravitationally, the electric field acts on charges, and can in principle be determined.Eq. (3.20) yields the prediction E(r) = [−rρ ′ tot (r)/2ε 0 ] 1/2 ; it is plotted in fig. 4 for A1689 and A1835.The central regions involve ∼ 1.5 and ∼ 0.5 kV/m, respectively.

T.M. Nieuwenhuizen
Solution of the dark matter riddle within standard model physics

Zero pressure EAE equation of state
On cosmological scales, one considers the spatial average of the mass density and pressure over a large cosmological volume V with many galaxies.
For a given galaxy (cluster), the electrons compensate the positively charged interior of the galaxy (cluster) beyond the crossover length R co , making ρ E decay faster than 1/r4 .The integral in eq.(3.17) vanishes in the limit r → ∞ , so as to keep the limits Likewise, consider the energy stored in the pressures, Eq. (3.22) yields the values A large cosmological volume V with many galaxies at random positions has a mass density and an isotropic pressure.The mass density can be obtained by smearing out the mass of each galaxy over V and summing.The 3 × 3 pressure matrix is, on the average, diagonal and isotropic.For each galaxy, the contribution to the pressure is obtained as 1/3V times the trace of the pressure matrix integrated over space.Eq. ( 8.3) implies that this trace vanishes for each galaxy and hence for their combination.This leads to the equation of state for the EAE p e = 0, w e ≡ p e ρ e = 0, (8.4) which coincides with the standard p dm = 0 for particle cold dark matter like ΛCDM.
In an alternative approach, we only use that the trace of the EM stress energy tensor vanishes.Averaging over a large volume V with many galaxies yields an isotropic pressure p E = 1 3 ρ E .With p λ = −ρ λ and ρ λ = 1 3 ρ E 4 , eq. ( 8.4) follows again from p e = p λ + p E .

The Hubble tension and the increasing amount of dark matter
The so-called Hubble tension is the fact that local measurements of the Hubble constant via supernovas yield H 0 ≈ 73 km/s Mpc (Brout et al., 2022), while the cosmic microwave background fixes it at ≈ 68 km/s Mpc (Aghanim et al., 2020), with smaller and smaller error bars bringing it to a ∼ 5σ discrepancy.It is often argued that late time physics cannot solve the Hubble tension, see, e.g., the "most general" scenario (Keeley and Shafieloo, 2023).The reason for this is simple, at the present low temperature of the cosmos, no new dark matter particles can be created, so there cannot be an increase of cold dark matter.
But it is also early-time new physics alone cannot solve the Hubble tension (Vagnozzi, 2023).Apparently, the situation is presently in a limbo.
While the EAE acts as a pressureless type of cold dark matter, its mass density can grow in time, since more and more ZPE/AE can be condense locally on BHs, galaxies, clusters, filaments, etc.This will lead to a growth equation for the cosmic dark matter fraction Ω dm ; a modelling is given in eq. ( 8.11) below.It fixes the "enslaved" global dark energy density Ω λ by energy conservation, see eq. (8.15) below.

Friedmannology for aether energy condensation
The Friedman equations for the cosmic scale factor a(t) with a(now) = 1, are Two of these three equations are independent; the last one expresses energy conservation.It is customary to divide out the critical density ρ c = 3H 2 0 /8πG, where the Hubble constant H 0 is the present value of H, and to split up in various components We will consider radiation (r), baryons (b), dark matter (dm), and aether energy (λ).
With w dm ≡ w e = 0, see eq. (8.4), the respective equation of state parameters are Energy conservation gets expressed as with the solution In the present epoch, Ω r (a) = Ω r and Ω b (a) = Ω b are constants, but they change when neutrinos become nonrelativistic due to their small, finite masses; in the past such changes happened during the freeze out of the various species, accompanied by a change in Ω λ .
Deviating from the standard assumption that also Ω dm and Ω λ are constants, we will consider an increasing Ω dm (a) with an appropriate, non-constant Ω λ (a).
The prediction of a growing Ω dm is supported by observational values for the cosmic matter fraction Ω m = Ω dm + Ω b .The early time value 0.315 ± 0.007 from the cosmic microwave background (CMB) observed by the Planck satellite (Aghanim et al., 2020), is smaller than the late time ("now") value 0.334± 0.018 deduced from supernovae in the nearby cosmos (Brout et al., 2022).
To deal properly with the problem, a unified approach covering these epochs is needed.Here we connect to ref. (Dainotti et al., 2021).By parting the Pantheon supernova data in redshift bins, a weak timedependence of Hubble constant H 0 is found and fit to the form H 0 (z) = H now 0 (1 + z) −α with α ≈ 0.01.Inspired by this, we consider the DM growth function (8.11)where δ e ∼ 0.02, deviating from δ e = 0 in ΛCDM.(For simplicity, we neglect a possible a-dependence of δ e ).This small value is compatible with the non-growth of dark matter halos of quasars since z ∼ 6 (Arita et al., 2023).
So, next to radiation (r), baryons (b), we consider EAE (e) dark matter, and a time dependent dark energy (λ).The respective equation of state parameters are Since Ω b and Ω e thus have negligible pressure, the total pressure reads The total energy content is now Like in other applications, the aether, here expressed in the component Ω λ , is enslaved.Eq.(8.9), expressing energy conservation, fixes it as with Ω Λ the present cosmological constant.The a-depencence exhibits that Ω λ (a) was much larger in the past than now, and will be smaller in the future: Throughout the history of the Universe, aether energy is turned into other forms of energy.In this setup, the smallness of Ω λ (1) = Ω Λ is not a result of fine-tuning, but of dynamics, set by the integration constant Ω Λ .
Eqs. (8.14) and (8.15) lead to which can be written in the effective form a 3 + ΩΛ , (8.17) with, as usual, the present-time sum rule (8.18) Eq. ( 8.17) has the familiar form, only modified by an increasing dark matter component.Unlike other rather ad hoc approaches, our ongoing EAE condensation Ansatz (8.11) fits in a bigger picture, the one where EAE also acts as the dark matter in galaxies and clusters, as ingredient of singularity-free black holes, and more.≈ 68 km/s Mpc, apparently solving the Hubble tension with late-time physics, an option unjustly ruled out by restricting in ref. (Keeley and Shafieloo, 2023) to the "most general" scenario.

A fit to CMB data
Of course, a more fundamental analysis is warranted.Employing the CLASS code (Blas et al., 2011;Di Dio et al., 2013), it is possible to find reasonable parameter fits to the Planck CMB TT, EE and TE spectra (Aghanim et al., 2020) up to angular index l = 2000.A specific case5 is depicted in fig. 1.While encouraging by its good fit, this is only indicative.In a proper approach one has to derive the theoretical CMB spectra for the EAE situation of (8.11), considering effects of the electric fields, and fit those predictions to the various data sets, such as CMB, baryon acoustic oscillations and supernovae.At the next level of description, one determines a practical shape for Ω dm (z) from Monte Carlo simulation or otherwise.These steps are beyond the aim of the present paper; we restrict ourselves to stating that the Hubble tension is eased and likely solvable in EAE.

The large Hubble constant and the age of the Universe
Though too seldomly stressed, a large Hubble constant leads to a small age of the Universe, bringing it close to physical lower bounds.The value H 0 = 73 km/s Mpc leads to Age of Universe ≈ 12.8 Gyr,(8.20)rather than the 13.8 Gyr in ΛCDM.This is slightly older than the age of the oldest cluster M4 based on main sequence stars 12.6 ± 1.1 Gyr, or based on the oldest white dwarfs 12.7 ± 0.7 Gyr (Rich, 2009).But it is younger than the estimated age 13.7 ± 0.7 Gyr for the Methusalah star HD 140283 (Creevey et al., 2015), and the accurate 13.535 ± 0.002 Gyr of the ultra-metal poor 2MASS J18082002-5104378 B (Schlaufman et al., 2018).
These and related estimates are based on fitting to simulations of ΛCDM, not to EAE theory.Its predicted early structure formation is suited for early objects.Anyhow, in no consistent theory an age larger than its age of the Universe should appear.

The lopsidedness of the cosmos and the axis of evil
Another challenge to the standard cosmological model concerns the cosmological principle: the expansion of the Universe is homogeneous and isotropic.The largest effect in the cosmic microwave background fluctuations is the so-called the dipole asymmetry.Is it simply due to the motion of the Galaxy through the cosmos, or is it due to a genuine asymmetry in the distribution of matter?Analysis of X ray galaxy clusters (Migkas et al., 2020) and radio galaxies and quasars (Secrest et al., 2022) suggests that the universe is lopsided in our frame.This correlates with the mysterious "Axis of evil" (Land and Magueijo, 2005), the fact that the plane of the Galaxy correlates with the alignment of the low-l (l = 2, 3, 4, 5) multipoles of the CMB: the "top half" of the cosmos is sllghtly cooler than the "bottom half", and the axes of the quadrupole and octopole correlate with it.
It is natural to imagine that the cosmic expansion occurred in an anisotropic way, and, consequently, that the cosmological constant was anisotropic.EAE theory can accommodate that, since it connects the cosmological constant (as vacuum energy) to matter, via the necessary electric fields carried by charge mismatches.

Beyond present
Except for special periods in the early Universe, Ω ′ r and Ω ′ b are zero, so that Ω r,b (a) = Ω r,b keep their present values in the future.(For simplicity, we neglect the fact that neutrinos have a small mass and that BHs radiate).If Ω e (a) were also constant, we would have the ΛCDM connection Ω λ (a) = Ω Λ .Since Ω e (a) = a 3 E dm (a) is the cosmic EAE fraction in a comoving volume, it can continue to increase by further condensation as DM.
Let us consider the far future where expansion leads to a scale factor a > 1 or even ≫ 1 and the integral in the expression for E(a) in (8.10) is negative.The EAE condensation will likely go on until the dissolution of EAE in galaxies and clusters possibly takes the overhand, which could be modeled by a parameter δ e < 0. Finally this leads to a "true" cosmological constant, (8.21) with finite or infinite a max .In the latter case, it results in The proposed shape Ω e (a) = Ω e a δe leads for δ e < 3 to (8.23) The case δ e = 1 connects to Ω e (a)/a 3 = Ω e (1)/a 2 which is commonly connected to curvature of space; here it is a special -and relatively large -parameter δ e .

Black holes and the big crunch
Ref. (Farrah et al., 2023a) considers that supermassive BHs in the redshift interval 0.7 < z < 2.5 have a mass M = M i (a/a i ) k growing with the scale factor as a k since time t i where a(t i ) = a i , where a = 1/(1 + z), and construct a histogram of the distribution p(k), normalized as dk p(k) = 1.Being centred around k = 3, it appears to have weight between k min = 0 and k max ≈ 6.These authors put forward that black holes are the source of dark energy, whereas we assume that dark energy can condense on galaxies and end up in black holes.
When supermassive black holes are relevant for the mass budget, there will be an extra term in (8.14).Assuming for simplicity that p(k) does not depend on a and holds also for a > 1, it takes the form (8.26) Black holes do not create pressure (Hawking radiation is negligible), so for w bh = 0 the energy conservation (8.9) yields an extra term to (8.10), (8.27) vanishing at a = 1, in accord with the definition Ω Λ = Ω λ (1).Eq. ( 8.16) now reads with the closure Ω r + Ω b + Ω e + Ω Λ + Ω bh = 1.For large a the integral behaves as ∼ −a kmax−3 ≪ −1, demonstrating that more energy is extracted from the aether than transferred to the black holes, due to the work cost to get it there.
Apparently, black holes have a tendency to make E negative, but since E = H 2 /H 2 0 , this is physically impossible.If E reaches zero at some a max , this implies ȧ = 0. Since values ä < 0 must then hold already, it describes a Big Crunch with ȧ < 0 from thereon.

Times near the Big Bang
Let us return to eq. (8.9) and add a term Ω φ (a) from an unspecified degree of freedom φ, Assuming that Ω r (a), Ω b (a), Ω e (a) and Ω φ (a) (but not Ω λ ) all start at zero at a = 0, and grow slowly enough, the quantity Ω dS , where dS stands for de Sitter, takes a finite value, Energy conservation (8.9) allows to express Ω λ with Ω λ (1) = Ω Λ as After partial integration, eq.(8.29) reads The decay of E(a) is expressed here by the growth of the integral.

High zero-point energy initial state: Automatic inflation
The product ρ c Ω dS is the aether energy density at the big bang (a ≈ 0).It may have the Planck value ∼ m 4 P , so that Ω dS ∼ 10 123 , which is commonly seen as a catastrophic mismatch between theory and observation.But, as mentioned in section 2.3, that refers to the bare ZPE, which is unphysical; here it refers to the physical energy content at the Big Bang.
Rather than doing away with the large ZPE, we make it a cornerstone.We now assume that ρ c Ω dS = ρ P ∼ m 4 P is the physical zero point energy density injected in the quantum aether during the Big Bang, where Ω dS ∼ 10 123 .This ZPE gets subsequently diluted by the expansion, by creating gravitational waves, by turning it into particles and electrostatic energy, and by participating in the dark matter.
The initial phase of the Universe is a de Sitter universe, with cosmological constant H dS The Friedman equation ȧ2 /a 2 = H 2 0 E(a) leads at early times, when only E ≈ Ω dS ≫ 1 matters, to exponential expansion, a(t) ≈ a P exp H dS t with H dS = √ Ω dS H 0 , the inverse of the de Sitter time τ dS = 1/ √ Ω dS H 0 ≈ 25t P for Ω P = 10 123 and 100t P for Ω P = 10 120 with t P = G/c 5 = 5.39 10 −44 s the Planck time.This behavior is called "inflation"; in EAE theory, it happens automatically; an inflaton field φ seems needed to end the inflation.In this classical approach, the initial time where a = 0 appears to be t i = −∞.
In the course of time, the subtraction terms in (8.31) and (8.32) grow in size and diminish E(a).We assume that, after a certain period, these integrals creep towards Ω dS , leaving a relatively small E(a), and essentially make an end to the period of inflation.Such an end is generally expected to be smooth and called "graceful exit" (gx).
A candidate for enforcing the gx may be primordial black holes.It happens around a = a gx where From then on, the two terms in (8.32) nearly cancel, and the expansion is more effectively expressed by the familiar form (8.14) with Ω φ (a) added, and with Ω λ from (8.31) and (8.30) expressed in only moderately large terms, In order not to "overshoot", i.e., not to make E negative, Ω φ (a) must have become relatively small near a gx .For black holes, this may occur by Hawking evaporation, which leads to particle creation and thus increase of their temperature, the so-called reheating.

Low entropy initial state
Roger Penrose has estimated the final entropy of the Universe by considering it as a huge black hole and applying the Bekenstein-Hawking formula.The result S ∼ 10 123 (Penrose, 1989) coincides with the above Ω dS ≈ m 4 P /ρ c .This is much larger than the entropy ∼ 10 88 in the CMB radiation.Due to the second law, the entropy kept on increasing in the past, so it must even have been much smaller during the Big Bang.The volume of phase space is V = exp S ∼ exp(10 123 ), hence the question arises how the Creator could select our low entropy Universe out of this enormous number of candidates (Penrose, 1989).
In the above EAE initial state, only EZP/VE is present but no particles, neither photons nor gluons nor black holes, so that all field modes lie in their quantum ground state.Though general quantum systems are described by a mixed state, this case can be described by a pure quantum state, roughly as a product of individual ground states, like the Hartle-Hawking state (Hartle and Hawking, 1983), but with modified individual zero point energies to code the ZPE injection.The fine grained entropy, in this situation the von Neumann entropy, vanishes in a pure state.Being conserved under quantum dynamics, it vanishes at all times.This puts many constraints on the ensuing dynamics, some of which are coded in known conservation laws.
Neglecting these correlations leads to a coarse grained entropy, that increases in time.Starting from zero, this solves Penrose's conundrum: in EAE theory a zero -or at least, a low -entropy initial state is a consistency property, for which neither fine-tuning nor selection out of a vast set of candidate universes is involved.

CONCLUSION
At the time of writing, there are two standard models.The first is the standard model of particle physics, formulated as a quantum field theory which is shown to be renormalizable by our teachers Gerard 't Hooft and Tiny Veltman.The second is ΛCDM, the standard model of cosmology, based on the assumptions of a cosmological constant and cold dark matter.Next to the no-show in multiple CDM searches, this approach suffers from the Hubble tension.
Here we put forward that such a new type of matter is neither wanted nor needed, and that standard models of cosmology and particle physics are actually one and the same.No dark matter particle, which would require an extension of the standard model, is involved; a new view on the zero point energy of the vacuum suffices to explain the main constituents of the Universe, the 95% fraction of dark matter and dark energy.Given that the Casimir effect for moving parallel conducting plates involves an inflow or outflow of zero point energy (ZPE), we consider the ZPE as a fluid that can partly act as dark matter.Interpretation of our analytical results leads to consider the energy of the vacuum itself as zero, or unmeasurable in any way, while energy added to it, or taken out from it, acts as "aether" energy (though not the historic aether ruled out by the Michelson-Morley experiment), a physical component, subject to the Einstein equations.
In this interpretation, we are led by the principle that one should first solve the mathematics and then provide a physical interpretation of the results, as applied to our approach to dynamics of quantum measurement and the ensuing statistical interpretation of quantum mechanics (Allahverdyan et al., 2013(Allahverdyan et al., , 2017(Allahverdyan et al., , 2024)).
The Einstein equations require that ZPE is assisted by an electric field, which can arise from a tiny mismatch between plus and minus charges in cosmic plasmas.The combination is termed electro-zeropoint energy (EZPE) or electro-aether-energy (EAE ), which aims to replace the popular cold dark matter.In fact, the connection to an electric field seems natural but is not compulsary; its role may be taken by any vector field producing the structure ρ E C µ ν in eq.(3.3).
Rather than invoking new physics, EAE theory takes a new view on the capacities of the zero point energy of the (quantum) vacuum or just the energy of the classical vacuum.We are led to view it not as a static, uniform entity but as a type of fluid, that can condense on mass concentrations.This application of the standard model appoints an indispensable, and even leading role for the aether, to function as the main actor in cosmic structures by providing, presently, 70% of the total mass/energy as dark energy and 6% involved as the ZPE part of the dark matter, combined with 19% electrostatic energy in the dark matter.Particles, in the form of normal matter, only play a secondary role, coming into existence later in the early Universe and forming presently some 5% of the total mass.EAE theory predicts that the dark matter present during the emission of cosmic background radiation arose from ZPE/AE condensation; this leads naturally to the assumption of primordial black holes.They may have grown by EAE condensation and merging.Black holes from stellar collapse can likewise grow by gentle inflow of EAE, filling "mass gaps" and triggering the growth of supermassive BHs not dominated by merging.Massive BHs may "steal" the EAE from small surrounding ones, the tidal field effect.
Next, the BHs organize the galaxy around them, by an interdependent propensity for the available zero point energy, that partly streams in from infinity and is partly taken out of the vacuum in the outskirts.To achieve this, the charge mismatch has to be optimal according to the Einstein-Coulomb equations.In a galaxy and in a cluster there is a dynamical connection with the baryons: in order to host more ZPE coming in from infinity, an adjustment of the local net charge mismatch has to take place.Flow of aether energy into black holes requires inflow of charges, assuring a dynamical connection between the central BH and the whole galaxy.
A dynamical instability is identified, which assists in a speedy buildup of galactic and cluster cores with constant DM density, supporting EAE theory and observations on the cusp-core problem.The reverse mechanism can explain the expansion and possible dissolution of galactic and cluster cores.
Hydrostatic equilibrium in galaxy clusters satisfied in EAE.
On cosmological scales the EAE acts as a pressureless type of cold dark matter.EAE theory goes even one step further: the "cosmological constant" measured from supernovae is merely the present value of the dynamical zero point energy, that may have started out at the Big Bang with the field theoretic value larger by some 123 orders of magnitude.In EAE theory there is no fine-tuning, the "cosmological constant" is small, since it is the present value of a decaying function.Despite Einstein's most famous equation E = mc 2 , EAE theory involves a discrepancy between mass and energy.Mass and matter are related to particles, including photons, while energy relates to a modified aether without further particles.The kinetic energy of particles remains included in the "mass".
The analysis of the present section shows that the final parsec problem is solved by mass accretion as it happens for EAE.Two supermassive BHs at parsec distance will finally merge by absorbing aether energy, which enhances the probability for observing gravitational waves from merging events by the future LISA system.

SUMMARY
Electro-aether energy (EAE ) relies on electric fields and the zero point energy of the quantum fields of the standard model of particle physics.Alternatively, it is just a property of the "classical" vacuum.The energy, often equated to a cosmological constant, actually gets depleted in its condensation as part of dark matter.These insights explain a cornucopia of phenomena.
After considering various aspects of ZPE in section 2, the EAE framework is laid out in sec.3.For spherically symmetric setups, it is shown how a non-uniform ZPE, combined with an electric field, is compatible with the Einstein equations.ZPE is absorbed from the environment, while subject to a reshuffling inside the galaxy or cluster; its density is positive inside a core region and negative in the halo region.There results a core with a net plus charge, surrounded by a halo with net minus charge; the total charge is zero.
A stability analysis is carried out.An inhomogeneous solution is connected to the formation of dark matter cores made up of EAE, and their later dissolution.
Estimates for various quantities are discussed in section 4. Particular attention is paid to the net charge fraction in the plasma.While standard estimates allow maximally a fraction of 10 −18 , EAE involves a fraction that can perhaps be 10 5 times larger.Analysis of the hydrostatic equilibrium shows that the strong Coulomb repulsion and attraction is counteracted by the negative casu quo positive gradient of the ZPE pressure.
Section 5 deals with black holes.They can grow by EAE accretion, which rules out the "mass gaps" from standard arguments, and is supported by some gravitational wave events.It is shown that the final parsec problem for black hole merging is overcome by ongoing mass accretion within EAE theory.
Section 6 considered the application to galaxies.It is postulated that results from Modified Newtonian Dynamics (MOND) can be modeled by EAE theory, and that the involved electric field and underlying charge density regulates a connection between the dark matter structure, the shape of the rotation curve and the central supermassive black hole.It is pointed out that constant-density dark matter cores, more than cusped ones, should be expected, and support for this is reviewed.
In the application to galaxy clusters of section 7, first an isothermal sphere-type of fit is worked out for strong lensing data of two fat clusters and the relation to the charge distributions and ZPE profiles is worked out.Special attention is payed to their hydrodynamic equilibrium puzzle, solved in EAE theory.
For the application to cosmology, section 8 first shows that at cosmological scales, the pressure connected to EAE dark matter vanishes, as desired.It is pointed out that ongoing ZPE condensation leads to a late-time increasing amount of dark matter.A fit of the ΛCDM theory for the Planck data for the cosmic microwave background already softens the Hubble tension between its value H 0 = 68 km/s Mpc and the late-time value H 0 = 73 km/s Mpc from supernovae.To investigate a full resolution of the problem, the CMB theory within EAE theory needs to be worked out and fitted.Next, the cosmologicallynear and far future is considered.A big crunch scenario is worked out involving a dominant role of black holes.Finally, attention is paid to the Big Bang period, where it is assumed that a large cosmological "constant" is inserted, leading to an initial state with large aether energy and zero entropy.Inflation occurs automatically.
Various further dark matter aspects in galaxies, clusters and cosmology seem to fall into place like pieces of a jigsaw puzzle, see Section 9 of ref. (Nieuwenhuizen, 2023b).

OUTLOOK
In EAE cosmology, there is no dark matter particle, as supported by the no-show in dark matter searches, but the theory gets ruled out when such a detection is made.
Simulations for the EAE paradigm are desired to test it on various observations, replacing the current ΛCDM simulations.Given the great expertise in the latter, the situation seems hopeful.Irrespective of our proposal, the recent James Webb Space Telescope observation of very early onset of massive galaxy formation (Labbé et al., 2022) already seems to demand a new understanding of structure formation.
With the arrival of a new standard model, many issues in cosmology may hope for explanation.We have mentioned the Hubble tension, softened already, the Lithium-7 problem and hinted at black holes for a big crunch and, perhaps, a gentle exit of inflation.
The predicted smaller age of the Universe of some 12.8 Gyr poses questions regarding the earliest stars and structures; these are not new, however, since they follow mainly from adopting the large value of the present Hubble constant.In this regard, early black hole growth by EAE accretion and early galaxy formation due to rupture of hydrogen clouds may emerge as a consistent picture, allowing vast polar structures of matter around them due to the EAE scaffold.Being charged locally, the expanding primordial hydrogen cloud will be subject to lightnings, after which the thunders may provide the onset of prolate and/or barred structure formation.
The question "why is the cosmological constant so small" gets the EAE answer: the cosmological "constant" depends on space and time; there is no fine-tuning, during the Big Bang a large zero point energy was inserted, which decreased.It was a cold Big Bang, that quickly became hot.The possibility for a recollapsing Universe may relate to a cyclic repetition of expansions and collapses.
In all these processes, aether energy (vacuum energy) is the ideal servant, an obedient, malleable agency, doing just the right thing at the right time.One may wonder whether it plays a similar role in standard, terrestrial electrostatic and perhaps magnetostatic problems.Progress in this direction will be reported elsewhere (Nieuwenhuizen, 2024).Lastly, one may wonder whether, as in black holes, also in cosmology, the rotation of structures can carry some of the burden of the net charges.

Figure 3 .
Figure 3. Cylindrical mass density as function of r in the galaxy clusters A1689 (upper) and A1835 (lower) with their fits to the truncated isothermal profile of Eq. (7.1) with index n co = 2.

Figure 4 .
Figure 4. Prediction for the static electric field (in V/m) as function of the radius in the galaxy clusters A1689 (upper, black) and A1835 (lower, blue).

Figure 5 .
Figure 5. Fits to the TT, EE and TE spectra of the cosmic microwave background observed by the Planck satellite for l up to 2000.The large value of the Hubble constant of 73 km/s Mpc is compatible with the data due to the dynamical nature of electro-zero-point energy.