How the Redshift of Gravitons Explains Dark Matter and Dark Energy ()
1. Introduction
This paper describes a theory of gravitons acting in the universe and it superscedes the theory expressed in [1] in regards to the expansion of the universe. We will also describe the action of gravitons in galaxies [2]. Assuming that gravitons are the agents of interaction in a gravitational field, then our goal is to describe how the gravitational redshift of gravitons shows up as what has traditionally been labeled as dark matter and dark energy. We will show how the redshift of gravitons can explain both of these phenomena using high precision data. Throughout this paper, all reference to sources of mass and test particles are assumed to be of baryonic nature, unless otherwise specified. When we speak of the graviton mass mg, it is a relativistic mass.
We assume gravitons are bosonic particles which travel at constant speed c in vacuum, where c equals the speed of light. Gravitons traveling in a gravitational field between a source mass M and a test mass m, modeled by the equivalence principle as an accelerating system, should experience an average energy loss of
due to motion in that field, over a short time period
, where the acceleration a at a point r in the field is given by
. The energy loss is expressed differentially as
(1)
where
is the average relativistic graviton mass, n is the number of gravitons, m is the test mass,
is the change in the free fall velocity of the system observed from an inertial reference frame, G is Newton’s gravitational constant, M is the baryonic mass of the field source, r is the distance between the center of the source and the location of the moving gravitons,
is the short travel time of the gravitons at speed c over distance
. The energy change is a loss (negative), because the velocity change
is in the same direction as the motion of the gravitons, so that for an inertial observer moving in the same direction as the graviton, the energy of the graviton is redshifted. We call this effect of energy loss a gravitational redshift which, as we have described it resembles a Doppler effect. Since gravitons are agents of the gravitational field, our model is essentially an attempt to describe the graviton-graviton interaction as in the self interaction of a quantum gravity theory [3].
Gravitons exist and travel in a gravitational field, which in principle is equivalent to an accelerating system. Assume that the total graviton energy for a system of two masses is expressed by,
(2)
where
is the total graviton mass associated with the test mass m, where mg is the average graviton mass and n is the number of gravitons. The total graviton energy decrease
due to its freefall in the gravitational field of mass M, when viewed from an inertial system, is expressed by
(3)
where, as in (1),
is the velocity increase in the accelerated reference frame equivalent, according to the principle of equivalence, to the gravitational field of mass M at the position r.
2. Gravitons in an Expanding Universe
Consider the universe as a sphere of interior baryonic mass M with a thin spherical shell of mass m. The masses M and m are constants. The thin shell has a radius
at time t. Only the mass interior to the shell has an effect on the shell. The total graviton energy
within the shell at time t is given by (2), where
. Assuming isotropic uniformity, the mass density
within the shell at time t is
(4)
At the present epoch of time
the baryon mass density is
, which is given by,
(5)
where
is the radius of the universe at the present epoch.
2.1. Energy Loss Due to Gravitational Redshift
Since the potential function
is defined as having zero energy at infinity and having a negative energy at position r, likewise we define the graviton energy loss to be zero at infinity and negative at position r. Applying (1) to the n gravitons in free fall in the expanding universe we have the energy loss,
(6)
where
is a coupling constant to be determined by observation and the extra minus sign accounts for an energy loss because the gravitons travel in the same direction as the free fall velocity change and appear redshifted to an inertial observer also moving in the same direction. Substituting for M from (5) and simplifying yields,
(7)
This component is the energy loss of the so called dark matter.
2.2. Energy Loss Due to Expansion
Gravitons traveling at speed c in the vacuum of the expanding universe undergo cosmological redshift in three dimensions on the way to interaction with the masses. We express this redshift by applying the 3-D velocity differential (
) to the total graviton energy
from (2), and applying (3) in the three spatial dimensions, given in the form,
(8)
where the negative sign is applied because the motion of the gravitons is in the same direction as the freefall in the field. We can convert the 3-D velocity differential to a ratio of 3-D volume differential by the construction,
(9)
where x, y and z are Cartesian co-ordinates,
,
and
are independent times and where
,
and
. Furthermore, we convert the volume differential in Cartesian co-ordianes to radial co-ordiates, in the form
(10)
Now, applying the transformations (9) and (10) to (8), while also moving the volume differential
to the left hand side of the equation we get,
(11)
The left hand side of (11) is a quadruple differential whilst the right hand side is a single differential. Integrating both sides of (11) yields,
(12)
where the expansion time is taken to be
(13)
where
is the present radius of the universe and
is a dimensionless coupling constant and where, in the final line we substituted for M from (5). Rearranging (12) we get the graviton energy loss due to the expansion of the universe, the so called dark energy component,
(14)
We will subsequently show that the graviton energy losses
and
can account for the expansion rate of the universe without dark matter or dark energy which is required by the Lamda Cold Dark Matter cosmological model [4].
2.3. Equation of the Expanding Universe
The total energy of the shell of mass m, where
, having kinetic energy, gravitational potential energy and energy loss due to the cosmological redshift of gravitons (14), is expressed by
(15)
where the term on the far right is the total energy, where k is a constant (curvature) with dimensionality [length]−2,
is the present radius of the universe and M is the mass of the universe. Multiplying (15) by
and simplifying, we get the expression for the expansion of the shell,
(16)
Note that the shell mass m can be made arbitrarily small compared to the universe mass M.
Define the distance r by
(17)
where the time varying scale factor a is dimensionless with
. Using (17), the velocity v takes the form
(18)
Substituting (18) into (16) gives us,
(19)
Putting (19) in terms of the mass densities we have,
(20)
where
(21)
is the mass density composed of baryons and graviton energy loss mass due to gravitational redshift (an apparent dark matter mass density),
(22)
is the graviton energy loss due to the expanding universe (an apparent dark energy mass density) and
(23)
is the curvature mass density.
Define the Hubble parameter
by
(24)
where, by (17),
. Equation (24) can also be written as
(25)
which is the Hubble law for
where
is the present epoch of cosmic time. Thus,
defined by (24) is the general form of Hubble’s law [5]. Substituting
for
in (20), with some manipulation, we get
(26)
where
is called the critical mass density at time t. Dividing (26) by
yields the parametric equation
(27)
where
,
and
. At the present epoch
, the mass density parameter is
(28)
where
is the baryon mass density parameter, the graviton expansion energy loss mass density parameter is
(29)
and the curvature density parameter is
(30)
3. General Relativity for an Expanding Universe with Graviton Interaction
The Friedmann-Lemaître-Robertson-Walker (FLRW) metric line element [6] [7] [8] [9] in terms of the scale factor
is
(31)
where the scale factor
and the curvature k has units of [length]−2 where
,
or
. The Einstein equations [10] [11] in trace reverse form is
(32)
where
is the Ricci tensor,
is the energy-momentum tensor, T is the contracted energy-momentum tensor,
is the metric tensor defined by (31) and
. Define the energy-momentum tensor
of a perfect fluid,
(33)
where, referring to (21) to (23) for the mass densities, the total mass density
is given by
(34)
where
(35)
(36)
(37)
and where the current baryon mass density is
(38)
where
is the baryon mass density parameter. We have assumed the equation of state for the relativistic particles
, with
. We are neglecting the radiation density (photon and neutrinos) and we will justify this when we fit the model to SNe Ia data.
Solving the Einstein Equations (32) given the mass-energy tensor (33), and simplifying the results, yields the equation for the rate of change of the scale factor,
(39)
Assuming
for no curvature, and the total mass density given by (34) - (38), the Hubble parameter
is
(40)
where
.
4. Fits to Type Ia Supernova Data and Comparison with the Standard Model
The scale factor a has the relation to the cosmological redshift z expressed by
(41)
An increment of comoving distance
defined in terms of the scale factor is
(42)
where we used the definition (24) of the Hubble parameter. In terms of the cosmological redshift
, the comoving distance is
(43)
where we used the fact that
to transform (42), dropping the minus sign.
The flux
from a distant light source at redshift z is defined in terms of the observed luminosity
, where L is the luminosity of the emitting source,
(44)
where
is the proper distance. The luminosity distance, from (44) is
(45)
The luminosity distance, from (43), is
(46)
where
is the co-moving distance and where
(47)
and we changed the negative sign to positive by inverting the limits of integration. The form (47) of the Hubble parameter
is identical to that of the standard model. Therefore, fitting to the SNe Ia data will be identical. The magnitude is defined, in the standard way,
(48)
where
is the source magnitude and
is an offset. Generally, the source magnitude is combined into
.
We applied (48) in a fit to 580 Type Ia supernovae (SNe Ia) magnitude data from the Super Nova Cosmology Project Union 2.1 data set [12]. With Hubble constant
, densities
,
[12] and offset
, the fit to the SNe Ia data set obtained a two parameter
. Figure 1 shows the fit to the data. Regarding neglecting the radiation density in our model, the radiation density parameter of photons, with
is given by [13],
. At the maximum redshift
Figure 1. Supernova cosmology project union 2.1 SNe Ia magnitude vs redshift data points with error bars. The solid line is the fit for the graviton model with
and
with a two parameter
. The dotted line is the fit for the LCDM model with the same parameters and the same
. The graviton model and the LCDM model fits are (obviously) identical.
of
for the SNe Ia data, the relative magnitude of the error
in comparing the radiation and matter densities is
(49)
which justifies ignoring the radiation density in the fitting.
Assuming that the curvature
, so that
, then from (27) we have for the present epoch
(50)
and from (29),
(51)
For example, for
and
which is in the big bang nucleosynthesis (BBN) allowable range [14] we get a value of the coupling constant
. Also, from (28) we get
. The expansion time, from (13) has the value
.
5. Gravitons in Galaxies
Integrating (1) up to radial distance r we obtain the average energy change per graviton
expressed by
(52)
Equation (52) describes the gravitational redshift of the energy of the average graviton as it travels from a lower, more negative potential to a position r of higher, less negative potential and is consistent with energy conservation.
Now, consider the energy for a galaxy of mass M interior of a small mass m in a circular orbit of radius r. The gravitons traversing the distance at lightspeed from the interior mass to the orbiting mass will experience a decrease in energy as described by (52),
, where n is the number of gravitons. Taking the energy loss of the gravitons into account, the total orbital energy of the orbiting mass m is
(53)
where v is the rotational velocity of the orbiting mass,
is a coupling coefficient, a constant for each galaxy, and the total energy
. Using (53) by expanding
using (52), with
, multiplying by 2/m and moving all terms except
to the right hand side, we obtain the expression for the orbital velocity,
(54)
As an approximation, we model the mass distribution of a spiral galaxy by a spherically symmetric distribution
, even though a mass density distribution consisting of a spherically symmetric central bulge surrounded by an axially symmetric thin disk would be more realistic. Then the mass
of the galaxy within the radial distance r from the galaxy center is given by,
(55)
Coupling Coefficient Kg
Under our assumption that the coupling coefficient is constant for each galaxy, (54) can be solved for
at the galaxy edge, where
and
, giving
(56)
where
is the total baryonic mass in the galaxy. By Newton’s gravitational law, at the galaxy edge, the final velocity
is related to the total mass
contained within the radial distance
by,
(57)
where
, where
is the total apparent mass due to the graviton energy loss, the mass of the so called “dark matter”. The total apparent dark matter in a galaxy is given at the galaxy edge by,
(58)
Using the results from analyzing galaxy rotational data we can estimate the dark matter by (58) as
, and substituting that and (57) into (56) we arrive at a formula for estimating
as,
(59)
6. The Einstein Equations for Galaxies with Gravitons
Assuming a spherically symmetric mass distribution for a galaxy we will use the Schwarzschild metric to describe the motions within the galaxy. The Einstein equations for the region outside of a spherical mass distribution is
(60)
where
is the Ricci tensor. The metric element we use is defined in spherical coordinates,
(61)
where
(62)
(63)
and where
is the mass within radius r. From the metric element (61) the metric tensor is
(64)
Since the metric (64) does not depend on time t and angles
and
, it follows from Hamilton’s principle and the Lagrange equations1 that
(65)
implying that the total energy E is
(66)
from which we get
(67)
And similarly for the
component, we have that
(68)
implying that the specific angular momentum is
(69)
from which we get,
(70)
For motion in the plane defined by r and
,
, where
is defined as
(71)
Equation of Motion
We derive the equation of motion using the metric element (61), where
where
is the differential of proper time
. Dividing the metric element by
and substituting from (66)-(71) with the metric element we get, dropping the
term since
,
(72)
Multiplying (72) by
and simplifying yields
(73)
Substituting for
from (63) into (73) and, since
from (62) and (63) so that
, we get after simplification,
(74)
We assume that the mass
within distance r from the galaxy center is given by the sum of the baryonic mass
within r and the mass
equivalent to the graviton gravitational redshift energy loss mass within r. The mass density of the baryons is given
. The total mass within r, which is the sum of baryonic mass and gravitonic energy loss mass, is expressed by
(75)
where the baryonic mass is
(76)
where
is the baryonic mass density and from (54) we express the gravitonic energy loss mass in terms of the baryonic mass density by,
(77)
where
is a constant coupling coefficient which is peculiar to each galaxy. For circular motion the specific angular momentum term in (74) is
(78)
Dropping the constant terms in (74) and substituting with (78), the velocity of a particle in orbit in the galaxy is expressed by
(79)
7. Results for SPARC Galaxies
We use the velocities from the Spitzer Photometry and Accurate Rotation Curves (SPARC) data base [15] [16], derived from near-infrared (NIR) surface photometry at 3.6 μm. We select the spiral galaxies NGC 2403, NGC 2841 and dwarf galaxies DDO 154 and NGC 2915. From the photometric data which has been reduced to the equivalent velocities for the galaxy bulge, disk and gas mass content, we approximate the baryonic mass as due to a spherically symmetric distribution, which is given by the Newtonian relation for the velocity to mass contained within the radius r from the galaxy center, expressed by
(80)
where
,
,
, N the number of radial distances observed, and the absolute values of the velocities are needed because they can sometimes be negative (Ref. [16], p. 5) . The velocities for the disk and gas from Table 2 of [15] are taken with
. In our Table 2,
and
are the
used in (80) to make the fits. Using the SPARC results for the mass at r in terms of the gas, disk and bulge velocities, where the mass internal to r is given by
of (80), the equivalent graviton energy loss mass is
(81)
The predicted velocity (79) is expressed in the form
(82)
where
,
.
We strived to obtain good fits to the velocity
at each radial distance
by minimizing the mean absolute error MAE between
and
while iterating
for the disk and, when available,
for the bulge. We constrained the mass to agree with the baryonic Tully-Fisher relation (BTFR) [17] [18] mass for each galaxy. Table 1 lists for each galaxy the baryonic mass
determined by the velocity profiles used in (80), the BTFR estimated galaxy mass, the mean data rotation velocity error
, the mean absolute fitted error MAE and the coefficient
.
Table 2 lists the average values for disk
and bulge
mass to light ratios which were used in making each galaxy fit. Also listed in the table are minimum and maximum
taken from Table 4 and Table 5 of [19] for
Table 1. Results of fits to SPARC galaxy data using the graviton model (82) with masses from (80) and (81). The
are the mean error of the reported stellar velocities. The MAE errors are the average absolute error for the fits.
† Using final velocity as flat velocity. ‡ Rotation curve does not flatten.
Table 2. Results of fits to SPARC galaxy data using the graviton model (82). The columns for
and
are the averages for the disk and bulge mass to light ratios determined by modeling the rotation velocity curve.
0 in the
column means the bulge velocity is zero in the data. † Min and max
values for high surface mass galaxies from [19]. NA means the result is not available.
comparison. The fitted spiral galaxy disk mass to light ratios are
except for NGC 2841 which fitted with
. A good reference for fits to SPARC data can also be found in [20], especially for spiral galaxy NGC 2841, where from our Table 2 we have an average of
and
, each of which is of the same order of magnitude as the fitted values
and
, respectively, from Figure 1 of [20]. For each plot of the galaxy rotation velocity, the Newtonian velocity curve is also displayed. Plots for the results of the SPARC galaxies can be found in Figure 2 and Figure 3.
SPARC Galaxies DDO 154 and NGC 2915
Two dwarf galaxies that have been difficult to understand in terms of missing mass are DDO 154 and NGC 2915. Figure 3 shows the results for the fits using the graviton model. The upper curve shows the mass density, the middle curve displays the M/L ratios and the lower plot shows the observed and predicted rotation velocities. With the graviton model for DDO 154 we obtained a total baryonic mass of
, which is less than 2× the detected luminous mass [21] of
. Significantly, it is only one-fifth of the estimated dark plus luminous matter of
. Our fitted disk mass to light ratios have an average of
which is 3× the published value of
.
For NGC 2915 we obtained a total baryonic mass of
, where the stellar mass and luminous HI mass [22] [23]
. This compares with the total dynamical mass
which is 7× the total baryonic mass derived by the graviton model.
8. Discussion
Regarding the baryon acoustic oscillations (BAO) in the primordial plasma and the resulting cosmic microwave background radiation (CMB), if the energy loss in the gravitational field due to the gravitational redshift of graviton energy is the cause of the apparent dark matter in the universe, implying that dark matter
particles do not exist, then in density perturbations of the primordial plasma just before recombination, the baryons, photons and gravitons are synchronized and there is not a central region of dark matter particles as theorized in the standard model [24], so that when recombination begins the photons dissociate from the electrons and the resultant hydrogen atoms begin to attract and form dense regions in the perturbation. Just how the graviton redshift energy loss affects the density formations in the primordial plasma and how this affects the determination of the Hubble constant from the BAO and CMB modeling [25] is an area which needs to be studied.
Regarding galaxy dynamics, with the modified Newtonian dynamics (MOND) [26] [27] there may be a way to compare it with the graviton gravitational redshift theory. MOND takes effect in galaxies at a distance where the central acceleration is around 1.2 × 10−10 m∙s−2 by effecting a transition in the gravitational attractive force from a
to a
form. On the other hand, the graviton gravitational redshift theory augments the Newtonian rotational velocity beginning from the galactic center, which allows to distinguish it from MOND. Referring to Figure 2, at a radius of
the spiral galaxy NGC 2841 has an observed velocity of
and we calculate the Newtonian velocity of
using the SPARC data. The radial acceleration at this distance is
which is at the beginning of the MOND regime and, assuming we have the correct baryon mass within that radius, MOND would predict an orbital velocity smoothly connected to the Newtonian value
but we see that it would be smaller than V by 106.5 km∙s−1.
Another interesting alternative approach to modified Newtonian dynamics is given in [28] which, by way of an additive inverse Yukawa-like term to Newtonian gravitation, purports to account for gravitational dynamics from solar systems to galaxies and galaxy clusters and to the large scale universe expansion. There is an analogy between (Ref. [29] Equation (19)) and our (54) which can to be further explored.
9. Conclusions
The graviton self interaction model describes the effect of gravitons in free fall in the gravitational field, losing energy by way of gravitational redshift and cosmological redshift without emitting any radiation. General Relativity was applied for both the universe expansion and for spiral galaxy rotation curves. By assuming a coupling coefficient
for the graviton redshift in free fall, the graviton model can account for the apparent dark matter in the universe being related to the baryon density. Likewise, the dark energy depends on the coupling constant
and the baryon density. Thus, the apparent dark matter and dark energy are replaced by two constants and the hypothesized redshift of graviton energy.
In the case of the rotational characteristics of spiral galaxies, the graviton theory well explains the greater than expected galaxy rotational velocities in the SPARC data with only the baryonic mass derived from the gas, disk and bulge velocity data with fitted
ratios, with a galaxy dependent coupling coefficient
and with the total baryon mass conforming to the BTFR.
It is apparent that the fundamental aspect of graviton redshift points to the need to include this in the General Relativity field equations.
NOTES
1Wikipedia: Schwarzshild geodesics, https://en.wikipedia.org/Schwarzschild_geodesics.