The LIGO observatories can potentially detect stochastic gravitational waves arising from phase transitions which happened in the early universe at temperatures around $T\sim 10^{8}$ GeV. This provides an extraordinary opportunity for discovering the phase transition associated with the breaking of the Peccei-Quinn symmetry, required in QCD axion models. Here we consider the simplest Peccei-Quinn models and study under which conditions a strong first-order phase transition can occur, analyzing its associated gravitational wave signal. To be detectable at LIGO, we show that some supercooling is needed, which can arise either in Coleman-Weinberg-type symmetry breaking or in strongly-coupled models. We also investigate phase transitions that interestingly proceed by first breaking the electroweak symmetry at large scales before tunneling to the Peccei-Quinn breaking vacuum. In this case, the associated gravitational wave signal is more likely to be probed at the proposed Einstein Telescope.

Oscillons are bound states sustained by self-interactions that appear in rather generic scalar models. They can be extremely long-lived and in the context of cosmology they have a built-in formation mechanism - parametric resonance instability. These features suggest that oscillons can affect the standard picture of scalar ultra-light dark matter (ULDM) models. We explore this idea along two directions. First, we investigate numerically oscillon lifetimes and their dependence on the shape of the potential. We find that scalar potentials that occur in well motivated axion-like models can lead to oscillons that live up to $10^8$ cycles or more. Second, we discuss the observational constraints on the ULDM models once the presence of oscillons is taken into account. For a wide range of axion masses, oscillons decay around or after matter-radiation equality and can thus act as early seeds for structure formation. We also discuss the possibility that oscillons survive up to today. In this case they can most easily play the role of dark matter.

We propose a mechanism to generate Primordial Black Holes (PBHs) which is independent of cosmological inflation and occurs slightly below the QCD phase transition. Our setup relies on the collapse of long-lived string-domain wall networks and is naturally realized in QCD axion models with domain wall number $N_{DW}>1$ and Peccei-Quinn symmetry broken after inflation. In our framework, dark matter is mostly composed of axions in the meV mass range along with a small fraction, $\Omega_{\text{PBH}} \gtrsim 10^{-6} \Omega_{\text{CDM}} $ of heavy $M \sim 10^4-10^7 M_\odot$ PBHs. The latter could play a role in alleviating some of the shortcomings of the $\Lambda$CDM model on sub-galactic scales. The scenario has distinct signatures in ongoing axion searches as well as gravitational wave observatories.

We present a standard hydrodynamical description for non-canonical scalar field theories with kinetic gravity braiding. In particular, this picture applies to the simplest galileons and k-essence. The fluid variables not only have a clear physical meaning but also drastically simplify the analysis of the system. The fluid carries charges corresponding to shifts in field space. This shift-charge current contains a spatial part responsible for diffusion of the charges. Moreover, in the incompressible limit, the equation of motion becomes the standard diffusion equation. The fluid is indeed imperfect because the energy flows neither along the field gradient nor along the shift current. The fluid has zero vorticity and is not dissipative: there is no entropy production, the energy-momentum is exactly conserved, the temperature vanishes and there is no shear viscosity. Still, in an expansion around a perfect fluid one can identify terms which correct the pressure in the manner of bulk viscosity. We close by formulating the non-trivial conditions for the thermodynamic equilibrium of this imperfect fluid.

We introduce a large class of scalar-tensor models with interactions containing the second derivatives of the scalar field but not leading to additional degrees of freedom. These models exhibit peculiar features, such as an essential mixing of scalar and tensor kinetic terms, which we have named kinetic braiding. This braiding causes the scalar stress tensor to deviate from the perfect-fluid form. Cosmology in these models possesses a rich phenomenology, even in the limit where the scalar is an exact Goldstone boson. Generically, there are attractor solutions where the scalar monitors the behaviour of external matter. Because of the kinetic braiding, the position of the attractor depends both on the form of the Lagrangian and on the external energy density. The late-time asymptotic of these cosmologies is a de Sitter state. The scalar can exhibit phantom behaviour and is able to cross the phantom divide with neither ghosts nor gradient instabilities. These features provide a new class of models for Dark Energy. As an example, we study in detail a simple one-parameter model. The possible observational signatures of this model include a sizeable Early Dark Energy and a specific equation of state evolving into the final de-Sitter state from a healthy phantom regime.

Horava's proposal for non-relativistic quantum gravity introduces a preferred time foliation of space-time which violates the local Lorentz invariance. The foliation is encoded in a dynamical scalar field which we call `khronon'. The dynamics of the khronon field is sensitive to the symmetries and other details of the particular implementations of the proposal. In this paper we examine several consistency issues present in three non-relativistic gravity theories: Horava's projectable theory, the healthy non-projectable extension, and a new extension related to ghost condensation. We find that the only model which is free from instabilities and strong coupling is the non-projectable one. We elaborate on the phenomenology of the latter model including a discussion of the couplings of the khronon to matter. In particular, we obtain the parameters of the post-Newtonian expansion in this model and show that they are compatible with current observations.

We address the consistency of Horava's proposal for a theory of quantum gravity from the low-energy perspective. We uncover the additional scalar degree of freedom arising from the explicit breaking of the general covariance and study its properties. The analysis is performed both in the original formulation of the theory and in the Stueckelberg picture. A peculiarity of the new mode is that it satisfies an equation of motion that is of first order in time derivatives. At linear level the mode is manifest only around spatially inhomogeneous and time-dependent backgrounds. We find two serious problems associated with this mode. First, the mode develops very fast exponential instabilities at short distances. Second, it becomes strongly coupled at an extremely low cutoff scale. We also discuss the "projectable" version of Horava's proposal and argue that this version can be understood as a certain limit of the ghost condensate model. The theory is still problematic since the additional field generically forms caustics and, again, has a very low strong coupling scale. We clarify some subtleties that arise in the application of the Stueckelberg formalism to Horava's model due to its non-relativistic nature.

We propose a natural extension of Horava's model for quantum gravity, which is free from the notorious pathologies of the original proposal. The new model endows the scalar graviton mode with a regular quadratic action and remains power-counting renormalizable. At low energies, it reduces to a Lorentz-violating scalar-tensor gravity theory. The deviations with respect to general relativity can be made weak by an appropriate choice of parameters.

We present a higher codimension generalization of the DGP scenario which, unlike previous attempts, is free of ghost instabilities. The 4D propagator is made regular by embedding our visible 3-brane within a 4-brane, each with their own induced gravity terms, in a flat 6D bulk. The model is ghost-free if the tension on the 3-brane is larger than a certain critical value, while the induced metric remains flat. The gravitational force law `cascades' from a 6D behavior at the largest distances followed by a 5D and finally a 4D regime at the shortest scales.

We revisit the issue of the stability in the Dvali-Gabadadze-Porrati model, by considering the nucleation of bubbles of the conventional branch within the self-accelerating branch. We construct an instanton describing this process in the thin wall approximation. On one side of the bubble wall, the bulk consists of the exterior of the brane while on the other side it is the interior. The solution requires the presence of a 2-brane (the bubble wall) which induces the transition. However, we show that this instanton cannot be realized as the thin wall limit of any smooth solution. Once the bubble thickness is resolved, the equations of motion do not allow O(4) symmetric solutions joining the two branches. We conclude that the thin wall instanton is unphysical, and that one cannot have processes connecting the two branches, unless negative tension bubble walls are introduced. This also suggests that the self-accelerating branch does not decay into the conventional branch nucleating bubbles. We comment on other kinds of bubbles that could interpolate between the two branches.

The Dvali-Gabadadze-Porrati model introduces a parameter, the cross-over scale $r_c$, setting the scale where higher dimensional effects are important. In order to agree with observations and to explain the current acceleration of the Universe, $r_c$ must be of the order of the present Hubble radius. We discuss a mechanism to generate a large $r_c$, assuming that it is determined by a dynamical field and exploiting the quantum effects of the graviton. For simplicity, we consider a scalar field $\Psi$ with a kinetic term on the brane instead of the full metric perturbations. We compute the Green function and the 1-loop expectation value of the stress tensor of $\Psi$ on the background defined by a flat bulk and an inflating brane (self-accelerated or not). We also include the flat brane limit. The quantum fluctuations of the bulk field $\Psi$ provide an effective potential for $r_c$. For a flat brane, the 1-loop effective potential is of the Coleman-Weinberg form, and admits a minimum for large $r_c$ without fine tuning. When we take into account the brane curvature, a sizeable contribution at the classical level changes this picture and the potential develops a (minimum) maximum for the (non-) self-accelerated branch.

We consider a class of 5-D brane-world solutions with a power-law warp factor $a(y)\propto y^{q}$, and bulk dilaton with profile $\phi \propto \ln y$, where $y$ is the proper distance in the extra dimension. This class includes the Heterotic M-theory brane-world of Refs. \cite{ovrut,ekpy} and the Randall-Sundrum (RS) model as a limiting case. In general, there are two moduli fields $y_{\pm}$, corresponding to the "positions" of two branes (which live at the fixed points of an orbifold compactification). Classically, the moduli are massless, due to a scaling symmetry of the action. However, in the absence of supersymmetry, they develop an effective potential at one loop. Local terms proportional to $K_\pm^4$, where $K_\pm=q/y_\pm$ is the local curvature scale at the location of the corresponding brane, are needed in order to remove the divergences in the effective potential. Such terms break the scaling symmetry and hence they may act as stabilizers for the moduli. When the branes are very close to each other, the effective potential induced by massless bulk fields behaves like $V\sim d^{-4}$, where $d$ is the separation between branes. When the branes are widely separated, the potentials for each one of the moduli generically develop a "Coleman-Weinberg"-type behaviour of the form $a^4(y_\pm)K_\pm^4 \ln(K_\pm/\mu_\pm)$, where $\mu_\pm$ are renormalization scales. In the RS case, the bulk geometry is $AdS$ and $K_\pm$ are equal to a constant, independent of the position of the branes, so these terms do not contribute to the mass of the moduli. However, for generic warp factor, they provide a simple stabilization mechanism. For $q\gtrsim 10$, the observed hierarchy can be naturally generated by this potential, giving the lightest modulus a mass of order $m_- \lesssim TeV$.

We show that in brane-world scenarios with warped extra dimensions, the Casimir force due to bulk matter fields may be sufficient to stabilize the radion field $\phi$. In particular, we calculate one loop effective potential for $\phi$ induced by bulk gravitons and other possible massless bulk fields in the Randall-Sundrum background. This potential has a local extremum, which can be a maximum or a minimum depending on the detailed bulk matter content. If the parameters of the background are chosen so that the hierarchy problem is solved geometrically, then the radion mass induced by Casimir corrections is hierarchycally smaller than the $TeV$. Hence, in this important case, we must invoke an alternative mechanism (classical or nonperturbative) which gives the radion a sizable mass, to make it compatible with observations.