Cosmological Constant, Fine Structure Constant and Beyond

In the present work, we consider the cosmological constant model $\Lambda\propto\alpha^{-6}$, which is well motivated from three independent approaches. As is well known, the hint of varying fine structure constant $\alpha$ was found in 1998. If $\Lambda\propto\alpha^{-6}$ is right, it means that the cosmological constant $\Lambda$ should also be varying. Here, we try to develop a suitable framework to model this varying cosmological constant $\Lambda\propto\alpha^{-6}$, in which we view it from an interacting vacuum energy perspective. Then, we consider the observational constraints on these models by using the 293 $\Delta\alpha/\alpha$ data from the absorption systems in the spectra of distant quasars. We find that the model parameters can be tightly constrained to the very narrow ranges of ${\cal O}(10^{-5})$ typically. On the other hand, we can also view the varying cosmological constant model $\Lambda\propto\alpha^{-6}$ from another perspective, namely it can be equivalent to a model containing"dark energy"and"warm dark matter", but there is no interaction between them. We find that this is also fully consistent with the observational constraints on warm dark matter.


I. INTRODUCTION
The cosmological constant has been one of the long-standing issues in physics and cosmology since it was introduced by Einstein in 1917 [1] for a static universe. However, Hubble discovered in 1929 [2] that the universe is expanding. Then, Einstein abandoned the cosmological constant as the "biggest blunder" of his life [3]. From 1929 until the early 1990s, most physicists and cosmologists assumed the cosmological constant to be zero. Since the vacuum energy is equivalent to the cosmological constant [4], an exactly zero cosmological constant requires that the bare cosmological constant should be exactly canceled by the vacuum energy. This is a difficult problem [5] (sometimes it is called the (old) cosmological constant problem in the literature). In 1998, the accelerated expansion of the universe was discovered [6], and since then dark energy has been one of the most active fields in cosmology [7,8]. So, the cosmological constant was revived again, since the simplest candidate of dark energy is a tiny positive cosmological constant. However, it is difficult to understand why the observable cosmological constant is about 120 orders of magnitude smaller than its natural expectation of the vacuum energy [5,7]. Now, the (new) cosmological constant problem becomes the question why the non-zero cosmological constant is so tiny. It means that a fine-tuning is necessary when the bare cosmological constant is canceled by the vacuum energy [5,7]. In fact, the cosmological constant is still an important topic in physics and cosmology by now.
It is commonly believed that the cosmological constant problem can only be solved ultimately in a unified theory of quantum gravity and the standard model of electroweak and strong interactions, which is still absent so far. Nevertheless, many attempts have been made in the literature. One of the interesting ideas is the so-called axiomatic approach to the cosmological constant [9]. In this approach, the cosmological constant is derived from four axioms, but the underlying physical origin (say, the theory of quantum gravity) is still unknown. It is proposed in close analogy to the Khinchin axioms in information theory. The well-known Khinchin axioms can uniquely derive the Shannon entropy, on which the entire mechanism of statistical mechanics is based (see the textbook e.g. [10]).
The Khinchin axioms in information theory [11] describe the most desirable properties an information measure I should have. Axiom K1 "fundamentality": an information measure I only depends on the probabilities p i (the fundamental quantities) of the events under consideration and nothing else. Axiom K2 "boundedness": there is a lower bound for the information measure I. Axiom K3 "simplicity": the information measure I should take the simplest description. Axiom K4 "invariance": there is a suitable scale transformation in the space of probabilities and information measures that leaves the physical contents invariant. These four Khinchin axioms look very natural and simple. However, from such natural and simple axioms, one can uniquely fix the functional form of the Shannon information which is extremely important for the statistical mechanics [10].
Inspired by the successful Khinchin axiomatic approach to the Shannon entropy in information theory, Beck [9] proposed four axioms in close analogy to the Khinchin axioms. Axiom B1 "fundamentality": the cosmological constant Λ only depends on fundamental constants of nature. Axiom B2 "boundedness": the cosmological constant is bounded from below, Λ > 0. Axiom B3 "simplicity": the cosmological constant Λ is given by the simplest possible formula consistent with the other axioms. Axiom B4 "invariance": the cosmological constant Λ formed with potentially different values of fundamental parameters leaves the large-scale physics of the universe scale invariant. These four axioms look also very natural and simple. From these natural and simple axioms, Beck [9] derived the explicit form of the cosmological constant, where α is the fine structure constant, G is the gravitational constant,h is the reduced Planck constant, and m e is the electron mass. Accordingly, the (observable) vacuum energy density is given by [9] where c is the speed of light. Numerically, this formula yields ρ Λ ≃ 4.0961 GeV/m 3 , which can pass the current observational constraints with flying colors. We refer to [9] for the detailed derivations.
In fact, Beck [9] is not the first and the only one who derived the cosmological constant given in Eq. (1). It was independently derived from other arguments in the literature. Using the generalized Buchdahl identity, Boehmer and Harko [12] argued that the existence of a non-negative Λ imposes a lower bound on the mass M and density ρ for general relativistic objects with radius R, On the other hand, Wesson [13] argued that the mass is quantized according to the rule m = (nh/c) Λ/3, and the minimum mass corresponding to the ground state n = 1 is given by which is indeed a very small mass. Boehmer and Harko [12,14] proposed to identify the minimum mass in Eq. (3) with the one in Eq. (4), and found that the radius corresponding to m P is given by where ℓ pl is the Planck length. Noting the radius R P is of the same order of magnitude as the classical radius of the electron r e = e 2 /(m e c 2 ) (where e is the electron charge), Boehmer and Harko [14] further proposed to formally equate R P with r e while the term of O(1) is neglected, and then they found that the cosmological constant is given by [14] in which we have used the definition α = e 2 /(hc). Clearly, the same result given in Eq. (1) has been independently derived from completely different arguments. The third independent approach to derive the cosmological constant Λ in Eq. (1) is the well-known Eddington-Dirac large number hypothesis [15][16][17]. Nottale in 1993 [18] (see also [9]) has written down a large number hypothesis connecting cosmological parameters with standard model parameters, where m pl is the Planck mass. It is easy to check that Eq. (7) is equivalent to Eqs. (1) and (6) in fact. We note that the cosmological constant Λ given in Eq. (1) is related to the fine structure constant α according to Λ ∝ α −6 . This is interesting for us. As is well known, in the same year 1998 when the accelerated cosmic expansion was discovered, the evidence for cosmological evolution of the fine structure constant α has also been found [19,20]. Using the absorption systems in the spectra of distant quasars, Webb et al. [19] found the first evidence for the time variation of α, namely ∆α/α ≡ (α − α 0 )/α 0 = (−1.1 ± 0.4) × 10 −5 over the redshift range 0.5 < z < 1.6, where α 0 is the present value of α. Three years later, they improved the evidence to 4σ, namely ∆α/α = (−0.72 ± 0.18) × 10 −5 over the redshift range 0.5 < z < 3.5 [20]. The fine structure constant α was smaller in the past, and it is not a true constant in fact. Nowadays, a time-varying α has been extensively discussed in the community. There are many works on this topic in the literature [21][22][23]60]. If the cosmological constant Λ given in Eq. (1) is right, it should also be time-varying, because Λ ∝ α −6 . In the literature (e.g. [24,25]), there exist some Λ(t) models already. However, most of them are written by hand, e.g. Λ ∝ H 2 , Λ ∝ä/a, Λ ∝ R sc , Λ ∝ ρ m , where H is the Hubble parameter, a is the scale factor, R sc is the scalar curvature, ρ m is the density of matter. Different from the Λ(t) models purely written by hand, the time-varying Λ ∝ α −6 given in Eq. (1) is well motivated, as is shown above.
In the present work, we are interested to study the varying cosmological constant Λ ∝ α −6 . In Sec. II, we try to develop a suitable framework to model the varying cosmological constant. In Sec. III, we consider the observational constraints on the varying Λ models. In Sec. IV, the possible connection between the varying cosmological constant and warm dark matter is discussed. In Sec. V, some brief concluding remarks are given.

II. VARYING COSMOLOGICAL CONSTANT AND FINE STRUCTURE CONSTANT
Here, we try to develop a suitable framework to model the varying cosmological constant Λ given in Eq. (1). For convenience, we instead use the vacuum energy density ρ Λ given in Eq. (2), which is equivalent to Λ in fact. Throughout this work, we use the terms "cosmological constant" and "vacuum energy" interchangeably. If the cosmological constant is varying, we haveρ Λ = −Q = 0, where a dot denotes the derivative with respect to cosmic time t. To preserve the total energy conservation equatioṅ ρ tot + 3H(ρ tot + p tot ) = 0, a coupling between the vacuum energy and the pressureless matter is necessary, and henceρ m + 3Hρ m = Q = 0, where ρ m is the density of pressureless matter, ρ tot = ρ Λ + ρ m , and p tot is the total pressure. Note that the equation-of-state parameter (EoS) of the cosmological constant w Λ = −1, and the EoS of the pressureless matter w m = 0. Throughout this work, we assume that only the fine structure "constant" α is varying, and all the other fundamental constantsh, G, c, m e are true constants, i.e. they do not vary indeed. Since α = e 2 /(hc), this means that only the electron charge e is varying. Therefore, we have ρ Λ ∝ Λ ∝ α −6 . It is easy to see thaṫ and then the total energy conservation equation can be preserved according tȯ The coupling term Q = 6ρ Λα /α = 0 if the fine structure "constant" α is varying. In this work, we consider a spatially flat Friedmann-Robertson-Walker (FRW) universe containing only the vacuum energy and the pressureless matter. H ≡ȧ/a is the Hubble parameter, and a = (1 + z) −1 is the scale factor (we have set a 0 = 1; the subscript "0" indicates the present value of corresponding quantity; z is the redshift). In this way, we can turn the varying cosmological constant model into an interacting vacuum energy model. The vacuum energy interacts with the pressureless matter by exchanging energy between them. Due to the interaction Q, the evolutions of ρ m and ρ Λ should deviate from the ones without interaction, namely ρ m ∝ a −3 and ρ Λ = const. If the coupling term Q is given, one can derive the evolutions of ρ m and ρ Λ . However, the logic can be reversed. If the deviated evolutions of ρ m and/or ρ Λ are given, we can find the corresponding interaction Q from Eqs. (9), (10), and then the evolution of α. Inspired by e.g. [24,[26][27][28][29], we consider two different types of models to characterize the deviated evolutions of ρ m and/or ρ Λ in the following two subsections, respectively.

A. Type I models
Inspired by e.g. [27][28][29], the type I models are characterized by where f (a) can be any function of scale factor a. If f (a) ∝ a 3 , it corresponds to ΛCDM model whose ρ Λ = const. and ρ m ∝ a −3 . From Eq. (11) and Friedmann equation H 2 = 8πG(ρ Λ + ρ m )/3, we have where Ω i ≡ 8πGρ i /(3H 2 ) (i = Λ, m) are the fractional energy densities of the vacuum energy and matter, respectively. Substituting ρ Λ = ρ m f (a) into Eq. (9) and usingρ m from Eq. (10), we find where a prime denotes the derivative with respect to a. From Eqs. (12) and (13), we obtaiṅ where E ≡ H/H 0 . If f ∝ a 3 , it is easy to see thatα = 0, namely α = const. On the other hand, one can recast the total energy conservation equationρ tot + 3Hρ tot (1 + w tot ) = 0 as in which we have used w tot = p tot /ρ tot = Ω Λ w Λ + Ω m w m and w Λ = −1, w m = 0. We can integrate Eq. (15) to get where const. is an integration constant. Using Eqs. (12), (16) and H 2 = 8πGρ tot /3, we find where Ω Λ and E 2 are given in Eqs. (12) and (17), respectively. It is easy to check that if f ∝ a 3 , we obtain ∆α/α = 0, namely α = const. In summary, if the function f (a) is given, one can get the cosmic expansion history from Eq. (17), and the variation of the fine structure "constant" α from Eqs. (18) or (14). Finally, it is worth noting that by definition (11), we obtain which is useful to fix one of the model parameters in f (a).

A. Observational data
In the literature, there are two kinds of observational data concerning the variation of the fine structure "constant", namely the data of ∆α/α and the data ofα/α. To the best of our knowledge, most of the observational data are given in terms of ∆α/α, and only a few of the observational data are given in terms ofα/α. On the other hand, comparingα/α in Eqs. (14) and (26) with ∆α/α in Eqs. (18) and (27), it is easy to see that the number of free parameters inα/α is always more than the one in ∆α/α, namely one more free parameter H 0 is required inα/α while ∆α/α need not. Due to the above two reasons, we only consider the observational data of ∆α/α in the present work.
Here, we consider the observational ∆α/α dataset given in [30], which consists of 293 ∆α/α data from the absorption systems in the spectra of distant quasars. This sample includes 154 quasar absorption systems from the Very Large Telescope (VLT) in Chile, and 141 quasar absorption systems from the Keck Observatory in Hawaii. The full numerical data of these 295 quasar absorption systems are available in [31] or [32]. According to [30] and the instructions of [31, 32], two outliers (J194454+770552 at z abs = 2.8433, and J000448-415728 at z abs = 1.5419) should be removed. Therefore, there are 293 usable data in the final dataset, over the absorption redshift range 0.2223 ≤ z abs ≤ 4.1798. Note that all these 293 ∆α/α data are of O(10 −5 ). The χ 2 from these 293 ∆α/α data is given by where σ 2 i = σ 2 stat,i + σ 2 rand,i (see Sec. 3.5.3 of [30] and the instructions of [31, 32] for the technical details of σ rand and the error budget). In fact, we have tested our two types of models with these 293 ∆α/α data, and found that these ∆α/α data can tightly constrain the model parameters in f (a) or ǫ(a), but the constraints on the model parameter Ω m0 are too loose. Therefore, other cosmological observations, for instance, type Ia supernovae (SNIa), cosmic microwave background (CMB), and baryon acoustic oscillation (BAO), are required to properly constrain the model parameter Ω m0 .
We further consider the Union2.1 SNIa dataset [33] consisting of 580 data points, which are given in terms of the distance modulus µ obs (z i ). The theoretical distance modulus is defined by whereμ 0 ≡ 42.3841 − 5 log 10 h (h is the Hubble constant H 0 in units of 100 km/s/Mpc), and in which E ≡ H/H 0 , and p denotes the model parameters. The χ 2 from 580 Union2.1 SNIa is given by The parameterμ 0 (equivalent to H 0 ) is a nuisance parameter, but it is independent of the data points. One can perform a uniform marginalization overμ 0 . However, there is an alternative way. Following [34], the minimization can be made by expanding χ 2 µ in Eq. (31) with respect toμ 0 as .
Eq. (32) has a minimum forμ 0 =B/C at Since χ 2 µ, min =χ 2 µ, min (up to a constant), we can instead minimizeχ 2 µ which is independent ofμ 0 . Since using the full data of CMB and BAO to perform a global fitting consumes a large amount of computation time and power, we instead use the shift parameter R from the observation of CMB, and the distance parameter A from the measurement of BAO, which are model-independent and contain the main information of the observations of CMB and BAO [35], respectively. The shift parameter R of CMB is defined by [35,36] where the redshift of recombination z * is determined to be 1089.90 by the Planck 2015 data [37]. On the other hand, the Planck 2015 data have also determined the observed value of shift parameter R obs to be 1.7382 ± 0.0088 [38]. The χ 2 from CMB is given by The distance parameter A of the measurement of the BAO peak in the distribution of SDSS luminous red galaxies [39] is given by where z b = 0.35. In [39], the value of A has been determined to be 0.469 (n s /0.98) −0.35 ± 0.017. Here the scalar spectral index n s is taken to be 0.9741 by the Planck 2015 data [38]. The χ 2 from BAO is given by The total χ 2 from the combined ∆α/α, SNIa, CMB and BAO data is given by The best-fit model parameters are determined by minimizing the total χ 2 . As in [34,40], the 68.3% confidence level is determined by ∆χ 2 ≡ χ 2 − χ 2 min ≤ 1.

C. Observational constraints on type II models
Let us turn to the observational constraints on type II models introduced in Sec. II B. Obviously, the simplest type II model is given by where ǫ = 3 is a constant. If ǫ = 0, it reduces to ΛCDM model and α = const. (nb. Eqs. (26) and (27)). For ǫ(a) = ǫ = 3, we find the analytical formulas for η(a) and E 2 in Eqs. (24) and (25), namely There are two free model parameters, namely Ω m0 and ǫ. By minimizing the corresponding total χ 2 in Eq. (36), we find the best-fit model parameters Ω m0 = 0.279 and ǫ = 0.430 × 10 −6 , while χ 2 min = 870.391 and χ 2 min /dof = 0.997. In Fig. 3, we also present the corresponding 68.3% and 95.4% confidence level contours in the Ω m0 − ǫ plane. The parameter ǫ is tightly constrained to a narrow range of O(10 −6 ), thanks to the 293 ∆α/α data of O(10 −5 ). From Fig. 3, we see that ǫ = 0 (corresponding to ΛCDM model and α = const.) is fully consistent with the observational data (in fact it is close to the best fit). So, Λ and α can be non-varying in the type II model characterized by Eq. (40).
Next, we consider another type II model characterized by a CPL-like ǫ(a), namely where ǫ 0 and ǫ 1 are constants. If ǫ 0 = ǫ 1 = 0, it reduces to ΛCDM model and α = const. (nb. Eqs. (26) and (27)). As mentioned above, this CPL-like ǫ(a) = ǫ 0 + ǫ 1 (1 − a) can be regarded as the Taylor series expansion of ǫ(a) with respect to the scale factor a up to first order (linear expansion), and hence it is well motivated. There are three free model parameters, namely Ω m0 , ǫ 0 and ǫ 1 . Note that there are no analytical formulas for η(a) and E 2 in this case, but we can get η(a) by using numerical integration in Eq. (24) and then E 2 in Eq. (25) Fig. 4, we also present the corresponding 68.3% and 95.4% confidence level contours in the ǫ 0 − ǫ 1 , Ω m0 − ǫ 0 and Ω m0 − ǫ 1 planes. Both the parameters ǫ 0 and ǫ 1 are tightly constrained to the narrow ranges of O(10 −5 ), thanks to the 293 ∆α/α data of O(10 −5 ). From Fig. 4, we note that ǫ 0 = ǫ 1 = 0 (corresponding to ΛCDM model and α = const.) deviate from the best fit far beyond 2σ. This indicates that the varying Λ and α are favored by the observational data.

IV. VARYING COSMOLOGICAL CONSTANT AND WARM DARK MATTER
In the previous sections, we turned the varying cosmological constant model into an interacting vacuum energy model. The vacuum energy interacts with the pressureless matter by exchanging energy between them. In this section, we would like to view this model from another perspective. As is shown in e.g. [42], an interacting dark energy model can be equivalent to a warm dark matter model without interaction between dark energy and dark matter, while these two different kinds of models can share both the same cosmic expansion history and growth history. To keep things simple, here we only consider the models from the side of expansion history. Although the cold dark matter (CDM) model is very successful in many fields, it has been seriously challenged recently. We refer to e.g. [43] for the detailed reviews on these challenges. Recently, warm dark matter (WDM) remarkably rose as an alternative of CDM. We refer to e.g. [44] for several comprehensive reviews. The leading WDM candidate is the keV scale sterile neutrino. In fact, the keV scale WDM is an intermediate case between the eV scale hot dark matter (HDM) and the GeV scale CDM. Unlike CDM which is challenged on the small/galactic scale, it is claimed that WDM can successfully reproduce the astronomical observations over all the scales (from small/galactic to large/cosmological scales) [44]. One of the key differences between WDM and CDM is their EoS. WDM has a fairly small but non-zero EoS, while the EoS of CDM is exactly zero. In the literature, many attempts have been made to determine the EoS of dark matter (see e.g. [45][46][47][48][49]), and it is found that the EoS of WDM are of O(10 −6 ), O(10 −5 ) or O(10 −3 ) (depending on the working assumptions and the observational data in use).
Let us come back to the starting point Eqs. (9) and (10), and view them from another perspective. In the form of Eqs. (9) and (10), the vacuum energy (whose EoS is w Λ = −1) interacts with the cold dark matter (whose EoS is w m = 0) through an interaction Q = 6ρ Λα /α = 0. Now, we recast them aṡ where In this new form, the varying cosmological constant model becomes a model containing "dark energy" (whose EoS w de = −1) and "warm dark matter" (whose EoS w dm = 0), but there is no interaction between them. Since the observational data concerning the time variation of α are of O(10 −5 ), it is natural to expect that w dm ∼ 1 + w de ∼ O(10 −5 ) or smaller. For type I models introduced in Sec. II A, substituting Eq. (13) into Eq. (45), we have where Ω Λ and Ω m are given in Eq. (12). Thus, it is easy to find the evolutions of w dm and w de if f (a) is given. In the top panels of Fig. 5, we plot w dm and 1 + w de as functions of the scale factor a for the type I models with f (a) = f 0 a ξ in Eq. (37) and f (a) = f 0 a ξ0+ξ1(1−a) in Eq. (39), while the corresponding best-fit model parameters obtained in Sec. III B are taken. As expected above, they are of order 10 −6 or 10 −5 . Thus, the effective "warm dark matter" from type I models of the varying cosmological constant Λ ∝ α −6 is fully consistent with the observational constraints on WDM (e.g. [45][46][47][48][49]58]). For type II models introduced in Sec. II B, substituting Eqs. (21), (20) and (23) into Eq. (45), we get where η(a) is given in Eq. (24). Thus, it is easy to find the evolutions of w dm and w de if ǫ(a) is given.

V. CONCLUDING REMARKS
In this work, we considered the cosmological constant model Λ ∝ α −6 , which is well motivated from three independent approaches as mentioned in Sec. I. As is well known, in the passed 18 years, the hint of varying fine structure constant α was found, and the observational data of varying α were accumulated. Nowadays, a time-varying α has been extensively discussed in the community. If Λ ∝ α −6 is right, it means that the cosmological constant Λ should also be varying. In this work, we tried to develop a suitable framework to model this varying cosmological constant Λ ∝ α −6 , in which we view it from an interacting vacuum energy perspective. We proposed two types of models to describe the evolutions of Λ and α. Then, we considered the observational constraints on these models, by using the 293 ∆α/α data from the absorption systems in the spectra of distant quasars, and the data of SNIa, CMB, BAO. We found that the model parameters can be tightly constrained to the narrow ranges of O(10 −5 ) typically, thanks to the 293 ∆α/α observational data of O(10 −5 ). In particular, 3 of 4 models considered in this work favor the varying Λ and α, while ΛCDM model and α = const. deviate from the best fit beyond 2σ or at least 1σ. On the other hand, we can also view the varying cosmological constant model Λ ∝ α −6 from another perspective, namely it can be equivalent to a model containing "dark energy" (whose EoS w de = −1) and "warm dark matter" (whose EoS w dm = 0), but there is no interaction between them. We derived the effective EoS of "warm dark matter" and "dark energy", and found that they are fully consistent with the observational constraints on warm dark matter. In summary, we consider that the varying cosmological constant model Λ ∝ α −6 is viable and deserves further studies.
Some remarks are in order. First, although the cosmological constant Λ ∝ α −6 is derived from three independent approaches as mentioned in Sec. I (see also e.g. [59]), the underlying fundamental theory for it is still unknown. Varying fundamental constants require new physics [21]. It is commonly believed that the cosmological constant problem can only be solved ultimately in a unified theory of quantum gravity and the standard model of electroweak and strong interactions, which is still absent so far. Nevertheless, we consider that the studies on such a cosmological constant Λ ∝ α −6 might shed new light on the possible ways to the unknown underlying theory.
Second, we have tightly constrained the model parameters besides Ω m0 , namely ξ, ξ 0 , ξ 1 , ǫ, ǫ 0 and ǫ 1 , to the narrow ranges of O(10 −5 ) typically, mainly by using the 293 ∆α/α observational data from the absorption systems in the spectra of distant quasars. In fact, these parameters were confronted with only the observations of SNIa, CMB and BAO in e.g. [29], and the corresponding constraints are of O(1). So, the significant leap from O(1) to O(10 −5 ) shows the great power of the 293 ∆α/α observational data. In fact, these 293 ∆α/α data have been used in many issues (see e.g. [50,51]). We advocate the further uses of these ∆α/α observational data in relevant studies.
Third, in addition to the well-known evidence of the time variation in the fine structure constant α, it was claimed that α is also spatially varying [30,52]. If Λ ∝ α −6 is right, the cosmological constant should be not only time-dependent but also space-dependent. This might bring about new features to this field, and deserve further detailed studies. For example, it is claimed that there exists a preferred direction in the CMB temperature map (known as the "Axis of Evil" in the literature) [53], the distribution of SNIa or gamma-ray bursts [54][55][56], and the quasar optical polarization data [57]. If the cosmological constant Λ ∝ α −6 is also space-dependent, it might be responsible for the possible anisotropy in the (accelerated) expansion of the universe. We leave this interesting issue to future work.
Fourth, besides the 293 ∆α/α data from the absorption systems in the spectra of distant quasars, there are more ∆α/α observational data from other types of observations, for examples, atomic clocks, Oklo natural nuclear reactor, meteorite dating, CMB, big bang nucleosynthesis. We refer to e.g. [21] for a comprehensive review. Thus, it is interesting to consider the constraints from these observational data. Since these data are subtle in some sense [21,23], we also leave this to future work.
Fifth, let us turn to the varying cosmological constant model itself. In the present work, we considered four particular parameterizations of the functions f (a) and ǫ(a). In fact, one can instead consider other parameterizations, for instance, f (a), ξ(a) or ǫ(a) characterized by c 0 + c 1 ln a or c 0 a c1 [29]. On the other hand, in the present work we assumed that only the fine structure constant α is varying, and the other fundamental constants G, c,h, m e do not vary indeed. However, the varying G, c,h, m e models do exist in the literature (see e.g. [21] for a comprehensive review). Since the cosmological constant Λ given in Eqs. (1) or (2) also depends on G, c,h and m e , there are diverse variants of the varying cosmological constant model in fact. These variants might bring about new features. Since this is beyond the scope of the present work, we again leave it to future work.
Sixth, it is worth noting that the data analysis based on Eq. (36) can only give a rough indication and cannot be used to infer that any of the models considered in the present work is better than ΛCDM (we thank Prof. Dominik Schwarz for pointing out this issue), because (a) the SNIa light curve fitters assumed ΛCDM. So, if one wants to fit a model other than ΛCDM, the light curve fitting procedure should be redone within the new model. (b) the use of the parameters R and A from CMB and BAO measurements is crude, since they cannot take some complicated effects (e.g. the integrated Sachs-Wolfe effect) into account. (c) Eq. (36) implicitly assumes that ∆α/α data, SNIa, CMB and BAO have equal statistical weight, but this is questionable. In addition, a full Markov Chain Monte Carlo (MCMC) analysis of the CMB data should be used to further test the models considered in this work (we thank Prof. Dominik Schwarz for pointing out this issue), and we leave it to future work.
Finally, it is important to clarify the two different (but equivalent) perspectives on the varying cosmological constant model considered in this work. The first perspective is to regard the varying cosmological constant as a fluid not interacting with dark matter. In fact, this is the case considered in Sec. IV. The conservation equation of this fluid is given by Eq. (44), in the form ofρ de + 3Hρ de (1 + w de ) = 0. In this case, we stress that the EoS of this fluid is not w de = −1. As is clearly shown in e.g. Fig. 5, the EoS of this fluid is time-dependent, rather than constant. Indeed, the EoS of this fluid w de = −1 in this case. So, one cannot say ρ de = const. We have not assumed w de = −1 in this perspective indeed, and we refer to Sec. IV for detailed discussions. On the other hand, the second perspective on the varying cosmological constant model is considered in Secs. II and III. In this case, one can regard the varying cosmological constant as a fluid interacting with dark matter (different from the first perspective considered in Sec. IV). Now, the conservation equation of this fluid is given by Eq. (9), in the form ofρ Λ + 3Hρ Λ (1 + w Λ ) = −Q = 0. Yes, we assumed the EoS of this fluid w Λ = −1 in the second perspective considered in Secs. II-III, and thenρ Λ = −Q = 0. However, due to the non-zero Q, again one cannot say ρ Λ = const. In fact, the two perspectives considered in Sec. IV and Secs. II-III are completely independent. There is no inconsistency in each independent perspective. One should not mix up these two different perspectives considered in this work, otherwise confusion and misunderstanding might arise.