Ultra-light Dark Matter (ULDM) models are suitable candidates for the cosmological Dark Matter that may leave characteristic imprints in many observables. Among other probes, signatures of ULDM can be searched for in pulsar timing data. In this work we describe the effects of spin-2 ULDM on pulsar timing arrays, extending previous results on lower spins. Spin-2 ULDM is universally coupled to standard matter with dimensionless strength $\alpha$. We estimate that current data could constrain this coupling in the mass range $m\lesssim4\times10^{-22}$ eV at the $10^{-5}$ to $10^{-6}$ level, which is the most competitive constraint in this mass range. A crucial feature of the spin-2 ULDM effect on pulsar timing is its anisotropic, quadrupolar shape. This feature can be instrumental in differentiating the effects sourced by spin-2 ULDM from, for instance, scalar ULDM, and the systematics of a PTA experiment.

Dark matter (DM) can consist of very light bosons behaving as a classical scalar field that experiences coherent oscillations. The presence of this DM field would perturb the dynamics of celestial bodies, either because the (oscillating) DM stress tensor modifies the gravitational potentials of the galaxy, or if DM is directly coupled to the constituents of the body. We study secular variations of the orbital parameters of binary systems induced by such perturbations. Two classes of effects are identified. Effects of the first class appear if the frequency of DM oscillations is in resonance with the orbital motion; these exist for general DM couplings including the case of purely gravitational interaction. Effects of the second class arise if DM is coupled quadratically to the masses of the binary system members and do not require any resonant condition. The exquisite precision of binary pulsar timing can be used to constrain these effects. Current observations are not sensitive to oscillations in the galactic gravitational field, though a discovery of pulsars in regions of high DM density may improve the situation. For DM with direct coupling to ordinary matter, the current timing data are already competitive with other existing constraints in the range of DM masses $\sim 10^{-22}-10^{-18}\,{\rm eV}$. Future observations are expected to increase the sensitivity and probe new regions of parameters.

We provide a general framework for self-interacting warm dark matter (WDM) in cosmological perturbations, by deriving from first principles a Boltzmann hierarchy which retains certain independence from a particular interaction Lagrangian. We consider elastic interactions among the massive particles, and obtain a hierarchy which is more general than the ones usually obtained for non-relativistic (as for cold DM) or for ultra-relativistic (as for neutrinos) approximations. The more general momentum-dependent kernel integrals in the Boltzmann collision terms, are explicitly calculated for different field-mediator models, including examples of a scalar field (either massive or massless) or a massive vector field. As an application, we study the evolution of the interaction rate per particle under the relaxation time approximation, and assess when a given self-interaction is relevant in comparison with the Hubble expansion rate. Our framework aims to be a useful starting point to evaluate DM self-interaction effects in the linear power spectrum, necessary to then study its evolution all the way to non-linear stages of structure formation, where certain DM interactions were proven to be relevant.

Binary pulsars can be excellent probes of ultra-light dark matter. We consider the scenario where the latter is represented by a spin-2 field. The coherent oscillations of the dark matter field perturb the dynamics of binary systems, leading to secular effects for masses that resonate with the binary systems. For the range $10^{-23}$eV $\lesssim m \lesssim 10^{-17}$eV we show that current timing data could potentially constrain the universal coupling strength of dark matter to ordinary matter at the level of $\alpha \simeq 10^{-5}$.

We study the secular effects that an oscillating background ultralight (fuzzy) cosmological vector field has on the dynamics of binary systems; such effects appear when the field and the binary are in resonance. We first consider the gravitational interaction between the field and the systems, and quantify the main differences with an oscillating background scalar field. If the energy density of such a field is sufficiently large, as required if it is supposed to be all of the dark matter, we show that the secular effects could yield potentially observable signatures in high precision time of arrival measurements of binary pulsars. We then analyse the secular effects that arise when the field is directly coupled to the bodies in the binary. We show that this study is particularly relevant for models where fuzzy dark matter mediates a baryonic force $B$ (or $B-L$, with $L$ the lepton number), due to the stellar amount of nucleons present in the stars. The constraints we obtain from current data are already competitive with (or even more constraining than) laboratory tests of the equivalence principle.

We consider the scenario where dark matter (DM) is represented by an ultra-light classical scalar field performing coherent periodic oscillations. We point out that such DM perturbs the dynamics of binary systems either through its gravitational field or via direct coupling to ordinary matter. This perturbation gets resonantly amplified if the frequency of DM oscillations is close to a (half-)integer multiple of the orbital frequency of the system and leads to a secular variation of the orbital period. We suggest to use binary pulsars as probes of this scenario and estimate their sensitivity. While the current accuracy of observations is not yet sufficient to probe the purely gravitational effect of DM, it already yields constraints on direct coupling that are competitive with other bounds. The sensitivity will increase with the upcoming radio observatories such as Square Kilometer Array.

We analyze the Extended Quasi-Dilaton Massive Gravity model around a Friedmann-Lema\^itre-Robertson-Walker cosmological background. We present a careful stability analysis of asymptotic fixed points. We find that the traditional fixed point cannot be approached dynamically, except from a perfectly fine-tuned initial condition involving both the quasi-dilaton and the Hubble parameter. A less-well examined fixed-point solution, where the time derivative of the 0-th St\"uckelberg field vanishes $\dot\phi^0=0$, encounters no such difficulty, and the fixed point is an attractor in some finite region of initial conditions. We examine the question of the presence of a Boulware-Deser ghost in the theory. We show that the additional constraint which generically allows for the elimination of the Boulware-Deser mode is only present under special initial conditions. We find that the only possibility corresponds to the traditional fixed point and the initial conditions are the same fine-tuned conditions that allow the fixed point to be approached dynamically.

We investigate the possible implications of the measured value of the scalar tilt $n_s$ for the tensor-to-scalar ratio $r$ in slow-roll, single-field inflationary models. The measured value of the tilt satisfies $n_s -1\sim 1/N_*$, where $N_* \sim 60$ is the number of $e$-folds for observationally relevant scales. If this is not a coincidence and the scaling holds for different values of $N$, it strongly suggests that either $r$ is as big as $10^{-1}$ (a possibility in tension with the recent data), or smaller than $10^{-2}$ and exponentially dependent on $n_s$. A large region of the ($n_s$,$r$) plane is not compatible with this scaling.

We update the forecasts for the measurement of the tensor-to-scalar ratio $r$ for various ground-based experiments (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), balloons (EBEX 10k and Spider) and satellites (CMBPol, COrE and LiteBIRD), taking into account the recent Planck data on polarized dust and using a component separation method. The forecasts do not change significantly with respect to previous estimates when at least three frequencies are available, provided foregrounds can be accurately described by few parameters. We argue that a theoretically motivated goal for future experiments is $r\sim2\times10^{-3}$, and that this is achievable if the noise is reduced to $\sim1\,\mu$K-arcmin and lensing is reduced to $10\%$ in power. We study the constraints experiments will be able to put on the frequency and $\ell$-dependence of the tensor signal as a check of its primordial origin. Futuristic ground-based and balloon experiments can have good constraints on these parameters, even for $r\sim2\times10^{-3}$. For the same value of $r$, satellites will marginally be able to detect the presence of the recombination bump, the most distinctive feature of the primordial signal.

We reanalyze the behavior of Friedmann-Lema\^itre-Robertson-Walker cosmologies in the recently proposed quasidilaton massive-gravity model, and discover that the background dynamics present hitherto unreported features that require unexpected fine-tuning of the additional fundamental parameters of the theory for an observationally consistent background cosmology. We also identify new allowed regions in the parameters space and exclude some of the previously considered ones. The evolution of the mass of gravitational waves reveals non-trivial behavior, exhibiting a mass squared that may be negative in the past, and that presently, while positive, is larger than the square of the Hubble parameter. These properties of the gravity-wave mass have the potential to lead to observational tests of the theory. While quasidilaton massive gravity is known to have issues with stability at short distances, the current analysis is a first step toward the investigation of the more stable extended quasidilaton massive gravity theory, with some expectation that both the fine-tuning of parameters and the interesting behavior of the gravity-wave mass will persist.

The simplest inflationary model $V=\frac12 m^2\phi^2$ represents the benchmark for future constraints. For a quadratic potential, the quantity $(n_s-1)+r/4+11 (n_s-1)^2/24$ vanishes (up to corrections which are cubic in slow roll) and can be used to parametrize small deviations from the minimal scenario. Future constraints on this quantity will be able to distinguish a quadratic potential from a pseudo-Nambu-Goldstone boson with $f \lesssim 30 M_{\rm pl}$ and set limits on the deviation from unity of the speed of sound $| c_s-1| \lesssim 3\times 10^{-2}$ (corresponding to an energy scale $\Lambda\gtrsim 2\times 10^{16}\, \mathrm{GeV}$), and on the contribution of a second field to perturbations ($\lesssim 6 \times 10^{-2}$). The limiting factor for these bounds will be the uncertainty on the spectral index. The error on the number of e-folds will be $\Delta N \simeq 0.4$, corresponding to an error on the reheating temperature $\Delta T_\mathrm{rh}/T_\mathrm{rh}\simeq 1.2$. We comment on the relevance of non-Gaussianity after BICEP2 results.

In the simplest inflationary model $V=\frac12 m^2\phi^2$, we provide a prediction accurate up to $1\%$ for the spectral index $n_s$ and the tensor-to-scalar ratio $r$ assuming instantaneous reheating and a standard thermal history: $n_s = 0.9668\pm0.0003$ and $r=0.131\pm 0.001$. This represents the simplest and most informative point in the $(n_s,r)$ plane. The result is independent of the details of reheating (or preheating) provided the conversion to radiation is sufficiently fast. A slower reheating or a modified post-inflationary evolution push towards smaller $n_s$ (and larger $r$), so that our prediction corresponds to the maximum $n_s$ (and minimum $r$) for the quadratic potential. We also derive similar results for a general $V \propto \phi^p$ potential.

We extend the resummation method of Anselmi & Pietroni (2012) to compute the total density power spectrum in models of quintessence characterized by a vanishing speed of sound. For standard $\Lambda$CDM cosmologies, this resummation scheme allows predictions with an accuracy at the few percent level beyond the range of scales where acoustic oscillations are present, therefore comparable to other, common numerical tools. In addition, our theoretical approach indicates an approximate but valuable and simple relation between the power spectra for standard quintessence models and models where scalar field perturbations appear at all scales. This, in turn, provides an educated guess for the prediction of nonlinear growth in models with generic speed of sound, particularly valuable since no numerical results are yet available.

We generalize the consistency condition for the three-point function in single field inflation to the case of dissipative, multi-field, single-clock models. We use the recently introduced extension of the effective field theory of inflation that accounts for dissipative effects, to provide an explicit proof to leading (non-trivial) order in the generalized slow roll parameters and mixing with gravity scales. Our results illustrate the conditions necessary for the validity of the consistency relation in situations with many degrees of freedom relevant during inflation, namely that there is a preferred clock. Departures from this condition in forthcoming experiments would rule out not only single field but also a large class of multi-field models.

We generalize the effective field theory of single clock inflation to include dissipative effects. Working in unitary gauge we couple a set of composite operators in the effective action which is constrained solely by invariance under time-dependent spatial diffeomorphisms. We restrict ourselves to situations where the degrees of freedom responsible for dissipation do no contribute to the density perturbations at late time. The dynamics of the perturbations is then modified by the appearance of `friction' and noise terms, and assuming certain locality properties for the Green's functions of these composite operators, we show that there is a regime characterized by a large friction term \gamma >> H in which the \zeta-correlators are dominated by the noise and the power spectrum can be significantly enhanced. We also compute the three point function <\zeta\zeta\zeta> for a wide class of models and discuss under which circumstances large friction leads to an increased level of non-Gaussianities. In particular, under our assumptions, we show that strong dissipation together with the required non-linear realization of the symmetries implies |f_NL| ~ \gamma/(c_s^2H) >> 1. As a paradigmatic example we work out a variation of the `trapped inflation' scenario with local response functions and perform the matching with our effective theory. A detection of the generic type of signatures that result from incorporating dissipative effects during inflation, as we describe here, would teach us about the dynamics of the early universe and also extend the parameter space of inflationary models.

We show how the quantum to classical transition of the cosmological fluctuations produced during inflation can be described by means of the influence functional and the master equation. We split the inflaton field into the system-field (long-wavelength modes), and the environment, represented by its own short-wavelength modes. We compute the decoherence times for the system-field modes and compare them with the other time scales of the model. We present the renormalized stochastic Langevin equation for an homogeneous system-field and then we analyze the influence of the environment on the power spectrum for some modes in the system.