JWST early Universe observations and {\Lambda}CDM cosmology

Deep space observations of the James Webb Space Telescope (JWST) have revealed that the structure and masses of very early Universe galaxies at high redshifts (z~15), existing at ~0.3 Gyr after the BigBang, maybe as evolved as the galaxies in existence for ~10 Gyr. The JWST findings are thus in strong tension with the ${\Lambda}$CDM cosmological model. While tired light (TL) models have been shown to comply with the JWST angular galaxy size data, they cannot satisfactorily explain isotropy of the cosmic microwave background (CMB) observations or fit the supernovae distance modulus vs. redshift data well. We have developed hybrid models that include the tired light concept in the expanding universe. The hybrid ${\Lambda}$CDM model fits the supernovae type 1a data well but not the JWST observations. We present a model with covarying coupling constants (CCC), starting from the modified FLRW metric and resulting Einstein and Friedmann equations, and a CCC+TL hybrid model. They fit the Pantheon+ data admirably, and the CCC+TL model is compliant with the JWST observations. It stretches the age of the universe to 26.7 Gyr with 5.8 Gyr at z=10 and 3.5 Gyr at z=20, giving enough time to form massive galaxies. It thus resolves the 'impossible early galaxy' problem without requiring the existence of primordial black hole seeds or modified power spectrum, rapid formation of massive population III stars, and super Eddington accretion rates. One could infer the CCC model as an extension of the ${\Lambda}$CDM model with a dynamic cosmological constant.


INTRODUCTION
Observations with the James Webb Space Telescope (JWST) have revealed the existence of massive, bright galaxies in the very young universe of age ~500 million years (e.g., Naidu et al. 2022aNaidu et al. , 2022b;;Labbe et al. 2023;Curtis-Lake et al. 2023;Hainline et al. 2023;Robertson et al. 2023), which is a small fraction of its current age estimate of 13.7 billion years according to the standard ΛCDM model.Angular diameters of many such galaxies are an order of magnitude smaller than expected from the ΛCDM model (e.g., Adams et al. 2022, Atek et al. 2022, Chen et al. 2022, Donnan et al. 2022, Finkelstein et al. 2022, Naidu et al. 2022a, 2022b, Ono et al. 2022, Tacchella et al. 2022, Wu et al. 2022, Yang et al. 2022, Austin et al. 2023, Baggen et al. 2023).Astronomers first identified such an 'impossible early galaxy' problem from observations with the Hubble Space Telescope (HST) at high redshifts ~10 (Melia 2014(Melia , 2020(Melia , 2023)).Dekel et al. (2023) stated 'JWST observations reveal a surprising excess of luminous galaxies at ~10 '.According to Boyett et al. (2023) 'JWST observations confirm the existence of galaxies as early as 300 Myr and at a higher number density than expected based on galaxy formation models and HST observations'.Looser et al. (2023) observed the existence of a quiescent galaxy when the Universe was only 700 Myr old (see also Long et al. 2023).Bunker et al. (2023) wrote 'Our NIRSpec spectroscopy confirms that GN-z11 is a remarkable galaxy with extreme properties seen 430 Myr after the Big Bang' (see also Tacchella et al. 2023); NIRSpec refers to 'near infrared spectrograph' and GN-z11 nomenclature means Good North survey of galaxy with  ≅ 11.Following spectroscopic confirmation of several photometric redshifts from JWST early galaxy observations, Haro et al. (2023) stated 'our results solidifies photometric evidence for a high space density of bright galaxies at  > 8 compared to theoretical model predictions'.Related to the massive quasar analysis observed at  > 6, Eilers et al. (2023) wrote '..this quasar hosts a ten billion solar mass black hole less than 1 Gyr after the Big Bang, which is challenging to explain with current black hole formation models'.'And the findings have been dazzling astronomers, revealing that stars and galaxies were forming and evolving much earlier than anyone had expected' as expressed by Alexandra Witze (2023) in a recent news article in Nature.How come the galaxies in the very early universe were as evolved as those with billions of years of evolution, some as early as less than ~300 million years after the big bang?These astonishing observations make the problem even more acute to resolve by tweaking wellestablished galaxy formation and cosmological models developed to satisfy lower redshifts observations (e.g., Haslbauer et al. 2022;Inayoshi et al. 2022;Kannan et al. 2022;Keller et al. 2022;Regan 2022;Yajima et al. 2022;Atek et al. 2023, Mason et al. 2023;Mirocha & Furlanetto 2023;Whitler et al. 2023aWhitler et al. , 2023b;;McCaffrey et al. 2023).Attempts have been made to compress time for the formation of population III stars and galaxies more and more, such as by considering the presence of primordial massive black hole seeds, and super-Eddington accretion rates in the early Universe (Ellis 2022, Bastian et al. 2023, Brummel-Smith 2023, Chantavat et al. 2023, Dolgov 2023, Larson et al. 2023, Maiolino et al. 2023).While analysing GN-Z11 JWST-NIRSpec data, Maiolino et al. (2023) concluded that the black hole seed of this exceptionally luminous galaxy at  = 10.6 must be accreting at an episodic rate of about five time the Eddington rate for 100 million years (since ~12 − 15) and is challenging for theoretical models (see also Schneider et al. 2023).Chen et al. (2023) considered the presence of massive dark matter halos (see also Mauerhofer and Dayal).As discussed elegantly by Melisa (2023), both are considered unrealistic based on the zero angular momentum argument and observations showing most of the distant quasars accreting below or at the Eddington limit (Melia 2023).
As stated by Wang and Liu (2023) 'JWST high redshift galaxy observations have a strong tension with Planck CMB measurements'.They could not resolve the tension using alternative cosmological models either, including dark matter-baryon interaction, f(R) gravity, and dynamical dark energy (see also Santini et al. 2023).Parashari and Laha (2023) have shown that a blue-tilted power spectrum could potentially alleviate the tension with a low to moderate star formation efficiency.Studying the Balmer breaks at highest redshift galaxies, Steinhardt et al (2023) proposed it as a test of the ΛCDM model.They emphasized that the existence of stronger Balmer breaks out to  ≳ 11, will either demonstrate the early galaxy formation templates are invalid at high redshift or imply new physics beyond the 'vanilla ' ΛCDM model. Lovyagin et al. (2022) have shown that the impossible early galaxy problem can be resolved amicably with the tired light theory in a steady-state Universe, first advocated by Zwicky (1929) to explain the early redshift observations by Hubble (1929).They have succinctly reviewed the current status of the problem and suggested some solutions.However, the tired light model cannot explain the extreme directional uniformity of the observed cosmic microwave background (CMB) radiation (Penzias and Wilson 1965).Additionally, the standard tired light model does not fit the supernovae type 1a data except at very low redshifts.An expanding Universe model can easily account for the observed redshift of distant galaxies and the CMB isotropy but has problems with the early Universe observations.However, the two approaches are not mutually exclusive.Gupta (2018aGupta ( , 2018b) ) has contemplated the existence of tired light in an expanding Universe.Our main objective in this paper is to explore if hybrid tired light and expanding Universe models can resolve the impossible early galaxy problem without conflicting with the relatively low redshift  ≤ 2.5 data from supernovae type 1a (SNe Ia) observations in distant galaxies such as the Pantheon+ data (Scolnic et al. 2022).Since the distance traveled by light reaching us is the same irrespective of whether in the tired light scenario or the expanding universe, it constrains the parameters in them; one does not need an extra parameter to include tired light in an expanding Universe model.We will use two expanding Universe models to develop hybrid models, compare them, and see how they fit the SNe Ia and the galaxy size data from HST and JWST observations.The first hybrid model combines the ΛCDM model with the tired light model and is named the ΛCDM+TL model.The second hybrid model incorporates a new model derived from covarying coupling constant (CCC) approach (Gupta 2022).We start by defining an appropriate metric for the CCC approach and develop Einstein and Friedmann equations for the new CCC cosmological model; the hybrid model is dubbed the CCC+TL model.
Section II of this paper develops the theory of the expanding Universe models and their hybrid counterparts, especially for the CCC model, and shows how it relates to the standard ΛCDM model.In Section III, we attempt to fit the Pantheon+ data (Scolnic et al. 2022, Brout et al. 2022) for the models developed in Section II to qualify them for testing further with the HST and JWST data.We discuss our findings in Section IV and present the conclusions in Section V.

II. THEORY
Basics: We will derive the Friedmann equations starting from the Einstein equations and see how they are modified when  and  are varying such that ~ 3 and the speed of light is used to measure distances (Gupta 2022).The basic Einstein equations (without the cosmological constant) are where   is the Einstein tensor and   is the stress-energy tensor.The FLRW metric takes the form where () = exp(( −  0 ) defines the variation of  and the distances are measured using the speed of light;  0 is the current time when ( 0 ) = 1,  is the scale factor, and  (−1, 0, 1) determines the geometry of the space.The most significant difference from the standard FLRW metric is that the  00 metric coefficient is time dependent (Gomide & Oehara 1981).Since  ̇ ⁄ =  and  ̈ =  2 ⁄ , and   is a diagonal matrix, the essential components of the Einstein tensor are 2 ), and (3) The stress-energy tensor for a homogeneous and isotropic universe considered a perfect fluid, has the form Here  is the energy density,  is the pressure,   are the metric elements of the FLRW metric, and the fourvelocity   are related through       =  2 .For a comoving observer,   = (, 0,0,0).Since   satisfies the contracted Bianchi identities and   obeys the local conservation laws, we have () 4   ) = 0, and ∇    = 0.

CCC Model:
Defining the Hubble expansion parameter as  = ̇ ⁄ , we may write Eq. ( 11) for a flat universe ( = 0) as In the matter-dominated universe of interest to us in this work, using Eq. ( 14), we may write Since the observations are made using redshift , we have to see how the scale factor  relates to  in the CCC model.Along the spatial geodesic ( and  constant) between the observer and the source at a fixed time  using the modified FLRW metric (Eq.2) Thus, the proper distance for commoving coordinate  [recall that ( 0 ) = 1 = ( 0 )] Since the light follows the null geodesic, Eq. ( 2) for a light emitted by a source at a time   and detected by the observer at a time  0 yields Thus, the expression for the proper distance is the same when using the standard FLRW metric.It can now be shown (Ryden 2017) that  = 1 (1 + ) ⁄ , i.e., the same as for the ΛCDM model.The next thing to consider is to transpose () to (), as it is the latter that we will require in calculating the proper distance.We may write Eq. ( 20) Its analytic solution (using WolframAlpha), with the boundary condition  = 1 at  =  0 , is It can be written as a cubic equation Its solution is Since the scale factor Thus, the above equation provides the function (,  0 , ).
Next, to determine the luminosity distance   of an object, we need to understand how the photon energy flux evolves with the redshift in the CCC model.The photon energy is reduced by a factor (1 + ).The effect of the interval between photo arrival time on the flux is an additional factor (1 + ).Since in CCC, the distance scales as the speed of light, their evolution cancels out: two emitted photons separated in time by   are separated in space by a distance   =     .This distance becomes  0 =  0   (1 + ) and the corresponding time interval  0 =  0   (1 + )  0 ⁄ =   (1 + ).Thus the luminosity distance scaling is the same as for the ΛCDM model, i.e.   =   ( 0 )(1 + ), and the expression for the distance modulus  also the same, viz.,  = 5 log 10 (  /1Mpc) + 5 log 10 (1 + ) + 25. (28) We will now focus on determining the expression for the proper distance   .Since  =  × / = /̇, we may write Eq. ( 23) and since  = 1 (1 + ) ⁄ ,  = −  (1 + ) 2 = − 2 ⁄ , we get, using Eq. ( 20) CDM model: The corresponding expression for the ΛCDM model (matter dominant, flat universe) is (Ryden 2017), where Ω ,0 is the relative matter energy density and (1 − Ω ,0 ) is the relative dark energy density.Notice that both the models have only two unknowns, i.e.,  0 and  for the CCC model, and  0 and Ω ,0 for the ΛCDM model.

Tired Light Model:
We would also like to consider the tired light model as it fits admirably well the high redshift galaxy size data (Lovyagin et al. 2022), especially from the James Webb Space Telescope (JWST).In this approach, one assumes that the photon energy loss, ℎ, is proportional to its energy, ℎ, and the distance traveled, , i.e.,  = , with  as the proportionality constant.Thus, Here   and  0 are the emitted and observed photon frequencies and   and  0 are the respective wavelengths.The redshift is defined as  = ( 0 −   )/  , i.e., (1 + ) =  0   ⁄ .Therefore, for the tired light It must reduce to Hubble's law  = / 0 in the limit of  ≪ 1, yielding  =  0 /.Thus, the proper distance in the tired light (TL) model is given by Hybrid Models: Since the distance traveled is the same in any model, the parameters of a hybrid tired light and expanding universe model are obtained by equating   for the two.Thus for the hybrid CCC+TL model, we equate Eqs. ( 30) and ( 34), and write (with subscript  for tired light and  for the CCC model) We will now determine   and   by considering the above expression in the limit of very low redshifts ( ≪ 1) and keeping the terms with up to second order in .Since the left-hand side involves an integral, we need to retain only up to first-order terms in the integrand.
Let us first consider ().It is not practical to determine it from Eq. ( 27) in the small redshift limit  = 0. We go back to Eq. ( 24) and rewrite it with  = 1 (1 + ) ⁄ : Since  −(1 2 ⁄ ) varies slowly as compared to (1 + ) 5 2 ⁄ , we can assume it is a constant equal to unity (i.e., its value at  = 0) in the above equation for solving it.With boundary condition  =  0 at  = 0, the solution in the limit of  ≪ 1 up to the first term in  can be easily determined: Since  = −  for the CCC case, the above expression for  reduces to The left-hand side of Eq. ( 35): Thus, by expanding the right-hand side, we may write Eq. ( 35) as ). (43) Let us now consider Eq. ( 35).With the above findings, we may write Given   and , this equation determines (1 +   ) for any (1 + ) value and, therefore, also (1 +   ).We can now use either of Eq. ( 31) or Eq. ( 34) to determine the proper distance   while replacing  0 and  in there with   and   or   and   as applicable.When determining the luminosity distance, only (1 +   ) is involved in time dilation, not (1 +   ).The luminosity distance in the hybrid model is then We can follow a similar approach for the ΛCDM+TL model (with subscript  for the ΛCDM model).Equating Eq. ( 31) and Eq. ( 34) In the limit of   ≪ 1 and   ≪ 1, and keeping the terms up to second order in redshift, we get the same results as for the CCC model by comparing the first-order terms but   = 1.5Ω ,0   by comparing the second-order terms, and Eq. ( 44) is replaced with

Fitting Pantheon+ Data:
We will now attempt fitting the Pantheon+ data (Scolnic et al. 2022, Brout et al. 2022) for supernovae type 1a (SNe1a).We have considered the following six models in matter dominated Euclidean Universe: 1.The ΛCDM model.
3. The hybrid model, comprising the ΛCDM model and the tired light model labeled ΛCDM+TL.
4. The hybrid model, comprising the CCC model and the tired light model labeled CCC+TL. 5.The tired light model TL, modified to include an unknown flux loss correction term, labeled TL+, added to the luminosity distance, i.e., 5 log(1 + )  , with  a data-fit parameter.6.The standard tired light model, TL, with luminosity distance corrected only for the photon energy loss since the time dilation is not relevant for it.We used the standard Curve Fitting tool in Matlab for this purpose by minimizing  2 (Gupta 2020).The first five have two free parameters; thus, the degrees of freedom are identical.They are all within the error bars of the data (Fig. 1; only the first four model curves are plotted).The  2 values of the first four models are within 0.1% of each other at 745.7 as its mean value (Table 1).Even the TL+ model's  2 value is within 0.6% of the mean value.However, as could be expected, the  2 value for the standard TL model is significantly larger at 2580.Thus, the Pantheon+ data rejects only the TL model but none of the others.
It should be mentioned that we have not taken into consideration the systematic uncertainties associated with the data.However, they would affect the fitting of all the models and  2 values similarly and thus are not expected to impact the comparative study of the models.The models' fit parameters and  2 values, along with the deceleration parameter and the age of the universe, are presented in Table 1.It is relevant at this point to refer to the work of Lopez-Corredoia and Calvo-Torel (2022) wherein they compare the fits of several alternative models without dark energy by fitting Pantheon data (Scolnic et al. 2018).Unlike in a previous paper (Gupta 2020), wherein the absolute magnitude   = −19.35was used to convert the apparent magnitudes reported in the Pantheon data (Scolnic 2018) to the distance moduli , Pantheon+ data (Scolnic 2022) includes distance modulus, apparently based on   = −19.253.Since  0 is degenerate with   , the choice of   affects the determination of  0 for different models.However, all the models yield about the same value of  0 .Thus, we do not expect the choice of absolute magnitude to affect the findings of this paper.

Angular-Diameter Angle and Distance:
These are essential for testing a model based on the data on the size of galaxies at high redshift, especially those recently observed with the James Webb Space Telescope and discussed extensively in the literature (e.g., Lovyagin et al. 2022 and references therein).
The angular diameter distance   is defined in terms of the physical size  of an object and its observed angular size  as   =   ⁄ .Using the metric (Eq.2), the object at a location (, ), i.e.,  = 0 and  = 0, at time  has a size given by  2 = () 2 () 2  2  2 ⇒  = ()() =  (48) Therefore, with  as the proper distance   , the angular diameter distance becomes We have to be cogent of the fact that the scale factor  in a hybrid model relates to the redshift   due to expanding universe only, i.e.,  = 1 (1 +   ) ⁄ , whereas an observer measures the total redshift .The function (), of course, has no significance for the non-CCC models.In addition, () is irrelevant for the tired light model.Thus, when we know   for a model, we can immediately compute   , and, therefore, the size of an object from its observed angular size.
Plots of angular-diameter distance for the five models are shown in Figure 2. Two of the plots jump out: i) the TL+ plot has no maximum, but that can be expected as no expansion of the universe is involved, and ii) the CCC+TL plot that has its peak value about four times higher than the ΛCDM value, and that too at  > 10.This dramatically affects the observed  and brightness of the distant galaxies.The plots are shown in Fig. 3 for 10 kpc objects against the backdrop of measured angular sizes of galaxies from multiple sources (the same as cited by Lovyagin et al. 2022; those points without errors in the original data are arbitrarily allowed a 5% error), including the latest JWST data.While all models are satisfactory at low redshifts, only TL+ and CCC+TL are acceptable for high redshifts.
It should be mentioned that the size evolution of galaxies, especially at z > 10, is not well measured.In particular, the measured angular sizes are subject to the distribution of stellar populations within galaxies and their light or mass profiles (e.g.outskirts of galaxies are likely too faint to be detected).Also, we have not addressed the uncertainties related to baryonic processes and feedback mechanisms that may potentially effective in shaping galaxy sizes of observed galaxies.Additionally, not all redshifts determined initially via photometric color selection, have been confirmed spectroscopically; some have been found spectroscopically at much lower redshift, e.g., a bright interloper at z spec = 4.91 that was claimed as a photometric candidate at z ∼ 16 (Harikane et al. 2023).

Age of The Universe and Redshift:
Galaxies at high redshifts, especially those observed by JWST, appear as evolved and as massive as those at lower redshifts.How can this be possible, considering that the age of the universe at very high redshifts was less than half a billion years as per the ΛCDM model?This subject has been extensively discussed in the literature (e.g., Adams et al. 2022, Atek et al. 2022, Chen et al. 2022, Donnan et al. 2022, Finkelstein et al. 2022, Naidu et al. 2022aand 2022b, Ono et al. 2022, Tacchella et al. 2022, Wu et al. 2022, Yang et al. 2022, Labbe et al. 2023), so we will not delve on it.Nevertheless, we will try to see how the Universe ages under different models.
Let us first consider the CCC model.As per Eq. ( 24), we may write Here   = 1 for the age of the universe today at   = 0. To find the age for a specific   , we replace 1 (1 +   ) ⁄ with   and integrate.For non-hybrid models,  replaces   .When working with a CCC+TL model, we have to find first the expanding universe component   of  as described under Eq.( 44), and then use Eq. ( 50) to determine the age corresponding to the observed .Similarly, we determine the ages for the ΛCDM and ΛCDM+TL models.The Universe age has no meaning for the TL and TL+ models.
Figure 4 shows how the age of the universe decreases with increasing redshift for the four models.Figure 5 shows the age increment with redshift for the two hybrid models compared to the currently expected age for the standard ΛCDM model.Both the hybrid models show significant increases, but the CCC+TL model provides a 10 to 20-fold increase at redshifts 10 to 20, giving enough time (5.8 Gyr at  = 10 and 3.5 Gyr at  = 20) for large, massive galaxies to form.In contrast, for the ΛCDM+TL model, corresponding age increases are 1.7 Gyr and 0.7 Gyr.

Tired Light Contribution to Redshift:
We were curious to know how the ratio of the tired light and expanding universe redshift evolves with the total redshift for the hybrid models.It is plotted in Figure 6.For the CCC+TL model, the tired light fraction starts at 22% at  ≈ 0 and declines slowly to 20% at  = 10, to 18% at  = 20, and to 3% at  = 1000.For the ΛCDM+TL model, the corresponding values are 20%, 11%, 7%, and 0.3%.Despite their relatively small contribution, the tired light modifies the observed redshift considerably since redshifts are multiplicative through the relation From the second Friedmann equation (Eq.12) in the matter-dominated universe, using the above equation Similarly,  0 for the ΛCDM model may be expressed (Ryden 2017), since Ω Λ,0 = 1 − Ω ,0 , as In the above expressions, we must use the expanding Universe component of the Hubble constant when applying them to the hybrid models.The deceleration parameters and the universe's ages for different models are shown in Table 1.

IV. DISCUSSION
The main objective of this paper is to explore if the hybrid models that include tired light cosmology can explain the deep space high redshift observation of JWST on the large-scale structures of the universe, that is, the structure and evolution of galaxies.These observations show that the very early Universe galaxies were almost as bright, massive, and structurally as evolved as the galaxies in the late universe but with rather small angular diameters.Under the standard ΛCDM model, the physical sizes of these galaxies turn out to be about 10% of their expected sizes.It is because, according to this model, the angular diameter distance of objects reaches a maximum of about ~1.7 Mpc at  = 1.6 and decreases at higher redshift leading to the angular diameter  of a standard size object increasing with increasing .Thus, if one observes a smaller , it is interpreted as a smaller size object.So, if a model has no angular diameter distance maximum, such as for the tired light models, or if the maximum is shifted to higher redshifts and thus has a higher value, then observed smaller  will translate into larger object sizes.Observed over density of galaxies in the region of ~10 (e.g.Whitler et al. 2023aWhitler et al. , 2023b) ) can also be explained accordingly.
However, before any model can be considered seriously, it must pass the most basic test fitting the supernovae type 1a data compiled as Pantheon+ by Scolnic et al. (2022).Additionally, it should, in principle, not conflict with the observed isotropy of the CMB radiation.The standard tired light model TL does not pass the basic test, as is evident from its  2 value of 2580 against all others less than 750 in Table 1.The modified tired light model TL+ with added unknown flux loss term 5 log(1 + )  in the distance modulus equation gives a good fit to Pantheon+ with a  2 value of 749.9 when  = 0.6418 ± 0.0196.The inclusion of the additional term indicates that the TL model is missing something and that missing something is possibly the expanding universe component.However, neither of the tired light models can explain the extreme isotropy of the CMB radiation.Nevertheless, we have included the TL+ model in some of the comparative figures for discussion purposes.All the remaining models provide similarly great fits to Pantheon+ data (Fig. 1 and Table 1) and involve expansion of the universe and thus are compliant with CMB isotropy.
One may be concerned about how the peak luminosities of supernovae type 1a in the Pantheon+ data will be affected due to the variation of  and other constants.The peak luminosity can be related to the Chandrasekhar mass  ℎ (Arnett 1982, Garcia-Berro et al 1999, Gaztañaga et al. 2001, Wright & Li 2018) with  ℎ ~(ℎ  ⁄ ) 3/2 .Since ℎ~ 2 , ~, and ~ 3 , we find that  ℎ ~0 .Therefore, the coupling constants' variation does not affect the peak supernovae luminosity in the CCC model1 .
Coming back to angular diameter distance, its peak value of 7.45 Gpc at  = 12.5 is the highest among the four models of interest, and its profile is closest to the TL+ profile (Fig. 2).Thus, a small observed angular size  of a galaxy does not translate into a smaller galaxy size.This is depicted in Fig. 3, which shows how one would observe  of a 10 kpc object with increasing redshift under different models; the curves are superimposed over the  of many galaxies from JWST observations and some pre-JWST observations.We notice that the CCC+TL model curve is very close to the TL+ curve, whereas the ΛCDM+TL model curve is not.This may mean that the ΛCDM+TL model (also, the CCC and ΛCDM models) do not faithfully represent the early universe evolution.
An alternative way to analyse the angular size data of the galaxies would be to see what each model yields as their physical size.Since the physical size is directly proportional to angular diameter distance, the size increase in different models compared to the ΛCDM model can be easily determined, as shown in Fig. 7.We see an insignificant difference using different models at redshifts less than 1.The difference becomes perceptible at higher and higher redshift.The multiplier for the CCC+Tl model becoming 5.6 at  = 5, 9.5 at  = 10, 12.8 at  = 15, and so on.Astronomers became concerned about the galaxy sizes using the ΛCDM model only at  ≥ 5, especially at  ≥ 10.Consequently, they suggested models and their improvements to create large structures of the universe in a shorter and shorter time (e.g., Adams et al. 2022, Atek et al. 2022, Chen et al. 2022, Finkelstein et al. 2022, Naidu et al. 2022, Ono et al. 2022, Tacchella et al. 2022, Wu et al. 2022, Yang et al. 2022, Labbe et al. 2023, Yung et al. 2023, Chantavat et al. 2023).Instead of attempting to compress the cosmic timeline for creating well-evolved massive galaxies in a very young universe to resolve the impossible early galaxy problem with the ΛCDM model (Melia 2023) and cosmological simulations (McCaffrey et al. 2023), it might be prudent to consider alternative models that stretch the timeline.An alternative scenario in support of the standard model has been proposed by Prada et al. (2023) by contending that during the early epochs of the universe the stellar mass-to-light ratio could not have reached the values reported by Labbe et al. (2023).However, one would need to explain the existence of very high  quasars, such as those considered by Latif et al. ( 2023), and remarkably high IR luminosity of massive galaxies at ~8 discovered by JWST (Akins et al. 2023).
Figure 7. Object size under different models as a multiple of its size using the ΛCDM model.
We have shown in Fig. 4 how the cosmic timeline is stretched in the two hybrid models.We show the age of the universe under the two hybrid models relative to the age under the ΛCDM model.While the stretch is modest (1.40 times for the ΛCDM+TL model vs. 1.94 times for the CCC+TL model) at the current epoch, it becomes very significant at high redshift (2.97,3.56,and 4.16 times at 𝑧 = 5,10,and 20,respectively,for the ΛCDM+TL model,and 7.76,12.3,and 19.3 times for the same redshifts, respectively, for the CCC+TL model).The cosmic time stretch is thus up to 5.8 Gyr at  = 10 and 3.5 Gyr at  = 20; there is no need to invent new physics for the rapid formation of galaxies.
It would appear that except for the CCC+TL model, none of the models are able to reasonably resolve the impossible early galaxy problem; the TL models are unrealistic as they cannot reproduce the isotropy of CMB, and others do not provide realistic sizes of the early galaxies.Thus, the only realistic option is the CCC+TL hybrid model.
The CCC model is derived from Dirac's hypothesis of varying physical constants (Dirac 1937).His analysis was based on his large number theory and predicted the variation of  and the fine structure constant .As discussed by Uzan (2011), if one constant varies, then others also must vary, and their variations correlated.From the local energy conservation laws, Gupta (2022) found that the variation of , , the Planck constant ℎ, and the Boltzmann constant  must follow ~ 3 ~ℎ1.5 ~1.5 when the distance is measured with .This leads to  ̇ ⁄ = 3 ̇ ⁄ = 1.5 ℎ ̇ℎ ⁄ = 1.5  ̇ ⁄ , and therefore their variation can be represented by a single dimensionless function (), i.e., if one constant varies, then all of them do.The standard procedure is to consider a power law function or a linear function.However, we found the most convenient form to use is () = exp[( −  0 )] with  determined by fitting the Pantheon+ data as  = −0.660 (= ̇ ⁄ );  may be considered to represent the strength of the coupling constants' variation.Since  ̇ ⁄ = 3 ̇ ⁄ , we get  ̇ ⁄ = −2.00 , about the same predicted by Dirac from his large number hypothesis.While there have been multiple attempts to constrain |  ̇ ⁄ | by various methods to several orders of magnitude lower values (e.g.,Teller 1948;Morrison 1973;Sisterna & Vucetich 1990;Degl'Innocenti et al. 1995;Thorsett 1996;Corsico et al. 2013;Sahini & Shtanov 2014;Ooba et al. 2017;Hofmann & Müller 2018;  Genova et al. 2018;Wright & Li 2018;Bellinger & Christensen-Dalsgaard 2019;Zhu et al. 2019;Alvey et al. 2020;Vijaykumar et al. 2021) they all consider other constants to be pegged to their current value.However, keeping any of the constants fixed automatically forces  = 0, and therefore  ̇ ⁄ = 0. We believe the two-parameter ( 0 , ) CCC model is possibly an extension of the ΛCDM model.Examining the Friedmann equations, one can imagine the CCC model as the ΛCDM model with a dynamic cosmological constant: (i) Λ → Λ() = −3( 2 + 2 ̇ ⁄ ), consisting of a static term and a dynamic term in the first Friedmann equation, and (ii) only the dynamic term in the second Friedmann equation.In the CCC model, recalling that  numerically turns out to be negative, it replaces Λ as the source of Universe expansion.One can even define energy density corresponding to  2 + 2 ̇ ⁄ .However, it remains to be seen how well the CCC model can fit the Planck data and explain other astrophysical and cosmological observations.
It should be mentioned that by treating the  00 metric coefficient time dependent in the FLRW metric (Eq.2), Gomide & Uehara (1981) were able to show that the local inertial effects are dependent on the overall structure of the Universe; cosmological models with positive curvature are Machian, whereas open ones are not.The CCC+TL model could, thus, be used to test the applicability of Mach's effect in the Universe.
It remains to be seen if the CCC+TL model can help address or explain other cosmology problems.Initially, we would want to determine if the new model is satisfactory for explaining CMB, BBN (big-bang nucleosynthesis), and BAO (baryonic acoustic oscillations) observations.

V. CONCLUSION
JWST is perhaps playing the same role as HST did in the 1990s -reinventing cosmology.HST put the ΛCDM model on the pedestal.JWST is challenging standard ΛCDM.In this paper, we have attempted to show that an extension of the ΛCDM model with deemed dynamical cosmological constant, when hybridized with the tired light concept and parameterized with Pantheon+ data, provides a model, dubbed CCC+TL, that is compliant with the deep space observation of JWST on the angular sizes of high redshift galaxies.It stretches the cosmic time, especially at high redshifts, to allow the formation of large galaxies.It eliminates the need for stretching and tuning existing models to produce such structures in the early universe, thus amicably resolving the impossible early galaxy problem.
=  0  −3 , and  =  0  −4 , (18) immediately reveals: (a) the cosmological constant Λ of the ΛCDM model is replaced with the constant  in the CCC model, (b) the continuity equation has an additional term involving  in the CCC model, and (c) energy density evolution has extra factors that must be considered for the CCC model.

Figure 1
Figure 1 Pantheon+ (Scolnic et al. 2022, Brout et al. 2022) data fit for the four models discussed in the text.

Figure 2 .
Figure 2. Angular diameter distance according to different models plotted against observed redshift.

Figure 3 .
Figure 3.The angular size of 10 kpc objects for the five models against the backdrop of measured angular sizes of galaxies from multiple sources, including the latest JWST data and some pre-JWST data labeled as NJWST (provided by Lovyagin et al. 2022).

Figure 4 .
Figure 4. Age as a function of the observed redshift for the four models.

Figure 5 .
Figure 5. Age advantage of the two hybrid models over the ΛCDM model

Deceleration Parameter :
We will now calculate the deceleration parameter for the models.It is defined by the relation  0 = −̈( 0 2 ⁄ ),  0 being the expanding Universe part of the Hubble constant at the current time.Let us first consider the CCC model.From the first Friedmann equation (Eq.11) for a flat Universe

TABLE 1 .
Pantheon+ data fit parameters,  2 values, the deceleration parameter, and the age of the universe.