“A Grand Design for Dark Matter” – Trinh Manh Do

Good morning! Today we are going to investigate Dark Energy or generally visible matter, radiation. Even many theoretical models describing dark matter have been persuasively demonstrated, dark matter has not been detected directly, making it become one of the greatest mysteries in our cosmos.

Dark matter is invisible. Based on the effect of gravitational lensing, a ring of dark matter has been inferred in this image of a galaxy cluster CL0024+17 and has been represented in blue color.

Observations have provided hints that the dark matter around one of the central four merging galaxies is not moving with the galaxy itself.

Strong gravitational lensing as observed by the Hubble Space Telescope in Abell 1689 indicates the presence of dark matter-enlarge the image to see the lensing arcs.

The Bullet Cluster: HST image with overlays. The total projected mass distribution reconstructed from strong and weak gravitational lensing is shown in blue, while the X-ray emitting hot has observed with Chandra is shown in red.

Dark Energy is a key explanation for accelerating universe expansion. You can obviously find the experimental and observational evidences of dark energy through many internet search engines. The most important point comes from the observation of the supernova in galaxies which are accelerating away from us. Suggesting that the universe expansion is itself accelerating. We can start our investigation from Einstein who was trying to understand the nature of the universe and he was being told by the cosmologists at the time that the universe was fixed and unchanging. On the other hand, it can be understood that the distances between galaxies in the universe do not change.

Albert Einstein 

Now this statement became an immediate problem for Einstein because he knew that the only long scale force that was operating that he was aware of was the gravitational force. As a result, all the galaxies would exert a gravitational force which was an attractive force on all the other galaxies. And they should be as it was attracting themselves and pulling all the galaxies in. So the universe ought to be contracting.

The gravitational waves are created causing the attracting force

The astronomers at the time was however quite adamant that the cosmos was fixed. Therefore Einstein concluded there must be another force which is called “cosmological constant”, which was equal and opposite to the force of gravity. This balance is essentially maintained the galaxies in space in a fixed position.

Now, before Einstein got very far with this analysis, a man named Edwin Hubble had come along and indicated that the universe was in fact expanding. We can easily understood by comparing this problem to the expansion of a rubber band.


This is a rubber band

Let’s spin the rubber band at the 5 different points A, B, C, D and E respectively. Assume that A is fixed.


With the scale factor, the distance between any two adjacent points is called a. Now, let’s stretch the rubber band, the coordinate points will move and a will obviously increase. It can be also noticed that as the point get further to A, the increasing amount of a increases.


Now we have a formula,

D(DE)=Δx.a(t)                                                                                                                   (1)

D(DE) be the distance between D and E. Δx is the difference in coordinate points (in this case, Δx(DE) would be equal to 4-1=3. a(t) is the function of time because a itself increases as the universe increases. Now let’s differentiate both sides of (1) with respect to time.

dD/dt = D’ = Δx.a'(t)                                                                                                     (2)

Δx remains unchanged, but a does change with time so we will need to take the differential of a with respect to time.

On the left hand side, dD/dt is simply a calculation of velocity V. Rewrite (2)

V =  Δx.a'(t)                                                                                                                         (3)

Now we can always multiply anything by one. Multiply (3) by a/a, which is simply equal to 1.

V =  Δx.a'(t).a/a                                                                                                                (4)

Rearrange (4)

V =  a'(t).Δx.a/a                                                                                                                (5)

It can be noticed that Δx.a=D from (1). Therefore

V =  D.a'(t)/a                                                                                                                      (6)

We can simply put (t) away from a’.

V =  D.a’/a                                                                                                                            (7)

Now a’/a is particularly called by Hubble constant H. Finally

V=HD                                                                                                                                    (8)

And that means if we know the velocity of a receding galaxy, you can calculate how far it is. It also implies that the further something is away from you, the faster it’s travelling. One important point to make the Hubble constant is that it’s constant everywhere in space today, but it does change with time. So it is a constant in space, not in time.

So now we have a new picture of universe where all the galaxies are moving away from each other. There is no sense in which there is a  center of the universe. That’d rather like dots on a balloon. You can take a balloon and mark dots on it. As the balloon expands, the dots moves away from each other.

Let A be the center of reference frame that we happened to choose and B be another galaxy. From here, d is considered to be distance between A and B. In this case, A is a fixed position and B therefore will be moving away from A, with a velocity V which was previously demonstrated to be equal to H.D.

Now the universe is generally regarded as being homogeneous and isotropic which basically means the same everywhere. In addition, the Newton’s theorem tells us that there will be a gravitational force acting on galaxy B. And the mass that will be acting on B, will be all the mass that is contained within a sphere of which B is at the H and A is at the center. Let M and m be the mass of all surrounding galaxies and B, respectively. So the formula to calculate the gravitational force acting on that galaxy would be

A____________________B————> (V)

F=GMm/(d²)                                                                                                                      (9)

Furthermore, the potential energy of that galaxy is equal to

PE = – GMm/d                                                                                                                (10)

Now the total energy of this galaxy will be equal to the sum of kinetic energy and potential energy. And the value of total energy must remain constant in time, which is named k.

Total E = KE + PE = 0.5mv² – GMm/d  = k                                                     (11)

Now, multiply both sides of (11) by 2

mv² – 2GMm/d = 2k                                                                                                   (12)

Divide m both sides of (12)

v² – 2GM/d = 2k/m                                                                                                      (13)

On the right hand side of (13), 2k/m must remain constant because m, which is known as the mass of galaxy B, is also a constant in time. 2k/m can be considered to be equal to another constant k.

v² – 2GM/d = k                                                                                                               (14)

Now, rewrite (1) and (3) below

D(DE)=Δx.a(t)                                                                                                                   (1)

V =  Δx.a'(t)                                                                                                                         (3)

Substitute (1) and (3) back into (14)

[Δx.a'(t)]² – 2GM/[Δx.a(t)] = k                                                                              (15)

Moreover, the value of mass m will be equal to volume V times density p. As the universe expands, the volume will increase, but the mass does not because the number of galaxy remain constant. The density therefore will theoretically decrease and vary with time.

M = V.p = [4πD³/3] . p(t)                                                                                         (16)

Now we substitute (16) back into M of (15)

[Δx.a'(t)]² – 2G . [4πD³/3] . p(t)  /  [Δx.a(t)] = k                                            (17)

Substitute (1) into (17)

[Δx.a'(t)]² – 2G . [4πΔx³.a(t)³/3] . p(t)  /  [Δx.a(t)] = k                               (18)

Simplify (18)

[Δx.a'(t)]² – 2G . [4πΔx².a(t)²/3] . p(t) = k                                                        (19)

In order to be dimensionally consistent, the equation (19) that must follow that K on the right hand side, has got a Δx in it, otherwise dimensionally it’s not going to work. Now we can cross out the term Δx² on the left hand side of the equation. It simply means we’re left with a different constant but it’s still a constant. Let’s rearrange (19)

a'(t)² – 8πGΔ.a(t)².p(t)/3 = k                                                                                  (19)

Finally, let’s multiply both side by a(t)²

[a'(t)²/a(t)²] – 8.π.G.p(t)/3 = k/a(t)²                                                                    (20)

Rearrange (20)

[a'(t)²/a(t)²] =  8.π.G.p(t)/3    –  k/a(t)²                                                               (21)

And (21) is known as a simple version of Friedman Robertson Walker formula. It is an exact solution of Einstein field equation of general relativity. The formula describes a homogeneous, isotropic expanding universe or contracting universe that may be simply connected or multiply connected.

Particularly, it can be noticed that 8.π.G.p(t)/3 is always positive. And if -k  is always positive, the whole term on the right hand side is positive. And the left hand side therefore must be also positive. It means the universe will constantly expand forever and that’s called an open universe. If -k was negative, there might come a point where the negative value of k/a(t)²  is larger the positive value of 8.π.G.p(t)/3 . And therefore [a'(t)²/a(t)²] will become negative. This situation would eventually stop the universe growing and start contracting. That would be called a closed universe. If k was 0, it would indicate that the kinetic energy and potential energy was balanced and this case is called a flat universe. On the other hand, the universe will continue to expand, but it will slow down and become an asymptotic value. It is only after an infinite amount of time.

A diagram is provided as follows

The time and space dimension represents x and y-axis, respectively. If we have an open universe, it will expand forever and never stop. If we have a closed universe, then the universe will expand, but it’s constantly being slowed by the gravitational forces, until the gravity eventually stops the universe expanding anymore and then causes it to collapse on itself. This process is named big crunch. The cosmos crunches back into a single singularity and can be understood as the reverse procedure of the Big Bang. In a flat universe, the universe expands based on the particular circumstances and being asymptotically towards a value, but never reaches it. It will, but only after an infinite amount of time.

Let’s consider a cube of universe. In this case, p, M, V and a will be the density, mass, volume and sides of this cubed universe, respectively.

Rewrite the relationship between mass, volume and density

p = M/V                                                                                                                             (22)

The volume of a cubed universe is considered to be equal to

V = a³                                                                                                                                   (23)

Substitute (23) back into (22)

p = M/a³                                                                                                                             (24)

Substitute (24) back into (21)

[a'(t)²/a(t)²] =  8.π.G.M/[3.a(t)³]    –    k/a(t)²                                                  (25)

If the universe is considered to be flat, it means k will be equal to 0.

[a'(t)²/a(t)²] =  8.π.G.M/[3.a(t)³]                                                                          (26)

Now we will discover how a(t) varies with time. The initial way to do that is to guess. Let’s a(t) be described by a new exponential equation.

a(t) = c.(t^p)                                                                                                                   (27)

Differentiate (27)

a'(t) = pc.t^(p-1)                                                                                                           (28)

Divide both sides of (28) by a(t)

[a'(t)/a(t)] = pc.t^(p-1)/a(t)                                                                                     (29)


[a'(t)/a(t)] = pc.t^(p-1)/[ct^p]                                                                               (30)

Simplify (30)

[a'(t)/a(t)] = pt^(p-1)/[t^p]                                                                                  

The right hand side of above equation is simply equal to

[a'(t)/a(t)] = p/t                                                                                                              (31)

Substitute (31) back into (26)

p²/t²8.π.G.M/[3.a(t)³]                                                                                           (32)

Substitute (27) back into (32)

p²/t² = 8.π.G.M/[3.c³.t^(3p)]                                                                                 (33)

Now, we have a t² term on the left hand side of the equation. And the only t term on the right hand side of the equation is t^(3p). In order to be dimensionally consistent, there must be

t² = t^(3p)                                                                                                                         (34)

Solve (34), the solution would be

p=2/3                                                                                                                                  (35)

Substitute (35) back into (27)

a = ct^(2/3)                                                                                                                     (36)

 Now we plot the value of a against the value of t.

At the beginning of time, in the first 10000 years, the universe would not have been matter dominated, but radiation dominated. When the Big Bang happened, the assumption is that a huge amount of matter was created. But all the matter that was created would have been equal to the amount of anti-matter created. Those matter and anti-matter particles would combine and annihilate and produce a huge amount of radiation in the form of photons. Surprisingly, the process also left a small surplus of matter, which is essentially everything we can see in the present universe. So now let’s consider our cubed universe which contains photons only.

Assume that a is the length of side of the cube. The energy of photon is Planck constant times the frequency.

E = h.v = hc/λ                                                                                                                  (37)

Now as the universe expands, a also expands and λ therefore increases.

a = k.λ                                                                                                                                 (38)

Consequently the energy of the photon must decrease when λ increases. The high energy photon, possibly gamma rays created through the annihilation of matter and anti-matter had their wavelengths stretched as the universe has expanded. And the wavelength is now in the microwave regions and that’s what we see in the sky called the cosmic microwave background radiation. This type of radiation comes from all directions. In general, E is inversely proportional to λ and a.

E = N/λ = N/a                                                                                                                 (39)

Where N is always a constant.

Next, what is the energy density? The value is calculated by dividing energy value by the volume.

E (density) = E/V = N/(a^4)                                                                                   (40)

Let’s p(R) and p(M) be the energy density for radiation and matter dominated universe, respectively. p(M) is the rewritten form of the equation (24).

p(R) = N/(a^4)                                                                                                                (41)

p(m) = M/(a^3)                                                                                                              (42)

Rewrite (26)

[a'(t)²/a(t)²] =  8.π.G.M/[3.a(t)³]                                                                           (26)

Substitute the energy density for radiation dominated universe (41) back into (26). Instead of having the mass density for matter dominated universe which is M/a(t)², we are now going to substitute the radiation density. Rewrite new equation

[a'(t)²/a(t)²] = 8.π.G.N/[3.a(t)^4]

And once again, let’s investigate how a increases with time. With matter dominated universe, a expands based on the power of 2/3, which is derived in (36).

Similarly, we have

a=c.(t^p)                                                                                                                           (37)

Rewrite (31)

[a'(t)/a(t)] = p/t     

Square both sides of (31)

[a'(t)/a(t)]² = p²/t² =  8.π.G.N/[3.a(t)^4]                                                          (38)

Substitute (37) back into (38)

[a'(t)/a(t)]² = p²/t² =  8.π.G.N/[3.(c^4).(t^4p)]                                             (39)

In order to be dimensionally consistent, it can be concluded that

t² = t^(4p)                                                                                                                         (40)

Solve (40)

p=0.5                                                                                                                                   (41)

Substitute (41) back into (37)

a =c.t^(0.5)                                                                                                                      (42)


 I have pointed out that at the beginning of the radiation dominated universe. And that radiation would have an energy which is equal to Planck constant times the frequency. Now if you have high energy, you also have high temperature as a result of being directly proportional to each other.

E=kT                                                                                                                                    (43)

Where E is the energy, k is the Boltzmann constant and T is the temperature. At around 3000K, the atoms are ionized. The electrons in those atoms have so much energy that they all disappear from the atoms and the atoms are simply left essentially as the nucleus. In the very early universe, when the temperature is above 3000K, there would have been no atoms. They would just have been as it were hydrogen nuclei which were protons and free electrons. A huge amount of photons that had been produced through the anti-matter matter annihilation, would not be able to get very far in the universe because they would constantly be interacting with electrons and swapping energy. And consequently, the universe in those days would not be observable because the photons couldn’t travel through it.

The temperature of the universe today is approximately 3K. At some points, when the universe cooled, it cooled enough for the electrons to be caught up by the protons to form hydrogen atoms. And when that happened, the photons would have a much easier direction to escape the current position. During this, the universe became transparent. On the other hand, it become observable.

Rewrite the relationship between E, T, λ and a.

E = k.T                                                                                                                                 (44)

E = hc/λ = m/λ = m/a                                                                                                 (45)

Where m is a constant. Based on (44) and (45), it can be noted that

T = m/a                                                                                                                               (46)

The temperature has fallen from 3000K to 3K by a factor of 1000. From here, it can be concluded that a must increase by a factor of 1000. So the scaling factor today is

[a-today/a-ion] = 1000                                                                                            (47)

In a matter dominated universe, a-today can be equal to c.t^(2/3), a-ion can be equal to c.(t-ion)^(2/3). t in this case is simply the age of the current universe.

[a-today/a-ion] = 1000 = [c.t^(2/3)] / [c.(t-ion)^(2/3)]                         (48)

Simplify (48)

[a-today/a-ion] = 1000 = [t^(2/3)] / [(t-ion)^(2/3)]                                (49)

As a result of that,

t/(t-ion) = 1000^(1.5)                                                                                               (50)

Rearrange the units

t-ion = t/[1000^(1.5)]                                                                                               (51)

Using a graphing calculator provides that 1000^1.5 will be equal to 31622.7766. And t means the current age of the universe, which is approximately 10^10 years.

t-ion = [10^10]/31622 = 316235 (years)                                                     (51)

Therefore 316235 years is the time when the atoms would be able to form after the Big Bang. Rewrite the Friedman Robertson Walker formula (21)

[a’²/a²] = 8.π.G.p/3

p is essentially density of the matter dominated universe which is equal to

p = M/a³

Consequently, the  acceleration has a scale factor of the universe that is proportional to t^(2/3). 

a = c.t^(2/3)

In a radiation dominated universe, the density p is

p = M/(a^4)

And the scaling factor is

a = c.t^(0.5)

Now let’s consider the pressure and thermodynamics. Firstly, we have a cubed matter dominated universe containing all galaxies inside. And the problem is what pressure do those galaxies exert on surface of the cubed universe?

The answer is 0 because those galaxies are bouncing backwards and forwards on the surface on the cubed universe. They are simply suspended in space as the cosmos expands. It will carry the galaxies with them, but there is no pressure on the sides. In general, in a matter dominated universe, the pressure on the sides of a cube is 0.

In a radiation dominated universe, let’s suppose we think about one dimensional issue. A photon is travelling backwards and forwards between these two points A and B, at the speed of light.

                 A     o→                                                                           B


So what is the pressure on this surface AB. We all know that

P = F/A                                                                                                                               (52)

Where P is pressure, F is force and A is area. And the force F in this case is the rate of change of momentum over time.

F = dp / dt                                                                                                                         (53)

Let’s suppose that we have a photon with momentum p, travelling towards B and coming back. So the total change of momentum will be 2p. And the photon takes an amount of time to complete the experiment.

Time = Distance/Speed = D/V = 2L/c                                                            (54) 

F = dp / dt = 2p / [2L/c] = 2.p.c / 2L = p.c / L                                                (55)

Based on research, it can be noticed that the energy of a photon is equal to the momentum times the velocity.

E = p.c                                                                                                                                 (56)

F = p/c / L = E / L                                                                                                           (57)

E / L is essentially the density term in one dimension. But we have to discover it in three dimensions so let’s now take our cube of photons which can travel in any direction.

P = F / dA = E/(L.dA) = E/V = p                                                                             (58)

So we conclude that for the radiation dominated universe, the pressure is apparently equal to density.

The pressure on one side of the cube then would be equal to

P = p/3                                                                                                                                (59)

In general, the total pressure formula is

P = w.p                                                                                                                                (60)

For w=0, P=0, then it’s the matter dominated universe

For w=1/3, P=p/3, then it’s the radiation dominated universe.

Now the work formula is provided as follows:

Work = Force x Distance = F.D = P.A.D = P.V or PdV                          (61)

With Force = Pressure x Area = P.A

With Volume = Area x Distance = D.A

The term “work” must depend on the term “energy” which comes from internal energy of molecules. So when the molecules lose the energy, and that means they lose temperature because energy and temperature relate. Now the change in energy is equal to

dE = – P.dV                                                                                                                       (62)

As volume increases, the energy decreases. We can also conclude that the total energy is equal to the energy density times the volume.

E = p.V      (63)

dE = p.dV + V.dp                                                                                                          (64)

Based on (62) and (64), we have

dE = p.dV + V.dp  =  – P.dV                                                                                     (65)

Rearrange (65)

V.dp = -(P + p).dV                                                                                                       (66)

Substitute (60) back into (66), it can be noticed that

V.dp = -(w.p + p).dV                                                                                                   (67)

Rearrange (67)

V.dp = -(w + 1).p.dV                                                                                                   (68)

 Rearrange (68)

[dp/p] = -(w + 1).p.[dV/V]                                                                                       (69)

 Rearrange (69)

d.log(p) = -(w + 1).d.log(V)                                                                                     (70)

 Rearrange (70)

log(p) = -(w + 1).log(V)                                                                                              (71)

 Rearrange (71)

log(p) = logV^[-(w + 1)]                                                                                            (72)

 Rearrange (72)

p = V^[-(w + 1)]                                                                                                            (73)

 Rearrange (73)

p = (c) / [V^(w + 1)]                                                                                                     (74)

Let’s think about a cube of space again, where a is length side and V is the volume. We also have

p = (c) / a^[3.(w + 1)]                                                                                                  (75)

What does it mean?

For w = 0, it’s the matter dominated universe, that gives us p=c/a³. And that’s exactly what we determined before.

For w = 1/3, it’s the radiation dominated universe, that gives us p=c/(a^4). And that’s exactly what we determined before.

By using thermodynamics, this problem has been successfully demonstrated. It means we discovered the same results. The equation (75) therefore is true.

Substitute (75) back into the Friedman Robertson Walker formula (21)

[a'(t)²/a(t)²] =  8.π.G.p/3  = [8.π.G.c] / 3a^[3.(w + 1)]                             (76)

For w = 0, it’s the matter dominated universe.

For w = 1/3, it’s the radiation dominated universe.

So what happens if w = (-1)?

Substitute  w = (-1)  into (76)

[a'(t)²/a(t)²] = [8.π.G/3] x p(0) = constant 

The result is a constant. On the other hand, it’s the cosmological constant. Now we have

[a’/a] = Hubble constant = H = √[(8.π.G.p(0)/3]                                        (77)

Rearrange (77)

a’ = a .√[(8.π.G.p(0)/3]

We can recognize that the condition for exponential equation where differential of a value is equal to the value itself. For example, dy/dt is proportional to y. In this case,

a = c.e^(u.t)

Where u = √[(8.π.G.p(0)/3]

So that means the scaling factor a exponentially increases with time. It’s accelerating with time and that’s how we think the universe is doing. And what is causing that acceleration?

=====>>>     DARK ENERGY

Trinh Manh Do

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s