The $\Lambda$-term in Einstein's equations is a fundamental building block of the `concordance' $\Lambda$CDM model of cosmology. Even though the model is not free of fundamental problems, they have not been circumvented by any alternative dark energy proposal either. Here we stick to the $\Lambda$-term, but we contend that it can be a `running quantity' in quantum field theory (QFT) in curved spacetime. A plethora of phenomenological works have shown that this option can be highly competitive with the $\Lambda$CDM with a rigid cosmological term. The, so-called, `running vacuum models' (RVM's) are characterized by the vacuum energy density, $\rho_{vac}$, being a series of (even) powers of the Hubble parameter and its time derivatives. Such theoretical form has been motivated by general renormalization group arguments, which look plausible. Here we dwell further upon the origin of the RVM structure within QFT in FLRW spacetime. We compute the renormalized energy-momentum tensor with the help of the adiabatic regularization procedure and find that it leads essentially to the RVM form. This means that $\rho_{vac}(H)$ evolves as a constant term plus dynamical components ${\cal O}(H^2)$ and ${\cal O}(H^4)$, the latter being relevant for the early universe only. However, the renormalized $\rho_{vac}(H)$ does not carry dangerous terms proportional to the quartic power of the masses ($\sim m^4$) of the fields, these terms being a well-known source of exceedingly large contributions. At present, $\rho_{vac}(H)$ is dominated by the additive constant term accompanied by a mild dynamical component $\sim \nu H^2$ ($|\nu|\ll1$), which mimics quintessence.

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.