In the forthcoming large volume galaxy surveys higher order statistics will provide complementary information to the usual two point statistics. Low variance estimators of the Three Point Correlation Function (3CPF) of discrete data count triangle configurations with vertices mixing data and random catalogues. Large density random catalogues are used to reduce the shot noise, which leads to a computational cost of one or two orders of magnitude more than the pure data histogram. In this paper, we explore time reductions of the isotropic 3PCF random sampling terms in periodic boxes without using random catalogues. In the first approach, based on Hamilton's construction of his famous two point estimator, we use an ad-hoc two point correlation term, while for the second procedure we construct the operators from a geometrical viewpoint, using two sides and their opening angle to describe the 3PCF triangle configurations. We map the last result to the three triangle side basis either numerically or analytically, and show that the latter approach performs best when applied to synthetic data. Moreover, we elaborate on going beyond periodic boxes, discuss other low variance n-point estimators and present useful 3PCF visualization schemes.

It is commonly believed that a successful period of inflation driven by a single or several scalar fields requires a specific hierarchy of masses given by $M_{inf}\ll H \ll M_{heavy}$, where $M_{inf}$ can correspond to several or a single light field and $M_{heavy}\,\,\,$ corresponds to any heavy field that might be integrated out if it satisfies suitable conditions. This is at the heart of the so called $\eta$-problem in inflation, since large contributions to the masses of the inflatons might spoil the slow-roll conditions required for inflation. We show that, while this is an unavoidable conclusion in single field inflation, in multifield inflation, heavy fields as defined above, may be fully responsible for a successful period of what we call fat slow-roll inflation. Moreover we show that in this scenario, the turning rate of the inflationary trajectory, $\Omega/H$, is larger than one. Thus, the $\eta$-problem is evaded with large turns in fat inflation. Depending on the perturbations' mass spectra, cosmological predictions will differ either slightly or largely with respect to those of the single field case. We illustrate this scenario in a concrete example in Type IIB string flux compactifications, where a probe D5-brane moving along the angular and radial directions in a warped throat drives fat D5-brane natural inflation. An instantaneous superplanckian decay constant can be defined, consistent with our low energy approximations, thanks to the strong warping of the geometry. We compute the cosmological observables, which differ from the single field case. We also discuss fat inflation in the context of recently proposed swampland de Sitter conjectures.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

Starting from the Gleyzes-Langlois-Piazza-Vernizzi action, we derive the most general effective theory that is invariant under internal shifts and a $\mathbb{Z}_2$ mirror symmetry in the scalar sector. Contrary to what one may think, this model presents a dark matter tracker previous to the dark energy domination. We show that, in an empty universe and to linear order in perturbations, the scalar mode clusters in exactly the same way as standard nonrelativistic cold dark matter. This also holds for the subsector of the theory where the speed of propagation of gravitational waves equals that of light, in agreement with the recent multimessenger observation. However, the inclusion of standard model particles introduces nontrivial couplings of the gravitational scalar mode to baryons, modifying their clustering properties. We argue that no arrangement of the parameters of the model can reduce the extra scalar to precisely behave as cold dark matter.

In a recent paper, Erik Verlinde has developed the interesting possibility that spacetime and gravity may emerge from the entangled structure of an underlying microscopic theory. In this picture, dark matter arises as a response to the standard model of particle physics from the delocalized degrees of freedom that build up the dark energy component of the Universe. Dark matter physics is then regulated by a characteristic acceleration scale $a_0$, identified with the radius of the (quasi)-de Sitter universe we inhabit. For a point particle matter source, or outside an extended spherically symmetric object, MOND's empirical fitting formula is recovered. However, Verlinde's theory critically departs from MOND when considering the inner structure of galaxies, differing by a factor of 2 at the centre of a regular massive body. For illustration, we use the eight classical dwarf spheroidal satellites of the Milky Way. These objects are perfect testbeds for the model given their approximate spherical symmetry, measured kinematics, and identified missing mass. We show that, without the assumption of a maximal deformation, Verlinde's theory can fit the velocity dispersion profile in dwarf spheroidals with no further need of an extra dark particle component. If a maximal deformation is considered, the theory leads to mass-to-light ratios that are marginally larger than expected from stellar population and formation history studies. We also compare our results with the recent phenomenological interpolating MOND function of McGaugh {\it et al}, and find a departure that, for these galaxies, is consistent with the scatter in current observations.

The direct detection of gravitational waves opens new perspectives for measuring properties of gravitationally bound compact objects. It is then important to investigate black holes and neutron stars in alternative theories of gravity, since they can have features that make them observationally distinguishable from their General Relativity (GR) counterparts. In this work, we examine a special case of vector Galileons, a vector-tensor theory of gravity with interesting cosmological properties, which consists of a one parameter modification of the Einstein-Maxwell action. Within this theory, we study configurations describing asymptotically flat, spherically symmetric black holes and neutron stars. The set of black hole solutions in this theory is surprisingly rich, generalising results found in GR or in related scalar-tensor theories. We investigate the properties and conserved charges of black holes, using both analytical and numerical techniques, highlighting configurations that are more compact than in GR. We then study properties of neutron stars, showing how the vector profile can influence the star internal structure. Depending on properties of matter and fields inside the star, neutron stars can be more massive than in GR, and they can be even more compact than Schwarzschild black holes, making these objects observationally interesting. We also comment on possible extensions of our configurations to magnetically charged or rotating configurations.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

In the context of a cubic Galileon model in which the Vainshtein mechanism suppresses the scalar field interactions with matter, we study low-density stars with slow rotation and static relativistic stars. We develop an expansion scheme to find approximated solutions inside the Vainshtein radius, and show that deviations from General Relativity (GR), while considering rotation, are also suppressed by the Vainshtein mechanism. In a quadratic coupling model, in which the scalarisation effect can significantly enhance deviations from GR in normal scalar tensor gravity, the Galileon term successfully suppresses the large deviations away from GR. Moreover, using a realistic equation of state, we construct solutions for a relativistic star, and show that deviations from GR are more suppressed for higher density objects. However, we found that the scalar field solution ceases to exist above a critical density, which roughly corresponds to the maximum mass of a neutron star. This indicates that, for a compact object described by a polytropic equation of state, the configuration that would collapse into a black hole cannot support a non-trivial scalar field.

We study static, spherically symmetric solutions in a recently proposed ghost-free model of non-linear massive gravity. We focus on a branch of solutions where the helicity-0 mode can be strongly coupled within certain radial regions, giving rise to the Vainshtein effect. We truncate the analysis to scales below the gravitational Compton wavelength, and consider the weak field limit for the gravitational potentials, while keeping all non-linearities of the helicity-0 mode. We determine analytically the number and properties of local solutions which exist asymptotically on large scales, and of local (inner) solutions which exist on small scales. We find two kinds of asymptotic solutions, one of which is asymptotically flat, while the other one is not, and also two types of inner solutions, one of which displays the Vainshtein mechanism, while the other exhibits a self-shielding behaviour of the gravitational field. We analyse in detail in which cases the solutions match in an intermediate region. The asymptotically flat solutions connect only to inner configurations displaying the Vainshtein mechanism, while the non asymptotically flat solutions can connect with both kinds of inner solutions. We show furthermore that there are some regions in the parameter space where global solutions do not exist, and characterise precisely in which regions of the phase space the Vainshtein mechanism takes place.

We study linear perturbations around time dependent spherically symmetric solutions in the Lambda_3 massive gravity theory, which self-accelerate in the vacuum. We find that the dynamics of the scalar perturbations depend on the coordinate choice for the background solutions. For particular choices of coordinates there is a symmetry enhancement, leaving no propagating scalar degrees of freedom at linear order in perturbations. In contrast, any other coordinate choice propagates a single scalar mode. We find that the Hamiltonian of this scalar mode is unbounded from below for all self-accelerating solutions, signalling an instability.

Starting from the general Horndeski action, we derive the most general effective theory for scalar perturbations around flat space that allows us to screen fifth forces via the Vainshtein mechanism. The effective theory is described by a generalization of the Galileon Lagrangian, which we use to study the stability of spherically symmetric configurations exhibiting the Vainshtein effect. In particular, we discuss the phenomenological consequences of a scalar-tensor coupling that is absent in the standard Galileon Lagrangian. This coupling controls the superluminality and stability of fluctuations inside the Vainshtein radius in a way that depends on the density profile of a matter source. Particularly we find that the vacuum solution is unstable due to this coupling.

Massive gravity is a good theoretical laboratory to study modifications of General Relativity. The theory offers a concrete set-up to study models of dark energy, since it admits cosmological self-accelerating solutions in the vacuum, in which the size of the acceleration depends on the graviton mass. Moreover, non-linear gravitational self-interactions, in the proximity of a matter source, manage to mimic the predictions of linearised General Relativity, hence agreeing with solar-system precision measurements. In this article, we review our work in the subject, classifying, on one hand, static solutions, and on the other hand, self-accelerating backgrounds. For what respects static solutions we exhibit black hole configurations, together with other solutions that recover General Relativity near a source via the Vainshtein mechanism. For the self-accelerating solutions we describe a wide class of cosmological backgrounds, including an analysis of their stability.

We investigate vector contributions to the Lagrangian of $\Lambda_3-$massive gravity in the decoupling limit, the less explored sector of this theory. The main purpose is to understand the stability of maximally symmetric %self-accelerating vacuum solutions. Around self-accelerating configurations, vector degrees of freedom become strongly coupled since their kinetic terms vanish, so their dynamics is controlled by higher order interactions. Even in the decoupling limit, the vector Lagrangian contains an infinite number of terms. We develop a systematic method to covariantly determine the vector Lagrangian at each order in perturbations, fully manifesting the symmetries of the system. We show that, around self-accelerating solutions, the structure of higher order $p$-form Galileons arise, avoiding the emergence of a sixth BD ghost mode. However, a careful analysis shows that there are directions along which the Hamiltonian is unbounded from below. This instability can be interpreted as one of the available fifth physical modes behaving as a ghost. Therefore, we conclude that self-accelerating configurations, in the decoupling limit of $\Lambda_3$-massive gravity, are generically unstable.

We expose a new symmetry for linear perturbations around a solution of non-linear Fierz-Pauli massive gravity plus a bare cosmological constant. The cosmological constant is chosen such that the background metric is flat while the Stuckelberg fields have a non-trivial profile. Around this background, at linear order the new symmetry reduces the propagating degrees of freedom to those of General Relativity, namely the massless helicity 2 modes only. We discuss the physical consequences and possible applications of these findings.

We investigate strong coupling effects in a covariant massive gravity model, which is a candidate for a ghost free non-linear completion of Fierz-Pauli. We analyse the conditions to recover general relativity via Vainshtein mechanism in the weak field limit, and find three main cases depending on the choice of parameters. In the first case, the potential is such that all non-linearities disappear and the vDVZ discontinuity cannot be avoided. In the second case, the Vainshtein mechanism allows to recover general relativity within a macroscopic radius from a source. In the last case, the strong coupling of the scalar graviton completely shields the massless graviton, and weakens gravity when approaching the source. In the second part of the paper, we explore new exact vacuum solutions, that asymptote de Sitter or anti de Sitter space depending on the parameter choice. The curvature of the space is proportional to the mass of the graviton, thus providing a cosmological background which may explain the present day acceleration in terms of the graviton mass. Moreover, by expressing the potential for non-linear massive gravity in a convenient form, we also suggest possible connections with a higher dimensional framework.

We explore the possibility of realising self-accelerated expansion of the Universe taking into account the vector components of a massive graviton. The effective action in the decoupling limit contains an infinite number of terms, once the vector degrees of freedom are included. These can be re-summed in physically interesting situations, which result in non-polynomial couplings between the scalar and vector modes. We show there are self-accelerating background solutions for this effective action, with the possibility of having a non-trivial profile for the vector fields. We then study fluctuations around these solutions and show that there is always a ghost, if a background vector field is present. When the background vector field is switched off, the ghost can be avoided, at the price of entering into a strong coupling regime, in which the vector fluctuations have vanishing kinetic terms. Finally we show that the inclusion of a bare cosmological constant does not change the previous conclusions and it does not lead to a ghost mode in the absence of a background vector field.

We study spherically symmetric solutions in a covariant massive gravity model, which is a candidate for a ghost-free non-linear completion of the Fierz-Pauli theory. There is a branch of solutions that exhibits the Vainshtein mechanism, recovering General Relativity below a Vainshtein radius given by $(r_g m^2)^{1/3}$, where $m$ is the graviton mass and $r_g$ is the Schwarzschild radius of a matter source. Another branch of exact solutions exists, corresponding to Schwarzschild-de Sitter spacetimes where the curvature scale of de Sitter space is proportional to the mass squared of the graviton.

We consider a magnetic Bianchi I braneworld, embedded in between two Schwarzschild-AdS spacetimes, boosted equal amounts in opposite directions and compare them to the analagous solution in four-dimensional General Relativity. The efficient dissipation of anisotropy on the brane is explicitly demonstrated, a process we dub braneworld isotropization. From the bulk point of view, we attribute this to anisotropic energy being carried into the bulk by hot gravitons leaving the brane. From the brane point of view this can be interpreted in terms of the production of particles in the dual CFT. We explain how this result enables us to gain a better understanding of the behaviour of anisotropic branes already studied in the literature. We also show how there is evidence of particles being over-produced, and comment on how this may ultimately provide a possible observational signature of braneworlds.