We have shown that the use of evolutionary gravitational constants G, the speed of light c, and cosmological constant Λ in the Einstein’s field equations leads to a very simple model that fits the supernovae Ia data with a single parameter as well as the standard ΛCDM model with two parameters, and has the predictive capability superior to the latter. We have shown that at the current time G and c both increase as dG/dt = 5.4GH₀ and dc/dt = 1.8cH₀ with H₀ as the Hubble constant, and Λ decreases as dΛ/dt = −1.2ΛH₀. Negative finding on the variation of G with time derived from the analysis of lunar laser ranging data in several studies has been attributed to the fact that these studies did not consider the variation of c simultaneously with G.

We have developed a cosmological model by allowing the speed of light c, gravitational constant G and cosmological constant Λ in the Einstein filed equation to vary in time, and solved them for Robertson-Walker metric. Assuming the universe is flat and matter dominant at present, we obtain a simple model that can fit the supernovae 1a data with a single parameter almost as well as the standard ΛCDM model with two parameters, and has the predictive capability superior to the latter. The model, together with the null results for the variation of G from the analysis of lunar laser ranging data determines that at the current time G and c both increase as dG/dt = 5.4GH0 and dc/dt = 1.8cH0 with H0 as the Hubble constant, and Λ decreases as dΛ/dt = -1.2ΛH0. This variation of G and c is all what is needed to account for the Pioneer anomaly, the anomalous secular increase of the Moon eccentricity, and the anomalous secular increase of the astronomical unit. We also show that the Planck’s constant ħ increases as dħ/dt = 1.8ħH0 and the ratio D of any Hubble unit to the corresponding Planck units increases as dD/dt = 1.5DH0. We have shown that it is essential to consider the variation of all the physical constants that may be involved directly or indirectly in a measurement rather than only the one whose variation is being considered.

We have shown that our varying physical constant model is consistent with the recently published variational approach wherein Einstein equations are modified to include the variation of the speed of light c, gravitational constant G and cosmological constant Λ using the Einstein-Hilbert action. The general constraint resulting from satisfying the local conservation laws and contracted Bianchi identities provides the freedom to choose the form of the variation of the constants as well as how their variations are related. When we choose dG/Gdt=3dc/cdt, ̇the same as in our quasi-phenomenological model, c=c₀ exp⁡(a^α-1), G=G₀ exp⁡[3(a^α-1)] and Λ=Λ₀ exp⁡[(a^-α-1)], where a is the scale factor and α=1.8, we are able to confirm the success of our the model in explaining three astrometric anomalies and the null results on the variation of G and the fine structure constant. We show that the model: (a) fits the supernovae 1a observational data better than the ΛCDM model; (b) determines the first peak in the power spectrum of the cosmic microwave background temperature anisotropies at multipole value of l=217.3; (c) calculates the age of the universe as 14.1 Gyr; and (d) finds the BAO acoustic scale to be 145.2 Mpc. These numbers are within less than 3% percent of the observed values and the values obtained by the ΛCDM model. Surprisingly we find that the dark-energy density is negative in a universe that has significant negative curvature and whose expansion is accelerating at a faster rate than predicted by the ΛCDM model.

We have shown that three astrometric solar-system anomalies can be explained satisfactorily by using evolutionary gravitational constant G and speed of light c in the Einstein’s field equation. These are: a) the Pioneer acceleration anomaly; b) the anomalous secular increase of Moon-orbit eccentricity; and c) the anomalous secular change in the astronomical unit AU. The gravitational constant G and the speed of light c both increase as dG/dt = 5.4GH₀ and dc/dt = 1.8cH₀ with H₀ as the Hubble constant. We also show that the Planck’s constant ħ increases as dħ/dt = 1.8ħH₀. Additionally, the new approach fits the supernovae Ia redshift vs distance modulus data as well as the standard ΛCDM model with just one adjustable parameter H₀.