Running vacuum in the Universe and the time variation of the fundamental constants of Nature

We compute the time variation of the fundamental constants (such as the ratio of the proton mass to the electron mass, the strong coupling constant, the fine structure constant and Newton's constant) within the context of the so-called running vacuum models (RVM's) of the cosmic evolution. Recently, compelling evidence has been provided showing that these models are able to fit the main cosmological data (SNIa+BAO+H(z)+LSS+BBN+CMB) significantly better than the concordance $\Lambda$CDM model. Specifically, the vacuum parameters of the RVM (i.e. those responsible for the dynamics of the vacuum energy) prove to be nonzero at a confidence level $\gtrsim3\sigma$. Here we use such remarkable status of the RVM's to make definite predictions on the cosmic time variation of the fundamental constants. It turns out that the predicted variations are close to the present observational limits. Furthermore, we find that the time variation of the dark matter particles should be crucially involved in the total mass variation of our Universe. A positive measurement of this kind of effects could be interpreted as strong support to the"micro and macro connection"(viz. the dynamical feedback between the evolution of the cosmological parameters and the time variation of the fundamental constants of the microscopic world), previously proposed by two of us (HF and JS).


Introduction
The possibility that the so-called "constants" of Nature (such as the particle masses and the couplings associated to their interactions) are actually not constants, but time evolving quantities following the slow pace of the current cosmological evolution, has been investigated in the literature since long time ago [1]. The history of these investigations traces back to early ideas in the thirties on the possibility of a time evolving gravitational constant G by Milne [2] and the suggestion by Dirac of the large number hypothesis [3], which led him also to propose in 1937 the time evolution of G. The same year Jordan speculated that the fine structure constant α em together with G could be both space and time dependent [4] -see also [5]. This, more field-theoretically oriented, point of view was retaken later on by Fierz [6] and, finally, in the sixties, G was formally associated to the existence of a dynamical scalar field φ ∼ 1/G coupled to the curvature. Such was the famous framework originally proposed by Brans and Dicke ("BD" for short) [7,8], which was the first historical attempt to extend General Relativity (GR) to accommodate variations in the Newtonian coupling G. To avoid conflict with the weak equivalence principle, the variability of φ in the BD-theory was originally restricted to a possible time evolution only, i.e. φ(t) with no spatial dependence [9]. A generalization of the BD approach in the seventies [10] subsequently led to a wide variety of scalar-tensor theories, which are still profusely discussed in the current literature. Furthermore, the subject of the time variation of the fine structure constant (cf. [11] for other early proposals) is particularly prolific. Numerous studies have been undertaken in our days from different perspectives, sometimes pointing to positive observational evidence [12] but often disputed by alternative observations [13] -see e.g. the reviews [14,15,16,17]. The possibility that the particle masses could also drift with the cosmic evolution has also become a favorite target of the modern astrophysical observations. Thus e.g. the dimensionless proton-to-electron mass ratio, µ ≡ m p /m e , has been carefully monitored using quasar absorption lines. Claims on significant time variationμ/µ at ∼ 4σ c.l. are available in the literature [18], although still unconfirmed by other observations [19]. Future high precision quantum optic experiments in the lab are also planned to test the possible time variation of these observables, and they will most likely be competitive [1]. Clearly, the time and space variation of the fundamental constants is a very active field of theoretical and experimental research that may eventually provide interesting surprises in the near future [20].
From a more modern perspective the dark energy (DE) theories of the cosmological evolution also predict the cosmic time variation of the fundamental constants. These include the stringdilaton models, Kaluza-Klein theories, chameleon models etc., in which the underlying dynamical fields are either massless or nearly massless -cf. e.g. [21] for some early attempts and [16] for a review of more recent proposals and the observational situation. The possibility that dark matter theories could also impinge on the time variation of fundamental constants has also been put forward [22]. Many other proposals are available in the literature, see e.g. [23] for more contemporary developments.
In this paper, we would like to focus our attention on a specific class of DE models leading to a time variation of the fundamental constants [24]. These models are particularly interesting since it has recently been shown that they prove capable of fitting the main cosmological data in a fully competitive way with the concordance ΛCDM model [25]. These are the so-called running vacuum models (RVM's) of the cosmological evolution, see e.g. [26] for a recent summarized discussion and [27,28,29] for a more extensive review and a comprehensive list of references. In these models there are no ultralight scalar fields and the time variation of the fundamental constants is effectively triggered via quantum effects induced by the cosmological renormalization group, whose flow is naturally set up by the expansion rate H, cf. [27] and references therein. The framework is characterized by a dynamical vacuum energy density, ρ Λ , which is a power series of H and its time derivatives. For the current Universe, it suffices to consider ρ Λ = ρ Λ (H,Ḣ) up to linear terms inḢ and quadratic in H [28,29,30,31,32]. However, extensions with higher powers have also been considered for describing inflation, see [28,33] and [34,35,36] for different kind of scenarios, including anomaly-induced inflation [37]. Because of the dynamical nature of the vacuum in these models, a natural feedback occurs between the cosmological parameters and the fundamental constants of the microscopic world, such as the particle masses and coupling constants. Remarkably, recent studies have shown that the class of RVM's can provide an excellent fit to the main cosmological data (SNIa+BAO+H(z)+LSS+BBN+CMB) which is highly competitive with that of the ΛCDM -see most particularly [38,39,40,41], and previous studies such as [30,31,42,43] and references therein. Therefore, there is plenty of motivation for further investigating these running vacuum models. In this paper, building upon the aforementioned works which single out the especial status of the RVM's, we wish to estimate the possible variation of the fundamental constants triggered by the dynamical interplay between the evolution of the vacuum energy density ρ Λ = ρ Λ (H) and the concomitant anomalous conservation law of matter (which may lead to time dependent particle masses) and/or the time evolution of the gravitational coupling G(H). Such feedback between the ultralarge scales of cosmology and the minute scales of the subatomic physics is what two of us have called elsewhere "the micro and macro connection" [44,24]. Recently it has been shown that it could also lead to a possible explanation for the origin of the Higgs potential in the gravitational framework [45].
This paper is organized as follows. After a preliminary historical discussion in the introduction, in sect. 2 we remind of the possibility of the cosmic time variation of the cosmological parameters in the context of GR without committing to any particular model. In sect. 3 we focus on the running vacuum models (RVM's) as a particularly appealing framework where to realize the time evolution of the fundamental constants. In sect. 4 we consider in detail the specific prediction of the RVM's concerning the time variation of the particle masses and couplings. In sect. 5 we briefly discuss alternative dynamical vacuum models. Finally, in sect. 6 we present our conclusions.

Cosmological models with time evolving parameters
We wish to explore the possibility that the fundamental 'constants' or parameters P of the subatomic world (such as masses and couplings) are actually slowly varying with (cosmic) time t and whether this feature might be related to the cosmic evolution of the fundamental constants of gravitation. If so the cosmic time evolution of P should be typically proportional to the rate of change of the scale factor of the cosmic expansion, i.e.Ṗ/P ∝ȧ/a ≡ H (the Hubble rate). From this point of view we should expect (at least in linear approximation) that the fractional cosmic drift rate of P(t) near our time is proportional to H 0 = 1.0227 h × 10 −10 yr −1 , with h ≃ 0.67, and hence very small in a natural way. SpecificallyṖ/P ∼ H 0 ∆P/P (∆P/P) 10 −10 yr −1 . From the existing bounds on the known particles, typically we expect that ∆P/P should vary between a few parts per million (ppm) to at most one part in a hundred thousand over a cosmological span of time, depending on the specific parameter (such as couplings, masses etc: P = G, m i , α em , α s , Λ QCD ...) However, the attributes of the dark matter (DM) particles (e.g. their masses and couplings) could vary faster. The possible cosmic time evolution of these parameters can equivalently be described as a redshift dependence:Ṗ where z = (1 − a)/a is the cosmic redshift, a(t) is the scale factor (normalized to a(t 0 ) = 1 at present) and P ′ (z) = dP/dz . It will prove useful to present most of our results in terms of the cosmological redshift.
The above micro-and-macro-connection scenario [24,44] can be implemented from Einstein's field equations in the presence of a time-evolving cosmological constant, Λ(t). In fact, without violating the Cosmological Principle in the context of the Friedmann-Lemaître-Robertson-Walker (FLRW) universe, nothing prevents Λ(t) and/or G = G(t) from being functions of the cosmic time 2 . The field equations can be written G µν = 8πGT µν , where G µν = R µν − 1/2g µν R and T µν = T µν + g µν ρ Λ stand respectively for the Einstein tensor and the full energy-momentum tensor, in which T µν is the ordinary part (involving only the matter fields) and g µν ρ Λ carries the time-evolving vacuum energy density ρ Λ (t) = Λ(t)/8πG(t).
The field equations for a variable gravitational "constant" and a dynamical vacuum energy density in the FLRW metric in flat space are derived in the standard way and are formally similar to the case with G and Λ constants. The corresponding generalization of the Friedmann equation reads where ρ m = ρ B + ρ DM represents the contribution of non-relativistic matter (baryons + dark matter) and ρ Λ (t) stands for the aforementioned dynamical vacuum energy density. An important consequence of the covariance of GR is the Bianchi identity involving the Einstein tensor: ∇ µ G µν = 0. Via the field equations, it implies ∇ µ (GT µν ) = 0. Using the FLRW metric and considering that the matter of the current Universe is made of pressureless dust, we find the explicit form of the Bianchi identity is:Ġ(ρ m + ρ Λ ) + G(ρ m +ρ Λ ) + 3GHρ m = 0. This equation plays the role of a generalized local conservation law. Using Eq. (1) it can be conveniently rewritten in terms of the redshift variable as follows: The above expression relates, in a fully general form, the evolution of the matter energy density and the vacuum energy density in the presence of a dynamical G. For G =const. and ρ Λ =const. (i.e. within the context of the ΛCDM) Eq. (3) integrates trivially and renders the canonical conservation law of non-relativistic matter: where ρ 0 m is the current matter energy density. Equations (2) and (3) cannot be solved beyond the ΛCDM unless some model or an ansatz is provided e.g. on the evolution of the vacuum energy density ρ Λ . In the next section we adopt the proposal for ρ Λ associated to the running vacuum models (RVM's). In such case the solution for the matter energy density will no longer be in general of the canonical form (4) and this, as we shall see, can be interpreted either as an anomalous matter conservation law or as a time evolution of the particle masses. Together with the cosmic evolution of G, we can see that this scenario may well lead us to a concrete realization of the time variation of the fundamental constants that can be consistent with GR.

Running vacuum in the expanding Universe
Even though we know that the standard matter conservation law (4) must be essentially correct, we wish to explore the possibility that small corrections could perhaps be accommodated. Let us first proceed phenomenologically and later on we will adduce some theoretical motivations supporting the obtained results. Suppose that the standard dilution law for matter in an expanding universe, ρ m ∝ a −3 , receives a small correction: ρ m ∝ a −3(1−νm) , with |ν m | ≪ 1 a dimensionless parameter. In terms of the redshift, we have where for ν m = 0 we retrieve of course the standard form (4). This possibility was first proposed and tested in the literature in Ref. [48], and later on it was addressed by other authors -see e.g. [49,50]. However, here we wish to reinterpret the anomalous conservation law (5) in a different way [24]. Rather than assuming that the presence of a nonvanishing ν m is related to an anomalous nonconservation of the number density of particle masses, we assume that it parameterizes the nonconservation of the particle masses themselves. In other words, we suppose that while there is a normal dilution of the particle number density with the expansion for all the particle species, i.e. n i = n 0 i a −3 = n 0 i (1 + z) 3 , the corresponding mass values m i are not conserved throughout the cosmic expansion: where m 0 i are the current values (z = 0) of the particle masses and the various ν i are the corresponding anomaly indices for the non-conservation of each species. This was the point of view adopted in [24]. Notice that ρ m i in Eq. (5), for the ith-species, is indeed equal to n i (z)m i (z), with m i (z) given by (6) and ρ 0 m i = n 0 i m 0 i . In the remaining of this section it will not be necessary to distinguish among the different values of ν i for each particle species, but we shall return to this point in the next section.

Running G and running ρ Λ
Let us assume that G is also a slowly varying cosmic variable. More specifically, we shall suppose that it varies logarithmically with the Hubble function as follows: where ν G is another small dimensionless parameter (|ν G | ≪ 1), different from ν m in general. A theoretical motivation for this expression within QFT in curved spacetime can be found in [37] see also [27]. There is also a practical reason for this form, as we will see in a moment. The above assumptions stand us in good stead to find out what is the corresponding dynamical evolution law for the vacuum energy density within GR. As indicated, we assume that the parameters are small since the ensuing variation of the fundamental constants must be small and the model cannot depart significantly from the standard cosmology. Having this in mind and substituting the expressions (5) and (7) in the generalized conservation law (3), as well as using the Friedmann equation (2), we can analytically integrate the evolution of ρ Λ as a function of the redshift. The final result reads as follows: Here ρ Λ (z = 0) = ρ 0 Λ is the current value of the vacuum energy density, and of course this value would stay constant for the entire cosmic history (as in the ΛCDM) if both ν m and ν G were zero. Therefore ν m + ν G can be interpreted as an effective coefficient for the running of ρ Λ . The obtained result (8) shows that the dynamical character of the vacuum energy can indeed be in interplay with both the dynamics of the gravitational coupling and the non-conservation of matter, in a manner perfectly compatible with the GR field equations. Amusingly, we note that if ν m = −ν G the vacuum energy density would remain constant as in the ΛCDM, although the model would not quite behave as the standard cosmological model because the anomalous matter conservation law (5) and the logarithmically evolving Newtonian coupling (7) both stay in force. In any case the departure from the standard cosmology in all cases will be small enough provided |ν m | ≪ 1 and |ν G | ≪ 1.
Defining the normalized Hubble function with respect to the current value, E = H/H 0 , explicit computation from Friedmann's equation, using (5) and (8), provides the following result: where G is given in our case by (7). Notice that for ν m = 0 (matter conservation) and ν G = 0 (G =const.) the Hubble function (9) boils down to the ΛCDM, as should be expected.
We may now understand why the forms (5) and (7), leading to (8), are particularly interesting. The obtained result for the vacuum energy density can be motivated from the integration of the renormalization group (RG) flow associated to the Hubble expansion, within the context of the running vacuum model (RVM) -see [27,28,30,31] and references therein. Let us note that the leading RG effects up to order ∼ H 4 take on the form of the power series [28] where the coefficients a i have different dimensionalities in natural units. Specifically, a 0 has dimension 4 since this is the dimension of ρ Λ ; a 1 and a 2 have both dimension 2; and, finally, a 3 , a 4 and a 5 are dimensionless. Notice that only the terms with an even number of derivatives of the scale factor are present in the above expression since the vacuum energy density is part of the effective action of QFT in curved spacetime and therefore must preserve general covariance. The O(H 4 ) terms can be important for the inflationary stage [28,33,34,35,36], but for the post-inflationary epoch and in particular for the current Universe they can be neglected. We are thus left with the O(H 2 ) terms only, namely H 2 andḢ. For simplicity we will focus on the former since the inclusion of the latter will not affect the main discussion in this study concerning the variation of the fundamental constants 3 . Therefore we concentrate here on the simplest form of the RVM density, which we can rewrite after an appropriate redefinition of a 0 and a 2 as follows: where we have arranged that for H = H 0 we recover the current value of the vacuum energy density: ρ Λ (H 0 ) = ρ 0 Λ . Here M P = 1/ √ G 0 is the Planck mass and ν is a small dimensionless coefficient which characterizes the dynamical evolution of ρ Λ (H); in fact ν plays the role of the β-function coefficients of the renormalization group equation (RGE) for the running vacuum 4 . In particular, for ν = 0 we recover the ΛCDM case. Notice that M 2 P H 2 is of order ρ 0 Λ and hence for ν = 0 the term ∼ ν M 2 P H 2 represents a small correction (for |ν| ≪ 1) to the constant value ρ 0 Λ and endows ρ Λ of a mild dynamical behavior which can be favorable to observations. Indeed, by fitting the parameter ν to the overall cosmological data one finds an improved fit as compared to the ΛCDM, provided |ν| ∼ 10 −3 -for details, see [38,39,30,31]. The connection between (11) and (8) can now be elucidated as follows [28,30]. To start with, take G =const. as this simplifies the structure of the generalized conservation law (3). Inserting (5) in that law we can solve for ρ Λ (z) and we find the expression (8) with ν G = 0. Knowing the matter and vacuum densities, Friedmann's equation immediately provides H as a function of the redshift and we arrive at Eq. (9) (with G = G 0 ). Combining this expression for H(z) with that of ρ Λ (z) we find (11).
One can also relate the running vacuum law (11) with a running gravitational coupling of the form (7) as follows. Assume now that the standard local matter conservation law (4) holds good, i.e. ν m = 0 in (5). In this way Eq. (3) boils down to (dG/G)(ρ m + ρ Λ ) + dρ Λ = 0. The latter can be solved using (11) with ν = ν G and noting also that ρ m + ρ Λ = 3H 2 /(8πG) by virtue of Friedmann's equation (2). The result is a simple differential equation for G(H): Integrating it with the boundary condition G(H 0 ) = G 0 = 1/M 2 P , the final result is indeed Eq. (7). In other words, if matter is conserved the running vacuum (11) leads to a logarithmic evolution of the gravitational coupling of the form (7). These two evolution laws take therefore the precise dynamical forms needed to fulfill the Bianchi identity of the field equations when matter is conserved.
In the previous considerations we have assumed either ν G = 0 or ν m = 0. But we can also determine the running of ρ Λ as a function of H when ν m and ν G are both nonvanishing. From (8) and (9) we find: (7), we find that Eq. (13) boils down to Eq. (11) at first order in the small parameters (ν m , ν G ). Recall that in all cases these parameters are assumed to satisfy |ν m , ν G | ≪ 1. It follows that Eq. (11) does correctly describe the running of the vacuum energy density in terms of H in all possible cases within this framework, i.e. even when matter is non-conserved and at the same time there is a running of the gravitational coupling. In other words, the running of ρ Λ as a function of H in (11) is controlled in all cases by ν = ν m + ν G , irrespective of whether one or none of these parameters is zero.

Free parameters
Let us summarize the situation with the free parameters in the class of RVM's considered here. The basic parameters are (ν m , ν G ), which are associated to the anomalous matter conservation law (5) and the time evolution of the gravitational coupling, Eq. (7). Given these two parameters the Bianchi identity (3) determines the evolution of the vacuum energy density, Eq. (8), through ν = ν m + ν G . However, as we have mentioned, there is no reason to expect that all of the particle masses should have the same anomaly index, and hence we expect that the anomaly mass index ν m can be expressed in terms of the different particle indices ν i defined in Eq. (6). For example, baryons can have the index ν B (assumed the same for all of them), but it could be different from the index for DM particles, X, which we call ν X . The relation between the overall anomaly index ν m and the specific indices (ν B , ν X ) will be discussed in sect. 4.
Furthermore, as indicated, if it turns that a 1 = 0 in Eq. (10), then an additional parameter is still possible for the running of the vacuum energy in the present universe, but we shall not take it into account since it is not necessary to illustrate the possible existence of the basic effects under study. In actual fact, in all phenomenological considerations we will assume a 1 = 0 together with one of the following two possibilities: either i) ν m = 0 with ν G = 0, or ii) ν G = 0 with ν m = 0. This will suffice to parametrize the time variation of the fundamental constants that we are considering here. In such context the Bianchi identity enforces the value of ν (the parameter that controls the running of the vacuum energy density in (11)) to be either ν m or ν G , depending on wether we assume either that G is fixed and the matter has some anomaly conservation law, or that the matter is strictly conserved and G has some evolution, but not both situations at the same time. While ν could perfectly receive simultaneous contributions from both kinds of effects -in the above mentioned form ν = ν m + ν G -at the moment it is not possible to individually disentangle them phenomenologically. Thus, in our numerical evaluations we will always assume the separate situations in which either ν = ν m or ν = ν G .
From the foregoing considerations we see that the RVM's offer several possibilities for the time variation of the fundamental "constants", all of them being connected through the evolution laws (5), (7), (8) and (11), in which there is a built-in principle for exchanging energy between matter, vacuum and the gravitational coupling in different combinations that are compatible with GR.

Time evolution of the fundamental constants in the RVM
In this section we wish to evaluate the specific impact of the running vacuum models (RVM's) on the time evolution of the fundamental constants of the standard model (SM) of particle physics and the fundamental constants of cosmology. Basically we will assess the predicted variation of the particle masses, most conspicuously the proton mass (through the proton-to-electron mass ratio µ ≡ m p /m e ), the QCD scale Λ QCD , the QED fine structure constant α em , the corresponding QCD coupling α s , the gravitational constant G and the cosmological constant Λ. We will also consider the possible implications from Grand Unified Theories (GUT's).

Time variation of masses and couplings in the SM
The framework outlined in sect. 3 suggests that basic quantities of the standard model (SM) of strong and electroweak interactions, such as the quark masses, the proton mass and the quantum chromodynamics (QCD) scale parameter, Λ QCD , might not be conserved in the course of the cosmological evolution [24,44]. Let us take, for example, the proton mass, which is given as follows: where m u,d.s are the quark masses and the last term represents the electromagnetic (em) contribution. Obviously the leading term is the first one, which is due to the strong binding energy of QCD. Thus, the nucleon mass can be expressed to within very good approximation as m p ≃ c QCD Λ QCD ≃ 938M eV , in which c QCD is a non-perturbative coefficient. The masses of the light quarks m u , m d and m s also contribute to the proton mass, although by less than 10% and can therefore be neglected for this purpose. It follows that time (or cosmic redshift) variations of the proton mass are essentially equivalent to time (redshift) variations of the QCD scale parameter: where in the last equality we have used Eq (1). On the other hand, the QCD scale parameter is related to the strong coupling constant α s = g 2 s /4π as follows (at 1-loop order): µ R being the renormalization point, n f the number of quark flavors and Λ QCD = 217 ± 25 MeV the measured value of the QCD scale paramete. However, if there is a crosstalk between the micro and macro world as suggested in the previous section, we expect that when we embed QCD in the context of a FLRW expanding universe the value of the proton mass, and hence of Λ QCD , need not remain constant anymore. The possible change of Λ QCD should be of course relatively small, and from the above considerations we envisage that its value may evolve with the rate of expansion of the universe, i.e. Λ QCD = Λ QCD (H). In such case the strong coupling constant α s becomes a function not only of the conventional renormalization scale µ R but also of the cosmic scale µ c ≡ H.
Since H = H(z) is a function of the cosmological redshift, we can write α s = α s (µ R , z). From Eq. (15) we find that the relative variations of the two QCD quantities with the Hubble rate are related (at one-loop) in the following manner: If the QCD coupling constant α s or the QCD scale parameter Λ QCD undergo a small cosmological shift, the nucleon masses and the masses of the atomic nuclei would also change along with Λ QCD . The cosmic dependence of the strong coupling α s (µ R , z) can be generalized to the electroweak couplings of the SM, including the fine structure constant in QED, α em , except that there is no electroweak equivalent for the Λ QCD scale. As a consequence there is no obvious connection of the variation of the particle masses with the variation of the electroweak coupling. This new sort of time variation would be possible only if the vacuum expectation value of the Higgs would be itself time-dependent, but we shall not address this option here. There is however an alternative link based on the framework of Grand Unified Theories, which will be explored in sect. 4.3.
In this paper we attribute the cosmic variation of the particle masses to the energy exchange with the cosmic vacuum according to the RVM framework outlined in the previous section, in which a possible additional variation of the gravitational constant may also concur. In order to estimate quantitatively these effects within the RVM, we take as a basis the numerical fit estimates obtained in [38,39,30,31] using the known data on SNIa+BAO+H(z)+LSS+BBN+CMB -see also [29] for a review. Among these observational sources (which involve several hundreds of data points indicated in these references) we use 36 data points on the Hubble rate H(z) at different reshifts in the range 0 < z ≤ 2.36, as compiled in [51,52], out of which 26 data points are inferred from the differential age method, whereas 10 correspond to measures obtained from the baryonic acoustic oscillation method (cf. Fig. 1). These data will play a significant role in our aim to constrain cosmological parameters because they are obtained from model-independent direct observations. In particular, we use this compilation for investigating a possible temporal evolution of the particle masses, both for baryons and dark matter. In Fig. 1 we plot the above mentioned data points H(z i ) (i = 1, 2, ..., 36) and at the same time we superimpose the Hubble functions H(z) for both the ΛCDM model (dashed line) and the RVM (solid one). The difference between the two is small, of course, because the RVM Hubble function (9) differs only mildly from the standard one owing to the parameters ν m , ν G being small in absolute value. However, the small differences are perfectly visible in Fig. 1 and are sufficient to improve significantly the overall fit to the SNIa+BAO+H(z)+LSS+BBN+CMB data points. According to [30], the data show a preference (at the level of 3σ) for a dynamical vacuum of the form (11) rather than the rigid vacuum (ν = 0) of the ΛCDM case, for which ρ Λ =const. The obtained fit values of the vacuum parameter ν stay in the ballpark of 10 −3 [38,39,29], and therefore we can use this order of magnitude determination as a characteristic input in our estimate of the variation of the fundamental constants. The predicted mass drift rates are indicated in Figs. 2-3, and in what follows we explain their origin.
To estimate the variation of the particle masses within the RVM, let us start by noting that the matter density of the universe can be approximated as: ρ m ≃ n p m p + n n m n + n X m X , where we neglect the leptonic contribution and the relativistic component (photons and neutrinos) [24]. Here n p , n n , n X (m p , m n , m X ) are the number densities (and corresponding masses) of protons, neutrons and dark matter (DM) particles X, respectively. Assuming that the mass non-conservation law in eq. (5) is to be attributed to the change of the mass of the particles -confer. Eq. (6) -the relative total time variation of the mass density associated to such mass anomaly can be estimated Here ν X = 10 −3 , which essentially saturates the fitted value of ν m (recall that ν B ≪ ν m , see the text).
as follows: where δρ m is obtained by differentiating ρ m with respect to time and subtracting the ordinary (i.e. fixed mass) time dilution of the number densities. In the above equation, Ω B and Ω DM represent as usual the fractional density of baryons and DM particles with respect to the critical density, respectively. Of course the total Ω m is the sum Ω B + Ω DM and in the numerical analysis we take as a current value Ω m = 0.30, according to the global fits obtained from the RVM models [38,39,40]. The approximation made on the r.h.s. of (17) uses the fact that ρ B /ρ DM = Ω B /Ω DM = O(10 −1 ), with ρ B = n p m p + n n m n and ρ DM = n X m X ≡ ρ X the density of DM particles. Equation (17) can be further expanded as follows. Let us take m n = m p ≡ m B so that ρ B = (n p + n n ) m B , and assumeṁ n =ṁ p ≡ṁ B . Since n n /n p is of order 10% after the primordial nucleosynthesis and Ω B /Ω DM is also of order 10%, we may neglect the product of these two terms or any higher power of them. In this way we are led to Note that the prefactor in the first term of the r.h.s. can be written Inserting (19) in (18) we see that the leading power of n p /n n appears at second order and hence can be neglected. The final result therefore reads: where we have introduced the anomaly indices ν B and ν X for the evolution of the baryon and DM masses. They define the corresponding mass drift rates for baryons and DM particles: The total drift rate from the time variation of the masses of all heavy and stable particles in the Universe (baryons + dark matter) reads The drift rates are of course functions of time and redshift. The corresponding relation with the variation of a particular mass m i , baryon or DM, within a cosmological span of time ∆t ∼ H −1 can be a complicated function of time, but it is usually approximated in a linear way, i.e. one assumes that on average the time variation of the mass was the same during the time interval ∆t. In this way we can writeṁ Thus, the anomaly indices ν i encode the typical mass variation of a given particle species (baryons or DM particles) in a cosmological span of time. Since, as mentioned, Ω B /Ω DM ≃ 0.1, we can neglect the square of this quantity and rewrite (21) in the more compact form This is the promised relation between the anomaly mass index in Eq. (5) and the specific baryonic and DM anomaly indices. Let us note that the presence of these anomaly indices for matter non-conservation can affect different aspects of the cosmic history, such as the precise moment of matter-radiation equality or the details of the growth of structure formation. These effects have been evaluated in different works, see e.g. [53,40,31], and can be important in the future when more precise observational data will be available. For the present work, it suffices to estimate the order of magnitude of the anomaly indices in (25) in order to transfer the possible impact on the time variation of the fundamental constants. For this reason we used here only the order of magnitude values of the fitting results obtained in [38,39,40]. As far as the lepton index ν L is concerned, it cannot have any significant effect in the above equation (25) since it is suppressed by the small lepton mass rate in the Universe as compared to baryons. See, however, Eq. (26) below for other effects from ν L that might be not so negligible. Equation (25) can actually be checked experimentally, for ν m is related to the running of the vacuum energy density (e.g. ν m = ν when ν G = 0, as explained in sect. 3) and ν can be fitted from cosmological observations based on SNIa+BAO+H(z)+LSS+BBN+CMB data [38,39,29], and it is found to be of order 10 −4 − 10 −3 , whereas ν B can be determined from astrophysical and lab experiments usually aimed at determining the time evolution of the ratio µ = m p /m e [16,23]. Thus, if the equation (21) -or equivalently (25) -must be fulfilled, we can indeed check if the DM part ν X (which, of course, cannot be measured individually) plays a significant role in it.
Observationally one finds from a rich variety of experimental situations both from astrophysical observations and direct lab measurements [16,17,54] (most of them compatible with a null result) that ∆µ/µ is at most in the ballpark of O(1 − 10) parts per million (ppm). Let us note that where we have used Eq. (24). The index ν B was applied to the proton as the only stable baryon, whereas ν L corresponds to the electron as the only stable lepton. It is usually assumed that ν B ≫ ν L and then ∆µ/µ ≃ ∆m p /m p . In this case the aforementioned limit on ∆µ/µ would imply ν B ∼ 10 −5 at most. However, a more symmetric option (which cannot be ruled out at present) is that the two indices ν B and ν L can be close to each other. In such case both could be of order 10 −4 and very similar; this case would still be compatible with the approximate bounds on ∆µ/µ of at most 10 ppm. We will keep in mind these two possibilities in our analysis 5 . Both of them, however, lead to ν X ∼ ν m ∼ 10 −3 via Eq. (25), what clearly points to the crucial role of the DM contribution to explain the bulk of the mass drift rate in the Universe (cf. Fig. 4).
Assuming that the anomaly indices for matter non-conservation are constant, we can integrate Eq. (22) and we find the evolution of the baryons and DM particle masses. We may most conveniently perform the integration in terms of the redshift using Eq. (1), and we find: Here we have defined ∆m i (z) = m i (z) − m 0 , with m 0 ≡ m(z = 0); and the index i = B, L, X runs over stable baryons, leptons and DM particles, respectively. The form (27) obtained for the mass variation of the particles with the redshift is indeed the one announced in (6). The total mass variation of the Universe is conceived here as a physical process connected with the variation of the particle masses themselves rather than the appearance or disappearance of new particles during the expansion process. In our framework the particle mass changes are possible thanks to the interaction with the dynamical vacuum and/or the evolution of the gravitational constant, as described in the previous section. In Fig. 4 we plot their evolution with the redshift, showing that the DM mass drift is the dominant one.

Time variation of G and Λ
The corresponding cosmic drift rates of the vacuum energy density and gravitational coupling ensue from (7) and (8). Using Eq. (1) we find, in leading order: where the Hubble function for the RVM is given by Eq. (9). It is convenient to trade the derivative of the Hubble function in the equation above in terms of the Hubble function itself. After some rearrangement we find: which is valid for any value of the redshift. While for small values z < 1 the previous expression behaves roughly asĠ G ≃ 3ν G Ω m (1 + z) 3 H 0 , for large values of the redshift we haveĠ G ≃ 3ν G H. The corresponding plots of G(z)/G 0 and of the cosmic drift rateĠ/G as a function of the redshift are depicted in Fig. 5. We can see from the plot on the left in that figure that the value of G decreases with the redshift and therefore G behaves as an asymptotically free coupling, that is to say G decreases in the past, which is the epoch where the Hubble rate (with natural dimension of energy) is higher. This is of course already obvious from Eq. (7) for ν G > 0. Moreover from the plot on the right in Fig. 5 we learn that the rhythm of variation of G slows down with the cosmic expansion, i.e. the rate of change is larger in the past. Following the fitting results to the overall cosmological observables obtained in [38,39,30] we can explore the parameter ν G in the typical range from 0.0005 to 0.001 and then evaluate the relative variation ∆G/G as a function of the redshift in that interval. The plot is depicted in Fig.  6. The larger is ν the smaller is G at any given redshift. The upper curve in Fig. 6 corresponds to ν G = 0.0005 whereas the lower one to ν G = 0.001. The shaded area gives the prediction of ∆G/G for the values of ν G comprised in that interval for each z. Worth noticing is that if we extend the domain of applicability of Eq. (7) up to the BBN epoch (z = z N ∼ 10 9 ), we find that the relative variation of G at the BBN time as compared to the present is (G(z N ) − G 0 )/G 0 ≃ −0.06, which is less than 10% in absolute value and hence compatible with the current BBN bounds [16].
Finally, in Fig. 7 we display the evolution of the vacuum energy density with the redshift according to the RVM formula (8), specifically the solid line corresponds to the relative correction with respect to the current value: where ν = ν m + ν G . We may use Eq. (11) and rewrite the evolution of the vacuum energy density as This formula displays the relative variation of the vacuum energy density with respect to the current value directly in terms of the Hubble function. In the RVM this expression is more fundamental than (31) because it is the solution of the RGE in terms of H. As we shall see in Sect. 5 , there are alternatives models that may lead to Eq. (31) using ad hoc assumptions on the interaction between DE and DM. By the same token the expression of G directly in terms of H is closer to the spirit of the RVM since it can be derived from the RG formalism of QFT in curved spacetime [27,37].

Time evolution of the fine-structure constant
Motivated by a possible indication of a variation (decrease) in the fine-structure constant at high redshift, as well as a possible spatial variation (see [55] and references therein, as well as the reviews [16,17])), we will address here this topic from the point of view of the implications of the running vacuum energy density throughout the cosmic history. We have mentioned before that in the electroweak sector of the SM is not possible to establish a connection between the cosmological evolution of the weak and electromagnetic couplings to the particle masses because there is no analogue in this sector of the QCD scale parameter Λ QCD . Notwithstanding, it is still possible to relate the electroweak couplings to Λ QCD itself in an indirect way if we use the hypothesis of Grand Unification of the SM couplings at a very high energy scale. We will focus here on the fine-structure constant α em = e 2 /4π and its correlated time-evolution with the strong coupling counterpart α s = g 2 s /4π, and ultimately with the time evolution of Λ QCD and µ = m p /m e .
In a Grand Unified Theory (GUT), the gauge couplings, and in particular the strong gauge coupling given in Eq. (15), can be made to converge to a unification point (at a high energy scale M X ) with the electroweak couplings. This is possible if more matter content is appropriately added (e.g. from supersymmetric particles), in which case M X ∼ 10 16 GeV [56]. This feature can be used as a theoretical argument to connect the possible time variation of the running coupling constants [24,57]. Let dα i /dz be the variation of α i with the cosmological redshift z. Such variation is possible if we have a consistent theoretical framework supporting this possibility, such as the RVM picture described in sect.3. Each of the couplings α i = g 2 i /4π is a function of the running scale µ R , and they follow the standard (1-loop) running laws to which we have appended the redshift variable to parameterize the cosmic evolution. Since the β-function coefficients b i of the running are constant in time and redshift, it follows that the it is a RG-invariant. Using this property and our ansatz concerning the cosmological evolution of the particle masses in the RVM, one can show that the running of the electromagnetic coupling α em is related to the corresponding cosmic running of the strong coupling α s as follows [24,57]: Combining this expression with Eq. (16) we can reexpress the cosmic running of α em in terms of the cosmic running of the QCD scale, and we find: At the Z-boson mass scale µ R = M Z , where both α em and α s are known with precision, one The above equation can now be nicely connected with our discussion of the cosmic running of the particle masses considered in sect. 4.1. Indeed, we have seen that the proton mass receives the bulk of its contribution from Λ QCD through m p ≃ c QCD Λ QCD (with a negligible contribution from the quark masses and the electromagnetism). Therefore, substituting this expression in Eq. (36) and integrating, and then inserting the redshift dependence of the proton mass through (27), we find: where α em (z) stands for the value of the fine structure constant at redshift z at a fixed value of µ R , and α 0 em (z) is its current value (z = 0). Since ν B is a small parameter, related to the fitted value ν m ∼ 10 −3 [38,39,40] through (25), we can estimate the relative variation of the electromagnetic fine structure constant with the redshift as follows: Defining ∆α em /α em = 3ν em by analogy with Eq. (24), we learn that the effective running index of the em coupling is some 30 times smaller than that of the baryonic index and with opposite sign, in other words ν em ≃ −0.03 ν B (up to logarithmic evolution with the redshift).
Let us now consider the observational situation concerning the measurements of the possible time variation of α em and the implied restrictions on the parameter ν B from Eq. (38). Table 1 shows recent measurements of α em in the redshift range 1.08 ≤ z ≤ 1.84. In Fig. 8 we have plotted these measurements. Let us also mention that apart from the astrophysical measurements of the time variation of α em , which give access to large look-back times of order of few billion years from now (corresponding to large redshifts of the order of those indicated in Table 1), there is a parallel research line of high precision laboratory based measurements whose look-back times are z ∆α/α(ppm) Ref. of course necessarily much more modest, but whose outstanding precision (based on state-of-art quantum optic techniques involving atomic clocks) can be highly competitive [1]. For instance, by comparing hyperfine transitions of different chemical elements with the cesium atomic clock, one can derive constrains on the time variation of the fundamental constants. Usually these measurements are correlated with the time dependence of the ratio µ = m p /m e through the the nuclear magnetic moment. In all these cases the typical result within errors (mostly compatible with zero) is |∆α em /α em | 10 −17 − 10 −16 yr −1 . These correspond once more to an upper bound on a relative variation |∆α em /α em | of order of 1 − 10 ppm, see [16,17], and therefore competitive with the astrophysical measurements. A particular lab experiment [64] employing narrow optical transitions in Hg + and Al + ions was able to directly measure changes in α em independent of other parameters. Comparing the transition frequencies over 12 months, the experiment renders a drift rate of ∆α em /α em = (−1.6 ± 2.3)10 −17 yr −1 , thus providing an upper bound on |∆α em /α em | of order of a few ppm. What is the possible impact of the RVM here? It is remarkable that the above mentioned results, whether from astrophysical or lab measurements, can be accounted for (in order of magnitude) within the RVM in combination with the GUT hypothesis. Indeed, we can see that the theoretical RVM prediction falls right within the order of magnitude of the typical measurements quoted in Table 1 and in Fig. 8, provided ν B lies in the range from 10 −4 to 10 −5 . This follows from Eq. (38), which, roughly speaking, says that the RVM prediction is of order ∆α em /α em ∼ −0.09ν B up to log corrections in the redshift. More precisely, in Fig. 8 we have superimposed the exact theoretical prediction α em (z) according to the formula (38). We can see that, notwithstanding the sizeable error bars, the trend of the measurements in Table 1 suggests a decrease of α em with the redshift (as there are more points compatible with ∆α em < 0 than points compatible with ∆α em > 0). This behavior has been previously noted in the literature [55] and is roughly in accordance with our theoretical curves in Fig. 8. But of course we need more precise measurements to confirm the real tendency of the data, as the errors are still too large and no firm conclusion is currently possible.
The following remarks should be emphasized. First, despite it is not possible to be more precise concerning the best fit value for ν B , in all cases the measurements in Table 1 indicate a maximum effect of 1−10 ppm, i.e. |∆α em /α em | at the level of 10 −6 to 10 −5 . This can be accommodated in the RVM framework since ν B describes the effect from only the baryonic component in Eq. (25). Such component should be naturally smaller than the value of the total index ν m ∼ 10 −3 fitted on the  Table 1 on the relative variation ∆α em /α em at different redshifts (in ppm). The solid and dashed lines correspond to the theoretical combined RVM-GUT prediction based on formula (38) for the values ν B = 10 −4 and 10 −5 , respectively. The tendency of the data to reflect smaller values of α em at large z is correctly described by the theoretical curves (which indicate ∆α em < 0), although the current observational errors are still too large.
basis of the overall analysis involving the SNIa+BAO+H(z)+LSS+BBN+CMB observables [38,39,30]. Second, the correct order of magnitude for ν B , which we have obtained from the direct ∆α em /α em observations (viz. ν B ∼ 10 −4 − 10 −5 ) does coincide with the result inferred from our previous considerations on the alternative observable ∆µ/µ in sect. 4.1. Put another way, we could have input the value of ν B needed to explain the typical measurements of the observable ∆µ/µ and we would have naturally predicted the typical range of values of ∆α em /α em derived from direct observations, and vice versa. As noted previously, this is because the RVM in combination with the GUT framework neatly predicts the relation ν em ≃ −0.03 ν B . In the light of the above results, the baryonic index ν B in Eq. (25) is definitely subdominant as compared to the dark matter one, |ν B | ≪ |ν X |, and hence ν X must be of order of the total matter index ν m ∼ 10 −3 fitted from the overall cosmological observations [38,39,30]. In other words, we find once more that it must be the DM component that provides the bulk of the contribution to the time variation of masses in the Universe. This fact was not obvious a priori, and is not necessarily related to the overwhelming abundance of DM as compared to baryons, for the large amounts of DM could simply remain passive and not evolve at all throughout the cosmic expansion. If the best fit value of the total mass variation index ν m would have been, say of order 10 −5 , equation (25) could have been naturally fulfilled with ν X << ν B ∼ 10 −4 and this would still be compatible with the measurements in Table 1. However, the fact that the value of ν m (obtained from the overall cosmological fit to the data within the RVM [38,39,40]) comes out significantly larger than the baryonic index ν B has nontrivial consequences and provides an independent hint of the need for (time-evolving) dark matter. Taking into account that we have been able to infer this same conclusion from the analysis of the two independent observables ∆µ/µ and ∆α em /α em , which become correlated in this theoretical framework, does reinforce the RVM scenario and places the contribution from the DM component to the forefront of our considerations concerning the total mass drift rate in the Universe [24].

Alternative dynamical vacuum models interacting with matter
Let us finally address some alternative scenarios for the non-conservation of the particle masses in an expanding universe. Our starting point is a phenomenological coupling between the vacuum energy density and the matter density, where for simplicity we now assume G =const. The background evolution is then encoded in the coupled system of local energy conservation equationṡ andρ where Q denotes the background energy source between dark matter and vacuum energy (or in general some dark energy source). From the previous equations we see that for Q > 0 the matter energy density increases whereas the vacuum energy density decreases, and hence the energy flows from vacuum to matter, and vice versa for Q featuring the opposite sign. In other words, for Q > 0 the vacuum is decaying into matter whereas for Q < 0 the matter decays into vacuum. Of these two options the naturally preferred one, at least from the point of view of the second principle of thermodynamics, should be the first one. Despite it is usually assumed that the vacuum decays only into DM, this effect has little quantitative implications for the present analysis since we have seen that the baryonic component is essentially conserved (ξ B << ξ X , see sect. 4.1). Therefore we shall not consider this correction here. For more details, see [29]. Many functions can be proposed for the interacting energy source Q, see e.g. [65,66,67]. For our purposes, it will suffice to consider two frequently discussed phenomenological expressions in the literature, namely a source proportional to ρ Λ in the form Q = q Λ Hρ Λ (hereafter called "q Λmodel"), and a source proportional to the matter density, Q = q m Hρ m ("q m -model"), where q Λ and q m are small dimensionless parameters (|q Λ |, |q m | ≪ 1). Considering the q Λ -model, we easily find from (40) the explicit form of the vacuum evolution law which can be substituted in (39) and upon integration we derive the corresponding matter density evolution in that model: Obviously equations (41) and (42) are very different from the corresponding ones from the RVM, see e.g. equations (5) and (8) for ν G = 0 (since we are now considering G =const.). With the help of these results and assuming that the number density of particles is conserved, i.e., n(a) = n 0 a −3 (cf. sect. 3), we can obtain the law for the cosmological evolution of masses. Notice that ρ m (a) = n(a)m(a) and ρ 0 m = n 0 m 0 . Phrased in terms of the redshift, the relative mass variation can be cast as follows: Notice that δm/δz = −q Λ m 0 (Ω Λ /Ω m ) (1 + z) q Λ −4 , so that the mass tends to decrease or increase with the redshift (equivalently, increase or decrease with the expansion) for q Λ > 0 or q Λ < 0 respectively (recall that |q Λ | ≪ 1). In the remote past (z >> 1), we have the initial value m i ≃ m 0 (1 − q Λ Ω Λ /3Ω m ), whereas in the remote future (z → −1, i.e. a → ∞) the asymptotic final mass value evolve as m f ≃ m 0 (q Λ Ω Λ /3Ω m ) a 3−q Λ . For q Λ < 0 we have Q < 0, which we already pointed out as being unfavored by the second law of thermodynamics since the energy flows from matter into vacuum. Even worse, for q Λ < 0 the asymptotic mass values m f become eventually negative, which would correspond to a rather unstable situation for the Universe. The lack of sympathy for this sign is explicitly confirmed by the analysis of the current observational data when confronted with the various types of dynamical vacuum models. One finds that q Λ < 0 gives a very bad fit as compared to q Λ > 0, see the study of Ref. [40]. Let us next assess the situation with the second alternative model for the dark energy source mentioned above, Q = q m Hρ m (the q m -model). We can easily integrate equations (39) and (40) anew, with the following results: and Clearly these expressions are formally similar to the corresponding ones in the RVM, i.e. equations (5) and (8), with the identification q m = 3ν m and taking into account that ν G = 0 here (i.e. G is constant). The mass evolution with the redshift is therefore as in Eq. (6), i.e. m(z) = m 0 (1 + z) −qm .
We find that m(z) decreases with z, and hence increases with the expansion, if q m > 0 (corresponding to a situation of decay of vacuum into matter); and m(z) decreases with the expansion if q m < 0 (when matter decays into vacuum). In contradistinction to the q Λ -model, we note that the masses remain now always positive irrespective of the sign of q m . We point out that in the two models considered in this section (the q Λ -model and the q m -model) there exists no special motivation for the particular form they have for the interacting source Q. In contrast, in the RVM case the form of Q is not taken as a mere phenomenological ansatz, it is theoretically determined. In fact, it can be derived from the dynamical vacuum equation (11), which leads to Eq. (8) when the model is explicitly solved. For example, take G =const. (ν G = 0) in which case the local conservation equation for the RVM can be put in the form (39) with Q = −ρ Λ . We may compute explicitly the time derivative of ρ Λ and reexpress the result in terms of the redshift using Eq. (1). We find: where use has been made of the expressions for ρ m (z) and ρ Λ (z) given in equations (5) and (8) with ν G = 0. Alternatively, we may use the generalized conservation law Eq. (3) with G ′ (z) = 0; and we find, once more, Q = ρ ′ Λ (z) (1 + z) H = (3ρ m (z) − ρ ′ m (z)(1 + z))H = +3 ν m ρ m H. In either way we arrive atρ This equation clearly shows that the RVM with G =const. behaves effectively (at least in the matter-dominated epoch) as a q m -model with q m = 3ν m . When radiation is included the two models present significant differences. These differences have been accounted for in [40,29] and they have implications for the overall fit to the data, especially for the high redshift data of course. The net outcome is that the RVM does better than the q m -model, but both of them do better than the ΛCDM, see [29] for details.
Recall that ultimately the behavior of the RVM stems from the dynamical vacuum equation (11), which follows from the RG-flow in QFT in curved spacetime [27,28]. As mentioned in sect. 3, Eq. (11) can be extended with higher powers of H so as to include inflation in a single unified theory describing the cosmic evolution from the early universe until the current one [28,33,34,35,36].
From the foregoing discussion it is clear that different dynamical vacuum energy models exist for describing the possible time variation of the fundamental constants in a framework which is consistent with GR. Some of these models are more phenomenological, but in the RVM case there is a more concrete theoretical motivation for the dynamical vacuum structure. It suggests that a slow change in the fundamental "constants" can be theoretically motivated and cannot be ruled out at present. The subject is therefore worthwhile being further investigated in the light of new data, as it can reveal new clues to fundamental physics.

Discussion and conclusions
In this paper we have addressed the running vacuum models (RVM) of the cosmic evolution and the possible implications they could have in explaining the reported hints of the time variation of the so-called fundamental constants of Nature, such as masses, coupling constants etc, including the gravitational coupling G and the Λ-term in Einstein's equations. The impact from the RVM on this issue stems from the cosmological energy exchange between vacuum, matter and the possible interplay with the Newtonian coupling G and the vacuum energy density ρ Λ = Λ/8πG. Because the possible cosmological running of these quantities is controlled by the Hubble parameter, the RVM predicts that the associated rhythm of change, i.e. the drift rate of the fundamental constants should naturally be as moderate as dictated by the expansion rate of the Universe at any given instant of the cosmic history. On this basis it is possible to compute the time evolution of the vacuum energy density and the corresponding change of the gravitational coupling and the particle masses. Combining this scenario with the GUT hypothesis we have obtained a prediction for the time evolution of the fine structure constant correlated with the time evolution of the proton mass, or more precisely the proton-to-electron mass ratio µ = m p /m e .
Taking into account that the small, but nonvanishing, rate of change of the vacuum energy density has been fitted to the cosmological data at a rather significant confidence level -see the recent studies [38,39,40,41] -and baring in mind that such vacuum rate of change impinges on the corresponding variation of the particle masses, we conclude that the mass variation must be essentially supported by the mass drift rate of the dark matter (DM) particles (since the time evolution of the baryonic component is found to be some two orders of magnitude smaller). This can be interpreted as an indirect alternative hint of the need for DM. If in the future the precision of these experiments further improves we might well find ourselves on the verge of measuring these subtle effects and perhaps be in position to check if they can be explained within the kind of theoretical running vacuum models that we have studied here.