We present a full-fledged analysis of Brans-Dicke cosmology with a cosmological constant and cold dark matter (BD-$\Lambda$CDM for short). We extend the scenarios where the current cosmological value of the BD-field is restricted by the local astrophysical domain to scenarios where that value is fixed only by the cosmological observations, which should be more natural in view of the possible existence of local screening mechanims. Our analysis includes both the background and perturbations equations in different gauges. We find that the BD-$\Lambda$CDM is favored by the overall cosmological data as compared to the concordance GR-$\Lambda$CDM model, namely data on distant supernovae, cosmic chronometers, local measurements of the Hubble parameter, baryonic acoustic oscillations, Large-Scale Structure formation and the cosmic microwave background under full Planck 2018 CMB likelihood. We also test the impact of Strong and Weak-Lensing data on our results, which can be significant. We find that the BD-$\Lambda$CDM can mimic effective quintessence with a significance of about $3-3.5\sigma$ c.l. (depending on the lensing datasets). The fact that the BD-$\Lambda$CDM behaves effectively as a Running Vacuum Model (RVM) when viewed from the GR perspective helps to alleviate some of the existing tensions with the data, such as the $\sigma_8$ excess predicted by GR-$\Lambda$CDM. On the other hand, the BD-$\Lambda$CDM model has a crucial bearing on the acute $H_0$-tension with the local measurements, which is rendered virtually harmless owing to the small increase of the effective value of the gravitational constant with the expansion. The simultaneous alleviation of the two tensions is a most remarkable feature of BD-gravity with a cosmological constant in the light of the current observations, and hence goes in support of BD-$\Lambda$CDM against GR-$\Lambda$CDM

We analyze Brans-Dicke gravity with a cosmological constant, $\Lambda$, and cold dark matter (BD-$\Lambda$CDM for short) in the light of the latest cosmological observations on distant supernovae, Hubble rate measurements at different redshifts, baryonic acoustic oscillations, large scale structure formation data, gravitational weak-lensing and the cosmic microwave background under full Planck 2015 CMB likelihood. Our analysis includes both the background and perturbations equations. We find that BD-$\Lambda$CDM is observationally favored as compared to the concordance $\Lambda$CDM model, which is traditionally defined within General Relativity (GR). In particular, some well-known persisting tensions of the $\Lambda$CDM with the data, such as the excess in the mass fluctuation amplitude $\sigma_8$ and specially the acute $H_0$-tension with the local measurements, essentially disappear in this context. Furthermore, viewed from the GR standpoint, BD-$\Lambda$CDM cosmology mimics quintessence at $\gtrsim3\sigma$ c.l. near our time.

Investigations on dark energy (DE) are currently inconclusive about its time evolution. Hints of this possibility do however glow now and then in the horizon. Herein we assess the current status of dynamical dark energy (DDE) in the light of a large body of updated $SNIa+H(z)+BAO+LSS+CMB$ observations, using the full Planck 2015 CMB likelihood. The performance of the $\Lambda$CDM model (with equation of state $w=-1$ for $\Lambda$) is confronted with that of the general XCDM and CPL parametrizations, as well as with the traditional $\Phi$CDM model based on the scalar field potential $V\sim \Phi^{-\alpha}$. In particular, we gauge the impact of the bispectrum in the LSS and BAO parts, and show that the subset of $CMB+BAO+LSS$ observations may contain the bulk of the DDE signal. The departure from $w=-1$ is significant: roughly $2.6\sigma$ for XCDM and $2.9\sigma$ for $\Phi$CDM. In both cases the full Bayesian evidence is found to be positive even for a prior range of the DDE parameters extending over several standard deviations from the mean when the bispectrum is taken into account. Positive signs follow as well from the preliminary results of Planck 2018 data using the compressed CMB likelihood. As a bonus we also find that the $\sigma_8$-tension becomes reduced with DDE.

We consider the impact of dynamical dark energy (DDE) in the possible solution of the existing tensions in the $\Lambda$CDM. We test both interacting and non-interacting DE models with dark matter (DM). Among the former, the running vacuum model (RVM) interacting with DM appears as a favored option. The non-interacting scalar field model based on the potential $V\sim \phi^{-\alpha}$, and the generic XCDM parametrization, also provide consistent signs of DDE. The important novelty of our analysis with respect to the existing ones in the literature is that we use the matter bispectrum, together with the power spectrum. Using a complete and updated set of cosmological observations on $SNIa+BAO+H(z)+LSS+CMB$, we find that the crucial triad $BAO+LSS+CMB$ (i.e. baryonic acoustic oscillations, large scale structure formation data and the cosmic microwave background) provide the bulk of the signal. The bispectrum data is instrumental to get hold of the DDE signal, as our analysis shows. If the bispectrum is not included, the DDE signal could not be currently perceived at a significant confidence level.

We study a generic class of time-evolving vacuum models which can provide a better phenomenological account of the overall cosmological observations as compared to the $\Lambda$CDM. Among these models, the running vacuum model (RVM) appears to be the most motivated and favored one, at a confidence level of $\sim 3\sigma$. We further support these results by computing the Akaike and Bayesian information criteria. Our analysis also shows that we can extract fair signals of dynamical dark energy (DDE) by confronting the same set of data to the generic XCDM and CPL parametrizations. In all cases we confirm that the combined triad of modern observations on Baryonic Acoustic Oscillations, Large Scale Structure formation, and the Cosmic Microwave Background, provide the bulk of the signal sustaining a possible vacuum dynamics. In the absence of any of these three crucial data sources, the DDE signal can not be perceived at a significant confidence level. Its possible existence could be a cure for some of the tensions existing in the $\Lambda$CDM when confronted to observations.

After about two decades of the first observational papers confirming the accelerated expansion of the universe, we are still facing the question whether the cause of it is a rigid cosmological constant $\Lambda$-term or a mildly evolving dynamical dark energy (DDE). While studies focusing mainly on CMB measurements do not perceive signs of physics beyond the $\Lambda$CDM, in this work we show that if we take a large string $SNIa+BAO+H(z)+LSS+CMB$ of modern cosmological observations, in which not only the CMB but also a rich sample of large scale structure formation data are included, one can extract $\sim 3.3\sigma$ signs of DDE using a simple XCDM parameterization. These signs can be enhanced up to near $3.8\sigma$ in the context of the running vacuum model (RVM), in which the vacuum energy density is in interaction with dark matter. Recently the RVM has been shown to provide an efficient and economical solution to the $\sigma_8$-tension, which is one of the intriguing phenomenological problems that has not been possible to solve within the $\Lambda$CDM so far. This fact contributes to strengthen the possibility that dynamical vacuum energy, or in general DDE, could be presently favored by the observations.

It has recently been shown that the class of running vacuum models (RVMs) has the capacity to fit the overall cosmological observations better than the concordance $\Lambda$CDM model, therefore supporting the possibility of dynamical dark energy (DE). Apart from the cosmic microwave background (CMB) anisotropies, the most crucial datasets involved are: i) baryonic acoustic oscillations (BAO), and ii) direct large scale structure (LSS) formation data. Analyses mainly focusing on CMB and with insufficient BAO+LSS input generally fail to capture the dynamical DE signature, whereas the few existing studies accounting for the wealth of known CMB+BAO+LSS data (see in particular Sol\`a, G\'omez-Valent \& de Cruz P\'erez 2015, 2017; and Zhao et al. 2017) do converge to the remarkable conclusion that dynamical DE might well be encoded in the current cosmological observations at a $3-4\sigma$ c.l. A decisive factor is the persistent $\sigma_8$-tension between the $\Lambda$CDM and the data. Because the issue is obviously pressing, we devote this work to explain how and why running vacuum in the expanding universe successfully relaxes the existing $\sigma_8$-tension and describes the LSS formation data significantly better than the $\Lambda$CDM.

Despite the outstanding achievements of modern cosmology, the classical dispute on the precise value of $H_0$, which is the first ever parameter of modern cosmology and one of the prime parameters in the field, still goes on and on after over half a century of measurements. Recently the dispute came to the spotlight with renewed strength owing to the significant tension (at $>3\sigma$ c.l.) between the latest Planck determination obtained from the CMB anisotropies and the local (distance ladder) measurement from the Hubble Space Telescope (HST), based on Cepheids. In this work, we investigate the impact of the running vacuum model (RVM) and related models on such a controversy. For the RVM, the vacuum energy density $\rho_{\Lambda}$ carries a mild dependence on the cosmic expansion rate, i.e. $\rho_{\Lambda}(H)$, which allows to ameliorate the fit quality to the overall $SNIa+BAO+H(z)+LSS+CMB$ cosmological observations as compared to the concordance $\Lambda$CDM model. By letting the RVM to deviate from the vacuum option, the equation of state $w=-1$ continues to be favored by the overall fit. Vacuum dynamics also predicts the following: i) the CMB range of values for $H_0$ is more favored than the local ones, and ii) smaller values for $\sigma_8(0)$. As a result, a better account for the LSS structure formation data is achieved as compared to the $\Lambda$CDM, which is based on a rigid (i.e. non-dynamical) $\Lambda$ term.

In this year, in which we celebrate 100 years of the cosmological term, $\Lambda$, in Einstein's gravitational field equations, we are still facing the crucial question whether $\Lambda$ is truly a fundamental constant or a mildly evolving dynamical variable. After many theoretical attempts to understand the meaning of $\Lambda$, and in view of the enhanced accuracy of the cosmological observations, it seems now mandatory that this issue should be first settled empirically before further theoretical speculations on its ultimate nature. In this work, we summarize the situation of some of these studies. Devoted analyses made recently show that the $\Lambda=$const. hypothesis, despite being the simplest, may well not be the most favored one. The overall fit to the cosmological observables $SNIa+BAO+H(z)+LSS+CMB$ singles out the class RVM of the "running" vacuum models, in which $\Lambda=\Lambda(H)$ is an affine power-law function of the Hubble rate. It turns out that the performance of the RVM as compared to the "concordance" $\Lambda$CDM model (with $\Lambda=$const.) is much better. The evidence in support of the RVM may reach $\sim 4\sigma$ c.l., and is bolstered with Akaike and Bayesian criteria providing strong evidence in favor of the RVM option. We also address the implications of this framework on the tension between the CMB and local measurements of the current Hubble parameter.

Recent analyses in the literature suggest that the concordance $\Lambda$CDM model with rigid cosmological term, $\Lambda=$const., may not be the best description of the cosmic acceleration. The class of "running vacuum models", in which $\Lambda=\Lambda(H)$ evolves with the Hubble rate, has been shown to fit the string of $SNIa+BAO+H(z)+LSS+CMB$ data significantly better than the $\Lambda$CDM. Here we provide further evidence on the time-evolving nature of the dark energy (DE) by fitting the same cosmological data in terms of scalar fields. As a representative model we use the original Peebles & Ratra potential, $V\propto\Phi^{-\alpha}$. We find clear signs of dynamical DE at $\sim 4\sigma$ c.l., thus reconfirming through a nontrivial scalar field approach the strong hints formerly found with other models and parametrizations.

Despite the fact that a rigid $\Lambda$-term is a fundamental building block of the concordance $\Lambda$CDM model, we show that a large class of cosmological scenarios with dynamical vacuum energy density $\rho_{\Lambda}$ and/or gravitational coupling $G$, together with a possible non-conservation of matter, are capable of seriously challenging the traditional phenomenological success of the $\Lambda$CDM. In this paper, we discuss these "running vacuum models" (RVM's), in which $\rho_{\Lambda}=\rho_{\Lambda}(H)$ consists of a nonvanishing constant term and a series of powers of the Hubble rate. Such generic structure is potentially linked to the quantum field theoretical description of the expanding Universe. By performing an overall fit to the cosmological observables $SNIa+BAO+H(z)+LSS+BBN+CMB$ (in which the WMAP9, Planck 2013 and Planck 2015 data are taken into account), we find that the class of RVM's appears significantly more favored than the $\Lambda$CDM, namely at an unprecedented level of $\gtrsim4.2\sigma$. Furthermore, the Akaike and Bayesian information criteria confirm that the dynamical RVM's are strongly preferred as compared to the conventional rigid $\Lambda$-picture of the cosmic evolution.

We determine the Hubble expansion and the general cosmic perturbations equations for a general system consisting of self-conserved matter and self-conserved dark energy (DE). While at the background level the two components are non-interacting, they do interact at the perturbations level. We show that the coupled system of matter and DE perturbations can be transformed into a single, third order, matter perturbation equation, which reduces to the (derivative of the) standard one in the case that the DE is just a cosmological constant. As a nontrivial application we analyze a class of dynamical models whose DE density $\rho_D$ consists of a constant term, $C_0$, and a series of powers of the Hubble rate. These models were previously analyzed from the point of view of dynamical vacuum models, but here we treat them as self-conserved DE models with a dynamical equation of state. We fit them to the wealth of expansion history and linear structure formation data and compare the obtained fit quality with that of the concordance $\Lambda$CDM model. Those with $C_0=0$ include the so-called "entropic-force" and "QCD-ghost" DE models, as well as the pure linear model $\rho_D\sim H$, all of which appear strongly disfavored. The models with $C_0\neq 0$, in contrast, emerge as promising dynamical DE candidates whose phenomenological performance is highly competitive with the rigid $\Lambda$-term inherent to the $\Lambda$CDM.

Recently there have been claims on model-independent evidence of dynamical dark energy. Herein we consider a fairly general class of cosmological models with a time-evolving cosmological term of the form $\Lambda(H)=C_0+C_H H^2+C_{\dot{H}} \dot{H}$, where $H$ is the Hubble rate. These models are well motivated from the theoretical point of view since they can be related to the general form of the effective action of quantum field theory in curved spacetime. Consistency with matter conservation can be achieved by letting the Newtonian coupling $G$ change very slowly with the expansion. We solve these dynamical vacuum models and fit them to the wealth of expansion history and linear structure formation data. The results of our analysis show a significantly better agreement as compared to the concordance $\Lambda$CDM model, thus supporting the possibility of a dynamical cosmic vacuum.

We focus on the class of cosmological models with a time-evolving vacuum energy density of the form $\rho_\Lambda=C_0+C_1 H+C_2 H^2$, where $H$ is the Hubble rate. Higher powers of $H$ could be important for the early inflationary epoch, but are irrelevant afterwards. We study these models at the background level and at the perturbations level, both at the linear and at the nonlinear regime. We find that those with $C_0=0$ are seriously hampered, as they are unable to fit simultaneously the current observational data on Hubble expansion and the linear growth rate of clustering. This is in contrast to the $C_0\neq 0$ models, including the concordance $\Lambda$CDM model. We also compute the redshift distribution of clusters predicted by all these models, in which the analysis of the nonlinear perturbations becomes crucial. The outcome is that the models with $C_0=0$ predict a number of counts with respect to the concordance model which is much larger, or much smaller, than the $\Lambda$CDM and the dynamical models with $C_0\neq 0$. The particular case $\rho_\Lambda\propto H$ (the pure lineal model), which in the past was repeatedly motivated by several authors from QCD arguments applied to cosmology, is also addressed and we assess in detail its phenomenological status. We conclude that the most favored models are those with $C_0\neq 0$, and we show how to discriminate them from the $\Lambda$CDM.

Perhaps the deepest mystery of our accelerating Universe in expansion is the existence of a tiny and rigid cosmological constant, $\Lambda$. Its size is many orders of magnitude below the expected one in the standard model of particle physics. However, an expanding Universe is not expected to have a static vacuum energy density. We should rather observe a mildly dynamical behavior $\delta\Lambda(t)\sim R\sim H^2(t)$ with the expansion rate $H$. At the same time, it is natural to think that the huge value of the primeval vacuum energy (presumably connected to some grand unified theory) was responsible for the initial inflationary phase. In the traditional inflaton models such phase is inserted by hand in the early epoch of the cosmic evolution, and it is assumed to match the concordance $\Lambda$CDM regime during the radiation epoch. Here, instead, we consider a class of dynamical vacuum models which incorporate into a single vacuum structure $\bar{\Lambda}(H)$ the rapid stage of inflation, followed by the radiation and cold matter epochs, until achieving our dark energy Universe. The early behavior of the model bares resemblance with Starobinsky's inflation and ptovides a solution to the large entropy problem. It is compatible with the latest cosmological data on Hubble expansion and structure formation, and presents distinctive observational features that can be tested in the near future.

Despite the many efforts, our theoretical understanding of the ultimate nature of the dark energy component of the universe still lags well behind the astounding experimental evidence achieved from the increasingly sophisticated observational tools at our disposal. While the canonical possibility is a strict cosmological constant, or rigid vacuum energy density $\rho_{\Lambda}=$const., the exceeding simplicity of this possibility lies also at the root of its unconvincing theoretical status, as there is no explanation for the existence of such constant for the entire cosmic history. Herein we explore general models of the vacuum energy density slowly evolving with the Hubble function $H$ and/or its time derivative, $\rho_{\Lambda}(H,\dot{H})$. Some of these models are actually well-motivated from the theoretical point of view and may provide a rich phenomenology that could be explored in future observations, whereas some others have more limitations. In this work, we put them to the test and elucidate which ones are still compatible with the present observations and which ones are already ruled out. We consider their implications on structure formation, in combination with data on type Ia supernovae, the Cosmic Microwave Background, the Baryonic Acoustic Oscillations, and the predicted redshift distribution of cluster-size collapsed structures. The relation of these vacuum models on possible evidence of dynamical dark energy recently pointed out in the literature is also briefly addressed.