Shift-symmetric Horndeski theories admit an interesting class of Schwarzschild black hole solutions exhibiting time-dependent scalar hair. By making use of Lema\^{i}tre coordinates, we analyze perturbations around these types of black holes, and demonstrate that scalar perturbations around black hole backgrounds inevitably have gradient instabilities. Taken together with previously established results, this newly-discovered instability rules out black holes with time-dependent scalar hair in Horndeski theories.

We propose a new measure for eternal inflation, based on search optimization and first-passage statistics. This work builds on the dynamical selection mechanism for vacua based on search optimization proposed recently by the author and Parrikar. The approach is motivated by the possibility that eternal inflation has unfolded for a finite time much shorter than the exponentially long mixing time for the landscape. The proposed accessibility measure assigns greater weight to vacua that are accessed efficiently under time evolution. It is the analogue of the closeness centrality index widely used in network science. The proposed measure enjoys a number of desirable properties. It is simultaneously time-reparametrization invariant, independent of initial conditions, and oblivious to physical vs comoving weighing of pocket universes. Importantly, the proposed measure makes concrete and testable predictions that are largely independent of anthropic reasoning. Firstly, it favors vacua residing in regions of the landscape with funnel-like topography, akin to the energy landscape of naturally-occurring proteins. Secondly, it favors regions of the landscape that are tuned at dynamical criticality, with vacua having an average lifetime of order the de Sitter Page time. Thus the predicted lifetime of our universe is of order its Page time, $\sim 10^{130}$ years, which is compatible with Standard Model estimates for electroweak metastability. Relatedly, the supersymmetry breaking scale should be high, at least $10^{10}$ GeV. The discovery of beyond-the-Standard Model particles at the Large Hadron Collider or future accelerators, including low-scale supersymmetry, would rule out the possibility that our vacuum lies in an optimal region of the landscape. The present framework suggests a correspondence between the near-criticality of our universe and dynamical critical phenomena on the string landscape.

The empirical scaling relations observed in disk galaxies remain challenging for models of galaxy formation. The most striking among these is the Mass Discrepancy-Acceleration Relation (MDAR), which encodes both a tight baryonic Tully-Fisher relation (BTFR) and the observed diversity of galaxy rotation curves through the central surface density relation (CSDR). Building on our earlier work, we propose here that the MDAR is the result of interactions between baryons and 'Baryon-Interacting Dark Matter' (BIDM), which heat up the dark matter. Following a bottom-up, hydrodynamical approach, we find that the MDAR follows if: $i)$ the BIDM equation of state approximates that of an ideal gas; $ii)$ the BIDM relaxation time is order the Jeans time; $iii)$ the heating rate is inversely proportional to the BIDM density. Remarkably, under these assumptions the set of hydrodynamical equations together with Poisson's equation enjoy an anisotropic scaling symmetry. In the BIDM-dominated regime, this gives rise to an enhanced symmetry which fully captures the low-acceleration limit of the MDAR. We then show that, assuming a cored pseudo-isothermal profile at equilibrium, this set of equations gives rise to parameters reproducing the MDAR. Specifically, in the flat part of the rotation curve the asymptotic rotational velocity matches the parametric dependence of the BTFR. Moreover, in the central region of high-surface brightness galaxies, the profile reproduces the CSDR. Finally, by studying the time-dependent approach to equilibrium, we derive a global combination of the BTFR and CSDR, which matches the expectations in low surface-brightness galaxies. The form of the heating rate also makes model-independent predictions for various cosmological observables. We argue that our scenario satisfies existing observational constraints, and, intriguingly, offers a possible explanation to the EDGES anomaly.

A striking feature of our universe is its near criticality. The cosmological constant and weak hierarchy problems, as well as the metastability of the electroweak vacuum, can all be understood as problems of criticality. This suggests a statistical physics approach, based on the string landscape. We present a dynamical selection mechanism for hospitable vacua based on search optimization. Instead of focusing on stationary probability distributions, we are interested in the approach to equilibrium. This is particularly relevant if cosmological evolution has occurred for a time shorter than the exponentially-long mixing time for the landscape. This imposes a strong selection pressure among hospitable vacua, favoring those in regions where the search algorithm is efficient. Specifically, the mean first passage time is minimized for hospitable vacua lying at the bottom of funnel-like regions, akin to the smooth folding funnels of proteins and the loss functions of well-trained deep neural networks. The optimality criterion is time-reparametrization invariant and defined by two competing requirements: search efficiency, which requires minimizing the mean first passage time, and sweeping exploration, which requires that random walks are recurrent. Optimal landscape regions reach a compromise by lying at the critical boundary between recurrence and transience, thereby achieving dynamical criticality. Remarkably, this implies that the optimal lifetime of vacua coincides with the de Sitter Page time. Our mechanism makes concrete phenomenological predictions: 1. The expected lifetime of our universe is $10^{130}$ years, which is ~$2\sigma$ from the Standard Model metastability estimate; 2. The SUSY breaking scale should be high, $>10^{10}$ GeV. The present framework suggests a correspondence between the near-criticality of our universe and non-equilibrium critical phenomena on the landscape.

We present a novel theory of a unified dark sector, where late-time cosmic acceleration emerges from the dark matter superfluid framework. The system is described by a superfluid mixture consisting of two distinguishable states with a small energy gap, such as the ground state and an excited state of dark matter. Given their contact in the superfluid, interaction between those states can happen, converting one state into the other. This long range interaction within the superfluid couples the two superfluid phonon species through a cosine potential motivated by Josephson/Rabi interactions. As a consequence of this potential, a new dynamics of late-time accelerated expansion emerges in this system, without the need of dark energy, coming from a universe containing only this two-state DM superfluid. Because the superfluid species are non-relativistic, their sound speeds remain suitably small throughout the evolution. We calculate the expansion history and growth of linear perturbations, and compare the results to $\Lambda$CDM cosmology. For the fiducial parameters studied here, the predicted expansion and growth function are close to those of $\Lambda$CDM, but the difference in the predicted growth rate is significant at late times. The present theory nicely complements the recent proposal of dark matter superfluidity to explain the empirical success of MOdified Newtonian Dynamics (MOND) on galactic scales, thus offering a unified framework for dark matter, dark energy, and MOND phenomenology.

We compute the dynamical friction on a small perturber moving through an inviscid fluid, i.e., a superfluid. Crucially, we account for the tachyonic gravitational mass for sound waves, reminiscent of the Jeans instability of the fluid, which results in non-zero dynamical friction even for subsonic velocities. Moreover, we illustrate that the standard leading order effective theory in the derivative expansion is in general inadequate for analysing supersonic processes. We show this in two ways: (i) with a fluid treatment, where we solve the linearized hydrodynamical equations coupled to Newtonian gravity; and (ii) with a quasiparticle description, where we study the energy dissipation of a moving perturber due to phonon radiation. Ordinarily a subsonic perturber moving through a superfluid is kinematically prohibited from losing energy, however the Jeans instability modifies the dispersion relation of the fluid which can result in a small but non-vanishing dynamical friction force. We also analyse the soft phonon bremsstrahlung by a subsonic perturber scattered off an external field.

We consider a massive complex scalar field with contact interactions with a source and show that, upon Bose-Einstein condensation, there is an emergent long-range interaction between sources. This interaction becomes long-range in the limit of vanishing self-interaction between Bose-Einstein constituents. More generally, the range is given by $\ell^{-1}\propto \sqrt{\lambda n/m}$, with $\lambda$ being the 2-body self-interaction coupling constant, $n$ the particle number density in the condensate, and $m$ the mass of the condensed particles. Naively this may sound surprising since in $\lambda\rightarrow 0$ limit gapless excitations of the condensate have dispersion relation $\omega_k=k^2/2m$, yet for the mediated force we have $F\propto 1/r^2$. The reason behind this seemingly counterintuitive result lies in the fact that the force is being mediated by the phonon, which happens to acquire a nontrivial derivative interaction with the source. We discuss the potential ramifications of this observation for dark matter models. In particular, we show that this force can compete with gravity on galactic scales for a wide range of dark matter mass, provided that the interaction with baryons allows the presence of an extended dark matter condensate core. The effect could be of particular interest in ultra-light dark matter models, such as Fuzzy Dark Matter.

We derive the finite-temperature equation of state of dark matter superfluids with 2-body and 3-body contact interactions. The latter case is relevant to a recently proposed model of dark matter superfluidity that unifies the collisionless aspects of dark matter with the empirical success of MOdified Newtonian Dynamics at fitting galactic rotation curves. The calculation uses a self-consistent mean-field approximation. It relies on the Hartree-Fock-Bogoliubov approximation and follows the Yukalov-Yukalova proposal to circumvent the well-known Hohenberg-Martin dilemma. The resulting equation of state is consistent with a gapless spectrum, and simultaneously satisfies the equation of motion for the condensate wavefunction. As an application, we derive the finite-temperature density profile for dark matter superfluids, assuming spherical symmetry and uniform temperature. The density profiles consist of a nearly homogeneous superfluid core, surrounded by an isothermal "atmosphere" of normal particles, with the transition taking place at the critical density.

We introduce a novel method to circumvent Weinberg's no-go theorem for self-tuning the cosmological vacuum energy: a Lorentz-violating finite-temperature superfluid can counter the effects of an arbitrarily large cosmological constant. Fluctuations of the superfluid result in the graviton acquiring a Lorentz-violating mass and we identify a unique class of theories that are pathology free, phenomenologically viable, and do not suffer from instantaneous modes. This new and hitherto unidentified phase of massive gravity propagates the same degrees of freedom as general relativity with an additional Lorentz-violating scalar that is introduced by higher-derivative operators in a UV insensitive manner. The superfluid is therefore a consistent infrared modification of gravity. We demonstrate how the superfluid can degravitate a cosmological constant and discuss its phenomenology.

The observed tightness of the mass discrepancy-acceleration relation (MDAR) poses a fine-tuning challenge to current models of galaxy formation. We propose that this relation could arise from collisional interactions between baryons and dark matter (DM) particles, without the need for modification of gravity or ad hoc feedback processes. We assume that these interactions satisfy the following three conditions: (i) the relaxation time of DM particles is comparable to the dynamical time in disk galaxies; (ii) DM exchanges energy with baryons due to elastic collisions; (iii) the product between the baryon-DM cross section and the typical energy exchanged in a collision is inversely proportional to the DM number density. We present an example of a particle physics model that gives a DM-baryon cross section with the desired density and velocity dependence. Direct detection constraints require our DM particles to be either very light ($m << m_b$) or very heavy ($m >> m_b$), corresponding respectively to heating and cooling of DM by baryons. In both cases, our mechanism applies and an equilibrium configuration can in principle be reached. Here, we focus on the heavy DM/cooling case as it is technically simpler. Under these assumptions, we find that rotationally-supported disk galaxies could naturally settle to equilibrium configurations satisfying a MDAR at all radii without invoking finely tuned feedback processes. We also discuss issues related to the small scale clumpiness of baryons, as well as predictions for pressure-supported systems. We argue in particular that galaxy clusters do not follow the MDAR despite being DM-dominated because they have not reached their equilibrium configuration. Finally, we revisit existing phenomenological, astrophysical and cosmological constraints on baryon-DM interactions in light of the unusual density dependence of the cross section.

Recently, a new form of dark matter has been suggested to naturally reproduce the empirically successful aspects of Milgrom's law in galaxies. The dark matter particle candidates are axion-like, with masses of order eV and strong self-interactions. They Bose-Einstein condense into a superfluid phase in the central regions of galaxy halos. The superfluid phonon excitations in turn couple to baryons and mediate an additional long-range force. For a suitable choice of the superfluid equation of state, this force can mimic Milgrom's law. In this paper we develop in detail some of the main phenomenological consequences of such a formalism, by revisiting the expected dark matter halo profile in the presence of an extended baryon distribution. In particular, we show how rotation curves of both high and low surface brightness galaxies can be reproduced, with a slightly rising rotation curve at large radii in massive high surface brightness galaxies, thus subtly different from Milgrom's law. We finally point out other expected differences with Milgrom's law, in particular in dwarf spheroidal satellite galaxies, tidal dwarf galaxies, and globular clusters, whose Milgromian or Newtonian behavior depends on the position with respect to the superfluid core of the host galaxy. We also expect ultra-diffuse galaxies within galaxy clusters to have velocities slightly above the baryonic Tully-Fisher relation. Finally, we note that, in this framework, photons and gravitons follow the same geodesics, and that galaxy-galaxy lensing, probing larger distances within galaxy halos than rotation curves, should follow predictions closer to the standard cosmological model than those of Milgrom's law.

It has recently been proposed, by assuming that dark matter is a superfluid, that MOND-like effects can be achieved on small scales whilst preserving the success of $\Lambda$CDM on large scales. Here we aim to provide the first set of spherical models of galaxy clusters in the context of superfluid dark matter. We first outline the theoretical structure of the superfluid core and the surrounding "normal phase" dark halo of quasi-particles in thermal equlibrium. The latter should encompass the largest part of galaxy clusters. Here, we set the SfDM transition at the radius where the density and pressure of the superfluid and normal phase coincides, neglecting the effect of phonons in the suprefluid core. We then apply the theory to a sample of galaxy clusters, and directly compare the SfDM predicted mass profiles to data. We find that the superfluid formulation can reproduce the X-ray dynamical mass profile of clusters. The SfDM fits however display slight under-predictions of the gravity in the central regions which might be partly related to our neglecting of the effect of phonons in these regions. We conclude that this superfluid formulation is successful in describing galaxy clusters, but further work will be needed to determine whether the parameter choice is consistent with galaxies. Our model could be made more realistic by exploring non-sphericity and the SfDM transition condition we impose.

We investigate the symmetry structure of inflation in 2+1 dimensions. In particular, we show that the asymptotic symmetries of three-dimensional de Sitter space are in one-to-one correspondence with cosmological adiabatic modes for the curvature perturbation. In 2+1 dimensions, the asymptotic symmetry algebra is infinite-dimensional, given by two copies of the Virasoro algebra, and can be traced to the conformal symmetries of the two-dimensional spatial slices of de Sitter. We study the consequences of this infinite-dimensional symmetry for inflationary correlation functions, finding new soft theorems that hold only in 2+1 dimensions. Expanding the correlation functions as a power series in the soft momentum $q$, these relations constrain the traceless part of the tensorial coefficient at each order in $q$ in terms of a lower-point function. As a check, we verify that the ${\cal O}(q^2)$ identity is satisfied by inflationary correlation functions in the limit of small sound speed.

Gravity is the weakest fundamental interaction and the only one that has not been measured at the particle level. Traditional experimental methods, from astronomical observations to torsion balances, use macroscopic masses to both source and probe gravitational fields. Matter wave interferometers have used neutrons, atoms and molecular clusters as microscopic test particles, but initially probed the field sourced by the entire earth. Later, the gravitational field arising from hundreds of kilograms of artificial source masses was measured with atom interferometry. Miniaturizing the source mass and moving it into the vacuum chamber could improve positioning accuracy, allow the use of monocrystalline source masses for improved gravitational measurements, and test new physics, such as beyond-standard-model ("fifth") forces of nature and non-classical effects of gravity. In this work, we detect the gravitational force between freely falling cesium atoms and an in-vacuum, centimeter-sized source mass using atom interferometry with state-of-the-art sensitivity. The ability to sense gravitational-strength coupling is conjectured to access a natural lower bound for fundamental forces, thereby representing an important milestone in searches for physics beyond the standard model. A local, in-vacuum source mass is particularly sensitive to a wide class of interactions whose effects would otherwise be suppressed beyond detectability in regions of high matter density. For example, our measurement strengthens limits on a number of cosmologically-motivated scalar field models, such as chameleon and symmetron fields, by over two orders of magnitude and paves the way toward novel measurements of Newton's gravitational constant G and the gravitational Aharonov-Bohm effect

Cosmic acceleration is widely believed to require either a source of negative pressure (i.e., dark energy), or a modification of gravity, which necessarily implies new degrees of freedom beyond those of Einstein gravity. In this paper we present a third possibility, using only dark matter and ordinary matter. The mechanism relies on the coupling between dark matter and ordinary matter through an effective metric. Dark matter couples to an Einstein-frame metric, and experiences a matter-dominated, decelerating cosmology up to the present time. Ordinary matter couples to an effective metric that depends also on the DM density, in such a way that it experiences late-time acceleration. Linear density perturbations are stable and propagate with arbitrarily small sound speed, at least in the case of `pressure' coupling. Assuming a simple parametrization of the effective metric, we show that our model can successfully match a set of basic cosmological observables, including luminosity distance, BAO measurements, angular-diameter distance to last scattering {\it etc.} For the growth history of density perturbations, we find an intriguing connection between the growth factor and the Hubble constant. To get a growth history similar to the $\Lambda$CDM prediction, our model predicts a higher $H_0$, closer to the value preferred by direct estimates. On the flip side, we tend to overpredict the growth of structures whenever $H_0$ is comparable to the Planck preferred value. The model also tends to predict larger redshift-space distortions at low redshift than $\Lambda$CDM.

In this talk I summarize a novel framework that unifies the stunning success of MOND on galactic scales with the triumph of the $\Lambda$CDM model on cosmological scales. This is achieved through the rich and well-studied physics of superfluidity. The dark matter and MOND components have a common origin, representing different phases of a single underlying substance. In galaxies, dark matter thermalizes and condenses to form a superfluid phase. The superfluid phonons couple to baryonic matter particles and mediate a MOND-like force. This framework naturally distinguishes between galaxies (where MOND is successful) and galaxy clusters (where MOND is not): dark matter has a higher temperature in clusters, and hence is in a mixture of superfluid and normal phase. The rich and well-studied physics of superfluidity leads to a number of striking observational signatures, which we briefly discuss. Remarkably the critical temperature and equation of state of the dark matter superfluid are similar to those of known cold atom systems. Identifying a precise cold atom analogue would give important insights on the microphysical interactions underlying DM superfluidity. Tantalizingly, it might open the possibility of simulating the properties and dynamics of galaxies in laboratory experiments.

In recent work we proposed a novel theory of dark matter (DM) superfluidity that matches the successes of the LambdaCDM model on cosmological scales while simultaneously reproducing MOdified Newtonian Dynamics (MOND) phenomenology on galactic scales. The agents responsible for mediating the MONDian force law are superfluid phonons that couple to ordinary (baryonic) matter. In this paper we propose an alternative way for the MOND phenomenon to emerge from DM superfluidity. The central idea is to use higher-gradient corrections in the superfluid effective theory. These next-to-leading order terms involve gradients of the gravitational potential, and therefore effectively modify the gravitational force law. In the process we discover a novel mechanism for generating the non-relativistic MOND action, starting from a theory that is fully analytic in all field variables. The idea, inspired by the symmetron mechanism, uses the spontaneous breaking of a discrete symmetry. For large acceleration, the symmetry is unbroken and the action reduces to Einstein gravity. For small acceleration, the symmetry is spontaneously broken and the action reduces to MONDian gravity. Cosmologically, however, the universe is always in the Einstein-gravity, symmetry-restoring phase. The expansion history and linear growth of density perturbations are therefore indistinguishable from LambdaCDM cosmology.

We propose a novel theory of dark matter (DM) superfluidity that matches the successes of the LambdaCDM model on cosmological scales while simultaneously reproducing the MOdified Newtonian Dynamics (MOND) phenomenology on galactic scales. The DM and MOND components have a common origin, representing different phases of a single underlying substance. DM consists of axion-like particles with mass of order eV and strong self-interactions. The condensate has a polytropic equation of state P~rho^3 giving rise to a superfluid core within galaxies. Instead of behaving as individual collisionless particles, the DM superfluid is more aptly described as collective excitations. Superfluid phonons, in particular, are assumed to be governed by a MOND-like effective action and mediate a MONDian acceleration between baryonic matter particles. Our framework naturally distinguishes between galaxies (where MOND is successful) and galaxy clusters (where MOND is not): due to the higher velocity dispersion in clusters, and correspondingly higher temperature, the DM in clusters is either in a mixture of superfluid and normal phase, or fully in the normal phase. The rich and well-studied physics of superfluidity leads to a number of observational signatures: array of low-density vortices in galaxies, merger dynamics that depend on the infall velocity vs phonon sound speed; distinct mass peaks in bullet-like cluster mergers, corresponding to superfluid and normal components; interference patters in super-critical mergers. Remarkably, the superfluid phonon effective theory is strikingly similar to that of the unitary Fermi gas, which has attracted much excitement in the cold atom community in recent years. The critical temperature for DM superfluidity is of order mK, comparable to known cold atom Bose-Einstein condensates.

Atom interferometry experiments are searching for evidence of chameleon scalar fields with ever-increasing precision. As experiments become more precise, so too must theoretical predictions. Previous work has made numerous approximations to simplify the calculation, which in general requires solving a 3-dimensional nonlinear partial differential equation (PDE). In this paper, we introduce a new technique for calculating the chameleonic force, using a numerical relaxation scheme on a uniform grid. This technique is more general than previous work, which assumed spherical symmetry to reduce the PDE to a 1-dimensional ordinary differential equation (ODE). We examine the effects of approximations made in previous efforts on this subject, and calculate the chameleonic force in a set-up that closely mimics the recent experiment of Hamilton et al. Specifically, we simulate the vacuum chamber as a cylinder with dimensions matching those of the experiment, taking into account the backreaction of the source mass, its offset from the center, and the effects of the chamber walls. Remarkably, the acceleration on a test atomic particle is found to differ by only 20% from the approximate analytical treatment. These results allow us to place rigorous constraints on the parameter space of chameleon field theories, although ultimately the constraint we find is the same as the one we reported in Hamilton et al. because we had slightly underestimated the size of the vacuum chamber. This new computational technique will continue to be useful as experiments become even more precise, and will also be a valuable tool in optimizing future searches for chameleon fields and related theories.

In this talk we present a novel framework that unifies the stunning success of MOND on galactic scales with the triumph of the LambdaCDM model on cosmological scales. This is achieved through the rich and well-studied physics of superfluidity. The dark matter and MOND components have a common origin, representing different phases of a single underlying substance. In galaxies, dark matter thermalizes and condenses to form a superfluid phase. The superfluid phonons couple to baryonic matter particles and mediate a MOND-like force. Our framework naturally distinguishes between galaxies (where MOND is successful) and galaxy clusters (where MOND is not): dark matter has a higher temperature in clusters, and hence is in a mixture of superfluid and normal phase. The rich and well-studied physics of superfluidity leads to a number of striking observational signatures, which we briefly discuss. Remarkably the critical temperature and equation of state of the dark matter superfluid are similar to those of known cold atom systems. Identifying a precise cold atom analogue would give important insights on the microphysical interactions underlying DM superfluidity. Tantalizingly, it might open the possibility of simulating the properties and dynamics of galaxies in laboratory experiments.

If dark energy --- which drives the accelerated expansion of the universe --- consists of a light scalar field, it might be detectable as a "fifth force" between normal-matter objects, in potential conflict with precision tests of gravity. Chameleon fields and other theories with screening mechanisms, however, can evade these tests by suppressing the forces in regions of high density, such as the laboratory. Using a cesium matter-wave interferometer near a spherical mass in an ultra-high vacuum chamber, we reduce the screening mechanism by probing the field with individual atoms rather than bulk matter. Thus, we constrain a wide class of dark energy theories, including a range of chameleon and other theories that reproduce the observed cosmic acceleration.

We propose a unified framework that reconciles the stunning success of MOND on galactic scales with the triumph of the LambdaCDM model on cosmological scales. This is achieved through the physics of superfluidity. Dark matter consists of self-interacting axion-like particles that thermalize and condense to form a superfluid in galaxies, with ~mK critical temperature. The superfluid phonons mediate a MOND acceleration on baryonic matter. Our framework naturally distinguishes between galaxies (where MOND is successful) and galaxy clusters (where MOND is not): dark matter has a higher temperature in clusters, and hence is in a mixture of superfluid and normal phase. The rich and well-studied physics of superfluidity leads to a number of striking observational signatures.

The cosmology of brane induced gravity in six infinite dimensions is investigated. It is shown that a brane with Friedmann-Robertson-Walker symmetries necessarily acts as a source of cylindrically symmetric gravitational waves, so called Einstein-Rosen waves. Their existence essentially distinguishes this model from its codimension-one counterpart and necessitates solving the nonlinear system of bulk and brane-matching equations. A numerical analysis is performed and two qualitatively different and dynamically separated classes of cosmologies are derived: degravitating solutions for which the Hubble parameter settles to zero despite the presence of a non-vanishing energy density on the brane and super-accelerating solutions for which Hubble grows unbounded. The parameter space of both the stable and unstable regime is derived and observational consequences are discussed: It is argued that the degravitating regime does not allow for a phenomenologically viable cosmology. On the other hand, the super-accelerating solutions are potentially viable, however, their unstable behavior questions their physical relevance.

After a decade and a half of research motivated by the accelerating universe, theory and experiment have a reached a certain level of maturity. The development of theoretical models beyond \Lambda, or smooth dark energy, often called modified gravity, has led to broader insights into a path forward, and a host of observational and experimental tests have been developed. In this review we present the current state of the field and describe a framework for anticipating developments in the next decade. We identify the guiding principles for rigorous and consistent modifications of the standard model, and discuss the prospects for empirical tests. We begin by reviewing attempts to consistently modify Einstein gravity in the infrared, focusing on the notion that additional degrees of freedom introduced by the modification must screen themselves from local tests of gravity. We categorize screening mechanisms into three broad classes: mechanisms which become active in regions of high Newtonian potential, those in which first derivatives become important, and those for which second derivatives are important. Examples of the first class, such as f(R) gravity, employ the familiar chameleon or symmetron mechanisms, whereas examples of the last class are galileon and massive gravity theories, employing the Vainshtein mechanism. In each case, we describe the theories as effective theories. We describe experimental tests, summarizing laboratory and solar system tests and describing in some detail astrophysical and cosmological tests. We discuss future tests which will be sensitive to different signatures of new physics in the gravitational sector. Parts that are more relevant to theorists vs. observers/experimentalists are clearly indicated, in the hope that this will serve as a useful reference for both audiences, as well as helping those interested in bridging the gap between them.

We propose an alternative to particle dark matter that borrows ingredients of MOdified Newtonian Dynamics (MOND) while adding new key components. The first new feature is a dark matter fluid, in the form of a scalar field with small equation of state and sound speed. This component is critical in reproducing the success of cold dark matter for the expansion history and the growth of linear perturbations, but does not cluster significantly on non-linear scales. Instead, the missing mass problem on non-linear scales is addressed by a modification of the gravitational force law. The force law approximates MOND at large and intermediate accelerations, and therefore reproduces the empirical success of MOND at fitting galactic rotation curves. At ultra-low accelerations, the force law reverts to an inverse-square-law, albeit with a larger Newton's constant. This latter regime is important in galaxy clusters and is consistent with their observed isothermal profiles, provided the characteristic acceleration scale of MOND is mildly varying with scale or mass, such that it is ~12 times higher in clusters than in galaxies. We present an explicit relativistic theory in terms of two scalar fields. The first scalar field is governed by a Dirac-Born-Infeld action and behaves as a dark matter fluid on large scales. The second scalar field also has single-derivative interactions and mediates a fifth force that modifies gravity on non-linear scales. Both scalars are coupled to matter via an effective metric that depends locally on the fields. The form of this effective metric implies the equality of the two scalar gravitational potentials, which ensures that lensing and dynamical mass estimates agree. Further work is needed in order to make both the acceleration scale of MOND and the fraction at which gravity reverts to an inverse-square law explicitly dynamical quantities, varying with scale or mass.

We derive novel identities satisfied by inflationary correlation functions in the limit where two external momenta are taken to be small. We derive these statements in two ways: using background-wave arguments and as Ward identities following from the fixed-time path integral. Interestingly, these identities allow us to constrain some of the O(q^2) components of the soft limit, in contrast to their single-soft analogues. We provide several nontrivial checks of our identities both in the context of resonant non-Gaussianities and in small sound speed models. Additionally, we extend the relation at lowest order in external momenta to arbitrarily many soft legs, and comment on the many-soft extension at higher orders in the soft momentum. Finally, we consider how higher soft limits lead to identities satisfied by correlation functions in large-scale structure.

We study the effect of the initial state on the consistency conditions for adiabatic perturbations. In order to be consistent with the constraints of General Relativity, the initial state must be diffeomorphism invariant. As a result, we show that initial wavefunctional/density matrix has to satisfy a Slavnov-Taylor identity similar to that of the action. We then investigate the precise ways in which modified initial states can lead to violations of the consistency relations. We find two independent sources of violations: i) the state can include initial non-Gaussianities; ii) even if the initial state is Gaussian, such as a Bogoliubov state, the modified 2-point function can modify the q->0 analyticity properties of the vertex functional and result in violations of the consistency relations.

We describe a positive energy theorem for Einstein gravity coupled to scalar fields with first-derivative interactions, so-called P(X,phi) theories. We offer two independent derivations of this result. The first method introduces an auxiliary field to map the theory to a lagrangian describing two canonical scalar fields, where one can apply a positive energy result of Boucher and Townsend. The second method works directly at the P(X,phi) level and uses spinorial arguments introduced by Witten. The latter approach follows that of arXiv:1310.1663, but the end result is less restrictive. We point to the technical step where our derivation deviates from theirs. One of the more interesting implications of our analysis is to show it is possible to have positive energy in cases where dispersion relations following from locality and S-Matrix analyticity are violated.

We study a new screening mechanism which is present in Dirac-Born-Infeld (DBI)-like theories. A scalar field with a DBI-like Lagrangian is minimally coupled to matter. In the vicinity of sufficiently dense sources, non-linearities in the scalar dominate and result in an approximately constant acceleration on a test particle, thereby suppressing the scalar force relative to gravity. Unlike generic P(X) theories, screening happens within the regime of validity of the effective field theory, thanks to the DBI symmetry. This symmetry also allows the removal of a constant field gradient, like in galileons. Not surprisingly, perturbations around the spherically-symmetry background propagate superluminally, but we argue for a chronology protection analogous to galileons. We derive constraints on the theory parameters from tests of gravity and discuss various extensions.

Single-field perturbations satisfy an infinite number of consistency relations constraining the squeezed limit of correlation functions at each order in the soft momentum. These can be understood as Ward identities for an infinite set of residual global symmetries, or equivalently as Slavnov-Taylor identities for spatial diffeomorphisms. In this paper, we perform a number of novel, non-trivial checks of the identities in the context of slow-roll single field inflationary models with arbitrary sound speed. We focus for concreteness on identities involving 3-point functions with a soft external mode, and consider all possible scalar and tensor combinations for the hard-momentum modes. In all these cases, we check the consistency relations up to and including cubic order in the soft momentum. For this purpose, we compute for the first time the 3-point functions involving 2 scalars and 1 tensor, as well as 2 tensors and 1 scalar, for arbitrary sound speed.

We show that the correlation functions of any single-field cosmological model with constant growing-modes are constrained by an infinite number of novel consistency relations, which relate (N+1)-point correlation functions with a soft-momentum scalar or tensor mode to a symmetry transformation on N-point correlation functions of hard-momentum modes. We derive these consistency relations from Ward identities for an infinite tower of non-linearly realized global symmetries governing scalar and tensor perturbations. These symmetries can be labeled by an integer n. At each order n, the consistency relations constrain - completely for n=0,1, and partially for n>= 2 - the q^n behavior of the soft limits. The identities at n=0 recover Maldacena's original consistency relations for a soft scalar and tensor mode, n=1 gives the recently-discovered conformal consistency relations, and the identities for n>= 2 are new. As a check, we verify directly that the n=2 identity is satisfied by known correlation functions in slow-roll inflation.

In these Lectures, I review various extensions of the Lambda-Cold Dark Matter model, characterized by additional light degrees of freedom in the dark sector. In order to reproduce the successful phenomenology of GR in the solar system, these fields must effectively decouple from matter on solar system/laboratory scales. This is achieved through screening mechanisms, which rely on the interplay between self-interactions and coupling to matter to suppress deviations from standard gravity. The manifestation of the new degrees of freedom depends sensitively on their environment, which in turn leads to striking experimental signatures.

We provide evidence that general relativity is the unique spatially covariant effective field theory of the transverse, traceless graviton degrees of freedom. The Lorentz covariance of general relativity, having not been assumed in our analysis, is thus plausibly interpreted as an accidental or emergent symmetry of the gravitational sector.

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

Chameleons are light scalar fields with remarkable properties. Through the interplay of self-interactions and coupling to matter, chameleon particles have a mass that depends on the ambient matter density. The manifestation of the fifth force mediated by chameleons therefore depends sensitively on their environment, which makes for a rich phenomenology. In this article, we review two recent results on chameleon phenomenology. The first result a pair of no-go theorems limiting the cosmological impact of chameleons and their generalizations: i) the range of the chameleon force at cosmological density today can be at most ~Mpc; ii) the conformal factor relating Einstein- and Jordan-frame scale factors is essentially constant over the last Hubble time. These theorems imply that chameleons have negligible effect on the linear growth of structure, and cannot account for the observed cosmic acceleration except as some form of dark energy. The second result pertains to the quantum stability of chameleon theories. We show how requiring that quantum corrections be small, so as to allow reliable predictions of fifth forces, leads to an upper bound of m < 0.0073 (\rho/ 10 g cm^{-3})^{1/3} eV for gravitational strength coupling, whereas fifth force experiments place a lower bound of m>0.0042 eV. An improvement of less than a factor of two in the range of fifth force experiments could test all classical chameleon field theories whose quantum corrections are well-controlled and couple to matter with nearly gravitational strength regardless of the specific form of the chameleon potential.

We systematically derive the consistency relations associated to the non-linearly realized symmetries of theories with spontaneously broken conformal symmetry but with a linearly-realized de Sitter subalgebra. These identities relate (N+1)-point correlation functions with a soft external Goldstone to N-point functions. These relations have direct implications for the recently proposed conformal mechanism for generating density perturbations in the early universe. We study the observational consequences, in particular a novel one-loop contribution to the four-point function, relevant for the stochastic scale-dependent bias and CMB mu-distortion.

We attempt to incorporate inflation into a string theory realization of the chameleon mechanism. Previously, it was found that the volume modulus, stabilized by the supersymmetric potential used by Kachru, Kallosh, Linde and Trivedi (KKLT) and with the right choice of parameters, can generically work as a chameleon. In this paper, we ask whether inflation can be realized in the same model. We find that we need a large extra dimensions set-up, as well as a semi-phenomenological deformation of the Kahler potential in the quantum region. We also find that an additional KKLT term is required so that there are now two pieces to the potential, one which drives inflation in the early universe, and one which is responsible for chameleon screening at late times. These two pieces of the potential are separated by a large flat desert in field space. The scalar field must dynamically traverse this desert between the end of inflation and today, and we find that this can indeed occur under the right conditions.

The pseudo-conformal universe is an alternative to inflation in which the early universe is described by a conformal field theory on approximately flat space-time. The fields develop time-dependent expectation values, spontaneously breaking the conformal symmetries to a de Sitter subalgebra, and fields of conformal weight zero acquire a scale invariant spectrum of perturbations. In this paper, we show that the pseudo-conformal scenario can be naturally realized within theories that would ordinarily be of interest for DBI inflation, such as the world-volume theory of a probe brane in an AdS bulk space-time. In this approach, the weight zero spectator field can be associated with a geometric flat direction in the bulk, and its scale invariance is protected by a shift symmetry.

The chameleon, or generalizations thereof, is a light scalar that couple to matter with gravitational strength, but whose manifestation depends on the ambient matter density. A key feature is that the screening mechanism suppressing its effects in high-density environments is determined by the local scalar field value. Under very general conditions, we prove two theorems limiting its cosmological impact: i) the Compton wavelength of such a scalar can be at most Mpc at present cosmic density, which restricts its impact to non-linear scales; ii) the conformal factor relating Einstein- and Jordan-frame scale factors is essentially constant over the last Hubble time, which precludes the possibility of self-acceleration. These results imply that chameleon-like scalar fields have a negligible effect on the linear-scale growth history; theories that invoke a chameleon-like scalar to explain cosmic acceleration rely on a form of dark energy rather than a genuine modified gravity effect. Our analysis applies to a broad class of chameleon, symmetron and dilaton theories.

Chameleon scalar fields are dark energy candidates which suppress fifth forces in high density regions of the universe by becoming massive. We consider chameleon models as effective field theories and estimate quantum corrections to their potentials. Requiring that quantum corrections be small, so as to allow reliable predictions of fifth forces, leads to an upper bound $m < 0.0073 (\rho / 10 {\rm g\, cm}^{-3})^{1/3}$eV for gravitational strength coupling whereas fifth force experiments place a lower bound of $m>0.0042$\,eV. An improvement of less than a factor of two in the range of fifth force experiments could test all classical chameleon field theories whose quantum corrections are well-controlled and couple to matter with nearly gravitational strength regardless of the specific form of the chameleon potential.

We present a novel theory of the very early universe which addresses the traditional horizon and flatness problems of big bang cosmology and predicts a scale invariant spectrum of perturbations. Unlike inflation, this scenario requires no exponential accelerated expansion of space-time. Instead, the early universe is described by a conformal field theory minimally coupled to gravity. The conformal fields develop a time-dependent expectation value which breaks the flat space so(4,2) conformal symmetry down to so(4,1), the symmetries of de Sitter, giving perturbations a scale invariant spectrum. The solution is an attractor, at least in the case of a single time-dependent field. Meanwhile, the metric background remains approximately flat but slowly contracts, which makes the universe increasingly flat, homogeneous and isotropic, akin to the smoothing mechanism of ekpyrotic cosmology. Our scenario is very general, requiring only a conformal field theory capable of developing the appropriate time-dependent expectation values, and encompasses existing incarnations of this idea, specifically the U(1) model of Rubakov and the Galileon Genesis scenario. Its essential features depend only on the symmetry breaking pattern and not on the details of the underlying lagrangian. It makes generic observational predictions that make it potentially distinguishable from standard inflation, in particular significant non-gaussianities and the absence of primordial gravitational waves.

We remark on the existence of non-linearly realized conformal symmetries for scalar adiabatic perturbations in cosmology. These conformal symmetries are present for any cosmological background, beyond any slow-roll or quasi-de Sitter approximation. The dilatation transformation shifts the curvature perturbation by a constant, and corresponds to the well-known symmetry under spatial rescaling. We argue that the scalar sector is also invariant under special conformal transformations, which shift the curvature perturbation by a term linear in the spatial coordinates. We discuss whether these conformal symmetries can be extended to include tensor perturbations. Tensor modes introduce their own set of non-linearly realized symmetries. We identify an infinite set of large gauge transformations which maintain the transverse, traceless gauge condition, while shifting the tensor mode non-trivially.

We offer predictions of symmetron modified gravity in the neighborhood of realistic dark matter halos. The predictions for the fifth force are obtained by solving the nonlinear symmetron equation of motion in the spherical NFW approximation. In addition, we compare the three major known screening mechanisms: Vainshtein, Chameleon, and Symmetron around such dark matter sources, emphasizing the significant differences between them and highlighting observational tests which exploit these differences. Finally, we demonstrate the host halo environmental screening effect ("blanket screening") on smaller satellite halos by solving for the modified forces around a density profile which is the sum of satellite and approximate host components.

General relativity is a covariant theory of two transverse, traceless graviton degrees of freedom. According to a theorem of Hojman, Kuchar, and Teitelboim, modifications of general relativity must either introduce new degrees of freedom or violate the principle of general covariance. In this paper, we explore modifications of general relativity that retain the same number of gravitational degrees of freedom, and therefore explicitly break general covariance. Motivated by cosmology, the modifications of interest maintain spatial covariance. Demanding consistency of the theory forces the physical Hamiltonian density to obey an analogue of the renormalization group equation. In this context, the equation encodes the invariance of the theory under flow through the space of conformally equivalent spatial metrics. This paper is dedicated to setting up the formalism of our approach and applying it to a realistic class of theories. Forthcoming work will apply the formalism more generally.

The symmetron is a scalar field associated with the dark sector whose coupling to matter depends on the ambient matter density. The symmetron is decoupled and screened in regions of high density, thereby satisfying local constraints from tests of gravity, but couples with gravitational strength in regions of low density, such as the cosmos. In this paper we derive the cosmological expansion history in the presence of a symmetron field, tracking the evolution through the inflationary, radiation- and matter-dominated epochs, using a combination of analytical approximations and numerical integration. For a broad range of initial conditions at the onset of inflation, the scalar field reaches its symmetry-breaking vacuum by the present epoch, as assumed in the local analysis of spherically-symmetric solutions and tests of gravity. For the simplest form of the potential, the energy scale is too small for the symmetron to act as dark energy, hence we must add a cosmological constant to drive late-time cosmic acceleration. We briefly discuss a class of generalized, non-renormalizable potentials that can have a greater impact on the late-time cosmology, though cosmic acceleration requires a delicate tuning of parameters in this case.

We derive all single-field cosmologies with unit sound speed that generate scale invariant curvature perturbations on a dynamical attractor background. We identify three distinct phases: slow-roll inflation; a slowly contracting adiabatic ekpyrotic phase, described by a rapidly-varying equation of state; and a novel adiabatic ekpyrotic phase on a slowly expanding background. All of these yield identical power spectra. The degeneracy is broken at the 3-point level: unlike the nearly gaussian spectrum of slow-roll inflation, adiabatic ekpyrosis predicts large non-gaussianities on small scales.

Galileon theories are of considerable interest since they allow for stable violations of the null energy condition. Since such violations could have occurred during a high-energy regime in the history of our universe, we are motivated to study supersymmetric extensions of these theories. This is carried out in this paper, where we construct generic classes of N=1 supersymmetric Galileon Lagrangians. They are shown to admit non-equivalent stress-energy tensors and, hence, vacua manifesting differing conditions for violating the null energy condition. The temporal and spatial fluctuations of all component fields of the supermultiplet are analyzed and shown to be stable on a large number of such backgrounds. In the process, we uncover a surprising connection between conformal Galileon and ghost condensate theories, allowing for a deeper understanding of both types of theories.

We show how to construct supersymmetric actions for higher-derivative scalar field theories of the form P(X,phi), within the context of d=4, N=1 supersymmetry. This construction is of general use, and is applied to write supersymmetric versions of the Dirac-Born-Infeld action. Our principal application of this formalism is to construct the supersymmetric extension of the ghost condensate. This allows us to study the interplay between supersymmetry, time-dependent backgrounds and violations of the null energy condition.

Chameleons are scalar fields that couple directly to ordinary matter with gravitational strength, but which nevertheless evade the stringent constraints on tests of gravity because of properties they acquire in the presence of high ambient matter density. Chameleon theories were originally constructed in a bottom-up, phenomenological fashion, with potentials and matter couplings designed to hide the scalar from experiments. In this paper, we attempt to embed the chameleon scenario within string compactifications, thus UV completing the scenario. We look for stabilized potentials that can realize a screening mechanism, and we find that the volume modulus rather generically works as a chameleon, and in fact the supersymmetric potential used by Kachru, Kallosh, Linde and Trivedi (KKLT) is an example of this type. We consider all constraints from tests of gravity, allowing us to put experimental constraints on the KKLT parameters.

We consider a cosmological scenario in which a scale-invariant spectrum of curvature perturbations is generated by a rapidly-evolving equation of state on a slowly expanding background. This scenario generalizes the "adiabatic ekpyrotic" mechanism proposed recently in arXiv:0910.2230. Whereas the original proposal assumed a slowly contracting background, the present work shows that the mechanism works equally well on an expanding background. This greatly expands the realm of broader cosmological scenarios in which this mechanism can be embedded. We present a phase space analysis and show that both the expanding and contracting versions of the scenario are dynamical attractors, with the expanding branch having a broader basin of attraction. In both cases, a finite range of scale invariant modes can be generated within the regime of validity of perturbation theory.

We show that curvature perturbations acquire a scale invariant spectrum for any constant equation of state, provided the fluid has a suitably time-dependent sound speed. In order for modes to exit the physical horizon, and in order to solve the usual problems of standard big bang cosmology, we argue that the only allowed possibilities are inflationary (albeit not necessarily slow-roll) expansion or ekpyrotic contraction. Non-Gaussianities offer many distinguish features. As usual with a small sound speed, non-Gaussianity can be relatively large, around current sensitivity levels. For DBI-like lagrangians, the amplitude is negative in the inflationary branch, and can be either negative or positive in the ekpyrotic branch. Unlike the power spectrum, the three-point amplitude displays a large tilt that, in the expanding case, peaks on smallest scales. While the shape is predominantly of the equilateral type in the inflationary branch, as in DBI inflation, it is of the local form in the ekpyrotic branch. The tensor spectrum is also generically far from scale invariant. In the contracting case, for instance, tensors are strongly blue tilted, resulting in an unmeasurably small gravity wave amplitude on cosmic microwave background scales.

We study a flat brane solution in an effective 5D action for cascading gravity and propose a mechanism to screen extrinsic curvature in the presence of a large tension on the brane. The screening mechanism leaves the bulk Riemann-flat, thus making it simpler to generalize large extra dimension dark energy models to higher codimensions. By studying an action with cubic interactions for the brane-bending scalar mode, we find that the perturbed action suffers from ghostlike instabilities for positive tension, whereas it can be made ghost-free for sufficiently small negative tension.

Despite the overwhelming evidence for the existence of dark energy and dark matter, their underlying fundamental physics remains unknown. This review article explores the tantalizing possibility that the dark sector includes new light degrees of freedom that mediate long-range forces on cosmological scales. To ensure consistency with laboratory and solar system tests of gravity, some screening mechanism is necessary to "hide" these degrees of freedom locally. I will focus on two broad classes of screening theories, chameleons and symmetrons, which rely respectively on the scalar field acquiring a large mass or weak coupling in the presence of large ambient matter density.

We perform a full perturbative stability analysis of the 6D cascading gravity model in the presence of 3-brane tension. We demonstrate that for sufficiently large tension on the (flat) 3-brane, there are no ghosts at the perturbative level, consistent with results that had previously only been obtained in a specific 5D decoupling limit. These results establish the cascading gravity framework as a consistent infrared modification of gravity.

We present a screening mechanism that allows a scalar field to mediate a long range (~Mpc) force of gravitational strength in the cosmos while satisfying local tests of gravity. The mechanism hinges on local symmetry restoration in the presence of matter. In regions of sufficiently high matter density, the field is drawn towards \phi = 0 where its coupling to matter vanishes and the \phi-> -\phi symmetry is restored. In regions of low density, however, the symmetry is spontaneously broken, and the field couples to matter with gravitational strength. We predict deviations from general relativity in the solar system that are within reach of next-generation experiments, as well as astrophysically observable violations of the equivalence principle. The model can be distinguished experimentally from Brans-Dicke gravity, chameleon theories and brane-world modifications of gravity.

The mounting evidence for anomalously large peculiar velocities in our Universe presents a challenge for the LCDM paradigm. The recent estimates of the large scale bulk flow by Watkins et al. are inconsistent at the nearly 3 sigma level with LCDM predictions. Meanwhile, Lee and Komatsu have recently estimated that the occurrence of high-velocity merging systems such as the Bullet Cluster (1E0657-57) is unlikely at a 6.5-5.8 sigma level, with an estimated probability between 3.3x10^{-11} and 3.6x10^{-9} in LCDM cosmology. We show that these anomalies are alleviated in a broad class of infrared-modifed gravity theories, called brane-induced gravity, in which gravity becomes higher-dimensional at ultra large distances. These theories include additional scalar forces that enhance gravitational attraction and therefore speed up structure formation at late times and on sufficiently large scales. The peculiar velocities are enhanced by 24-34% compared to standard gravity, with the maximal enhancement nearly consistent at the 2 sigma level with bulk flow observations. The occurrence of the Bullet Cluster in these theories is 10^4 times more probable than in LCDM cosmology.

Modifications of general relativity provide an alternative explanation to dark energy for the observed acceleration of the universe. We review recent developments in modified gravity theories, focusing on higher dimensional approaches and chameleon/f(R) theories. We classify these models in terms of the screening mechanisms that enable such theories to approach general relativity on small scales (and thus satisfy solar system constraints). We describe general features of the modified Friedman equation in such theories. The second half of this review describes experimental tests of gravity in light of the new theoretical approaches. We summarize the high precision tests of gravity on laboratory and solar system scales. We describe in some detail tests on astrophysical scales ranging from ~kpc (galaxy scales) to ~Gpc (large-scale structure). These tests rely on the growth and inter-relationship of perturbations in the metric potentials, density and velocity fields which can be measured using gravitational lensing, galaxy cluster abundances, galaxy clustering and the Integrated Sachs-Wolfe effect. A robust way to interpret observations is by constraining effective parameters, such as the ratio of the two metric potentials. Currently tests of gravity on astrophysical scales are in the early stages --- we summarize these tests and discuss the interesting prospects for new tests in the coming decade.

We develop a fully covariant, well-posed 5D effective action for the 6D cascading gravity brane-world model, and use this to study cosmological solutions. We obtain this effective action through the 6D decoupling limit, in which an additional scalar degree mode, \pi, called the brane-bending mode, determines the bulk-brane gravitational interaction. The 5D action obtained this way inherits from the sixth dimension an extra \pi self-interaction kinetic term. We compute appropriate boundary terms, to supplement the 5D action, and hence derive fully covariant junction conditions and the 5D Einstein field equations. Using these, we derive the cosmological evolution induced on a 3-brane moving in a static bulk. We study the strong- and weak-coupling regimes analytically in this static ansatz, and perform a complete numerical analysis of our solution. Although the cascading model can generate an accelerating solution in which the \pi field comes to dominate at late times, the presence of a critical singularity prevents the \pi field from dominating entirely. Our results open up the interesting possibility that a more general treatment of degravitation in a time-dependent bulk, or taking into account finite brane-thickness effects, may lead to an accelerating universe without a cosmological constant.

The universe can be made flat and smooth by undergoing a phase of ultra-slow (ekpyrotic) contraction, a condition achievable with a single, canonical scalar field and conventional general relativity. It has been argued, though, that generating scale-invariant density perturbations, requires at least two scalar fields and a two-step process that first produces entropy fluctuations and then converts them to curvature perturbations. In this paper, we identify a loophole in the argument and introduce an ekpyrotic model based on a single, canonical scalar field that generates nearly scale-invariant curvature fluctuations through a purely "adiabatic mechanism" in which the background evolution is a dynamical attractor. The resulting spectrum can be slightly red with distinctive non-gaussian fluctuations.

In the Cascading Gravity brane-world scenario, our 3-brane lies within a succession of lower-codimension branes, each with their own induced gravity term, embedded into each other in a higher-dimensional space-time. In the 6+1-dimensional version of this scenario, we show that a 3-brane with tension remains flat, at least for sufficiently small tension that the weak-field approximation is valid. The bulk solution is nowhere singular and remains in the perturbative regime everywhere.

We perform N-body simulations of theories with infinite-volume extra dimensions, such as the Dvali-Gabadadze-Porrati (DGP) model and its higher-dimensional generalizations, where 4D gravity is mediated by massive gravitons. The longitudinal mode of these gravitons mediates an extra scalar force, which we model as a density-dependent modification to the Poisson equation. This enhances gravitational clustering, particularly on scales that have undergone mild nonlinear processing. While the standard non-linear fitting algorithm of Smith et al. overestimates this power enhancement on non-linear scales, we present a modified fitting formula that offers a remarkably good fit to our power spectra. Due to the uncertainty in galaxy bias, our results are consistent with precision power spectrum determinations from galaxy redshift surveys, even for graviton Compton wavelengths as small as 300 Mpc. Our model is sufficiently general that we expect it to capture the phenomenology of a wide class of related higher-dimensional gravity scenarios.

We study the cosmology of a galileon scalar-tensor theory, obtained by covariantizing the decoupling lagrangian of the Dvali-Gabadadze-Poratti (DGP) model. Despite being local in 3+1 dimensions, the resulting cosmological evolution is remarkably similar to that of the full 4+1-dimensional DGP framework, both for the expansion history and the evolution of density perturbations. As in the DGP model, the covariant galileon theory yields two branches of solutions, depending on the sign of the galileon velocity. Perturbations are stable on one branch and ghost-like on the other. An interesting effect uncovered in our analysis is a cosmological version of the Vainshtein screening mechanism: at early times, the galileon dynamics are dominated by self-interaction terms, resulting in its energy density being suppressed compared to matter or radiation; once the matter density has redshifted sufficiently, the galileon becomes an important component of the energy density and contributes to dark energy. We estimate conservatively that the resulting expansion history is consistent with the observed late-time cosmology, provided that the scale of modification satisfies r_c > 15 Gpc.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

We present a case that current observations may already indicate new gravitational physics on cosmological scales. The excess of power seen in the Lyman-alpha forest and small-scale CMB experiments, the anomalously large bulk flows seen both in peculiar velocity surveys and in kinetic SZ, and the higher ISW cross-correlation all indicate that structure may be more evolved than expected from LCDM. We argue that these observations find a natural explanation in models with infinite-volume (or, at least, cosmological-size) extra dimensions, where the graviton is a resonance with a tiny width. The longitudinal mode of the graviton mediates an extra scalar force which speeds up structure formation at late times, thereby accounting for the above anomalies. The required graviton Compton wavelength is relatively small compared to the present Hubble radius, of order 300-600 Mpc. Moreover, with certain assumptions about the behavior of the longitudinal mode on super-Hubble scales, our modified gravity framework can also alleviate the tension with the low quadrupole and the peculiar vanishing of the CMB correlation function on large angular scales, seen both in COBE and WMAP. This relies on a novel mechanism that cancels a late-time ISW contribution against the primordial Sachs-Wolfe amplitude.

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.

In this paper, we present a new scenario of the early Universe that contains a pre big bang Ekpyrotic phase. By combining this with a ghost condensate, the theory explicitly violates the null energy condition without developing any ghost-like instabilities. Thus the contracting universe goes through a non-singular bounce and evolves smoothly into the expanding post big bang phase. The curvature perturbation acquires a scale-invariant spectrum well before the bounce in this scenario. It is sourced by the scale-invariant entropy perturbation engendered by two ekpyrotic scalar fields, a mechanism recently proposed by Lehners et al. Since the background geometry is non-singular at all times, the curvature perturbation remains nearly constant on super horizon scales. It emerges from the bounce unscathed and imprints a scale-invariant spectrum of density fluctuations in the matter-radiation fluid at the onset of the hot big bang phase. The ekpyrotic potential can be chosen so that the spectrum has a ``red'' tilt, in accordance with the recent data from WMAP. As in the original Ekpyrotic scenario, the model predicts a negligible gravity wave signal on all observable scales. As such ``New Ekpyrotic Cosmology" provides a consistent and distinguishable alternative to inflation to account for the origin of the seeds of large scale structure.

New Ekpyrotic Cosmology is an alternative scenario of early universe cosmology in which the universe existed before the big bang. The simplest model relies on two scalar fields, whose entropy perturbation leads to a scale-invariant spectrum of density fluctuations. The ekpyrotic solution has a tachyonic instability along the entropy field direction which, a priori, appears to require fine-tuning of the initial conditions. In this paper, we show that these can be achieved naturally by adding a small positive mass term for the tachyonic field and coupling to light fermions. Then, for a wide range of initial conditions, the tachyonic field gets stabilized with the appropriate values well before the onset of the ekpyrotic phase. Furthermore, we show that ekpyrotic theory is successful in solving the flatness, horizon and homogeneity problems of standard big bang cosmology. Motivated by the analysis of the tachyonic instability, we propose a simplification of the model in terms of new field variables. Instead of requiring two exponential potentials, one for each scalar field, it suffices to consider a single nearly exponential potential for one of the fields and a tachyonic mass term along the orthogonal direction in field space. All other terms in the potential are essentially arbitrary. This greatly widens the class of ekpyrotic potentials and allows substantial freedom in determining the spectral index and possible non-Gaussianity. We present a generalized expression for the spectral index which is easily consistent with the observed red tilt. We also argue that for a wide range of potentials non-Gaussianity can be substantial, within the reach of current observations.

We study the possibility of decoupling gravity from the vacuum energy. This is effectively equivalent to promoting Newton's constant to a high-pass filter that degravitates sources of characteristic wavelength larger than a certain macroscopic (super) horizon scale L. We study the underlying physics and the consistency of this phenomenon. In particular, the absence of ghosts, already at the linear level, implies that in any such theory the graviton should either have a mass 1/L, or be a resonance of similar width. This has profound physical implications for the degravitation idea.

Mach's principle is the concept that inertial frames are determined by matter. We propose and implement a precise formulation of Mach's principle in which matter and geometry are in one-to-one correspondence. Einstein's equations are not modified and no selection principle is applied to their solutions; Mach's principle is realized wholly within Einstein's general theory of relativity. The key insight is the observation that, in addition to bulk matter, one can also add boundary matter. Specification of both boundary and bulk stress tensors uniquely specifies the geometry and thereby the inertial frames. Our framework is similar to that of the black hole membrane paradigm and, in asymptotically AdS space-times, is consistent with holographic duality.

String gas cosmology is rewritten in the Einstein frame. In an effective theory in which a gas of closed strings is coupled to a dilaton gravity background without any potential for the dilaton, the Hagedorn phase which is quasi-static in the string frame corresponds to an expanding, non-accelerating phase from the point of view of the Einstein frame. The Einstein frame curvature singularity which appears in this toy model is related to the blowing up of the dilaton in the string frame. However, for large values of the dilaton, the toy model clearly is inapplicable. Thus, there must be a new string phase which is likely to be static with frozen dilaton. With such a phase, the horizon problem can be successfully addressed in string gas cosmology. The generation of cosmological perturbations in the Hagedorn phase seeded by a gas of long strings in thermal equilibrium is reconsidered, both from the point of view of the string frame (in which it is easier to understand the generation of fluctuations) and the Einstein frame (in which the evolution equations are well known). It is shown that fixing the dilaton at some early stage is important in order to obtain a scale-invariant spectrum of cosmological fluctuations in string gas cosmology.

Scalar self interactions are known to weaken considerably the current constraints on scalar-mediated fifth forces. We consider a scalar field with a quartic self interaction and gravitation-strength Yukawa couplings to matter particles. After discussing the phenomenology of this scalar field, we assess the ability of ongoing and planned experiments to detect the fifth force mediated by such a field. Assuming that the quartic and matter couplings are of order unity, the current-generation Eot-Wash experiment at the University of Washington will be able to explore an interesting subset of parameter space. The next-generation Eot-Wash experiment is expected to be able to detect, or to rule out, the fifth force due to such a scalar with unit quartic and matter couplings at the 3 sigma confidence level.

We show that an interaction between dark matter and dark energy generically results in an effective dark energy equation of state of w<-1. This arises because the interaction alters the redshift-dependence of the matter density. An observer who fits the data treating the dark matter as non-interacting will infer an effective dark energy fluid with w<-1. We argue that the model is consistent with all current observations, the tightest constraint coming from estimates of the matter density at different redshifts. Comparing the luminosity and angular-diameter distance relations with LambdaCDM and phantom models, we find that the three models are degenerate within current uncertainties but likely distinguishable by the next generation of dark energy experiments.

Large future galaxy cluster surveys, combined with cosmic microwave background observations, can achieve a high sensitivity to the masses of cosmologically important neutrinos. We show that a weak lensing selected sample of ~100,000 clusters could tighten the current upper bound on the sum of masses of neutrino species by an order of magnitude, to a level of 0.03 eV. Since this statistical sensitivity is below the best existing lower limit on the mass of at least one neutrino species, a future detection is likely, provided that systematic errors can be controlled to a similar level.

We show that the chameleon scalar field can drive the current phase of cosmic acceleration for a large class of scalar potentials that are also consistent with local tests of gravity. These provide explicit realizations of a quintessence model where the quintessence scalar field couples directly to baryons and dark matter with gravitational strength. We analyze the cosmological evolution of the chameleon field and show the existence of an attractor solution with the chameleon following the minimum of its effective potential. For a wide range of initial conditions, spanning many orders of magnitude in initial chameleon energy density, the attractor is reached before nucleosynthesis. Surprisingly, the range of allowed initial conditions leading to a successful cosmology is wider than in normal quintessence. We discuss applications to the cyclic model of the universe and show how the chameleon mechanism weakens some of the constraints on cyclic potentials.

We show that the abundance and redshift distribution ($dN/dz$) of galaxy clusters in future high--yield cluster surveys, combined with the spatial power spectrum ($P_c(k)$) of the same clusters, can place significant constraints on the evolution of the dark energy equation of state, $w=w(a)$. We evaluate the expected errors on $w_a=-dw/da$ and other cosmological parameters using a Fisher matrix approach, and simultaneously including cluster structure evolution parameters in our analysis. We study three different types of forthcoming surveys that will identify clusters based on their X-ray emission (such as DUO, the Dark Universe Observatory), their Sunyaev--Zel'dovich (SZ) decrement (such as SPT, the South Pole Telescope), or their weak lensing (WL) shear (such as LSST, the Large Synoptic Survey Telescope). We find that combining the cluster abundance and power spectrum significantly enhances constraints from either method alone. We show that the weak-lensing survey can deliver a constraint as tight as $\Delta w_a \sim 0.1$ on the evolution of the dark energy equation of state, and that the X-ray and SZ surveys each yield $\Delta w_a \sim 0.4$ separately, or $\Delta w_a \sim 0.2$ when these two surveys are combined. For the X-ray and SZ surveys, constraints on dark energy parameters are improved by a factor of two by combining the cluster data with cosmic microwave background (CMB) anisotropy measurements by Planck, but degrade by a factor of two if the survey is required to solve simultaneously for cosmological and cluster structure evolution parameters. The constraint on $w_a$ from the weak lensing survey is improved by $\sim 25$% with the addition of Planck data.

Chameleons are scalar fields whose mass depends on the environment, specifically on the ambient matter density. While nearly massless in the cosmos, where the matter density is tiny, their mass is of order of an inverse millimeter on Earth, where the density is high. In this note, we review how chameleons can satisfy current experimental constraints on deviations from General Relativity (GR). Moreover, we study the cosmological evolution with a chameleon field and show the existence of an attractor solution, akin to the tracker solution in quintessence models. We discuss how chameleons can naturally drive the observed acceleration of the universe

We present a novel scenario where a scalar field acquires a mass which depends on the local matter density: the field is massive on Earth, where the density is high, but is essentially free in the solar system, where the density is low. All existing tests of gravity are satisfied. We predict that near-future satellite experiments could measure an effective Newton's constant in space different than that on Earth, as well as violations of the equivalence principle stronger than currently allowed by laboratory experiments.

The mass of a scalar field mediating a fifth force is tightly constrained by experiments. We show, however, that adding a quartic self-interaction for such a scalar makes most tests much less constraining: the non-linear equation of motion masks the coupling of the scalar to matter through the chameleon mechanism. We discuss consequences for fifth force experiments. In particular, we find that, with quartic coupling of order unity, a gravitational strength interaction with matter is allowed by current constraints. We show that our chameleon scalar field results in experimental signatures that could be detected through modest improvements of current laboratory set-ups.

This is a short overview of the ekpyrotic/cyclic model of the universe, an alternative to the standard big bang inflationary paradigm.

We analyze the cosmological consequences of locked inflation, a model recently proposed by Dvali and Kachru that can produce significant amounts of inflation without requiring slow-roll. We pay particular attention to the end of inflation in this model, showing that a secondary phase of saddle inflation can follow the locked inflationary era. However, this subsequent period of inflation results in a strongly scale dependent spectrum that can lead to massive black hole formation in the primordial universe. Avoiding this disastrous outcome puts strong constraints on the parameter space open to models of locked inflation.

Recent advances in understanding the propagation of perturbations through the transition from big crunch to big bang (esp. Tolley et al. hep-th/0306109) make it possible for the first time to consider the full set of phenomenological constraints on the scalar field potential in cyclic models of the universe. We show that cyclic models require a comparable degree of tuning to that needed for inflationary models. The constraints are reduced to a set of simple design rules including "fast-roll" parameters analogous to the "slow-roll" parameters in inflation.

The evidence for the accelerated expansion of the universe and the time-dependence of the fine-structure constant suggests the existence of at least one scalar field with a mass of order H_0. If such a field exists, then it is generally assumed that its coupling to matter must be tuned to unnaturally small values in order to satisfy the tests of the Equivalence Principle (EP). In this paper, we present an alternative explanation which allows scalar fields to evolve cosmologically while having couplings to matter of order unity. In our scenario, the mass of the fields depends on the local matter density: the interaction range is typically of order 1 mm on Earth (where the density is high) and of order 10-10^4 AU in the solar system (where the density is low). All current bounds from tests of General Relativity are satisfied. Nevertheless, we predict that near-future experiments that will test gravity in space will measure an effective Newton's constant different by order unity from that on Earth, as well as EP violations stronger than currently allowed by laboratory experiments. Such outcomes would constitute a smoking gun for our scenario.

We present a simple, nearly model-independent estimate that yields the predictions of the simplest inflationary and ekpyrotic/cyclic models for the spectral tilt of the primordial density inhomogeneities. Remarkably, we find that the simplest models yield an identical result: n_s is approximately 0.95. For inflation, the same estimate predicts a ratio of tensor to scalar contributions to the low multipoles of the microwave background anisotropy of T/S=20%; the tensor contribution is negligible for ekpyrotic/cyclic models.

We analyze the general conditions on the equation of state $w$ required for quantum fluctuations of a scalar field to produce a scale-invariant spectrum of density perturbations, including models which (in the four dimensional effective description) bounce from a contracting to an expanding phase. We show that there are only two robust cases: $w\approx -1$ (inflation) and $w \gg 1$ (the ekpyrotic/cyclic scenario). All other cases, including the $w \approx 0$ case considered by some authors, require extreme fine-tuning of initial conditions and/or the effective potential. For the ekpyrotic/cyclic ($w \gg 1$) case, we also analyze the small deviations from scale invariance.

The Friedmann law on the brane generically depends quadratically on the brane energy density and involves a ``dark radiation'' term due to the bulk Weyl tensor. Despite its unfamiliar form, we show how it can be derived from a standard four-dimensional Brans-Dicke theory at low energy. In particular, the dark radiation term is found to depend linearly on the brane energy densities. For any equation of state on the branes, the radion evolves such as to generate radiation-dominated cosmology. The radiation-dominated era is conventional and consistent with nucleosynthesis.

We study the generation of density perturbations in the ekpyrotic scenario for the early universe, including gravitational backreaction. We expose interesting subtleties that apply to both inflationary and ekpyrotic models. Our analysis includes a detailed proposal of how the perturbations generated in a contracting phase may be matched across a `bounce' to those in an expanding hot big bang phase. For the physical conditions relevant to the ekpyrotic scenario, we re-obtain our earlier result of a nearly scale-invariant spectrum of energy density perturbations. We find that the perturbation amplitude is typically small, as desired to match observation.

We consider conditions under which a universe contracting towards a big crunch can make a transition to an expanding big bang universe. A promising example is 11-dimensional M-theory in which the eleventh dimension collapses, bounces, and re-expands. At the bounce, the model can reduce to a weakly coupled heterotic string theory and, we conjecture, it may be possible to follow the transition from contraction to expansion. The possibility opens the door to new classes of cosmological models. For example, we discuss how it suggests a major simplification and modification of the recently proposed ekpyrotic scenario.

We propose a cosmological scenario in which the hot big bang universe is produced by the collision of a brane in the bulk space with a bounding orbifold plane, beginning from an otherwise cold, vacuous, static universe. The model addresses the cosmological horizon, flatness and monopole problems and generates a nearly scale-invariant spectrum of density perturbations without invoking superluminal expansion (inflation). The scenario relies, instead, on physical phenomena that arise naturally in theories based on extra dimensions and branes. As an example, we present our scenario predominantly within the context of heterotic M-theory. A prediction that distinguishes this scenario from standard inflationary cosmology is a strongly blue gravitational wave spectrum, which has consequences for microwave background polarization experiments and gravitational wave detectors.

We respond to the criticisms by Kallosh, Kofman and Linde concerning our proposal of the ekpyrotic universe scenario. We point out a number of errors in their considerations and argue that, at this stage, the ekpyrotic model is a possible alternative to inflationary cosmology as a description of the very early universe.

It is shown that there exist large classes of BPS vacua in heterotic M-theory which have negative tension on the visible orbifold plane, positive tension on the hidden plane and positive tension, physical five-branes in the bulk space. Explicit examples of such vacua are presented. Furthermore, it is demonstrated that the ratio, beta/|alpha|, of the bulk five-brane tension to the visible plane tension can, for several large classes of such vacua, be made arbitrarily small. Hence, it is straightforward to find vacua with the properties required in the examples of the Ekpyrotic theory of cosmology - a visible brane with negative tension and beta/|alpha| small. This contradicts recent claims in the literature.

We consider the cosmology of a pair of domain walls bounding a five-dimensional bulk space-time with negative cosmological constant, in which the distance between the branes is not fixed in time. Although there are strong arguments to suggest that this distance should be stabilized in the present epoch, no such constraints exist for the early universe and thus non-static solutions might provide relevant inflationary scenarios. We find the general solution for the standard ansatz where the bulk is foliated by planar-symmetric hypersurfaces. We show that in all cases the bulk geometry is that of anti-de Sitter (AdS_5). We then present a geometrical interpretation for the solutions as embeddings of two de Sitter (dS_4) surfaces in AdS_5, which provide a simple interpretation of the physical properties of the solutions. A notable feature explained in the analysis is that two-way communication between branes expanding away from one another is possible for a finite amount of time, after which communication can proceed in one direction only. The geometrical picture also shows that our class of solutions (and related solutions in the literature) are not completely general, contrary to some claims. We then derive the most general solution for two walls in AdS_5. This includes novel cosmologies where the brane tensions are not constrained to have opposite signs. The construction naturally generalizes to arbitrary FRW cosmologies on the branes.

We present a study of the scaling properties of bipolar outflows, generic to the circumstellar environment of young stellar objects, and which are generally believed to be magnetohydrodynamically driven. Our work consists of a multifractal analysis of near-infrared light scattered by dust grains in the lobes of seven bipolar nebulae, namely GGD 18, V380 Orionis, LkH\alpha 101/NGC 1579, LkH\alpha 233, PV Cephei, V645 Cygni (GL 2789), and V633 Cassiopeiae (LkH\alpha 198). We argue that these objects pertain to a multifractal universality class, with the corresponding universal parameters averaged over the ensemble yielding alpha= 1.96 +/- 0.02, C_1=0.04 +/- 0.02, and H=0.7 +/- 0.2; these results might suggest that the dynamics of the observed objects are similar. We then proceed to investigate the existence of anisotropy in the scaling of GGD 18 using the Generalized Scale Invariance (GSI) formalism. It is found that differential rotation and stratification are involved in the outflow mechanism. The stratification is believed to be dynamical in nature i.e., gravity is apparently not the dominant stratifying force.