We have used the varying physical constant approach to resolve the flatness problem in cosmology. Friedmann equations are modified to include variability of speed of light, gravitational constant, cosmological constant and the curvature constant. The continuity equation obtained with such modifications includes scale factor dependent cosmological term as well as the curvature term along with the standard energy-momentum term. The result is that as the scale factor tends to zero (i.e. as the big-bang is approached) the universe becomes strongly curved rather than flatter and flatter in the standard cosmology. We have used the supernovae 1a redshift versus distance modulus data to determine the curvature variation parameter of the new model, which yields a better fit to the data than the standard ΛCDM model. The universe is found to be open type with radius of curvature R_c = 1.64 (1+z)^-3.3 c₀/H₀, where z is the redshift, c₀ is the current speed of light and H₀ is the Hubble constant.

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₀.

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.

Many models have been proposed to explain the intergalactic redshift using different observational data and different criteria for the goodness-of-fit of a model to the data. The purpose of this paper is to examine several suggested models using the same supernovae Ia data and gamma-ray burst (GRB) data with the same goodness-of-fit criterion and weigh them against the standard $\text{Λ}$ CDM model. We have used the redshift – distance modulus ( $z-\mu$ ) data for 580 supernovae Ia with $0.015\le z\le 1.414$ to determine the parameters for each model, and then use these model parameter to see how each model fits the sole SNe Ia data at $z=1.914$ and the GRB data up to $z=8.1$ . For the goodness-of-fit criterion, we have used the chi-square probability determined from the weighted least square sum through non-linear regression fit to the data relative to the values predicted by each model. We find that the standard ΛCDM model gives the highest chi-square probability in all cases albeit with a rather small margin over the next best model - the recently introduced nonadiabatic Einstein de Sitter model. We have made ( $z-\mu$ ) projections up to $z=1096$ for the best four models. The best two models differ in $\mu$ only by 0.328 at $z=1096$ , a tiny fraction of the measurement errors that are in the high redshift datasets.

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.