We study perturbation theory for large-scale structure in the most general scalar-tensor theories propagating a single scalar degree of freedom, which include Horndeski theories and beyond. We model the parameter space using the effective field theory of dark energy. For Horndeski theories, the gravitational field and fluid equations are invariant under a combination of time-dependent transformations of the coordinates and fields. This symmetry allows one to construct a physical adiabatic mode which fixes the perturbation-theory kernels in the squeezed limit and ensures that the well-known consistency relations for large-scale structure, originally derived in general relativity, hold in modified gravity as well. For theories beyond Horndeski, instead, one generally cannot construct such an adiabatic mode. Because of this, the perturbation-theory kernels are modified in the squeezed limit and the consistency relations for large-scale structure do not hold. We show, however, that the modification of the squeezed limit depends only on the linear theory. We investigate the observational consequences of this violation by computing the matter bispectrum. In the squeezed limit, the largest effect is expected when considering the cross-correlation between different tracers. Moreover, the individual contributions to the 1-loop matter power spectrum do not cancel in the infrared limit of the momentum integral, modifying the power spectrum on non-linear scales.

We develop an approach to compute observables beyond the linear regime of dark matter perturbations for general dark energy and modified gravity models. We do so by combining the Effective Field Theory of Dark Energy and Effective Field Theory of Large-Scale Structure approaches. In particular, we parametrize the linear and nonlinear effects of dark energy on dark matter clustering in terms of the Lagrangian terms introduced in a companion paper, focusing on Horndeski theories and assuming the quasi-static approximation. The Euler equation for dark matter is sourced, via the Newtonian potential, by new nonlinear vertices due to modified gravity and, as in the pure dark matter case, by the effects of short-scale physics in the form of the divergence of an effective stress tensor. The effective fluid introduces a counterterm in the solution to the matter continuity and Euler equations, which allows a controlled expansion of clustering statistics on mildly nonlinear scales. We use this setup to compute the one-loop dark-matter power spectrum.

We develop an approach to parametrize cosmological perturbations beyond linear order for general dark energy and modified gravity models characterized by a single scalar degree of freedom. We derive the full nonlinear action, focusing on Horndeski theories. In the quasi-static, non-relativistic limit, there are a total of six independent relevant operators, three of which start at nonlinear order. The new nonlinear couplings modify, beyond linear order, the generalized Poisson equation relating the Newtonian potential to the matter density contrast. We derive this equation up to cubic order in perturbations and, in a companion article, we apply it to compute the one-loop matter power spectrum. Within this approach, we also discuss the Vainshtein regime around spherical sources and the relation between the Vainshtein scale and the nonlinear scale for structure formation.

We study infrared effects in perturbation theory for large-scale structure coupled to the effective field theory of dark energy, focusing on, in particular, Degenerate Higher-Order Scalar-Tensor (DHOST) theories. In the subhorizon, Newtonian limit, DHOST theories introduce an extra large-scale velocity $v^i_\pi$ which is in general different from the matter velocity $v^i$. Contrary to the case in Horndeski theories, the presence of this extra large-scale velocity means that one cannot eliminate the long-wavelength effects of both $v^i$ and $v^i_\pi$ with a single coordinate transformation, and thus the standard $\Lambda$CDM consistency relations for large-scale structure are violated by terms proportional to the relative velocity $v^i - v^i_\pi$. We show, however, that in non-linear quantities this violation is determined by the linear equations and the symmetries of the fluid system. We find that the size of the baryon acoustic oscillations in the squeezed limit of the bispectrum is modified, that the bias expansion contains extra terms which contribute to the squeezed limit of the galaxy bispectrum, that infrared modes in the one-loop power spectrum no longer cancel, and that the equal-time double soft limit of the tree-level trispectrum is non-vanishing. In addition, we give explicit expressions for how these violations depend on the relative velocity. Many of our computations are also relevant for perturbation theory in $\Lambda$CDM with exact time dependence.

We study the screening mechanism in the most general scalar-tensor theories that leave gravitational waves unaffected and are thus compatible with recent LIGO/Virgo observations. Using the effective field theory of dark energy approach, we consider the general action for perturbations beyond linear order, focussing on the quasi-static limit. When restricting to the subclass of theories that satisfy the gravitational wave constraints, the fully nonlinear effective Lagrangian contains only three independent parameters. One of these, $\beta_1$, is uniquely present in degenerate higher-order theories. We compute the two gravitational potentials for a spherically symmetric matter source and we find that for $\beta_1 \ge 0$ they decrease as the inverse of the distance, as in standard gravity, while the case $\beta_1 < 0$ is ruled out. For $\beta_1 > 0$, the two potentials differ and their gravitational constants are not the same on the inside and outside of the body. Generically, the bound on anomalous light bending in the Solar System constrains $\beta_1 \lesssim 10^{-5}$. Standard gravity can be recovered outside the body by tuning the parameters of the model, in which case $\beta_1 \lesssim 10^{-2}$ from the Hulse-Taylor pulsar.

We study the decay of gravitational waves into dark energy fluctuations $\pi$, through the processes $\gamma \to \pi\pi$ and $\gamma \to \gamma \pi$, made possible by the spontaneous breaking of Lorentz invariance. Within the EFT of Dark Energy (or Horndeski/beyond Horndeski theories) the first process is large for the operator $\frac12 \tilde m_4^2(t) \, \delta g^{00}\, \left( {}^{(3)}\! R + \delta K_\mu^\nu \delta K^\mu_\nu -\delta K^2 \right)$, so that the recent observations force $\tilde m_4 =0$ (or equivalently $\alpha_{\rm H}=0$). This constraint, together with the requirement that gravitational waves travel at the speed of light, rules out all quartic and quintic GLPV theories. Additionally, we study how the same couplings affect the propagation of gravitons at loop order. The operator proportional to $\tilde m_4^2$ generates a calculable, non-Lorentz invariant higher-derivative correction to the graviton propagation. The modification of the dispersion relation provides a bound on $\tilde m_4^2$ comparable to the one of the decay. Conversely, operators up to cubic Horndeski do not generate sizeable higher-derivative corrections.

We develop an analytic method for implementing the IR-resummation of arXiv:1404.5954, which allows one to correctly and consistently describe the imprint of baryon acoustic oscillations (BAO) on statistical observables in large-scale structure. We show that the final IR-resummed correlation function can be computed analytically without relying on numerical integration, thus allowing for an efficient and accurate use of these predictions on real data in cosmological parameter fitting. In this work we focus on the one-loop correlation function and the BAO peak. We show that, compared with the standard numerical integration method of IR-resummation, the new method is accurate to better than 0.2 %, and is quite easily improvable. We also give an approximate resummation scheme which is based on using the linear displacements of a fixed fiducial cosmology, which when combined with the method described above, is about six times faster than the standard numerical integration. Finally, we show that this analytic method is generalizable to higher loop computations.

We compare analytical computations with numerical simulations for dark-matter clustering, in general relativity and in the normal branch of DGP gravity (nDGP). Our analytical frameword is the Effective Field Theory of Large-Scale Structure (EFTofLSS), which we use to compute the one-loop dark-matter power spectrum, including the resummation of infrared bulk displacement effects. We compare this to a set of 20 COLA simulations at redshifts $z = 0$, $z=0.5$, and $z =1$, and fit the free parameter of the EFTofLSS, called the speed of sound, in both $\Lambda$CDM and nDGP at each redshift. At one-loop at $z = 0$, the reach of the EFTofLSS is $k_{\rm reach}\approx 0.14 \, h { \rm Mpc^{-1}}$ for both $\Lambda$CDM and nDGP. Along the way, we compare two different infrared resummation schemes and two different treatments of the time dependence of the perturbative expansion, concluding that they agree to approximately $1\%$ over the scales of interest. Finally, we use the ratio of the COLA power spectra to make a precision measurement of the difference between the speeds of sound in $\Lambda$CDM and nDGP, and verify that this is proportional to the modification of the linear coupling constant of the Poisson equation.

We further develop the description of redshift space distortions within the Effective Field Theory of Large Scale Structures. First, we generalize the counterterms to include the effect of baryonic physics and primordial non-Gaussianity. Second, we evaluate the IR-resummation of the dark matter power spectrum in redshift space. This requires us to identify a controlled approximation that makes the numerical evaluation straightforward and efficient. Third, we compare the predictions of the theory at one loop with the power spectrum from numerical simulations up to $\ell=6$. We find that the IR-resummation allows us to correctly reproduce the BAO peak. The $k$-reach, or equivalently the precision for a given $k$, depends on additional counterterms that need to be matched to simulations. Since the non-linear scale for the velocity is expected to be longer than the one for the overdensity, we consider a minimal and a non-minimal set of counterterms. The quality of our numerical data makes it hard to firmly establish the performance of the theory at high wavenumbers. Within this limitation, we find that the theory at redshift $z=0.56$ and up to $\ell=2$ matches the data to percent level approximately up to $k \sim 0.13 \, h { \rm Mpc^{-1}}$ or $k \sim 0.18 \, h { \rm Mpc^{-1}}$, depending on the number of counterterms used, with potentially large improvement over former analytical techniques.

In the next few years, we are going to probe the low-redshift universe with unprecedented accuracy. Among the various fruits that this will bear, it will greatly improve our knowledge of the dynamics of dark energy, though for this there is a strong theoretical preference for a cosmological constant. We assume that dark energy is described by the so-called Effective Field Theory of Dark Energy, which assumes that dark energy is the Goldstone boson of time translations. Such a formalism makes it easy to ensure that our signatures are consistent with well-established principles of physics. Since most of the information resides at high wavenumbers, it is important to be able to make predictions at the highest wavenumber that is possible. The Effective Field Theory of Large-Scale Structure (EFTofLSS) is a theoretical framework that has allowed us to make accurate predictions in the mildly non-linear regime. In this paper, we derive the non-linear equations that extend the EFTofLSS to include the effect of dark energy both on the matter fields and on the biased tracers. For the specific case of clustering quintessence, we then perturbatively solve to cubic order the resulting non-linear equations and construct the one-loop power spectrum of the total density contrast.

Because large-scale structure surveys may very well be the next leading sources of cosmological information, it is important to have a precise understanding of the cosmological observables; for this reason, the Effective Field Theory of Large-Scale Structure (EFTofLSS) was developed. So far, most results in the EFTofLSS have used the so-called Einstein-de Sitter approximation, an approximation of the time dependence which is known to be accurate to better than one percent. However, in order to reach even higher accuracy, the full time dependence must be used. The computation with exact time dependence is sensitive to both infrared (IR) and ultraviolet (UV) effects in the loop integrands, and while these effects must cancel because of diffeomorphism invariance, they make numerical computation much less efficient. We provide a formulation of the one-loop, equal-time exact-time-dependence power spectrum of density perturbations which is manifestly free of these spurious IR and UV divergences at the level of the integrand. We extend our results to the total matter mode with clustering quintessence, show that IR and UV divergences cancel, and provide the associated IR- and UV-safe integrand. This also establishes that the consistency conditions are satisfied in this system. We then use our one-loop result to do an improved precision comparison of the two-loop dark-matter power spectrum with the Dark Sky N-body simulation.

The large scale structures of the universe will likely be the next leading source of cosmological information. It is therefore crucial to understand their behavior. The Effective Field Theory of Large Scale Structures provides a consistent way to perturbatively predict the clustering of dark matter at large distances. The fact that baryons move distances comparable to dark matter allows us to infer that baryons at large distances can be described in a similar formalism: the backreaction of short-distance non-linearities and of star-formation physics at long distances can be encapsulated in an effective stress tensor, characterized by a few parameters. The functional form of baryonic effects can therefore be predicted. In the power spectrum the leading contribution goes as $\propto k^2 P(k)$, with $P(k)$ being the linear power spectrum and with the numerical prefactor depending on the details of the star-formation physics. We also perform the resummation of the contribution of the long-wavelength displacements, allowing us to consistently predict the effect of the relative motion of baryons and dark matter. We compare our predictions with simulations that contain several implementations of baryonic physics, finding percent agreement up to relatively high wavenumbers such as $k\simeq 0.3\,h\, Mpc^{-1}$ or $k\simeq 0.6\, h\, Mpc^{-1}$, depending on the order of the calculation. Our results open a novel way to understand baryonic effects analytically, as well as to interface with simulations.

We make an extension to recent calculations of the probability density \rho(V) for the volume of the universe after inflation. Previous results have been accurate to leading order in the slow roll parameters \epsilon=\dot{H}/H^2 and \eta=\ddot{\phi}/(\dot{\phi} H), and 1/N_c, where H is the Hubble parameter and N_c is the classical number of e-foldings. Here, we present a modification which captures effects of order \epsilon N_c, which amounts to letting the parameters of inflation H and \dot{\phi} depend on the value of the inflaton \phi. The phase of slow roll eternal inflation can be defined as when the probability to have an infinite volume is greater than zero. Using this definition, we study the Laplace transform of \rho(V) numerically to determine the condition that triggers the transition to eternal inflation. We also study the average volume <V> analytically and show that it satisfies the universal volume bound. This bound states that, in any realization of inflation which ends with a finite volume, an initial volume must grow by less than a factor of exp(S_{dS}/2), where S_{dS} is the de Sitter (dS) entropy.

We analyze the signatures of inflationary models that are coupled to strongly interacting field theories, a basic class of multifield models also motivated by their role in providing dynamically small scales. Near the squeezed limit of the bispectrum, we find a simple scaling behavior determined by operator dimensions, which are constrained by the appropriate unitarity bounds. Specifically, we analyze two simple and calculable classes of examples: conformal field theories (CFTs), and large-N CFTs deformed by relevant time-dependent double-trace operators. Together these two classes of examples exhibit a wide range of scalings and shapes of the bispectrum, including nearly equilateral, orthogonal and local non-Gaussianity in different regimes. Along the way, we compare and contrast the shape and amplitude with previous results on weakly coupled fields coupled to inflation. This signature provides a precision test for strongly coupled sectors coupled to inflation via irrelevant operators suppressed by a high mass scale up to 1000 times the inflationary Hubble scale.

We study k-defects - topological defects in theories with more than two derivatives and second-order equations of motion - and describe some striking ways in which these defects both resemble and differ from their analogues in canonical scalar field theories. We show that, for some models, the homotopy structure of the vacuum manifold is insufficient to establish the existence of k-defects, in contrast to the canonical case. These results also constrain certain families of DBI instanton solutions in the 4-dimensional effective theory. We then describe a class of k-defect solutions, which we dub doppelgangers, that precisely match the field profile and energy density of their canonical scalar field theory counterparts. We give a complete characterization of Lagrangians which admit doppelganger domain walls. By numerically computing the fluctuation eigenmodes about domain wall solutions, we find different spectra for doppelgangers and canonical walls, allowing us to distinguish between k-defects and the canonical walls they mimic. We search for doppelgangers for cosmic strings by numerically constructing solutions of DBI and canonical scalar field theories. Despite investigating several examples, we are unable to find doppelganger cosmic strings, hence the existence of doppelgangers for defects with codimension >1 remains an open question.

We present photometric and spectral observation for four novae: V2362 Cyg, V2467 Cyg, V458 Vul, V2491 Cyg. All objects belongs to the "fast novae" class. For these stars we observed different departures from a typical behavior in the light curve and spectrum.

We present the results of unfiltered and UBVR_cI_c band CCD photometry fast nova V2467Cyg. Our analysis of the data gives two distinct frequencies corresponding to periods of ~3.8h and ~35 min. The observed light curve of V2467 Cyg is typical for an intermediate polar.