We present a new isocurvature mode present in scalar-tensor theories of gravity that corresponds to a regular growing solution in which the energy of the relativistic degrees of freedom and the scalar field that regulates the gravitational strength compensate during the radiation dominated epoch on scales much larger than the Hubble radius. We study this isocurvature mode and its impact on anisotropies of the cosmic microwave background for the simplest scalar-tensor theory, i.e. the extended Jordan-Brans-Dicke gravity, in which the scalar field also drives the acceleration of the Universe. We use Planck data to constrain the amplitude of this isocurvature mode in the case of fixed correlation with the adiabatic mode and we show how this mode could be generated in a simple two field inflation model.

We use $Planck$ 2018 data to constrain the simplest models of scalar-tensor theories characterized by a coupling to the Ricci scalar of the type $F(\sigma) R$ with $F(\sigma) = N_{pl}^2 + \xi \sigma^2$. We update our results with previous $Planck$ and BAO data releases obtaining the tightest constraints to date on the coupling parameters, that is $\xi < 5.5 \times 10^{-4}$ for $N_{pl}=0$ (induced gravity or equivalently extended Jordan-Brans-Dicke) and $(N_{pl} \sqrt{8 \pi G})-1 < 1.8 \times 10^{-5}$ for $\xi = -1/6$ (conformal coupling), both at 95% CL. Because of a modified expansion history after radiation-matter equality compared to the $\Lambda$CDM model, all these dynamical models accommodate a higher value for $H_0$ and therefore alleviate the tension between $Planck$/BAO and distance-ladder measurement from SNe Ia data from $4.4\sigma$ at best to $2.3\sigma$. We show that all these results are robust to changes in the neutrino physics. In comparison to the $\Lambda$CDM model, partial degeneracies between neutrino physics and the coupling to the Ricci scalar allow for smaller values $N_{\rm eff} \sim 2.8$, $1\sigma$ lower compared to the standard $N_{\rm eff} = 3.046$, and relax the upper limit on the neutrino mass up to 40%.

Branches of cosmic inflationary models, such as slow-roll inflation, predict a background of primordial gravitational waves that imprints a unique odd-parity B-mode pattern in the Cosmic Microwave Background (CMB) at amplitudes that are within experimental reach. The BICEP/Keck (BK) experiment targets this primordial signature, the amplitude of which is parameterized by the tensor-to-scalar ratio r, by observing the polarized microwave sky through the exceptionally clean and stable atmosphere at the South Pole. B-mode measurements require an instrument with exquisite sensitivity, tight control of systematics, and wide frequency coverage to disentangle the primordial signal from the Galactic foregrounds. BICEP Array represents the most recent stage of the BK program, and comprises four BICEP3-class receivers observing at 30/40, 95, 150 and 220/270 GHz. The 30/40 GHz receiver will be deployed at the South Pole during the 2019/2020 austral summer. After 3 full years of observations with 30,000+ detectors, BICEP Array will measure primordial gravitational waves to a precision $\sigma(r)$ between 0.002 and 0.004, depending on foreground complexity and the degree of lensing removal. In this paper we give an overview of the instrument, highlighting the design features in terms of cryogenics, magnetic shielding, detectors and readout architecture as well as reporting on the integration and tests that are ongoing with the first receiver at 30/40 GHz.

The BICEP/Keck (BK) experiment aims to detect the imprint of primordial gravitational waves in the Cosmic Microwave Background polarization, which would be direct evidence of the inflation theory. While the tensor-to-scalar ratio has been constrained to be r_0.05 < 0.06 at 95% c.l., further improvements on this upper limit are hindered by polarized Galactic foreground emissions and removal of gravitational lensing polarization. The 30/40 GHz receiver of the BICEP Array (BA) will deploy at the end of 2019 and will constrain the synchrotron foreground with unprecedented accuracy within the BK sky patch. We will show the design of the 30/40 GHz detectors and test results summarizing its performance. The low optical and atmospheric loading at these frequencies requires our TES detectors to have low saturation power in order to be photon-noise dominated. To realize the low thermal conductivity required from a 250 mK base temperature, we developed new bolometer leg designs. We will present the relevant measured detector parameters: G, Tc, Rn, Psat , and spectral bands, and noise spectra. We achieved a per bolometer NEP including all noise components of 2.07E-17 W/sqrt(Hz), including an anticipated photon noise level 1.54E-17 W/sqrt(Hz).

We study the cosmological constraints on the variation of the Newton's constant and on post-Newtonian parameters for simple models of scalar-tensor theory of gravity beyond the extended Jordan-Brans-Dicke theory. We restrict ourselves to an effectively massless scalar field with a potential $V \propto F^2$, where $F(\sigma)=N_{pl}^2+\xi\sigma^2$ is the coupling to the Ricci scalar considered. We derive the theoretical predictions for cosmic microwave background (CMB) anisotropies and matter power spectra by requiring that the effective gravitational strength at present is compatible with the one measured in a Cavendish-like experiment and by assuming adiabatic initial condition for scalar fluctuations. When comparing these models with $Planck$ 2015 and a compilation of baryonic acoustic oscilation (BAO) data, all these models accomodate a marginalized value for $H_0$ higher than in $\Lambda$CDM. We find no evidence for a statistically significant deviation from Einstein's general relativity. We find $\xi < 0.064$ ($|\xi| < 0.011$) at 95 % CL for $\xi > 0$ (for $\xi < 0$, $\xi \ne -1/6$). In terms of post-Newtonian parameters, we find $0.995 < \gamma_{\rm PN} < 1$ and $0.99987 < \beta_{\rm PN} < 1$ ($0.997 < \gamma_{\rm PN} < 1$ and $1 < \beta_{\rm PN} < 1.000011$) for $\xi >0$ (for $\xi < 0$). For the particular case of the conformal coupling, i.e. $\xi=-1/6$, we find constraints on the post-Newtonian parameters of similar precision to those within the Solar System.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

The extended Jordan-Brans-Dicke (eJBD) theory of gravity is constrained by a host of astrophysical and cosmological observations spanning a wide range of scales. The current cosmological constraints on the first post-Newtonian parameter in these simplest eJBD models in which the recent acceleration of the Universe is connected with the variation of the effective gravitational strength are consistent, but approximately two order of magnitude larger than the time-delay test within the Solar System. We forecast the capabilities of future galaxy surveys in combination with current and future CMB anisotropies measurements to further constrain the simplest dark energy models within eJBD theory of gravity. By considering two cases of a monomial potential (a quartic potential or a cosmological constant), we show how Euclid-like galaxy clustering and weak lensing data in combination with BOSS and future CMB observations have the potential to reach constraints on the first post-Newtonian parameter $\gamma_{\rm PN}$ comparable to those from the Solar System.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

We investigate the extent to which foreground cleaned CMB maps can be used to estimate the cosmological parameters at small scales. We use the SMICA method, a blind separation technique which works directly at the spectral level. In this work we focus on the small scales of the CMB angular power spectrum, which are chiefly affected by noise and extragalactic foregrounds, such as point sources. We adapt SMICA to use only cross-spectra between data maps, thus avoiding the noise bias. In this study, performed both on simulations and on Planck 2015 data, we fit for extragalactic point sources by modeling them as shot noise of two independent populations. In simulations we correctly recover the point source emission law, and obtain a CMB angular power spectrum that has an average foreground residual of one fifth of the CMB power at $\ell \geq$ 2200. On Planck data, the recovered point source emission law corresponds to external estimates, with some offsets at the highest and lowest frequencies, possibly due to frequency decoherence of point sources. The CMB angular power spectrum residuals are consistent with what we find in simulations. The cosmological parameters obtained from the simulations and the data show offsets up to 1$\sigma$ on average from their expected values. Biases on cosmological parameters in simulations represent the expected level of bias in Planck data. The results on cosmological parameters depend on the detail of the foreground residual contamination in the spectrum, and need a tailored modeling of the likelihood foreground model.

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

We study induced gravity dark energy models coupled with a simple monomial potential $\propto \sigma^n$ and a positive exponent $n$. These simple potentials lead to viable dark energy models with a weak dependence on the exponent, which characterizes the accelerated expansion of the cosmological model in the asymptotic attractor, when ordinary matter becomes negligible. We use recent cosmological data to constrain the coupling $\gamma$ to the Ricci curvature, under the assumptions that the scalar field starts at rest deep in the radiation era and that the gravitational constant in the Einstein equations is compatible with the one measured in a Cavendish-like experiment. By using $Planck$ 2015 data only, we obtain the 95 % CL bound $\gamma < 0.0017$ for $n=4$, which is further tightened to $\gamma < 0.00075$ by adding Baryonic Acoustic Oscillations (BAO) data. This latter bound improves by $\sim 30$ % the limit obtained with the $Planck$ 2013 data and the same compilation of BAO data. We discuss the dependence of the $\gamma$ and $\dot G_N/G_N (z=0)$ on $n$.

We study the predictions for structure formation in an induced gravity dark energy model with a quartic potential. By developing a dedicated Einstein-Boltzmann code, we study self-consistently the dynamics of homogeneous cosmology and of linear perturbations without using any parametrization. By evolving linear perturbations with initial conditions in the radiation era, we accurately recover the quasi-static analytic approximation in the matter dominated era. We use Planck 2013 data and a compilation of baryonic acoustic oscillation (BAO) data to constrain the coupling $\gamma$ to the Ricci curvature and the other cosmological parameters. By connecting the gravitational constant in the Einstein equation to the one measured in a Cavendish-like experiment, we find $\gamma < 0.0012$ at 95% CL with Planck 2013 and BAO data. This is the tightest cosmological constraint on $\gamma$ and on the corresponding derived post-Newtonian parameters. Because of a degeneracy between $\gamma$ and the Hubble constant $H_0$, we show how larger values for $\gamma$ are allowed, but not preferred at a significant statistical level, when local measurements of $H_0$ are combined in the analysis with Planck 2013 data.