A Probabilistic Method to Determine Whether the Speed of Light Is Constant
Author(s) Motohisa Osaka
ABSTRACT
Although the formula of mass-energy equivalence was derived from the hypothesis that the speed of light in free space is constant, conversely, the purpose of this research is to show that a method of probabilistically determining whether the speed of light is constant is derived from this formula. By considering the formula of mass-energy equivalence to be a function of the energy of an object moving at speed V, the probability density function (PDF) of the energy can be obtained using the inverse function of this formula, if the speed of light obeys a probability distribution. The main result is that the PDF of the energy diverges to infinity at a certain energy value regardless of the PDF of the speed of light. Thus, when the speed calculated from this value enters a certain range of the speed of light as V increases stepwise from below 299,792,458 m/s, the PDF of the energy should increase abruptly. If not, then the speed of light is constant. This is the method of probabilistically determining whether the speed of light is constant. An experimental method is proposed to confirm this.

1. Introduction

Albert Einstein published the theory of special relativity in 1905  . The theory is on the relationship between space and time. One of its results is mass-energy equivalence: E = mc2, where E is the energy of an object when it is moving, m is its mass while moving and c is the speed of light. This is derived from two hypotheses. One is that the speed of light in free space is constant for all observers, regardless of their relative motion or of the motion of the light source. This hypothesis is generally considered verified by the Michelson-Morley experiment, which shows the differences between the speed of light in the direction of motion of the earth and that in different directions are within experimental errors  . This was supported by similar experiments with higher resolutions   . However, these experiments do not show that the speed of light in free space is constant in all inertial systems. A team at University of Glasgow reported that photon group velocity was reduced using time-correlated photon pairs and that the delay was several micrometers over a propagation distance of the order of 1 m  . Although they showed that adding spatial structure to an optical beam of single photons reduced the speed of light, the significance of their study was considered to be limited. The findings do not affect the formula of energy-mass equivalence because c in the formula is still regarded as the maximum speed of all objects in free space. However, the question is whether the maximum speed is constant. This study presents a probabilistic method derived from E = mc2 to determine whether the speed of light is constant.

2. Mathematical Steps

It assumes that the speed of light obeys a probability distribution.

The formula E = mc2 is also expressed as

$E=\frac{{m}_{0}}{\sqrt{1-{\left(\frac{V}{c}\right)}^{2}}}{c}^{2},$ (1)

where m0 is the rest mass of the object, and V is its speed; V < c.

The assumptions are:

1) m0 and V are constant.

2) c obeys a probability distribution between ca and cb; ca < cb. Two probability distributions are adopted: a uniform distribution and a triangular distribution.

Then, the probability density function (PDF) of E is calculated by obtaining the inverse function of E as follows.

Step 1: Determining the inverse function of E

As 0 < V/c < 1, V/c is defined by sinθ (0 < θ < π/2).

From (1),

$E=\frac{{m}_{0}{V}^{2}}{\mathrm{cos}\theta -{\mathrm{cos}}^{3}\theta }.$ (2)

Then, cosθ is represented as y:

$c=\frac{V}{\sqrt{1-{y}^{2}}}.$ (3)

Representing m0V2 as P0, the following third-order equation for y is obtained from (2):

$f\left(y\right)={y}^{3}-y+\frac{{P}_{0}}{E}=0.$ (4)

A positive value of E results from any value of θ on the interval 0 < θ < π/2. In other words, at least one of the three roots of f(y) is between 0 and 1. As the roots of f(y) are intersections of $g\left(y\right)={y}^{3}-y=y\left(1-y\right)\left(1+y\right)$ and h(y) = −P0/E, all three roots are real. Two of them are between 0 and 1 and the remaining root is negative. The two roots between 0 and 1 are denoted as y1 and y2; y1 ≤ y2. As 0 < y < 1, the negative root is neglected. Since y1 and y2 are determined by E, it is necessary to find out whether the inverse function of E, c º h(E), is a two-valued or single-valued function.

Step 2: Determining whether E is a two-valued or singled-valued function

From (1), E has the minimum $\frac{3\sqrt{3}}{2}{P}_{0}$ (P0 = m0V2) at ${c}_{m}\equiv \frac{\sqrt{6}}{2}V$ . Three cases are examined on the basis of whether cm is between ca and cb.

Case 1: ${c}_{m}\le {c}_{a}$

Then $V\le \frac{\sqrt{6}}{3}{C}_{a}$ . As E increases monotonically between ca and cb,

$\frac{V}{\sqrt{1-{y}_{1}^{2}}}<{c}_{a}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{c}_{a}\le \frac{V}{\sqrt{1-{y}_{2}^{2}}}\le {c}_{b}.$ (5)

Therefore, only y2 is accepted. In this case c = h(E) is a single-valued function (Figure 1(a)).

Figure 1. Case of ${c}_{m}\le {c}_{a}$ . The object speed and limits of c are set to V = 0.5, ca = 2.2 and cb = 4.2. Then, cm ≈ 0.6124. (a) Inverse function of E for c; (b) Probability distribution of E when c obeys the uniform distribution; (c) Probability distribution of E for three cases: the c-value of the vertex of the triangular probability distribution of c is 2.7 (solid line), 3.2 (dashed line) or 3.7 (dotted line). The probability distribution of E is triangle-like with the vertex moving rightward with the probability distribution of c.

Case 2: ${c}_{a}<{c}_{m}<{c}_{b}$

Then $\frac{\sqrt{6}}{3}{C}_{a} . As E is parabolic and downward convex between ca and cb, the following relationship occurs:

${c}_{a}<\frac{V}{\sqrt{1-{y}_{1}^{2}}}<\frac{V}{\sqrt{1-{y}_{2}^{2}}}<{c}_{b}.$ (6)

Both y1 and y2 are accepted and c = h(E) is a two-valued function (Figure 2(a)). Otherwise, the following occurs:

$\frac{V}{\sqrt{1-{y}_{1}^{2}}}<{c}_{a}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{c}_{b}<\frac{V}{\sqrt{1-{y}_{2}^{2}}}.$ (7)

Then c = h(E) is single-valued function (only y2 is accepted in Figure 2(a)).

Figure 2. Case of ${c}_{a}<{c}_{m}<{c}_{b}$ . The object speed and limits of c are set to V = 2, ca = 2.2 and cb = 4.2. Then, cm ≈ 2.4995. (a) The inverse function of E is a two-valued function between minimum E and the value marked with an arrow and a one-valued function between the marked value and maximum E; (b) Probability distribution of E when c obeys the uniform distribution; (c) Probability distribution of E for two cases: the c-value of the vertex of the triangular probability distribution of c is less than cm, (2.3 < cm) (solid line), and larger than cm, (3.7 > cm), (dashed line). In the former, the PDF of E decreases monotonically. In the latter, it has a vertex. If ${c}_{a}<{c}_{m}<{c}_{b}$ , the probability density of E is maximum at minimum E in either distribution of c.

Case 3: ${c}_{m}\ge {c}_{b}$

Then $V\ge \frac{\sqrt{6}}{3}{C}_{b}$ . As E decreases monotonically between ca and cb,

${c}_{a}\le \frac{V}{\sqrt{1-{y}_{1}^{2}}}\le {c}_{b}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{c}_{b}<\frac{V}{\sqrt{1-{y}_{2}^{2}}}.$ (8)

Therefore, only y1 is accepted. In this case c = h(E) is also a single-valued function (Figure 3(a)).

Step 3: Calculation of PDF of E

From the PDF of c, fc(c), and the inverse function of E, c = h(E), the PDF of the random variable E, fE(E), can be obtained as

${f}_{E}\left(E\right)={f}_{c}\left[h\left(E\right)\right]|\frac{\text{d}}{\text{d}E}h\left(E\right)|$ (9)

Figure 3. Case of ${c}_{m}\ge {c}_{b}$ . The object speed and limits of c are set to V = 2, ca = 2.2 and cb = 2.4. Then, cm ≈ 2.4995. (a) Inverse function of E for c; (b) Probability distribution of E when c obeys the uniform distribution; (c) Probability distribution of E for three cases: the c-value of the vertex of the triangular probability distribution of c is 2.25 (solid line), 2.3 (dashed line) or 2.35 (dotted line). The probability distribution of E changes from trapezoid-like shape to triangle-like shape as the vertex of the probability distribution of c moves rightward.

where $|\frac{\text{d}}{\text{d}E}h\left(E\right)|$ is the Jacobian of the transformation. The absolute value should be taken since the PDF must be positive:

$|\frac{\text{d}c}{\text{d}E}|=|\frac{\text{d}}{\text{d}E}h\left(E\right)|=|\frac{\text{d}c}{\text{d}\theta }||\frac{\text{d}\theta }{\text{d}E}|$ (10)

From c = V/sinθ and y = cosθ,

$|\frac{\text{d}c}{\text{d}\theta }|=\sqrt{\frac{{P}_{0}}{{m}_{0}}}|-\frac{\mathrm{cos}\theta }{1-{\mathrm{cos}}^{2}\theta }|=\sqrt{\frac{{P}_{0}}{{m}_{0}}}\frac{y}{1-{y}^{2}}.$ (11)

From (2) and y = cosθ,

$|\frac{\text{d}\theta }{\text{d}E}|=\frac{{P}_{0}}{{E}^{2}}\frac{1}{\mathrm{sin}\theta |1-3{\mathrm{cos}}^{2}\theta |}=\frac{{P}_{0}}{{E}^{2}}\frac{1}{\sqrt{1-{y}^{2}}|1-3{y}^{2}|}$ . (12)

From (10), (11) and (12),

$|\frac{\text{d}c}{\text{d}E}|=|\frac{\text{d}}{\text{d}E}h\left(E\right)|=\frac{{P}_{0}\sqrt{{P}_{0}}}{\sqrt{{m}_{0}}}\frac{y}{\left(1-{y}^{2}\right)\sqrt{1-{y}^{2}}|1-3{y}^{2}|}\frac{1}{{E}^{2}}$ . (13)

From (9) and (13), the PDF of E can be obtained:

$\begin{array}{c}{f}_{E}\left(E\right)={f}_{c}\left[h\left(E\right)\right]|\frac{\text{d}}{\text{d}E}h\left(E\right)|\\ ={f}_{c}\left[h\left(E\right)\right]\frac{{P}_{0}\sqrt{{P}_{0}}}{\sqrt{{m}_{0}}}\frac{y}{\left(1-{y}^{2}\right)\sqrt{1-{y}^{2}}|1-3{y}^{2}|}\frac{1}{{E}^{2}}\end{array}$ (14)

As 0 < y < 1, fE(E) diverges to infinity at $y=\frac{1}{\sqrt{3}}$ . From (3), $c=\frac{\sqrt{6}}{2}V$ . This value of c is equal to cm. Then $E=\sqrt{3}{m}_{0}{c}_{m}^{2}$ .

Step 4: Setting of parameters

As the value of mo has no qualitative effect on the relationship between E and c, mo is set equal to 1. The values of V, ca and cb are set arbitrarily, depending on whether cm is within the range of c: ${c}_{a}\le c\le {c}_{b}$ .

3. Results

Case 1: ${c}_{m}\le {c}_{a}$

The object speed and limits of c are set to V = 0.5, ca = 2.2 and cb = 4.2. Then cm ≈ 0.6124. Figure 1(a) shows the relationship between E and c. The inverse function of E is a single-valued monotonically increasing function. Figure 1(b) shows the PDF of E, which decreases monotonically when c obeys the uniform distribution. Figure 1(c) shows three cases: the c-value of the vertex of the triangular probability distribution of c is 2.7, 3.2 or 3.7.

Case 2: ${c}_{a}<{c}_{m}<{c}_{b}$

The object speed and limits of c are set to V = 2, ca = 2.2 and cb = 4.2. Then cm ≈ 2.4495. Figure 2(a) shows the relationship between E and c. The inverse function of E is a two-valued function between minimum E and a certain value marked with an arrow and a single-valued function between the marked value and maximum E. In Figure 2(b), the PDF of E decreases monotonically when c has the uniform probability distribution. Figure 2(c) shows two cases: the c-value of the vertex of the triangular probability distribution of c is less than cm, (2.3 < cm), and higher than cm, (3.7 > cm).

Case 3: ${c}_{m}\ge {c}_{b}$

The object speed and limits of c are set to V = 2, ca = 2.2 and cb = 2.4. Then cm ≈ 2.4495. In Figure 3(a), c is a single-valued monotonically decreasing function. In Figure 3(b) the PDF of E again decreases monotonically when c has the uniform distribution. Figure 3(c) shows three cases: the c-value of the vertex of the triangular probability distribution of c is 2.25, 2.3 or 2.35.

4. Discussion

When the speed of light is assumed to be variable, its probability distribution is unknown. Although this study only examines two probability distributions for it, the probability distribution of E has certain characteristics. If the distribution of c is uniform, the probability density of E is always maximum at minimum E and decreases monotonically regardless of whether cm is within the range of c or not. If the distribution of c is triangular and ${c}_{m}\le {c}_{a}$ , the distribution of E is also triangle-like with the vertex moving rightward with c. If the distribution of c is triangular and ${c}_{m}\ge {c}_{b}$ , the probability distribution of E changes from trapezoid-like shape to triangle-like shape. In both cases the vertex of the probability distribution of E moves together with the probability distribution of c. In contrast with these cases, if ${c}_{a}<{c}_{m}<{c}_{b}$ , the PDF of E is always maximum at minimum E in either distribution of c. This is because the PDF of E diverges to infinity at minimum E in either distribution of c from Equation (14). In practice, the PDF of E at minimum E increases even faster as the calculation step becomes smaller. This equation shows that if ${c}_{a}<{c}_{m}<{c}_{b}$ , the PDF of E is always maximum at minimum E: $E=\sqrt{3}{m}_{0}{c}_{m}^{2}$ regardless of the distribution of c (that is, diverges to infinity). This suggests that even if the distribution of c is unknown, E will rapidly increase as soon as cm enters a certain range of c as the speed V of an object increases.

On this basis, the following method is proposed to detect any range of c. As the speed of light is defined as 299,792,458 m/s º cL, ${c}_{a}\le {c}_{L}\le {c}_{b}$ .

If

${c}_{m}=\frac{\sqrt{6}}{2}V<{c}_{a}\le {c}_{L}\le {c}_{b}$ (15)

then

$V<\frac{\sqrt{6}}{3}{c}_{L}$ . (16)

For example, the speed of one thousand electrons or protons is increased stepwise by an accelerator from below $\frac{\sqrt{6}}{3}{c}_{L}$ to cL. For each step, a frequency distribution of E will be obtained. As V is increased from below $\frac{\sqrt{6}}{3}{c}_{L}$ , cm will

enter the range of c at the critical value of V. Then the probability density of minimum E will increase sharply. Since the speed of light has been measured with very fine precision  , the range of c would be very narrow. Then the speed will need to be more finely increased bit by bit (in steps of 100 m/s if possible). As V increases after cm exceeds cb, the probability density of minimum E will decrease abruptly. If these phenomena are observed, then c is variable. If not, then c is constant.

5. Conclusion

If it is possible that the speed of light in free space is variable, then a probabilistic method to detect the variability is applicable. This assumes that c obeys a probability distribution. From mass-energy equivalence, the PDF of E can be obtained

using the inverse function of E. The energy is minimum at ${c}_{m}=\frac{\sqrt{6}}{2}V$ , and the PDF of E diverges to infinity at $E=\sqrt{3}{m}_{0}{c}_{m}^{2}$ . Thus, when cm enters the range of c as V is increased stepwise from below $\frac{\sqrt{6}}{3}{c}_{L}$ , the PDF of E increases abruptly

regardless of the PDF of c. If this is observed by accelerating a beam of electrons or photons, it will show that c is variable; otherwise, c is constant.

Acknowledgements

Mark Kurban, M.Sc., from Edanz Group (http://www.edanzediting.com/ac) edited a draft of this manuscript.

Cite this paper
Osaka, M. (2019) A Probabilistic Method to Determine Whether the Speed of Light Is Constant. Applied Mathematics, 10, 51-59. doi: 10.4236/am.2019.102005.
References
   Einstein, A. (1905) Ist die Trägheit eines Körpers von seinem Energieinhalt abhängig? Annalen der Physik, 323, 639-641.
https://doi.org/10.1002/andp.19053231314

   Michelson, A.A. and Morley, E.W. (1887) On the Relative Motion of the Earth and the Luminiferous Ether. American Journal of Science, 34, 333-345.
https://doi.org/10.2475/ajs.s3-34.203.333

   Evenson, K.M., Wells, J.S., Petersen, F.R., Danielson, B.L. and Day, G.W. (1973) Accurate Frequencies of Molecular Transitions Used in Laser Stabilization: The 3.39-μm Transition in CH4 and the 9.33 - and 10.18-μm Transitions in CO2. Applied Physics Letter, 22, 192-195.
https://doi.org/10.1063/1.1654607

   Herrmann, S., Senger, A., Möhle, K., Nagel, M., Kovalchuk, E.V. and Peters, A. (2009) Rotating Optical Cavity Experiment Testing Lorentz Invariance at the 10 - 17 Level. Physical Review D, 80, Article ID: 105011.
https://doi.org/10.1103/PhysRevD.80.105011

   Giovannini1, D., Romero1, J., Potoček, V., Ferenczi, G., Speirits, F., Barnett, S.M., Faccio, D. and Padgett, M.J. (2015) Spatially Structured Photons That Travel in Free Space Slower than the Speed of Light. Science, 347, 857-860.
https://doi.org/10.1126/science.aaa3035

Top