We study the cosmology of a specific class of nonlocal model of modified gravity, the so-called Deser-Woodard (DW) model, modifying the Einstein-Hilbert action by a term $\sim R f(\Box^{-1}R)$, where $f$ is a free function. Choosing $f$ so as to reproduce the $\Lambda{\rm CDM}$ cosmological background expansion history within the nonlocal model, we implement the model in a cosmological linear Einstein--Boltzmann solver and study the deviations to GR the model induces in the scalar and tensor perturbations. We observe that the DW nonlocal model describes a modified propagation for the gravitational waves, as well as a lower linear growth rate and a stronger lensing power as compared to $\Lambda{\rm CDM}$, up to several percents. Such prominent growth and lensing features lead to the inference of a significantly smaller value of $\sigma_8$ with respect to the one in $\Lambda{\rm CDM}$, given \textit{Planck} CMB+lensing data. The prediction for the linear growth rate $f \sigma_8$ within the DW model is therefore significantly smaller than the one in $\Lambda{\rm CDM}$ and the addition of growth rate data $f \sigma_8$ from Redshift-space distortion measurements to \textit{Planck} CMB+lensing, opens a (dominant) tension between Redshift-space distortion data and the reconstructed \textit{Planck} CMB lensing potential. However, model selection issues only result in "weak" evidences for $\Lambda{\rm CDM}$ against the DW model given the data. Such a fact shows that the datasets we consider are not constraining enough for distinguishing between the models. As we discuss, the addition of galaxy WL data or cosmological constraints from future galaxy clustering, weak lensing surveys, but also third generation gravitational wave interferometers, prove to be useful for discriminating modified gravity models such as the DW one from $\Lambda{\rm CDM}$, within the close future.

The recent report that the gravitational wave speed equals the light speed puts strong constraints on the anisotropic stress parameter of many modified gravity models, a quantity that is directly observable through large-scale structure. We show here that models without a mass scale completely escape these constraints. We discuss a few relevant cases in detail: Brans-Dicke theory, nonlocal models, and Galileon Lagrangian.

In this work, we extend previous analyses of the structure formation in the $f(\Box^{-1}R)$ model of nonlocal gravity proposed by Deser and Woodard (DW), which reproduces the background expansion of $\Lambda$CDM with no need of a cosmological constant nor of any dimensional constant beside Newton's one. A previous analysis based on redshift-space distortions (RSD) data concluded that the model was ruled out. In this work we revisit the issue and find that, when recast in a localized model, the DW model is not ruled out and actually gives a better fit to RSD data than $\Lambda$CDM. At the same time, the model predicts a slightly lower value of $\sigma_{8}$ than $\Lambda$CDM, in agreement with recent estimates based on lensing. We also produce analytical approximations of the two modified gravity functions and of $f\sigma_{8}(z)$ as a function of redshift. Finally, we also show how much the fit depends on initial conditions when these are generalized with respect to a standard matter-dominated era.

We discuss the cosmological implications of nonlocal modifications of general relativity containing tensorial structures. Assuming the presence of standard radiation- and matter-dominated eras, we show that, except in very particular cases, the nonlocal terms contribute a rapidly growing energy density. These models therefore generically do not have a stable cosmological evolution.

Many theories of modified gravity with higher order derivatives are usually ignored because of serious problems that appear due to an additional ghost degree of freedom. Most dangerously, it causes an immediate decay of the vacuum. However, breaking Lorentz invariance can cure such abominable behavior. By analyzing a model that describes a massive graviton together with a remaining Boulware-Deser ghost mode we show that even ghostly theories of modified gravity can yield models that are viable at both classical and quantum levels and, therefore, they should not generally be ruled out. Furthermore, we identify the most dangerous quantum scattering process that has the main impact on the decay time and find differences to simple theories that only describe an ordinary scalar field and a ghost. Additionally, constraints on the parameters of the theory including some upper bounds on the Lorentz-breaking cutoff scale are presented. In particular, for a simple theory of massive gravity we find that a breaking of Lorentz invariance is allowed to happen even at scales above the Planck mass. Finally, we discuss the relevance to other theories of modified gravity.

We discuss the cosmological implications of the $R~\Box^{-2}R$ nonlocal modification to standard gravity. We relax the assumption of special initial conditions in the local formulation of the theory, perform a full phase-space analysis of the system, and show that the late-time cosmology of the model exhibits two distinct evolution paths, on which a large range of values for the present equation of state can be reached. We then compare the general solutions to supernovae data and place constraints on the parameters of the model. In particular, we find that the mass parameter of the theory should be smaller than 1.2 in Hubble units.

Massive bigravity models are interesting alternatives to standard cosmology. In most cases, however, these models have been studied for a simplified scenario in which both metrics take homogeneous and isotropic forms [Friedmann-Lema\^itre-Robertson-Walker (FLRW)] with the same spatial curvatures. The interest to consider more general geometries arises, in particular, in view of the difficulty so far encountered in building stable cosmological solutions with homogeneous and isotropic metrics. Here we consider a number of cases in which the two metrics take more general forms, namely FLRW with different spatial curvatures---Lema\^itre, Lema\^itre-Tolman-Bondi (LTB), and Bianchi I---as well as cases where only one metric is linearly perturbed. We discuss possible consistent combinations and find that only some special cases of FLRW--Lema\^itre, LTB--LTB, and FLRW--Bianchi I combinations give consistent, nontrivial solutions.

Two-point functions, the mean field squared and the vacuum expectation value (VEV) of the energy-momentum tensor are investigated for the electromagnetic field in the geometry of parallel plates on background of $(D+1)$% -dimensional dS spacetime. We assume that the field is prepared in the Bunch-Davies vacuum state and on the plates a boundary condition is imposed that is a generalization of the perfectly conducting boundary condition for an arbitrary number of spatial dimensions. It is shown that for $D\geq 4$ the background gravitational field essentially changes the behavior of the VEVs at separations between the plates larger than the curvature radius of dS spacetime. At large separations, the Casimir forces are proportional to the inverse fourth power of the distance for all values of spatial dimension $D\geq 3$. For $D\geq 4$ this behavior is in sharp contrast with the case of plates in Minkowski bulk where the force decays as the inverse $(D+1)$th power of the distance.

We evaluate the two-point function for the electromagnetic field tensor in $(D+1)$-dimensional de Sitter spacetime assuming that the field is prepared in Bunch-Davies vacuum state. This two-point function is used for the investigation of the vacuum expectation values (VEVs) of the field squared and the energy-momentum tensor in the presence of a conducting plate. The VEVs are decomposed into the boundary-free and plate-induced parts. For the latter, closed form analytical expressions are given in terms of the hypergeometric function. For $3\leqslant D\leqslant 8$ the plate-induced part in the VEV of the electric field squared is positive everywhere, whereas for $D\geqslant 9$ it is positive near the plate and negative at large distances. The VEV of the energy-momentum tensor, in addition to the diagonal components, contains an off-diagonal component which corresponds to the energy flux along the direction normal to the plate. Simple asymptotic expressions are provided at small and large distances from the plate compared with the de Sitter curvature scale. For $D\geqslant 4$, all the diagonal components of the plate-induced vacuum energy-momentum tensor are negative and the energy flux is directed from the plate.

Vacuum expectation values (VEVs) of the current densities for charged scalar and Dirac spinor fields are investigated in (D+1)-dimensional de Sitter (dS) spacetime with toroidally compactified spatial dimensions. Along compact dimensions we impose quasiperiodicity conditions with arbitrary phases. In addition, the presence of a classical constant gauge field is assumed. The VEVs of the charge density and of the components for the current density along noncompact dimensions vanish. The gauge field leads to Aharonov-Bohm-like oscillations of the components along compact dimensions as functions of the magnetic flux. For small values of the comoving length of a compact dimension, compared with the dS curvature scale, the current density is related to the corresponding current in the Minkowski spacetime by a conformal relation. For large values of the comoving length and for a scalar field, depending on the mass of the field, two different regimes are realized with the monotonic and oscillatory damping of the current density. For a massive spinor field, the decay of the current density is always oscillatory. In supersymmetric models on the background of Minkowski spacetime with equal number of scalar and fermionic degrees of freedom and with the same phases in the periodicity conditions, the total current density vanishes due to the cancellation between the scalar and fermionic parts. The background gravitational field modifies the current densities for scalar and fermionic fields in different ways and for massive fields there is no cancellation in the dS spacetime.

The scattering and transformation of the waves propagating in magnetized plasma on a heavy stationary charged particle located at a plane plasma-vacuum boundary is considered. The scattering (transformation) occurs due to the nonlinear coupling of the incident wave with the polarization (shielding) cloud surrounding the particle. It is shown that the problem is reduced to the determination of the nonlinear (three index) dielectric tensor of magnetized plasma. The angular distribution and the cross section for scattering (transformation) of high-frequency ordinary and extraordinary waves and low-frequency upper-hybrid, low-hybrid, and magnetosonic waves are investigated within a cold plasma (hydrodynamic) model.