The expansion history of the Universe reconstructed from a combination of recent data indicates a preference for a changing Dark Energy (DE) density. Moreover, the DE density appears to be increasing with cosmic time, with its equation of state being below -1 on average, and possibly crossing the so-called phantom divide. Scalar-tensor theories, in which the scalar field mediates a force between matter particles, offer a natural framework in which the effective DE equation of state can be less than -1 and cross the phantom barrier. We consider the generalized Brans-Dicke (GBD) class of scalar-tensor theories and reconstruct their Lagrangian given the effective DE density extracted from recent data. Then, given the reconstructed Lagrangian, we solve for the linear perturbations and investigate the characteristic signatures of these reconstructed GBD in the cosmological observables, such as the cosmic microwave background (CMB) anisotropy, the galaxy number counts, and their cross-correlations. In particular, we demonstrate that the Integrated Sachs-Wolfe (ISW) effect probed by the cross-correlation of CMB with the matter distribution can rule out scalar-tensor theories as the explanation of the observed DE dynamics independently from the laboratory and solar system fifth force constraints.

Constant-roll inflation was recently introduced by Motohashi, Starobinsky and Yokoyama as a phenomenological way to parametrize deviations from the slow-roll scenarios. In this paper, we investigate the dynamics of both the background and the perturbations in this model, without making any slow-roll assumptions. The perturbation spectra are computed with an efficient and accurate novel method that allowed us to quickly scan the parameter space of constant-roll inflation. We derive the constraints on the model parameters from the cosmic microwave background anisotropy measurements provided by the joint analysis of the Planck Collaboration and the BICEP2/Keck Array data.

We present a major upgrade of MGCAMB, a patch for the Einstein-Boltzmann solver CAMB used for phenomenological tests of general relativity against cosmological datasets. This new version is compatible with the latest CosmoMC code and includes a consistent implementation of massive neutrinos and dynamical dark energy. The code has been restructured to make it easier to modify to fit the custom needs of specific users. To illustrate the capabilities of the code, we present joint constraints on the modified growth, massive neutrinos and the dark energy equation of state from the latest cosmological observations, including the recent galaxy counts and weak lensing measurements from the Dark Energy Survey, and find a good consistency with the {\Lambda}CDM model.

We reconstruct evolution of the dark energy (DE) density using a nonparametric Bayesian approach from a combination of latest observational data. We caution against parameterizing DE in terms of its equation of state as it can be singular in modified gravity models, and using it introduces a bias preventing negative effective DE densities. We find a $3.7\sigma$ preference for an evolving effective DE density with interesting features. For example, it oscillates around the $\Lambda$CDM prediction at $z\lesssim0.7$, and could be negative at $z\gtrsim2.3$; dark energy can be pressure-less at multiple redshifts, and a short period of cosmic deceleration is allowed at $0.1 \lesssim z\lesssim 0.2$. We perform the reconstruction for several choices of the prior, as well as a evidence-weighted reconstruction. We find that some of the dynamical features, such as the oscillatory behaviour of the DE density, are supported by the Bayesian evidence, which is a first detection of a dynamical DE with a positive Bayesian evidence. The evidence-weighted reconstruction prefers a dynamical DE at a $(2.5\pm0.06)\sigma$ significance level.

Was the primordial universe magnetized? The answer to this question would help explain the origin of micro-Gauss strength magnetic fields observed in galaxies. It is also of fundamental importance in developing a complete theory of the early universe. While there can be other signatures of cosmological magnetic fields, a signature in the cosmic microwave background (CMB) would prove their primordial origin. The B-mode polarization of CMB is particularly promising in this regard because there are relatively few other sources of B-modes, and because the vortical modes sourced by the primordial magnetic field (PMF) survive diffusion damping up to a small fraction of the Silk length. At present, the Planck temperature and polarization spectra combined with the B-mode spectrum measured by the South Pole Telescope (SPT) constrain the PMF strength to be no more than $\sim 1$ nano-Gauss (nG). Because of the quartic scaling of the CMB anisotropy spectra with the PMF strength, this bound will not change by much even with the significantly better measurements of the B-mode spectrum by the Stage III and Stage IV CMB experiments. On the other hand, tightening the bound well below the $1$ nG threshold would rule out the purely primordial origin (requiring no dynamo action) of galactic fields. Considering Faraday rotation, which converts some of the E-modes into B-modes and scales linearly with the field strength, will help to achieve this goal. As we demonstrate, the upcoming experiments, such as SPT-3G and the Simons Observatory, will be sensitive to fields of $\sim 0.5$ nG strength thanks to the mode-coupling signature induced by Faraday rotation. A future Stage IV ground based experiment or a Space Probe will be capable of probing fields below $0.1$ nG, and would detect a scale-invariant PMF of $0.2$ nG strength without de-lensing or subtracting the galactic rotation measure.

A primordial magnetic field (PMF) present before recombination can leave specific signatures on the cosmic microwave background (CMB) fluctuations. Of particular importance is its contribution to the B-mode polarization power spectrum. Indeed, vortical modes sourced by the PMF can dominate the B-mode power spectrum on small scales, as they survive damping up to a small fraction of the Silk length. Therefore, measurements of the B-mode polarization at high-$\ell$ , such as the one recently performed by the South Pole Telescope (SPT), have the potential to provide stringent constraints on the PMF. We use the publicly released SPT B-mode polarization spectrum, along with the temperature and polarization data from the Planck satellite, to derive constraints on the magnitude, the spectral index and the energy scale at which the PMF was generated. We find that, while Planck data constrains the magnetic amplitude to $B_{1 \, \text{Mpc}} < 3.3$ nG at 95\% confidence level (CL), the SPT measurement improves the constraint to $B_{1 \, \text{Mpc}} < 1.5$ nG. The magnetic spectral index, $n_B$, and the time of the generation of the PMF are unconstrained. For a nearly scale-invariant PMF, predicted by simplest inflationary magnetogenesis models, the bound from Planck+SPT is $B_{1 \, \text{Mpc}} < 1.2$ nG at 95% CL. For PMF with $n_B=2$, expected for fields generated in post-inflationary phase transitions, the 95% CL bound is $B_{1 \, \text{Mpc}} < 0.002$ nG, corresponding to the magnetic fraction of the radiation density $\Omega_{B\gamma} < 10^{-3}$ or the effective field $B_{\rm eff} < 100$ nG. The patches for the Boltzmann code CAMB and the Markov Chain Monte Carlo engine CosmoMC, incorporating the PMF effects on CMB, are made publicly available.

Modified gravity theories often contain a scalar field of gravitational strength which interacts with matter. We examine constraints on the range and the coupling strength of a scalar gravitational degree of freedom using a subset of current data that can be safely analyzed within the linear perturbation theory. Using a model-independent implementation of scalar-tensor theories in MGCAMB in terms of two functions of the scale factor describing the mass and the coupling of the scalar degree of freedom, we derive constraints on the $f(R)$, generalized chameleon, Symmetron and Dilaton models. Since most of the large scale structure data available today is from relatively low redshifts, only a limited range of observed scales is in the linear regime, leading to relatively weak constraints. We then perform a forecast for a future large scale structure survey, such as Large Synoptic Survey Telescope (LSST), which will map a significant volume at higher redshifts, and show that it will produce much stronger constraints on scalar interactions in specific models. We also perform a principal component analysis and find that future surveys should be able to provide tight constraints on several eigenmodes of the scalar mass evolution.