With the advent of gravitational wave astronomy, searching for gravitational wave echoes from black holes (BHs) is becoming an interesting probe of their quantum nature near their horizons. Newborn BHs may be strong emitters of echoes, as they accompany large perturbations in the surrounding spacetime upon formation. Utilizing the Quantum Black Hole Seismology framework (Oshita et al. 2020), we study the expected echoes upon BH formation resulting from neutron star mergers and failed supernovae. For BH remnants from neutron star mergers, we evaluate the consistency of these models with the recent claim on the existence of echoes following the neutron star merger event GW170817. We find that the claimed echoes in GW170817, if real, suggest that overtones contribute a significant amount of energy in the ringdown of the remnant BH. We finally discuss the detectability of echoes from failed supernovae by second and third-generation gravitational wave detectors, and find that current (future) detectors constrain physical reflectivity models for events occurring within a few Mpc (a few x 10 Mpc). Detecting such echo signals may significantly constrain the maximum mass and equation of state of neutron stars.

Searches for gravitational wave echoes in the aftermath of mergers and/or formation of astrophysical black holes have recently opened a novel and surprising window into the quantum nature of their horizons. Similar to astro- and helioseismology, study of the spectrum of quantum black holes provides a promising method to understand their inner structure, what we call $\textit{quantum black hole seismology}$. We provide a detailed numerical and analytic description of this spectrum in terms of the properties of the Kerr spacetime and quantum black hole horizons, showing that it drastically differs from their classical counterparts. Our most significant findings are the following: (1) If the temperature of quantum black hole is $\lesssim 2 \times$ Hawking temperature, then it will not suffer from ergoregion instability (although the bound is looser at smaller spins). (2) We find how quantum black hole spectra pinpoint the microscopic properties of quantum structure. For example, the detailed spacing of spectral lines can distinguish whether quantum effects appear through compactness (i.e., exotic compact objects) or frequency (i.e., modified dispersion relation). (3) We find out that the overtone quasinormal modes may strongly enhance the amplitude of echo in the low-frequency region. (4) We show the invariance of the spectrum under the generalized Darboux transformation of linear perturbations, showing that it is a genuine covariant observable.

One of the most powerful tools in the arsenal of theoretical cosmologists is the halo model of large scale structure, which provides a phenomenological description of nonlinear structure in our universe. However, it is well known that there is no simple way to impose conservation laws in the halo model. This can severely impair the predictions on large scales for observables such as weak lensing or the kinematic Sunyaev-Zel'dovich effect, which should satisfy mass and momentum conservations, respectively. For example, the standard halo model overpredicts weak lensing power spectrum by $> 8\%$ on scales $>20$ degrees. To address this problem, we present an $Amended ~ Halo ~ Model$, explicitly separating the linear perturbations from $compensated$ halo profiles. This is guaranteed to respect conservation laws, as well as linear theory predictions on large scales. We then provide a simple fitting function for the compensated halo profiles, and discuss the modified predictions for 1-halo and 2-halo terms, as well as other cosmological observations such weak lensing power spectrum. Furthermore, we argue that the amended halo model provides a more accurate framework to capture physical effects that happen in the process of cosmological structure formation.

Characterizing the physical properties of the stochastic gravitational waves background (SGWB) is a key step towards identifying the nature of its possible origin. We focus our analysis on SGWB anisotropies. The existence of a non-trivial primordial scalar-tensor-tensor (STT) correlation in the squeezed configuration may be inferred from the effect that a long wavelength scalar mode has on the gravitational wave power spectrum: an anisotropic contribution. Crucially, such contribution is correlated with temperature anisotropies in the cosmic microwave background (CMB). We show that, for inflationary models that generate suitably large STT non-Gaussianities, cross-correlating the CMB with the stochastic background of gravitational waves is a very effective probe of early universe physics. The resulting signal can be a smoking-gun for primordial SGWB anisotropies.

The cold dark matter model predicts that dark matter haloes are connected by filaments. Direct measurements of the masses and structure of these filaments are difficult, but recently several studies have detected these dark-matter-dominated filaments using weak lensing. Here we study the efficiency of galaxy formation within the filaments by measuring their total mass-to-light ratios and stellar mass fractions. Specifically, we stack pairs of Luminous Red Galaxies (LRGs) with a typical separation on the sky of $8 h^{-1}$ Mpc. We stack background galaxy shapes around pairs to obtain mass maps through weak lensing, and we stack galaxies from the Sloan Digital Sky Survey (SDSS) to obtain maps of light and stellar mass. To isolate the signal from the filament, we construct two matched catalogues of physical and non-physical (projected) LRG pairs, with the same distributions of redshift and separation. We then subtract the two stacked maps. Using LRG pair samples from the BOSS survey at two different redshifts, we find that the evolution of the mass in filament is consistent with the predictions from perturbation theory. The filaments are not entirely dark: their mass-to-light ratios ($M/L = 351\pm87$ in solar units in the $r$-band) and stellar mass fractions ($M_{\rm stellar}/M = 0.0073\pm0.0020$) are consistent with the cosmic values (and with their redshift evolutions)

We study the reflectivity of quantum black hole (BH) horizons using detailed balance and fluctuation-dissipation theorem, finding a universal flux reflectivity given by the Boltzmann factor ${\mathcal R} = \exp\left(-{\hbar |\omega| \over k T_{\rm H}}\right)$, where $\omega$ is frequency in the horizon frame and $T_{\rm H}$ is the Hawking temperature. This implies CP-symmetry (or $\bf{RP}^3$ topology) of the extended BH spacetime. We then briefly discuss related physical implications: We predict echoes in the ringdown of Kerr BHs, but they do not exhibit ergoregion instability. The viscosity in the membrane paradigm is modified to $\eta = \frac{c^3}{16\pi G}\tanh\left({\hbar |\omega| \over 4 k T_{\rm H}}\right)$, only approaching General Relativistic value at high frequencies.

Black Holes are possibly the most enigmatic objects in our Universe. From their detection in gravitational waves upon their mergers, to their snapshot eating at the centres of galaxies, black hole astrophysics has undergone an observational renaissance in the past 4 years. Nevertheless, they remain active playgrounds for strong gravity and quantum effects, where novel aspects of the elusive theory of quantum gravity may be hard at work. In this review article, we provide an overview of the strong motivations for why "Quantum Black Holes" may be radically different from their classical counterparts in Einstein's General Relativity. We then discuss the observational signatures of quantum black holes, focusing on gravitational wave echoes as smoking guns for quantum horizons (or exotic compact objects), which have led to significant recent excitement and activity. We review the theoretical underpinning of gravitational wave echoes and critically examine the seemingly contradictory observational claims regarding their (non-)existence. Finally, we discuss the future theoretical and observational landscape for unraveling the "Quantum Black Holes in the Sky".

The first direct observation of a binary neutron star (BNS) merger was a watershed moment in multi-messenger astronomy. However, gravitational waves from GW170817 have only been observed prior to the BNS merger, but electromagnetic observations all follow the merger event. While post-merger gravitational wave signal in general relativity is too faint (given current detector sensitivities), here we present the first tentative detection of post-merger gravitational wave "echoes" from a highly spinning "black hole" remnant. The echoes may be expected in different models of quantum black holes that replace event horizons by exotic Planck-scale structure and tentative evidence for them has been found in binary black hole merger events. The fact that the echo frequency is suppressed by $\log M$ (in Planck units) puts it squarely in the LIGO sensitivity window, allowing us to build an optimal model-agnostic search strategy via cross-correlating the two detectors in frequency/time. We find a tentative detection of echoes at $f_{\rm echo} \simeq 72$ Hz, around 1.0 sec after the BNS merger, consistent with a 2.6-2.7 $M_\odot$ "black hole" remnant with dimensionless spin $0.84-0.87$. Accounting for all the "look-elsewhere" effects, we find a significance of $4.2 \sigma$, or a false alarm probability of $1.6\times 10^{-5}$, i.e. a similar cross-correlation within the expected frequency/time window after the merger cannot be found more than 4 times in 3 days. If confirmed, this finding will have significant consequences for both physics of quantum black holes and astrophysics of binary neutron star mergers [Note added: This result is independently confirmed by arXiv:1901.04138, who use the electromagnetic observations to infer $t_{\rm coll}=0.98^{+0.31}_{-0.26}$ sec for black hole formation].

Gravitational wave echoes provide our most direct and surprising observational window into quantum nature of black holes. Three years ago, the first search for echoes from Planck-scale modifications of general relativity near black hole event horizons led to tentative evidence at false detection probability of 1\% arXiv:1612.00266 . The study introduced a naive phenomenological model and used the public data release by the Advanced LIGO gravitational wave observatory for the first observing run O1 (GW150914, GW151226, and LVT151012, now GW151012). Here, we provide a status update on various observational searches for echoes by independent groups, and argue that they can all be consistent if echoes are most prominent at lower frequencies and/or in binary mergers of more extreme mass ratio. We also point out that the only reported "detection" of echoes (with $>4\sigma$ confidence) at 1.0 second after the binary neutron star merger GW170817 arXiv:1803.10454 is coincident with the formation time of the black hole inferred from electromagnetic observations.

One of the most ubiquitous features of quantum theories is the existence of zero-point fluctuations in their ground states. For massive quantum fields, these fluctuations decouple from infrared observables in ordinary field theories. However, there is no "decoupling theorem" in Quantum Gravity, and we recently showed that the vacuum stress fluctuations of massive quantum fields source a red spectrum of metric fluctuations given by $\sim$ mass$^5$/frequency in Planck units. I show that this signal is consistent with the reported unattributed persistent noise, or "mystery" noise, in the Laser Interferometer Gravitational-Wave Observatory (LIGO), for the Standard Model of Particle Physics. If this interpretation is correct, then it implies that: 1) This will be a fundamental irreducible noise for all gravitational wave interferometers, and 2) There is no fundamental weakly-coupled massive particle heavier than those in the Standard Model.

Most measurements of mass in Astronomy that use kinematics of stars or gas rely on assumptions of equilibrium that are often hard to verify. Instead, we develop a novel idea that uses the clustering in action space, as a probe of underlying gravitational potential: the correct potential should maximize small-scale clustering in the action space. We provide a first-principle derivation of likelihood using the two-point correlation function in action space, and test it against simulations of stellar streams. We then apply this method to the 2nd data release of Gaia, and use it to measure the radial force fraction $f_h$ and logarithmic slope $\alpha$ of dark matter halo profile. We investigate stars within 9-11 kpc and 11.5-15 kpc from Galactic centre, and find $(f_h,\alpha)= (0.391\pm 0.009, 1.835\pm 0.092) $ and $(0.351\pm 0.012,1.687\pm 0.079)$, respectively. We also confirm that the set of parameters that maximize the likelihood function do correspond to the most clustering in the action space. The best-fit circular velocity curve for Milky Way potential is consistent with past measurements (although it is $\sim$ 5-10\% lower than previous methods that use masers or globular clusters). Our work provides a clear demonstration of the full statistical power that lies in the full phase space information, relieving the need for {\it ad hoc} assumptions such as virial equilibrium, circular motion, or steam-finding algorithms.

It has been suggested that the cosmic history might repeat in cycles, with an infinite series of similar aeons in the past and the future. Here, we instead propose that the cosmic history repeats itself exactly, constructing a universe on a periodic temporal history, which we call Periodic Time Cosmology. In particular, the primordial power spectrum, convolved with the transfer function throughout the cosmic history, would form the next aeon's primordial power spectrum. By matching the big bang to the infinite future using a conformal rescaling (a la Penrose), we uniquely determine the primordial power spectrum, in terms of the transfer function up to two free parameters. While nearly scale invariant with a red tilt on large scales, using Planck and Baryonic Acoustic Oscillation observations, we find the minimal model is disfavoured compared to a power-law power spectrum at $5.1\sigma$. However, extensions of $\Lambda$CDM cosmic history change the large scale transfer function and can provide better relative fits to the data. For example, the best fit seven parameter model for our Periodic Time Cosmology, with $w=-1.024$ for dark energy equation of state, is only disfavoured relative to a power-law power spectrum (with the same number of parameters) at $1.8\sigma$ level. Therefore, consistency between cosmic history and initial conditions provides a viable description of cosmological observations in the context of Periodic Time Cosmology.

We propose a cosmological model, \"u$\Lambda$CDM, based on {\it \"uber-gravity}, which is a canonical ensemble average of many theories of gravity. In this model, we have a sharp transition from (a purely) $\Lambda$CDM era to a phase in which the Ricci scalar is a constant. This transition occurs when the Ricci scalar reaches a critical scale or alternatively at a transition redshift, $z_\oplus$. We use the observations of baryonic acoustic oscillations (BAO) and Supernovae Ia (SNe), as well as the cosmic microwave background (CMB) data to constrain \"u$\Lambda$CDM. This yields $H_0=70.6_{-1.3}^{+1.1}$ km/s/Mpc, $\Omega_m=0.2861 \pm{0.0092} $ and $z_\oplus=0.537_{-0.375}^{+0.277}$, providing a marginally better fit with a Akaike information criterion of 0.8. %After adding in the data from BAO in the Lyman-$\alpha$ forest in quasar spectra, we get $\Delta\chi^2=-2.14$, where the best fit occurs at $H_0=71.2_{-1.5}^{+1.8}$ km/s/Mpc, $\Omega_m=0.280\pm 0.013$ and $z_{\oplus}= 0.41^{+0.31}_{-0.29}$. Therefore, \"u$\Lambda$CDM can ease the $H_0$-tension, albeit marginally, with one additional free parameter. We also provide a preliminary study of the linear perturbation theory in \"u$\Lambda$CDM which points to interesting potential {\it smoking guns} in the observations of large scales structure at $z < z_{\oplus}$.

The discrepancy between expected and observed cooling rates of X-ray emitting gas has led to the {\it cooling flow problem} at the cores of clusters of galaxies. A variety of models have been proposed to model the observed X-ray spectra and resolve the cooling flow problem, which involves heating the cold gas through different mechanisms. As a result, realistic models of X-ray spectra of galaxy clusters need to involve both heating {\it and} cooling mechanisms. In this paper, we argue that the heating time-scale is set by the magnetohydrodynamic (MHD) turbulent viscous heating for the Intracluster plasma, parametrised by the Shakura-Sunyaev viscosity parameter, $\alpha$. Using a cooling+heating flow model, we show that a value of $\alpha\simeq 0.05$ (with 10\% scatter) provides improved fits to the X-ray spectra of cooling flow, while at the same time, predicting reasonable cooling efficiency, $\epsilon_{cool} = 0.33^{+0.63}_{-0.15}$. Our inferred values for $\alpha$ based on X-ray spectra are also in line with direct measurements of turbulent pressure in simulations and observations of galaxy clusters. This simple picture unifies astrophysical accretion, as a balance of MHD turbulent heating and cooling, across more than 16 orders of magnitudes in scale, from neutron stars to galaxy clusters.

One of the most triumphant predictions of the theory if general relativity was the recent LIGO-Virgo detection of gravitational wave (GW) signals produced in binary black hole (BH) mergers. However, it is suggested that exotic compact objects, proposed in quantum gravity models of BHs, may produce similar classical GW waveforms, followed by delayed repeating "echoes". In a companion paper [1], we have presented different arguments for a universal Boltzmann reflectivity of quantum BH horizons. Here, we investigate the resulting echoes from this prescription. We derive corresponding quasi-normal modes (QNMs) for quantum BHs analytically, and show how their initial conditions can be related to the QNMs of classical BHs. Ergoregion instability is suppressed by the imperfect reflectivity. We then compare the analytic and numerical predictions for echoes in real time, verifying their consistency. In particular, we find that the amplitudes of the first ~20 echoes decay inversely with time, while the subsequent echoes decay exponentially. Finally, we present predictions for the signal-to-noise ratio of echoes for spinning BHs, which should be imminently detectable for massive remnants, subject to the uncertainty in the nonlinear initial conditions of the BH merger.

We investigate a Bose-Einstein condensate (BEC) as a gravitational wave detector, and study its sensitivity by optimizing the properties of the condensate and the measurement duration. We show that detecting kilohertz gravitational waves is limited by current experimental techniques in squeezing BEC phonons. Future improvements in technology to squeeze BEC states can make them competitive detectors for gravitational waves of astrophysical and/or cosmological origin.

Quantum nature of black hole horizons has been a subject of recent interest and scrutiny. In particular, a near-horizon quantum violation of the equivalence principle has been proposed as a resolution of the black hole information paradox. Such a violation may lead to a modified dispersion relation at high energies, which could become relevant due to the intense gravitational blueshift experienced by ingoing gravitational waves. We investigate the ringdown for a perturbed black hole with such a modified dispersion relation and find that infalling gravitational waves are partially reflected near the horizon. This results in the appearance of late-time ${\it echoes}$ in the ringdown phase of black hole merger events, with similar properties to those (arguably) seen in the Advanced LIGO observations. Current measurements suggest a Lorentz-violation scale of $10^{13 \pm 2}$ GeV for gravitational waves, with comparable dissipation and dispersion. Therefore, if confirmed, black hole ringdown echoes probe the microstructure of horizons and thus can test Lorentz-violating UV completions, such as in Ho$\v{r}$ava-Lifshitz gravity.

Individual extragalactic dark matter halos, such as those associated with nearby galaxies and galaxy clusters, are promising targets for searches for gamma-rays from dark matter annihilation. We review the predictions for the annihilation flux from individual halos, focusing on the effect of current uncertainties in the concentration-mass relation and the contribution from halo substructure, and also estimating the intrinsic halo-to-halo scatter expected. After careful consideration of recent simulation results, we conclude that the concentrations of the smallest halos, while well determined at high redshift, are still uncertain by a factor of 4-6 when extrapolated to low redshift. This in turn produces up to two orders of magnitude uncertainty in the predicted annihilation flux for any halo mass above this scale. Substructure evolution, the small-scale cutoff to the power spectrum, cosmology, and baryonic effects all introduce smaller, though cumulative, uncertainties. We then consider intrinsic variations from halo to halo. These arise from variations in concentration and substructure, leading to a scatter of $\sim$ 2.5 in the predicted annihilation luminosity. Finally, we consider the problem of detecting gamma-rays from annihilation, given the expected contributions from other sources. We estimate the signal-to-noise ratio for gamma-ray detection as a function of halo mass, assuming that cosmic rays from star formation are the main noise source in the detection. This calculation suggests that group-scale halos, individually or in stacks, may be a particularly interesting target for the next generation of annihilation searches.

We devise an optimal method to measure the temporal power spectrum of the lensing and intrinsic fluctuations of multiply-imaged strongly lensed quasar light curves, along with the associated time delays. The method is based on a Monte-Carlo Markov Chain (MCMC) sampling of a putative gaussian likelihood, and accurately recovers the input properties of simulated light curves, as well as the "Time Delay Challenge". We apply this method to constrain the dimensionless cosmological (non-linear) matter power spectrum on milliparsec scales (comparable to the size of the solar system), to $\Delta_{\rm NL}^2< 4 \times 10^7$ at $k_{\rm NL} \sim 10^3 {\rm pc}^{-1}$. Using a semi-analytic nonlinear clustering model which is calibrated to simulations, the corresponding constraint on the primordial (linear) scalar power spectrum is ${\cal P}_{\cal R} < 3 \times 10^{-9}$ at $k_{\rm L}\sim$ 3 pc$^{-1}$. This is the strongest constraint on primordial power spectrum at these scales, and is within an order of magnitude from the standard $\Lambda$CDM prediction. We also report measurements of temporal spectra for intrinsic variabilities of quasar light curves, which can be used to constrain the size of the emitting region in accretion disks. Future cadenced optical imaging surveys, such as LSST, should increase the number of observed strongly lensed quasars by 3 orders of magnitude and significantly improve these measurements, even though improvements in modelling quasar accretion and stellar microlensing are necessary.

While recent detections of gravitational waves from the mergers of binary black holes match well with the predictions of General Relativity (GR), they cannot directly confirm the existence of event horizons. Exotic compact objects (ECOs) are motivated by quantum models of black holes, and can have exotic structure (or a "wall") just outside the (would-be) horizon. ECOs produce similar ringdown waveforms to the GR black holes, but they are followed by delayed "echoes". By solving linearized Einstein equations we can model these echoes and provide analytic templates that can be used to compare to observations. For concreteness, we consider GW150914 event, detected by the LIGO/Virgo collaboration, and study the model dependence of its echo properties. We find that echoes are reasonably approximated by complex gaussians, with amplitudes that decay as a power law in time, while their width in time (frequency) grows (shrinks) over subsequent echoes. We also show that trapped modes between a perfectly reflecting wall and angular momentum barrier in Kerr metric can exhibit superradiant instability over long times, as expected.

Everpresent $\Lambda$ is a cosmological scenario in which the observed cosmological "constant" $\Lambda$ fluctuates between positive and negative values with a vanishing mean, and with a magnitude comparable to the critical density at any epoch. In accord with a longstanding heuristic prediction of causal set theory, it postulates that $\Lambda$ is a stochastic function of cosmic time that will vary from one realization of the scenario to another. Herein, we consider two models of "dark energy" that exhibit these features. Via Monte Carlo Markov chains, we explore the space of cosmological parameters and the set of stochastic realizations of these models, finding that Everpresent $\Lambda$ can fit the current cosmological observations as well as the $\Lambda$CDM model does. Furthermore, it removes observational tensions with $\Lambda$CDM, for low redshift measurements of Hubble constant, and the Baryonic Acoustic Oscillations (BAO) in Lyman-$\alpha$ forest at $z\sim 2-3$. However, we also find that Everpresent $\Lambda$ does not significantly help with the growth of ultramassive black holes at high redshift, and the Lithium problem in Big Bang Nuclesynthesis. Future measurements of "dark energy" at high redshifts will further test the viability of Everpresent $\Lambda$ as an alternative to the $\Lambda$CDM cosmology.

In a recent publication (1612.00266), we demonstrated that the events in the first observing run of the Advanced LIGO gravitational wave observatory (aLIGO O1) showed tentative evidence for repeating "echoes from the abyss" caused by Planck-scale structure near black hole horizons. By considering a phenomenological echo model, we showed that the pure noise hypothesis is disfavored with a p-value of 1%, i.e. higher amplitude for echoes than those in aLIGO O1 events are only recovered in 1% of random noise realizations. A recent preprint by Westerweck, et al. (1712.09966) provides a careful re-evaluation of our analysis which claims "a reduced statistical significance ... entirely consistent with noise". It is a mystery to us why the authors make such a statement, while they also find a p-value of 2 $\pm$ 1% (given the Poisson error in their estimate) for the same model and dataset. This is p-erfectly consistent with our results, which would be commonly considered as disfavoring the null hypothesis, or "moderate to significant" evidence for "echoes". Westerweck, et al. also point to diversity of the observed echo properties as evidence for statistical fluke, but such a diversity is neither unique nor surprising for complex physical phenomena in nature.

In classical General Relativity (GR), an observer falling into an astrophysical black hole is not expected to experience anything dramatic as she crosses the event horizon. However, tentative resolutions to problems in quantum gravity, such as the cosmological constant problem, or the black hole information paradox, invoke significant departures from classicality in the vicinity of the horizon. It was recently pointed out that such near-horizon structures can lead to late-time echoes in the black hole merger gravitational wave signals that are otherwise indistinguishable from GR. We search for observational signatures of these echoes in the gravitational wave data released by advanced Laser Interferometer Gravitational-Wave Observatory (LIGO), following the three black hole merger events GW150914, GW151226, and LVT151012. In particular, we look for repeating damped echoes with time-delays of $8 M \log M$ (+spin corrections, in Planck units), corresponding to Planck-scale departures from GR near their respective horizons. Accounting for the "look elsewhere" effect due to uncertainty in the echo template, we find tentative evidence for Planck-scale structure near black hole horizons at false detection probability of $1\%$ (corresponding to $2.5\sigma$ significance level). Future observations from interferometric detectors at higher sensitivity, along with more physical echo templates, will be able to confirm (or rule out) this finding, providing possible empirical evidence for alternatives to classical black holes, such as in ${\it firewall}$ or ${\it fuzzball}$ paradigms.

Standard accretion disk theory (Shakura & Sunyaev 1973) predicts that the total pressure in disks at typical (sub-)Eddington accretion rates becomes radiation pressure dominated. However, radiation pressure dominated disks are thermally unstable. Since these disks are observed in approximate steady state over the instability time-scale, our accretion models in the radiation pressure dominated regime (i.e. inner disk) need to be modified. Here, we present a modification to the SS model, where radiation pressure is in equipartition with gas pressure in the inner region. We call these flows Accretion in Radiative Equipartition (AiRE) Disks. We introduce the basic features of AiRE disks and show how they modify disk properties such as the Toomre parameter and central temperature. We then show that the accretion rate of AiRE disks is limited from above and below, by Toomre and nodal sonic point instabilities, respectively. The former leads to a strict upper limit on the mass of supermassive black holes as a function of cosmic time (and spin), while the latter could explain the transition between hard and soft states of X-ray binaries.

We introduce a new family of primordial cosmological perturbations that are not described by traditional power spectra. At the linear level, these perturbations live in the kernel of the spatial Laplacian operator, and thus we call them cosmological zero modes. We compute the cosmic microwave background (CMB) temperature and polarization anisotropy induced by these modes, and forecast their detection sensitivity using a cosmic-variance limited experiment. In particular, we consider two configurations for the zero modes: The first configuration consists of stochastic metric perturbations described by white noise on a "holographic screen" located at our cosmological horizon. The amplitude of the power spectrum of this white noise can be constrained to be $\lesssim 9 \times 10^{-14}$. The second configuration is a primordial monopole beyond our cosmological horizon. We show that such a monopole, with "charge" $Q$, can be detected in the CMB sky up to a distance of $11.6 ~ Q^{1/4}\times$ horizon radius (or $160~ Q^{1/4}$ Gpc). More generally, observational probes of cosmological zero modes can shed light on non-perturbative phenomena in the primordial universe, beyond our observable horizon.

Holographic cosmology offers a novel framework for describing the very early Universe in which cosmological predictions are expressed in terms of the observables of a three dimensional quantum field theory (QFT). This framework includes conventional slow-roll inflation, which is described in terms of a strongly coupled QFT, but it also allows for qualitatively new models for the very early Universe, where the dual QFT may be weakly coupled. The new models describe a universe which is non-geometric at early times. While standard slow-roll inflation leads to a (near-)power-law primordial power spectrum, perturbative superrenormalizable QFT's yield a new holographic spectral shape. Here, we compare the two predictions against cosmological observations. We use CosmoMC to determine the best fit parameters, and MultiNest for Bayesian Evidence, comparing the likelihoods. We find that the dual QFT should be non-perturbative at the very low multipoles ($l \lesssim 30$), while for higher multipoles ($l \gtrsim 30$) the new holographic model, based on perturbative QFT, fits the data just as well as the standard power-law spectrum assumed in $\Lambda$CDM cosmology. This finding opens the door to applications of non-perturbative QFT techniques, such as lattice simulations, to observational cosmology on gigaparsec scales and beyond.

Standard big bang cosmology predicts a cosmic neutrino background at $T_\nu \simeq 1.95~$K. Given the current neutrino oscillation measurements, we know most neutrinos move at large, but non-relativistic, velocities. Therefore, dark matter haloes moving in the sea of primordial neutrinos form a neutrino wake behind them, which would slow them down, due to the effect of {\it dynamical friction}. In this paper, we quantify this effect for realistic haloes, in the context of the halo model of structure formation, and show that it scales as $m_\nu^4$ $\times$ relative velocity, and monotonically grows with the halo mass. Galaxy redshift surveys can be sensitive to this effect (at $>3\sigma$ confidence level, depending on survey properties, neutrino mass and hierarchy) through redshift space distortions (RSD) of distinct galaxy populations.

We argue that massive quantum fields source low-frequency long-wavelength metric fluctuations through the quantum fluctuations of their stress-energy, given reasonable assumptions about the analytic structure of its correlators. This can be traced back to the non-local nature of the gauge symmetry in General Relativity, which prevents an efficient screening of UV scales (what we call the cosmological non-constant problem). We define a covariant and gauge-invariant observable which probes line-of-sight spacetime curvature fluctuations on an observer's past lightcone, and show that current pulsar timing data constrains any massive particle to $m\lesssim 600$ GeV. This astrophysical bound severely limits the possibilities for physics beyond the standard model below the scale of quantum gravity.

The 5D Holographic Big Bang is a novel model for the emergence of the early universe out of a 5D collapsing star (an apparent white hole), in the context of Dvali-Gabadadze-Porrati (DGP) cosmology. The model does not have a big bang singularity, and yet can address cosmological puzzles that are traditionally solved within inflationary cosmology. In this paper, we compute the exact power spectrum of cosmological curvature perturbations due to the effect of a thin atmosphere accreting into our 3-brane. The spectrum is scale-invariant on small scales and red on intermediate scales, but becomes blue on scales larger than the height of the atmosphere. While this behaviour is broadly consistent with the non-parametric measurements of the primordial scalar power spectrum, it is marginally disfavoured relative to a simple power law (at 2.7$\sigma$ level). Furthermore, we find that the best fit nucleation temperature of our 3-brane is at least 3 orders of magnitude larger than the 5D Planck mass, suggesting an origin in a 5D quantum gravity phase.

In a recent paper (arXiv:1612.00266), we reported the results of the first search for echoes from Planck-scale modifications of general relativity near black hole event horizons using the public data release by the Advanced LIGO gravitational wave observatory. While we found tentative evidence (at $\simeq 3 \sigma$ level) for the presence of these echoes, our statistical methodology was challenged by Ashton, et al. (arXiv:1612.05625), just in time for the holidays! In this short note, we briefly address these criticisms, arguing that they either do not affect our conclusion or change its significance by $\lesssim 0.3\sigma$. The real test will be whether our finding can be reproduced by independent groups using independent methodologies (and ultimately more data).

We test a class of holographic models for the very early universe against cosmological observations and find that they are competitive to the standard $\Lambda$CDM model of cosmology. These models are based on three dimensional perturbative super-renormalizable Quantum Field Theory (QFT), and while they predict a different power spectrum from the standard power-law used in $\Lambda$CDM, they still provide an excellent fit to data (within their regime of validity). By comparing the Bayesian evidence for the models, we find that $\Lambda$CDM does a better job globally, while the holographic models provide a (marginally) better fit to data without very low multipoles (i.e. $l\lesssim 30$), where the dual QFT becomes non-perturbative. Observations can be used to exclude some QFT models, while we also find models satisfying all phenomenological constraints: the data rules out the dual theory being Yang-Mills theory coupled to fermions only, but allows for Yang-Mills theory coupled to non-minimal scalars with quartic interactions. Lattice simulations of 3d QFT's can provide non-perturbative predictions for large-angle statistics of the cosmic microwave background, and potentially explain its apparent anomalies.

One of the most fundamental features of a black hole in general relativity is its event horizon: a boundary from which nothing can escape. There has been a recent surge of interest in the nature of these event horizons and their local neighbourhoods. In an attempt to resolve black hole information paradox(es), and more generally, to better understand the path towards quantum gravity, "firewalls" have been proposed as an alternative to black hole event horizons. In this paper, we explore the phenomenological implications of black holes possessing a surface or "firewall", and predict a potentially detectable signature of these firewalls in the form of a high energy astrophysical neutrino flux. We compute the spectrum of this neutrino flux in different models and show that it is a possible candidate for the source of the PeV neutrinos recently detected by IceCube. This opens up a new area of research, bridging the non-perturbative physics of quantum gravity with the observational black hole and high energy astrophysics.

We explore the space of scalar-tensor theories containing two non-conformalmetrics, and find a discontinuity pointing to a "critical" cosmological solution. Due to the different maximal speeds of propagation for matter and gravity,the cosmological fluctuations start off inside the horizon even without inflation, and will more naturally have a thermal origin (since there is never vacuum domination). The critical model makes an unambiguous, non-tuned prediction for the spectral index of the scalar fluctuations: $n_S= 0.96478(64)$. Considering also that no gravitational waves are produced, we have unveiled the most predictive model on offer. The model has a simple geometrical interpretation as a probe 3-brane embedded in an $EAdS_2\times E_3$ geometry.

We study the influence of the fluctuations of a Lorentz invariant and conserved vacuum on cosmological metric perturbations, and show that they generically blow up in the IR. We compute this effect using the K\"all\'en-Lehmann spectral representation of stress correlators in generic quantum field theories, as well as the holographic bound on their entanglement entropy, both leading to an IR cut-off that scales as the fifth power of the highest UV scale (in Planck units). One may view this as analogous to the Heisenberg uncertainty principle, which is imposed on the phase space of gravitational theories by the Einstein constraint equations. The leading effect on cosmological observables come from anisotropic vacuum stresses which imply: i) any extension of the standard model of particle physics can only have masses (or resonances) $\lesssim$ 24 TeV, and ii) perturbative quantum field theory or quantum gravity becomes strongly coupled beyond a cut-off scale of $\Lambda\lesssim1$ PeV. Such a low cut-off is independently motivated by the Higgs hierarchy problem. This result, which we dub the cosmological non-constant problem, can be viewed as an extension of the cosmological constant (CC) problem, demonstrating the non-trivial UV-IR coupling and (yet another) limitation of effective field theory in gravity. However, it is more severe than the old CC problem, as vacuum fluctuations cannot be tuned to cancel due to the positivity of spectral densities or entropy. We thus predict that future advances in cosmological observations and collider technology will sandwich from above and below, and eventually discover, new (non-perturbative) physics beyond the Standard Model within the TeV-PeV energy range.

We study a novel proposal for the origin of cosmological cold dark matter (CDM) which is rooted in the quantum nature of spacetime. In this model, off-shell modes of quantum fields can exist in asymptotic states as a result of spacetime nonlocality (expected in generic theories of quantum gravity), and play the role of CDM, which we dub off-shell dark matter (OfDM). However, their rate of production is suppressed by the scale of non-locality (e.g. Planck length). As a result, we show that OfDM is only produced in the first moments of big bang, and then effectively decouples (except through its gravitational interactions). We examine the observational predictions of this model: In the context of cosmic inflation, we show that this proposal relates the reheating temperature to the inflaton mass, which narrows down the uncertainty in the number of e-foldings of specific inflationary scenarios. We also demonstrate that OfDM is indeed cold, and discuss potentially observable signatures on small scale matter power spectrum.

We have recently introduced a novel statistical measure of dark matter clustering in phase space, the particle phase space average density ($P^2SAD$). In a two-paper series, we studied the structure of $P^2SAD$ in the Milky-Way-size Aquarius haloes, constructed a physically motivated model to describe it, and illustrated its potential as a powerful tool to predict signals sensitive to the nanostructure of dark matter haloes. In this letter, we report a remarkable universality of the clustering of dark matter in phase space as measured by $P^2SAD$ within the subhaloes of host haloes across different environments covering a range from dwarf-size to cluster-size haloes ($10^{10}-10^{15}$ M$_\odot$). Simulations show that the universality of $P^2SAD$ holds for more than 7 orders of magnitude, over a 2D phase space, covering over 3 orders of magnitude in distance/velocity, with a simple functional form that can be described by our model. Invoking the universality of $P^2SAD$, we can accurately predict the non-linear power spectrum of dark matter at small scales all the way down to the decoupling mass limit of cold dark matter particles. As an application, we compute the subhalo boost to the annihilation of dark matter in a wide range of host halo masses.

The smallest dark matter haloes are the first objects to form in the hierarchical structure formation of cold dark matter (CDM) cosmology and are expected to be the densest and most fundamental building blocks of CDM structures in our universe. Nevertheless, the physical characteristics of these haloes have stayed illusive, as they remain well beyond the current resolution of N-body simulations (at redshift zero). However, they dominate the predictions (and uncertainty) in expected dark matter annihilation signal, amongst other astrophysical observables. Using the conservation of total energy and the ellipsoidal collapse framework, we can analytically find the mean and scatter of concentration $c$ and 1-D velocity dispersion $\sigma_{\rm 1d}$ for haloes of different virial mass $M_{200}$. Both $c$ and $\sigma_{\rm 1d}/M_{200}^{1/3}$ are in good agreement with numerical results within the regime probed by simulations -- slowly decreasing functions of mass that approach constant values at large masses. In particular, the predictions for the 1-D velocity dispersion of cluster mass haloes are surprisingly robust as the inverse heat capacity of cosmological haloes crosses zero at $M_{200} \sim 10^{14} M_\odot$. However, we find that current extrapolations from simulations to smallest CDM haloes dramatically depend on the assumed profile (e.g. NFW vs. Einasto) and fitting function, which is why theoretical considerations, such as the one presented here, can significantly constrain the range of feasible predictions.

Is Cosmic Censorship special to General Relativity, or can it survive a violation of equivalence principle? Recent studies have shown that singularities in Lorentz violating Einstein-Aether (or Horava-Lifhsitz) theories can lie behind a universal horizon in simple black hole spacetimes. Even infinitely fast signals cannot escape these universal horizons. We extend this result, for an incompressible aether, to 3+1d dynamical or spinning spacetimes which possess inner killing horizons, and show that a universal horizon always forms in between the outer and (would-be) inner horizons. This finding suggests a notion of Cosmic Censorship, given that geometry in these theories never evolves beyond the universal horizon (avoiding potentially singular inner killing horizons). A surprising result is that there are 3 distinct possible stationary universal horizons for a spinning black hole, only one of which matches the dynamical spherical solution. This motivates dynamical studies of collapse in Einstein-Aether theories beyond spherical symmetry, which may reveal instabilities around the spherical solution.

Despite the tremendous empirical success of equivalence principle, there are several theoretical motivations for existence of a preferred reference frame (or aether) in a consistent theory of quantum gravity. However, if quantum gravity had a preferred reference frame, why would high energy processes enjoy such a high degree of Lorentz symmetry? While this is often considered as an argument against aether, here I provide three independent arguments for why perturbative unitarity (or weak coupling) of the Lorentz-violating effective field theories put stringent constraints on possible observable violations of Lorentz symmetry at high energies. In particular, the interaction with the scalar graviton in a consistent low-energy theory of gravity and a (radiatively and dynamically) stable cosmological framework, leads to these constraints. The violation (quantified by the relative difference in maximum speed of propagation) is limited to $\lesssim 10^{-10} E({\rm eV})^{-4}$ (superseding all current empirical bounds), or the theory will be strongly coupled beyond meV scale. The latter happens in extended Horava-Lifshitz gravities, as a result of a previously ignored quantum anomaly. Finally, given that all cosmologically viable theories with significant Lorentz violation appear to be strongly coupled beyond meV scale, we conjecture that, similar to color confinement in QCD, or Vainshetin screening for massive gravity, high energy theories (that interact with gravity) are shielded from Lorentz violation (at least, up to the scale where gravity is UV-completed). In contrast, microwave or radio photons, cosmic background neutrinos, or gravitational waves may provide more promising candidates for discovery of violations of Lorentz symmetry.

A century after the advent of Quantum Mechanics and General Relativity, both theories enjoy incredible empirical success, constituting the cornerstones of modern physics. Yet, paradoxically, they suffer from deep-rooted, so-far intractable, conflicts. Motivations for violations of the notion of relativistic locality include the Bell's inequalities for hidden variable theories, the cosmological horizon problem, and Lorentz-violating approaches to quantum geometrodynamics, such as Horava-Lifshitz gravity. Here, we explore a recent proposal for a "real ensemble" non-local description of quantum mechanics, in which "particles" can copy each others' observable values AND phases, independent of their spatial separation. We first specify the exact theory, ensuring that it is consistent and has (ordinary) quantum mechanics as a fixed point, where all particles with the same values for a given observable have the same phases. We then study the stability of this fixed point numerically, and analytically, for simple models. We provide evidence that most systems (in our study) are locally stable to small deviations from quantum mechanics, and furthermore, the phase variance per value of the observable, as well as systematic deviations from quantum mechanics, decay as $\sim$ (Energy$\times$Time)$^{-2n}$, where $n \geq 1$. Interestingly, this convergence is controlled by the absolute value of energy (and not energy difference), suggesting a possible connection to gravitational physics. Finally, we discuss different issues related to this theory, as well as potential novel applications for the spectrum of primordial cosmological perturbations and the cosmological constant problem.

It has been clear for some time now that super-critical surface magnetic fields, exceeding 4 x 10^13 G, exist on a subset of neutron stars. These magnetars may harbor interior fields many orders of magnitude larger, potentially reaching equipartition values. However, the impact of these strong fields on stellar structure has been largely ignored, potentially complicating attempts to infer the high density nuclear equation of state. Here we assess the effect of these strong magnetic fields on the mass-radius relationship of neutron stars. We employ an effective field theory model for the nuclear equation of state that includes the impact of hyperons, anomalous magnetic moments, and the physics of the crust. We consider two magnetic field geometries, bounding the likely magnitude of the impact of magnetic fields: a statistically isotropic, tangled field and a force-free configuration. In both cases even equipartition fields have at most a 30% impact on the maximum mass. However, the direction of the effect of the magnetic field depends on the geometry employed - force-free fields leading to reductions in the maximum neutron star mass and radius while tangled fields increase both - challenging the common intuition in the literature on the impact of magnetic fields.

We investigate pair creation by an electric field in four-dimensional de Sitter space. The expectation value of the induced current is computed, using the method of adiabatic regularization. Under strong electric fields the behavior of the current is similar to that in flat space, while under weak electric fields the current becomes inversely proportional to the mass squared of the charged field. Thus we find that the de Sitter space obtains a large conductivity under weak electric fields in the presence of a charged field with a tiny mass. We then apply the results to constrain electromagnetic fields in the early universe. In particular, we study cosmological scenarios for generating large-scale magnetic fields during the inflationary era. Electric fields generated along with the magnetic fields can induce sufficiently large conductivity to terminate the phase of magnetogenesis. For inflationary magnetogenesis models with a modified Maxwell kinetic term, the generated magnetic fields cannot exceed 10^{-30} G on Mpc scales in the present epoch, when a charged field carrying an elementary charge with mass of order the Hubble scale or smaller exists in the Lagrangian. Similar constraints from the Schwinger effect apply for other magnetogenesis mechanisms.

We show that the generalized McVittie spacetime, which represents a black hole with time-dependent mass in an expanding universe, is an exact solution of a subclass of the Horndeski family of actions. The heat-flow term responsible for the energy transfer between the black hole and the cosmological background is generated by the higher-order kinetic gravity braiding term, which generalizes the cuscuton action that yields McVittie with constant mass as a solution. Finally, we show that this generalization can be understood in terms of a duality realized by a disformal transformation, connecting the cuscuton field theory to an extension of the Horndeski action which does not propagate any scalar degrees of freedom. Our finding opens a novel window into studies of non-trivial interactions between dark energy/modified gravity theories and astrophysical black holes.

An important and long-standing puzzle in the history of modern physics is the gross inconsistency between theoretical expectations and cosmological observations of the vacuum energy density, by at least 60 orders of magnitude, otherwise known as the \textit{cosmological constant problem}. A characteristic feature of vacuum energy is that it has a pressure with the same amplitude, but opposite sign to its energy density, while all the precision tests of General Relativity are either in vacuum, or for media with negligible pressure. Therefore, one may wonder whether an anomalous coupling to pressure might be responsible for decoupling vacuum from gravity. We test this possibility in the context of the \textit{Gravitational Aether} proposal, using current cosmological observations, which probe the gravity of relativistic pressure in the radiation era. Interestingly, we find that the best fit for anomalous pressure coupling is about half-way between General Relativity (GR), and Gravitational Aether (GA), if we include \textit{Planck} together with \textit{WMAP} and BICEP2 polarization cosmic microwave background (CMB) observations. Taken at face value, this data combination excludes both GR and GA at around the 3-sigma level. However, including higher resolution CMB observations ("highL") or baryonic acoustic oscillations (BAO) pushes the best fit closer to GR, excluding the Gravitational Aether solution to the cosmological constant problem at the 4--5-sigma level. This constraint effectively places a limit on the anomalous coupling to pressure in the parametrized post-Newtonian (PPN) expansion, $\zeta_4 = 0.105 \pm 0.049$ (+highL CMB), or $\zeta_4 = 0.066 \pm 0.039$ (+BAO). These represent the most precise measurement of this parameter to date, indicating a mild tension with GR (for $\Lambda$CDM including tensors, with $\zeta_4=0$), and also among different data sets.

Motivations for the existence of a fundamental preferred frame range from pure phenomenology to attempts to solve the non-renormalizability of quantum gravity, the problem of time (and scale), and the cosmological constant problem(s). In many explicit constructions, such as Einstein-Aether or Gravitational Aether theories, K-essence, Cuscuton theory, Shape Dynamics, or (non-projectable) Horava-Lifshitz gravity, the low energy theory contains a fluid (which defines a preferred frame) with superluminal or incompressible excitations. We study here the formation of black holes in the presence of such a fluid. In particular, we focus on the incompressible limit of the fluid (or Constant Mean Curvature foliation) in the space-time of a spherically collapsing shell within an asymptotically cosmological space-time. In this case, ignoring the fluid back reaction, we can analytically show that an observer inside 3/4 of the Schwarzschild radius cannot send a signal outside, after a stage in collapse, even using signals that propagate infinitely fast in the preferred frame. This confirms the dynamical emergence of universal horizons that have been previously found in static solutions. We argue that this universal horizon should be considered as the future boundary of the classical space-time.

We present a novel perspective on the clustering of dark matter in phase space by defining the particle phase space average density ($P^2SAD$) as a two-dimensional extension of the two-point correlation function averaged within a certain volume in phase space. This statistics is a sensitive measure of small scale (sub-)structure of dark matter haloes. By analysing the structure of $P^2SAD$ in Milky-Way-size haloes using the Aquarius simulations, we find it to be nearly universal at small scales, i.e. small separations in phase space, where substructures dominate. This remarkable universality occurs across time and in regions of substantially different ambient densities (by nearly four orders of magnitude), with typical variations in $P^2SAD$ of a factor of a few. The maximum variations occur in regions where substructures have been strongly disrupted. The universality is also preserved across haloes of similar mass but diverse mass accretion histories and subhalo distributions. The universality is also broken at large scales, where the smooth dark matter distribution in the halo dominates. Although at small scales the structure of $P^2SAD$ is roughly described by a subhalo model, we argue that the simulation data is better fitted by a family of superellipse contours. This functional shape is inspired by a model that extends the stable clustering hypothesis into phase space. In a companion paper, we refine this model and show its advantages as a method to obtain predictions for non-gravitational signatures of dark matter.

We present a model for the structure of the particle phase space average density ($P^2SAD$) in galactic haloes, introduced recently as a novel measure of the clustering of dark matter. Our model is based on the stable clustering hypothesis in phase space, the spherical collapse model, and tidal disruption of substructures, which is calibrated against the Aquarius simulations. Using this model, we can predict the behaviour of $P^2SAD$ in the numerically unresolved regime, down to the decoupling mass limit of generic WIMP models. This prediction can be used to estimate signals sensitive to the small scale structure of dark matter. For example, the dark matter annihilation rate can be estimated for arbitrary velocity-dependent cross sections in a convenient way using a limit of $P^2SAD$ to zero separation in physical space. We illustrate our method by computing the global and local subhalo annihilation boost to that of the smooth dark matter distribution in a Milky-Way-size halo. Two cases are considered, one where the cross section is velocity independent and one that approximates Sommerfeld-enhanced models. We find that the global boost is $\sim10-30$, which is at the low end of current estimates (weakening expectations of large extragalactic signals), while the boost at the solar radius is below the percent level. We make our code to compute $P^2SAD$ publicly available, which can be used to estimate various observables that probe the nanostructure of dark matter haloes.

An intriguing possibility that can address pathologies in both early universe cosmology (i.e. the horizon problem) and quantum gravity (i.e. non-renormalizability), is that particles at very high energies and/or temperatures could propagate arbitrarily fast. A concrete realization of this possibility for the early universe is the Tachyacoustic (or Speedy Sound) cosmology, which could also produce a scale-invariant spectrum for scalar cosmological perturbations. Here, we study Thermal Tachyacoustic Cosmology (TTC), i.e. this scenario with thermal initial conditions. We find that a phase transition in the early universe, around the scale of Grand Unified Theories (GUT scale; $T\sim 10^{15}$ GeV), during which the speed of sound drops by $25$ orders of magnitude within a Hubble time, can fit current CMB observations. We further discuss how production of primordial black holes constrains the cosmological acoustic history, while coupling TTC to Horava-Lifshitz gravity leads to a lower limit on the amplitude of tensor modes ($r \gtrsim 10^{-3}$), that are detectable by CMBPol (and might have already been seen by the BICEP-Keck collaboration).

Both cosmological expansion and black holes are ubiquitous features of our observable Universe, yet exact solutions connecting the two have remained elusive. To this end, we study self-gravitating classical fields within dynamical spherically symmetric solutions that can describe black holes in an expanding universe. After attempting a perturbative approach of a known black-hole solution with scalar hair, we show by exact methods that the unique scalar field action with first-order derivatives that can source shear-free expansion around a black hole requires noncanonical kinetic terms. The resulting action is an incompressible limit of k-essence, otherwise known as the cuscuton theory, and the spacetime it describes is the McVittie metric. We further show that this solution is an exact solution to the vacuum Ho\v{r}ava-Lifshitz gravity with anisotropic Weyl symmetry.

We propose a novel solution for the endpoint of gravitational collapse, in which spacetime ends (and is orbifolded) at a microscopic distance from black hole event horizons. This model is motivated by the emergence of singular event horizons in the gravitational aether theory, a semi-classical solution to the cosmological constant problem(s), and thus suggests a catastrophic breakdown of general relativity close to black hole event horizons. A similar picture emerges in fuzzball models of black holes in string theory, as well as the recent firewall proposal to resolve the information paradox. We then demonstrate that positing a surface fluid in thermal equilibrium with Hawking radiation, with vanishing energy density (but non-vanishing pressure) at the new boundary of spacetime, which is required by Israel junction conditions, yields a thermodynamic entropy that is identical to the Bekenstein-Hawking area law, $S_{BH}$, for charged rotating black holes. To our knowledge, this is the first derivation of black hole entropy which only employs local thermodynamics. Furthermore, a model for the microscopic degrees of freedom of the surface fluid (which constitute the micro-states of the black hole) is suggested, which has a finite, but Lorentz-violating, quantum field theory. Finally, we comment on the effects of physical boundary on Hawking radiation, and show that relaxing the assumption of equilibrium with Hawking radiation sets $S_{BH}$ as an upper limit for Black Hole entropy.

We study the time variation of the apparent flux of cosmological point sources due to the transient weak lensing by dark matter microhaloes. Assuming a transverse motion of microhaloes with respect to our line of sight, we derive the correspondence between the temporal power spectrum of the weak lensing magnification, and the spatial power spectrum of density on small scales. Considering different approximations for the small scale structure of dark matter, we predict the apparent magnitude of cosmological point sources to vary by as much as $10^{-4}-10^{-3}$, due to this effect, within a period of a few months. This red photometric noise has an almost perfect gaussian statistics, to one part in $\sim 10^4$. We also compare the transient weak lensing power spectrum with the background effects such as the stellar microlensing on cosmological scales. A quasar lensed by a galaxy or cluster like SDSSJ1004+4112 strong lensing system, with multiple images, is a suitable system for this study as: (i) using the time-delay method between different images, we can remove the intrinsic variations of the quasar, and (ii) strong lensing enhances signals from the transient weak lensing. We also require the images to form at large angular separations from the center of the lensing structure, in order to minimize contamination by the stellar microlensing. With long-term monitoring of quasar strong lensing systems with a 10-meter class telescope, we can examine the existence of dark microhaloes as the building blocks of dark matter structures. Failure to detect this signal may either be caused by a breakdown of cold dark matter (CDM) hierarchy on small scales, or rather interpreted as evidence against CDM paradigm, e.g. in favor of modified gravity models.

While most of the singularities of General Relativity are expected to be safely hidden behind event horizons by the cosmic censorship conjecture, we happen to live in the causal future of the classical big bang singularity, whose resolution constitutes the active field of early universe cosmology. Could the big bang be also hidden behind a causal horizon, making us immune to the decadent impacts of a naked singularity? We describe a braneworld description of cosmology with both 4d induced and 5d bulk gravity (otherwise known as Dvali-Gabadadze-Porati, or DGP model), which exhibits this feature: The universe emerges as a spherical 3-brane out of the formation of a 5d Schwarzschild black hole. In particular, we show that a pressure singularity of the holographic fluid, discovered earlier, happens inside the white hole horizon, and thus need not be real or imply any pathology. Furthermore, we outline a novel mechanism through which any thermal atmosphere for the brane, with comoving temperature of 20% of the 5D Planck mass can induce scale-invariant primordial curvature perturbations on the brane, circumventing the need for a separate process (such as cosmic inflation) to explain current cosmological observations. Finally, we note that 5D space-time is asymptotically flat, and thus potentially allows an S-matrix or (after minor modifications) AdS/CFT description of the cosmological big bang.

We present a statistically-optimal and model-independent method to extract the pressure profile of hot gas in the intracluster medium (ICM). Using the thermal Sunyaev-Zeldovich effect, we constrain the mean pressure profile of the ICM by appropriately considering all primary cosmic microwave background (CMB) and instrumental noise correlations, while using the maximum resolution and sensitivity of all frequency channels. As a first application, we analyze CMB maps of WMAP 9-year data through a study of the Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC). We constrain the universal pressure profile out to 4R_500 with 15-sigma confidence, though our measurements are only significant out to R_200. Using a temperature profile constrained from X-ray observations, we measure the mean gas mass fraction out to R_200. Within statistical and systematic uncertainties, our constraints are compatible with the cosmic baryon fraction and the expected gas fraction in halos. While Planck multi-frequency CMB data are expected to reduce statistical uncertainties by a factor of 20, we argue that systematic errors in determining mass of clusters dominate the uncertainty in gas mass fraction measurements at the level of 20 percent.

Peculiar velocities in the nearby Universe can be measured via the kinetic Sunyaev-Zel'dovich (kSZ) effect. Using a statistical method based on an optimised cross-correlation with nearby galaxies, we extract the kSZ signal generated by plasma halo of galaxies from the Cosmic Microwave Background (CMB) temperature anisotropies observed by the Wilkinson Microwave Anisotropy Probe (WMAP). Marginalising over the thermal Sunyaev-Zel'dovich contribution from clusters of galaxies, possible unresolved point source contamination, and Galactic foregrounds (dust, synchrotron and free-free emission), we report a kSZ bulk flow signal present at the 90% confidence level in the seven-year WMAP data. When only galaxies within 50 Mpc/h are included in the kSZ template we find a bulk flow in the CMB frame of |V|=533 +/- 263 km/s, in the direction l=324 +/- 27, b=-7 +/- 17, consistent with bulk flow measurements on a similar scale using classical distance indicators. We show how this comparison constrains the (ionised) baryonic budget in the local universe. On very large (~ 500 Mpc/h) scales, we find a 95% upper limit of 470 km/s, inconsistent with some analyses of bulk flow of clusters from the kSZ. We estimate that the significance of the bulk flow signal may increase to 3-5 sigma using data from the PLANCK probe.

We define a distinguished "ground state" or "vacuum" for a free scalar quantum field in a globally hyperbolic region of an arbitrarily curved spacetime. Our prescription is motivated by the recent construction of a quantum field theory on a background causal set using only knowledge of the retarded Green's function. We generalize that construction to continuum spacetimes and find that it yields a distinguished vacuum or ground state for a non-interacting, massive or massless scalar field. This state is defined for all compact regions and for many noncompact ones. In a static spacetime we find that our vacuum coincides with the usual ground state. We determine it also for a radiation-filled, spatially homogeneous and isotropic cosmos, and show that the super-horizon correlations are approximately the same as those of a thermal state. Finally, we illustrate the inherent non-locality of our prescription with the example of a spacetime which sandwiches a region with curvature in-between flat initial and final regions.

The existence of an extended hot gaseous corona surrounding clusters, groups and massive galaxies is well established by observational evidence and predicted by current theories of galaxy formation. When a small galaxy collides with a larger one, their coronae are the first to interact, producing disturbances that remove gas from the smaller system and settle it into the corona of the larger one. For a Milky-Way-size galaxy merging into a low-mass group, ram pressure stripping and the Kelvin-Helmholtz instability are the most relevant of these disturbances. We argue that the turbulence generated by the latter mixes the material of both coronae in the wake of the orbiting satellite creating a "warm phase" mixture with a cooling time a factor of several shorter than that of the ambient intragroup gas. We reach this conclusion using analytic estimates, as well as adiabatic and dissipative high resolution numerical simulations of a spherical corona subject to the ablation process of a constant velocity wind with uniform density and temperature. Although this is a preliminary analysis, our results are promising and we speculate that the mixture could potentially trigger in situ star formation and/or be accreted into the central galaxy as a cold gas flow resulting in a new mode of star formation in galaxy groups and clusters.

The observed amount of lithium for low metallicity population II stars (known as the Spite plateau) is a factor of $\sim 3-5$ lower than the predictions of the standard cosmology. Since the observations are limited to the local Universe (halo stars, globular clusters and satellites of the Milky Way) it is possible that certain physical processes may have led to the spatial separation of lithium and local reduction of [Li/H]. We study the question of lithium diffusion after the cosmological recombination in sub-Jeans dark matter haloes, taking into account that more than 95% of lithium remains in the singly-ionized state at all times. Large scattering cross sections on the rest of the ionized gas leads to strong coupling of lithium to protons and its initial direction of diffusion coincides with that of H$^+$. In the rest frame of the neutral gas this leads to the diffusion of H$^+$ and Li$^+$ out of overdensities with the trend of reducing [Li/H] in the minima of gravitational wells relative to the primordial value. We quantify this process and argue that, with certain qualifications, it may have played a significant role in creating local lithium deficiency within the primordial dark matter haloes, comparable to those observed along the Spite plateau.

We show that in cases of marginal detections (~ 3\sigma), such as that of Baryonic Acoustic Oscillations (BAO) in cosmology, the often-used Gaussian approximation to the full likelihood is very poor, especially beyond ~3\sigma. This can radically alter confidence intervals on parameters and implies that one cannot naively extrapolate 1\sigma-errorbars to 3\sigma, and beyond. We propose a simple fitting formula which corrects for this effect in posterior probabilities arising from marginal detections. Alternatively the full likelihood should be used for parameter estimation rather than the Gaussian approximation of a just mean and an error.

Einstein's theory of General Relativity is the benchmark example for empirical success and mathematical elegance in theoretical physics. However, in spite of being the most successfully tested theory in physics, there are strong theoretical and observational arguments for why General Relativity should fail. It is not a question of if, but rather a question of where and when! I start by recounting the tremendous success in observational cosmology over the past three decades, that has led to the era of precision cosmology. I will then summarize the pathologies in Einstein's theory of gravity, as the cornerstone of standard cosmological model. Attempts to address these pathologies are either inspired by mathematical elegance, or empirical falsifiability. Here, I provide different arguments for why a falsifiable solution should violate Lorentz symmetry, or revive "gravitational aether". Deviations from Einstein's gravity are then expected in: 1) cosmological matter-radiation transition, 2) neutron stars, 3) gravitomagnetic effect, 4) astrophysical black holes, and their potential connection to dark energy, and 5) early Universe, where the predictions are ranked by their degree of robustness and falsifiability.

One of the deepest and most long-standing mysteries in physics has been the huge discrepancy between the observed vacuum density and our expectations from theories of high energy physics, which has been dubbed the Old Cosmological Constant problem. One proposal to address this puzzle at the semi-classical level is to decouple quantum vacuum from space-time geometry via a modification of gravity that includes an incompressible fluid, known as Gravitational Aether. In this paper, we discuss classical predictions of this theory along with its compatibility with cosmological and experimental tests of gravity. We argue that deviations from General Relativity (GR) in this theory are sourced by pressure or vorticity. In particular, the theory predicts that the gravitational constant for radiation is 33% larger than that of non-relativistic matter, which is preferred by (most) cosmic microwave background (CMB), Lyman-Alpha forest, and Lithium-7 primordial abundance observations, while being consistent with other cosmological tests at ~2-sigma level. It is further shown that all Parametrized Post-Newtonian (PPN) parameters have the standard GR values aside from the anomalous coupling to pressure, which has not been directly measured. A more subtle prediction of this model (assuming irrotational aether) is that the (intrinsic) gravitomagnetic effect is 33% larger than GR prediction. This is consistent with current limits from LAGEOS and Gravity Probe B at ~2-sigma level.

The gravitational aether theory is a modification of general relativity that decouples vacuum energy from gravity, and thus can potentially address the cosmological constant problem. The classical theory is distinguishable from general relativity only in the presence of relativistic pressure (or vorticity). Since the interior of neutron stars has high pressure and as their mass and radius can be measured observationally, they are the perfect laboratory for testing the validity of the aether theory. In this paper, we solve the equations of stellar structure for the gravitational aether theory and find the predicted mass-radius relation of non-rotating neutron stars using two different realistic proposals for the equation of state of nuclear matter. We find that the maximum neutron star mass predicted by the aether theory is 12% - 16% less than the maximum mass predicted by general relativity assuming these two equations of state. We also show that the effect of aether is similar to modifying the equation of state in general relativity. The effective pressure of the neutron star given by the aether theory at a fiducial density differs from the values given by the two nuclear equations of state to an extent that can be constrained using future gravitational wave observations of neutron stars in compact systems. This is a promising way to test the aether theory if further progress is made in constraining the equation of state of nuclear matter in densities above the nuclear saturation density.

One of the open questions of modern cosmology is the nature and properties of the Dark Matter halo and its substructures. In this work we study the gravitational effect of dark matter substructures on pulsar timing observations. Since millisecond pulsars are stable and accurate emitters, they have been proposed as plausible astrophysical tools to probe the gravitational effects of dark matter structures. We study this effect on pulsar timing through Shapiro time delay (or Integrated Sachs-Wolfe (ISW) effect) and Doppler effects statistically, showing that the latter dominates the signal. For this task, we relate the power spectrum of pulsar frequency change to the matter power spectrum on small scales, which we compute using the stable clustering hypothesis. We compare this power spectrum with the reach of current and future observations of pulsar timing designed for gravitational wave (GW) detection. Our results show that while current observations are unable to detect these signals, the sensitivity of the upcoming Square Kilometer Array (SKA) is only a factor of few weaker than our optimistic predictions.

Light scalar fields are expected to arise in theories of high energy physics (such as string theory), and find phenomenological motivations in dark energy, dark matter, or neutrino physics. However, the coupling of light scalar fields to ordinary (or dark) matter is strongly constrained from laboratory, solar system, and astrophysical tests of fifth force. One way to evade these constraints in dense environments is through the chameleon mechanism, where the field's mass steeply increases with ambient density. Consequently, the chameleonic force is only sourced by a thin shell near the surface of dense objects, which significantly reduces its magnitude. In this paper, we argue that thin-shell conditions are equivalent to "conducting" boundary conditions in electrostatics. As an application, we use the analogue of the method of images to calculate the back-reaction (or self-force) of an object around a spherical gravitational source. Using this method, we can explicitly compute the violation of equivalence principle in the outskirts of galactic haloes (assuming an NFW dark matter profile): Intermediate mass satellites can be slower than their larger/smaller counterparts by as much as 10% close to a thin shell.

We present a simple inflationary scenario that can produce arbitrarily large spherical underdense or overdense regions embedded in a standard Lambda cold dark matter paradigm, which we refer to as bubbles. We analyze the effect such bubbles would have on the Cosmic Microwave Background (CMB). For super-horizon sized bubble in the vicinity of the last scattering surface, a signal is imprinted onto CMB via a combination of Sach-Wolfe and an early integrated Sach-Wolfe (ISW) effects. Smaller, sub-horizon sized bubbles at lower redshifts (during matter domination and later) can imprint secondary anisotropies on the CMB via Rees-Sciama, late-time ISW and Ostriker-Vishniac effects. Our scenario, and arguably most similar inflationary models, produce bubbles which are over/underdense in potential: in density such bubbles are characterized by having a distinct wall with the interior staying at the cosmic mean density. We show that such models can potentially, with only moderate fine tuning, explain the \emph{cold spot}, a non-Gaussian feature identified in the Wilkinson Microwave Anisotropy Probe (WMAP) data by several authors. However, more detailed comparisons with current and future CMB data are necessary to confirm (or rule out) this scenario.

Structure formation provides a strong test of any cosmic acceleration model because a successful dark energy model must not inhibit {\black or overpredict} the development of observed large-scale structures. Traditional approaches to studies of structure formation in the presence of dark energy or a modified gravity implement a modified Press-Schechter formalism, which relates the linear overdensities to the abundance of dark matter haloes it at the same time. We critically examine the universality of the Press-Schechter formalism for different cosmologies, and show that the halo abundance is best correlated with spherical linear overdensity at 94% of collapse (or observation) time. We then extend this argument to ellipsoidal collapse (which decreases the fractional time of best correlation for small haloes, and show that our results agree with deviations from modified Press-Schechter formalism seen in simulated mass functions. This provides a novel universal prescription to measure linear density evolution, based on current and future observations of cluster (or dark matter) halo mass function. In particular, even observations of cluster abundance in a single epoch will constrain the entire history of linear growth of cosmological of perturbations.

Einstein's theory of general relativity describes gravity as the interaction of particles with space-time geometry, as opposed to interacting with a physical fluid, as in the old gravitational aether theories. Moreover, any theoretical physicist would tell you that, despite its counter-intuitive structure, general relativity is one of the simplest, most beautiful, and successful theories in physics, that has withstood a diverse battery of precision tests over the past century. So, is there any motivation to relax its fundamental principle, and re-introduce a gravitational aether? Here, I give a short and non-technical account of why quantum gravity and cosmological constant problems provide this motivation.

The two apparently distinct phenomena of dark energy (or late-time cosmic acceleration) and quantum gravity dominate physics on extremely low, and extremely high energies, but do not seem to have any apparent empirical connection. Nevertheless, the two have a theoretical connection, through the cosmological constant problem. I argue that the finite temperature quantum gravitational corrections to black hole entropy yields a pressure for the gravitational vacuum (or gravitational aether). Assuming that the relative corrections are linear in horizon temperature (i.e. are suppressed by one power of Planck energy), the pressure is comparable to that of dark energy for astrophysical black holes. This implies that the observation of late-time cosmic acceleration may have provided us with the first precision measurement of quantum gravity, i.e. that of black hole entropy.

We develop a theoretical framework for describing the hierarchical structure of the phase space of cold dark matter haloes, due to gravitationally bound substructures. Because it includes the full hierarchy of the cold dark matter initial conditions and is hence complementary to the halo model, the stable clustering hypothesis is applied for the first time here to the small-scale phase space structure. As an application, we show that the particle dark matter annihilation signal could be up to two orders of magnitude larger than that of the smooth halo within the Galactic virial radius. The local boost is inversely proportional to the smooth halo density, and thus is O(1) within the solar radius, which could translate into interesting signatures for dark matter direct detection experiments: The temporal correlation of dark matter detection can change by a factor of 2 in the span of 10 years, while there will be significant correlations in the velocity space of dark matter particles. This can introduce O(1) uncertainty in the direction of local dark matter wind, which was believed to be a benchmark of directional dark matter searches or the annual modulation signal.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

A proposal for a power-counting renormalizable theory of quantum gravity at a Lifshitz point was recently put forth by Horava (arXiv:0901.3775), and has been since dubbed as Horava-Lifshitz gravity. The theory explicitly breaks Lorentz invariance, which introduces an apparent extra scalar degree of freedom. In this note, we show that the low energy limit of (non-projectible) Horava-Lifshitz gravity is uniquely given by the quadratic cuscuton model: a covariant scalar field theory with an infinite speed of sound and a quadratic potential, which is minimally coupled to Einstein gravity. This implies that the extra scalar is non-dynamical to all orders in perturbation theory. Using current cosmological constraints on the quadratic cuscuton model, we can constrain the low energy Lorentz breaking parameter of Horava-Lifshitz theory (which leads to a running of Planck mass on the Hubble scale) to |lambda-1| < 0.014 (at 95% confidence level). We also point out that, with reasonable boundary conditions, the spatial hypersurfaces in this theory are constant mean curvature (CMC) or uniform expansion surfaces at low energies, and introduce geometrical symmetries that can protect the non-dynamical nature of these theories from quantum corrections. We also notice that the theory with lambda < 1 might suffer from a non-perturbative UV instability.

The discovery of cosmic acceleration has presented a unique challenge for cosmologists. As observational cosmology forges ahead, theorists have struggled to make sense of a standard model that requires extreme fine tuning. This challenge is known as the cosmological constant problem. The theory of gravitational aether is an alternative to general relativity that does not suffer from this fine-tuning problem, as it decouples the quantum field theory vacuum from geometry, while remaining consistent with other tests of gravity. In this paper, we study static black hole solutions in this theory and show that it manifests a UV-IR coupling: Aether couples the spacetime metric close to the black hole horizon, to metric at infinity. We then show that using the Trans-Planckian ansatz (as a quantum gravity effect) close to the black hole horizon, leads to an accelerating cosmological solution, far from the horizon. Interestingly, this acceleration matches current observations for stellar mass black holes. Based on our current understanding of the black hole accretion history in the Universe, we then make a prediction for how the effective dark energy density should evolve with redshift, which can be tested with future dark energy probes.

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 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.

Most of the mass content of dark matter haloes is expected to be in the form of tidal debris. The density of debris is not constant, but rather can grow due to formation of caustics at the apocenters and pericenters of the orbit, or decay as a result of phase mixing. In the phase space, the debris assemble in a hierarchy which is truncated by the primordial temperature of dark matter. Understanding this phase structure can be of significant importance for the interpretation of many astrophysical observations and in particular dark matter detection experiments. With this purpose in mind, we develop a general theoretical framework to describe the hierarchical structure of the phase space of cold dark matter haloes. We do not make any assumption of spherical symmetry and/or smooth and continuous accretion. Instead, working with correlation functions in the action-angle space, we can fully account for the hierarchical structure (predicting a two-point correlation function ~ (\Delta J)^{-1.6} in the action space), as well as the primordial discreteness of the phase space. As an application, we estimate the boost to the dark matter annihilation signal due to the structure of the phase space within virial radius: the boost due to the hierarchical tidal debris is of order unity, whereas the primordial discreteness of the phase structure can boost the total annihilation signal by up to an order of magnitude. The latter is dominated by the regions beyond 20% of the virial radius, and is largest for the recently formed haloes with the least degree of phase mixing.

We develop a systematic derivation for the Limber approximation to the angular cross-power spectrum of two random fields, as a series expansion in 1/(\ell+1/2). This extended Limber approximation can be used to test the accuracy of the Limber approximation and to improve the rate of convergence at large \ell's. We show that the error in ordinary Limber approximation is O(1/\ell^2). We also provide a simple expression for the second order correction to the Limber formula, which improves the accuracy to O(1/\ell^4). This correction can be especially useful for narrow redshift bins, or samples with small redshift overlap, for which the zeroth order Limber formula has a large error. We also point out that using \ell instead of (\ell+1/2), as is often done in the literature, spoils the accuracy of the approximation to O(1/\ell).

Cosmological constant problem (in its various versions) is arguably the deepest gap in our understanding of theoretical physics, the solution to which may very likely require revisiting the Einstein theory of gravity. In this letter, I argue that the simplest consistent way to decouple gravity from the vacuum energy (and hence solve the problem) is through the introduction of an incompressible gravitational aether fluid. The theory then predicts that gravitational constant for radiation is 33% larger than that of non-relativistic matter, which is preferred by most cosmological observations (with the exception of light element abundances), but is not probed by current precision tests of gravity. I also show that slow-roll inflation can happen in this theory, with only minor modifications. Finally, interpreting gravitational aether as a thermodynamic description of gravity, I propose a finite-temperature correction to the equation of state of gravity, which would explain the present-day acceleration of the cosmic expansion as a consequence of the formation of stellar mass black holes.

Any hint of non-gaussianity in the cosmological initial conditions will provide us with a unique window into the physics of early universe. We show that the impact of a small local primordial non-gaussianity (generated on super-horizon scales) on the statistics of collapsed objects (such as galaxies or clusters) can be approximated by using slightly modified, but gaussian, initial conditions, which we describe through simple analytic expressions. Given that numerical simulations with gaussian initial conditions are relatively well-studied, this equivalence provides us with a simple tool to predict signatures of primordial non-gaussianity in the statistics of collapsed objects. In particular, we describe the predictions for non-gaussian mass function, and also confirm the recent discovery of a non-local bias on large scales (arXiv:0710.4560, arXiv:0801.4826) as a signature of primordial non-gaussianity. We then study the potential of galaxy surveys to constrain non-gaussianity using their auto-correlation and cross-correlation with the CMB (due to the Integrated Sachs-Wolfe effect), as a function of survey characteristics, and predict that they will eventually yield an accuracy of Delta f_{NL} ~ 0.1 and 3 respectively, which will be better than or competitive with (but independent of) the best predicted constraints from the CMB. Interestingly, the cross-correlation of CMB and NVSS galaxy survey already shows a hint of a large local primordial non-gaussianity: f_{NL} = 236 +/- 127.

Sunyaev-Zel'dovich (SZ) cluster surveys are considered among the most promising methods for probing dark energy up to large redshifts. However, their premise is hinged upon an accurate mass-observable relationship, which could be affected by the (rather poorly understood) physics of the intracluster gas. In this letter, using a semi-analytic model of the intracluster gas that accommodates various theoretical uncertainties, I develop a Fundamental Plane relationship between the observed size, thermal energy, and mass of galaxy clusters. In particular, I find that M ~ (Y_{SZ}/R_{SZ,2})^{3/4}, where M is the mass, Y_{SZ} is the total SZ flux or thermal energy, and R_{SZ,2} is the SZ half-light radius of the cluster. I first show that, within this model, using the Fundamental Plane relationship reduces the (systematic+random) errors in mass estimates to 14%, from 22% for a simple mass-flux relationship. Since measurement of the cluster sizes is an inevitable part of observing the SZ clusters, the Fundamental Plane relationship can be used to reduce the error of the cluster mass estimates by ~ 34%, improving the accuracy of the resulting cosmological constraints without any extra cost. I then argue why our Fundamental Plane is distinctly different from the virial relationship that one may naively expect between the cluster parameters. Finally, I argue that while including more details of the observed SZ profile cannot significantly improve the accuracy of mass estimates, a better understanding of the impact of non-gravitational heating/cooling processes on the outskirts of the intracluster medium (apart from external calibrations) might be the best way to reduce these errors.

Einstein's theory of general relativity, which contains a universal value of the Planck mass, has been so far successfully invoked to explain gravitational dynamics from sub-millimeter scales to the scale of the cosmological horizon. However, one may envisage that in alternative theories of gravity, the effective value of the Planck mass (or Newton's constant), which quantifies the coupling of matter to metric perturbations, can run on the cosmological horizon scale. In this letter, we study the consequences of a glitch in the Planck mass from sub-horizon to super-horizon scales. We first give three examples of models that naturally exhibit this feature, and then show that current cosmological observations severely constrain this glitch to less than 1.2%. This is the strongest constraint to date, on natural (i.e. non-fine-tuned) deviations from Einstein gravity on the cosmological horizon scale.

In a companion paper (hep-th/0609150), we have introduced a model of scalar field dark energy, Cuscuton, which can be realized as the incompressible (or infinite speed of sound) limit of a k-essence fluid. In this paper, we study how Cuscuton modifies the constraint sector of Einstein gravity. In particular, we study Cuscuton cosmology and show that even though Cuscuton can have an arbitrary equation of state, or time dependence, and is thus inhomogeneous; its perturbations do not introduce any additional dynamical degree of freedom and only satisfy a constraint equation, amounting to an effective modification of gravity on large scales. Therefore, Cuscuton can be considered to be a minimal theory of evolving dark energy, or a minimal modification of a cosmological constant, as it has no internal dynamics. Moreover, this is the only modification of Einstein gravity to our knowledge that does not introduce any additional degrees freedom (and is not conformally equivalent to the Einstein gravity). We then study two simple Cuscuton models, with quadratic and exponential potentials. The quadratic model has the exact same expansion history as LCDM, and yet contains an early dark energy component with constant energy fraction, which is constrained to < 2%, mainly from WMAP Cosmic Microwave Background (CMB) and SDSS Lyman-alpha forest observations. The exponential model has the same expansion history as the DGP self-accelerating braneworld model, but generates a much smaller Integrated Sachs-Wolfe (ISW) effect, and is thus consistent with the CMB observations. Finally, we show that the evolution is local on super-horizon scales, implying that there is no gross violation of causality, despite Cuscuton's infinite speed of sound.

We introduce a model of scalar field dark energy, Cuscuton, which can be realized as the incompressible (or infinite speed of sound) limit of a scalar field theory with a non-canonical kinetic term (or k-essence). Even though perturbations of Cuscuton propagate superluminally, we show that they have a locally degenerate phase space volume (or zero entropy), implying that they cannot carry any microscopic information, and thus the theory is causal. Even coupling to ordinary scalar fields cannot lead to superluminal signal propagation. Furthermore, we show that the family of constant field hypersurfaces are the family of Constant Mean Curvature (CMC) hypersurfaces, which are the analogs of soap films (or soap bubbles) in a Euclidian space. This enables us to find the most general solution in 1+1 dimensions, whose properties motivate conjectures for global degeneracy of the phase space in higher dimensions. Finally, we show that the Cuscuton action can model the continuum limit of the evolution of a field with discrete degrees of freedom and argue why it is protected against quantum corrections at low energies. While this paper mainly focuses on interesting features of Cuscuton in a Minkowski spacetime, a companion paper (astro-ph/0702002) examines cosmology with Cuscuton dark energy.

The Sunyaev-Zel'dovich (SZ) effect is a direct probe of thermal energy content of the Universe, induced in the cosmic microwave background (CMB) sky through scattering of CMB photons off hot electrons in the intracluster medium (ICM). We report a 9-sigma detection of the SZ signal in the CMB maps of Wilkinson Microwave Anisotropy Probe (WMAP) 3yr data, through study of a sample of 193 massive galaxy clusters with observed X-ray temperatures greater than 3 keV. For the first time, we make a model-independent measurement of the pressure profile in the outskirts of the ICM, and show that it closely follows the profiles obtained by X-ray observations and numerical simulations. We find that our measurements of the SZ effect would account for only half of the thermal energy of the cluster, if all the cluster baryons were in the hot ICM phase. Our measurements indicate that a significant fraction (35 +/- 8 %) of baryonic mass is missing from the hot ICM, and thus must have cooled to form galaxies, intracluster stars, or an unknown cold phase of the ICM. There does not seem to be enough mass in the form of stars or cold gas in the cluster galaxies or intracluster space, signaling the need for a yet-unknown baryonic component (at 3-sigma level), or otherwise new astrophysical processes in the ICM.

There has been some recent activity in trying to understand the dark matter clustering properties in the quasilinear regime, through resummation of perturbative terms, otherwise known as the renormalized perturbation theory (astro-ph/0509418), or the renormalization group method (astro-ph/0606028). While it is not always clear why such methods should work so well, there is no reason for them to capture non-perturbative events such as shell-crossing. In order to estimate the magnitude of non-perturbative effects, we introduce a (hypothetical) model of sticky dark matter, which only differs from collisionless dark matter in the shell-crossing regime. This enables us to show that the level of non-perturbative effects in the dark matter power spectrum at k ~ 0.1 Mpc^{-1}, which is relevant for baryonic acoustic oscillations, is about a percent, but rises to order unity at k ~ 1 Mpc^{-1}.

Wilkinson Microwave Anisotropy Probe (WMAP) has provided us with the yet highest resolution all-sky maps of the Cosmic Microwave Background. As a result of thermal Sunyaev-Zel'dovich effect, clusters of galaxies are imprinted as tiny, poorly resolved dips on top of primary CMB anisotropies in these maps. Here, I describe different efforts to extract the physics of Intracluster Medium (ICM) from the sea of primary CMB, through combining WMAP with low-redshift galaxy or X-ray cluster surveys. This finally culminates at a mean (universal) ICM pressure profile, which is for the first time directly constrained from WMAP 3yr maps, and leads to interesting constraints on the ICM baryonic budget.

Recently, new arguments (astro-ph/0501152, hep-th/0503117) for how corrections from super-Hubble modes can explain the present-day acceleration of the universe have appeared in the literature. However, in this letter, we argue that, to second order in spatial gradients, these corrections only amount to a renormalization of local spatial curvature, and thus cannot account for the negative deceleration. Moreover, cosmological observations already put severe bounds on such corrections, at the level of a few percent, while in the context of inflationary models, these corrections are typically limited to ~ 10^{-5}. Currently there is no general constraint on the possible correction from higher order gradient terms, but we argue that such corrections are even more constrained in the context of inflationary models.

We study the possible origin of hydrodynamic turbulence in cold accretion disks such as those in star-forming systems and quiescent cataclysmic variables. As these systems are expected to have neutral gas, the turbulent viscosity is likely to be hydrodynamic in origin, not magnetohydrodynamic. Therefore MRI will be sluggish or even absent in such disks. Although there are no exponentially growing eigenmodes in a hydrodynamic disk, because of the non-normal nature of the eigenmodes, a large transient growth in the energy is still possible, which may enable the system to switch to a turbulent state. For a Keplerian disk, we estimate that the energy will grow by a factor of 1000 for a Reynolds number close to a million.

An interesting dynamical model for dark energy which does not require extremely light scalar fields such as quintessence, and at the same time explains the (near-) coincidence between the neutrino and dark energy densities is the model of dark energy coupled to mass varying neutrinos (MaVaNs). Despite the attractions of this model, we show that, generically, this model contains a catastrophic instability which occurs when neutrinos become non-relativistic. As a result of this instability, as neutrinos become non-relativistic, they condense into neutrino nuggets which redshift away similar to cold dark matter, and thus cease to act as dark energy. Any stable MaVaNs dark energy model is extremely contrived, and is virtually indistinguishable from a cosmological constant.

We devise a Monte-Carlo based, optimized filter match method to extract the thermal Sunyaev-Zel'dovich (SZ) signature of a catalog of 116 low-redshift X-ray clusters from the first year data release of the Wilkinson Microwave Anisotropy Probe (WMAP). We detect an over-all amplitude for the SZ signal at the ~ 8-sigma level, yielding a combined constraint of f_{gas}h = 0.08 +/- 0.01 (ran) +/- 0.01 (sys) on the gas mass fraction of the Intra-Cluster Medium. We also compile X-ray estimated gas fractions from the literature for our sample, and find that they are consistent with the SZ estimates at the 2-sigma level, while both show an increasing trend with X-ray temperature. Nevertheless, our SZ estimated gas fraction is 30-40% smaller than the concordance LCDM cosmic average. We also express our observations in terms of the SZ flux-temperature relation, and compare it with other observations, as well as numerical studies. Based on its spectral and spatial signature, we can also extract the microwave point source signal of the clusters at the 3-sigma level, which puts the average microwave luminosity (at ~ 41 GHz) of bright cluster members (M_K < -21) at (2.4 +/- 0.8) x 10^{27} h^{-2} erg/s/Hz. Furthermore, we can constrain the average dark matter halo concentration parameter to c_{vir}=3.4+0.6-0.9, for clusters with T_x > 5 kev. Our work serves as an example for how correlation of SZ surveys with cluster surveys in other frequencies can significantly increase our physical understanding of the intra-cluster medium.

Despite observational evidence for cold neutral astrophysical accretion disks, the viscous process which may drive the accretion in such systems is not yet understood. While molecular viscosity is too small to explain the observed accretion efficiencies by more than ten orders of magnitude, the absence of any linear instability in Keplerian accretion flows is often used to rule out the possibility of turbulent viscosity. Recently, the fact that some fine tuned disturbances of any inviscid shear flow can reach arbitrarily large transient growth has been proposed as an alternative route to turbulence in these systems. We present an analytic study of this process for 3D plane wave disturbances of a general rotating shear flow in Lagrangian coordinates, and demonstrate that large transient growth is the generic feature of non-axisymmetric disturbances with near radial leading wave vectors. The maximum energy growth is slower than quadratic, but faster than linear in time. The fastest growth occurs for two dimensional perturbations, and is only limited by viscosity, and ultimately by the disk vertical thickness. After including viscosity and vertical structure, we find that, as a function of the Reynolds number, R, the maximum energy growth is approximately 0.4 (R/log R)^{2/3}, and put forth a heuristic argument for why R > 10^4 is required to sustain turbulence in Keplerian disks. Therefore, assuming that there exists a non-linear feedback process to replenish the seeds for transient growth, astrophysical accretion disks must be well within the turbulent regime. However, large 3D numerical simulations running for many orbital times, and/or with fine tuned initial conditions, are required to confirm Keplerian hydrodynamic turbulence on the computer.

Cold accretion disks such as those in star-forming systems, quiescent cataclysmic variables, and some active galactic nuclei, are expected to have neutral gas which does not couple well to magnetic fields. The turbulent viscosity in such disks must be hydrodynamic in origin, not magnetohydrodynamic. We investigate the growth of hydrodynamic perturbations in a linear shear flow sandwiched between two parallel walls. The unperturbed flow is similar to plane Couette flow but with a Coriolis force included. Although there are no exponentially growing eigenmodes in this system, nevertheless, because of the non-normal nature of the eigenmodes, it is possible to have a large transient growth in the energy of perturbations. For a constant angular momentum disk, we find that the perturbation with maximum growth has a wave-vector in the vertical direction. The energy grows by more than a factor of 100 for a Reynolds number R=300 and more than a factor of 1000 for R=1000. Turbulence can be easily excited in such a disk, as found in previous numerical simulations. For a Keplerian disk, on the other hand, similar vertical perturbations grow by no more than a factor of 4, explaining why the same simulations did not find turbulence in this system. However, certain other two-dimensional perturbations with no vertical structure do exhibit modest growth. For the optimum two-dimensional perturbation, the energy grows by a factor of ~100 for R~10^4.5 and by a factor of 1000 for R~10^6. It is conceivable that these two-dimensional disturbances might lead to self-sustained turbulence. The Reynolds numbers of cold astrophysical disks are much larger even than 10^6, therefore, hydrodynamic turbulence may be possible in disks.

Turbulent viscosity in cold accretion disks is likely to be hydrodynamic in origin. We investigate the growth of hydrodynamic perturbations in a small region of a disk, which we model as a linear shear flow with Coriolis force, between two parallel walls. Although there are no exponentially growing eigenmodes in this system, because of the non-normal nature of the modes, it is possible to have a large transient growth in the energy of certain perturbations. For a constant angular momentum disk, the energy grows by more than a factor of 1000 for a Reynolds number of only 1000, and so turbulence is easily excited. For a Keplerian disk, the growth is more modest, and energy growth by a factor of 1000 requires a Reynolds number of nearly a million. Accretion disks have even larger Reynolds numbers than this. Therefore, transient growth of perturbations could seed turbulence in such disks.

[Abridged] The polarization of the CMB is widely recognized as a potential source of information about primordial gravitational waves. The gravitational wave contribution can be separated from the dominant CMB polarization created by density perturbations because it generates both E and B polarization modes, whereas the density perturbations create only E polarization. The limits of our ability to measure gravitational waves are thus determined by statistical and systematic errors, foregrounds, and nonlinear evolution effects such as lensing of the CMB. Usually it is assumed that most foregrounds can be removed because of their frequency dependence, however Thomson scattering of the CMB quadrupole by electrons in the Galaxy or nearby structures shares the blackbody frequency dependence of the CMB. If the optical depth from these nearby electrons is anisotropic, the polarization generated can include B modes even without tensor perturbations. We estimate this effect for the Galactic disk and nearby extragalactic structures, and find that it contributes to the B polarization at the level of ~(1--2)x10^-4\mu K per logarithmic interval in multipole L for L<30. This is well below the detectability level even for a future CMB polarization satellite. Depending on its structure and extent, the Galactic corona may be a source of B-modes comparable to the residual large-scale lensing B-mode after the latter has been cleaned using lensing reconstruction techniques. For an extremely ambitious post-Planck CMB experiment, Thomson scattering in the Galactic corona is thus a potential contaminant of the gravitational wave signal; conversely, if the other foregrounds can be cleaned out, such an experiment might be able to constrain models of the corona.

The Integrated Sachs-Wolfe (ISW) effect is a direct signature of the presence of dark energy in the universe, in the absence of spatial curvature. A powerful method for observing the ISW effect is through cross-correlation of the Cosmic Microwave Background (CMB) with a tracer of the matter in the low redshift universe. In this paper, we describe the dependence of the obtained cross-correlation signal on the geometry and other properties of a survey of the low redshift universe. We show that an all-sky survey with about 10 million galaxies, almost uniformly distributed within 0<z<1 should yield a near optimal ISW detection, at ~ 5\sigma level. In order to achieve this level of signal-to-noise, the systematic anisotropies in the survey must be below ~ 0.1 %, on the scale of ~ 10 degrees on the sky, while the systematic error in redshift estimates must be less than 0.05. Then, we argue that, while an ISW detection will not be a good way of constraining the conventional properties of dark energy, it could be a valuable means of testing alternative theories of gravity on large physical scales.

We introduce a self-consistent stochastic coarse-graining method, which includes both metric and scalar field fluctuations, to investigate the back reaction of long wavelength perturbations in single-scalar driven inflation, up to the second (one loop) order. We demonstrate that, although back reaction cannot be significant during the last 70 e-foldings of inflation with a smooth potential, there exist non-smooth inflaton potentials which allow significant back reaction, and are also consistent with cosmological observations. Such non-smooth potentials may lead to the generation of massive primordial black holes, which could be further used to constrain/verify these models.

We cross-correlate the Cosmic Microwave Background (CMB) temperature anisotropies observed by the Wilkinson Microwave Anisotropy Probe (WMAP) with the projected distribution of extended sources in the Two Micron All Sky Survey (2MASS). By modelling the theoretical expectation for this signal, we extract the signatures of dark energy (Integrated Sachs-Wolfe effect;ISW), hot gas (thermal Sunyaev-Zeldovich effect;thermal SZ), and microwave point sources in the cross-correlation. Our strongest signal is the thermal SZ, at the 3.1-3.7 \sigma level, which is consistent with the theoretical prediction based on observations of X-ray clusters. We also see the ISW signal at the 2.5 \sigma level, which is consistent with the expected value for the concordance LCDM cosmology, and is an independent signature of the presence of dark energy in the universe. Finally, we see the signature of microwave point sources at the 2.7 \sigma level.

A numerical model of a steady state, thin accretion disk with a constant effective speed of sound is presented. We demonstrate that `zero torque' inner boundary condition is a reasonable approximation provided that the disk thickness, including the thickness of the torquing magnetic fields, is small everywhere. It is likely that this conclusion is correct also for non-steady disks, as long as the total thickness at the sonic point, H_c, is much smaller than the radius there, r_c ~ r_{ms}. The very existence of thin disks is not proved or disproved in this work, but such disks are believed to exist for moderate accretion rates. Within our model there is a small torque at r_{ms}, which may increase disk luminosity by several percent. An important result of our analysis is that the physically acceptable steady state solutions in our toy model exist only for \alpha < 0.14 (100 v_s/c)^{1/3}. A significant torque may be applied to a thin disk if there is a large scale magnetic field, like in a modified Blandford-Znajek mechanism.

We present measurements of the angular cross-correlation between luminous red galaxies from the Sloan Digital Sky Survey and the cosmic microwave background temperature maps from the Wilkinson Microwave Anisotropy Probe. We find a statistically significant achromatic positive correlation between these two data sets, which is consistent with the expected signal from the late Integrated Sachs-Wolfe (ISW) effect. We do not detect any anti-correlation on small angular scales as would be produced from a large Sunyaev-Zel'dovich (SZ) effect, although we do see evidence for some SZ effect for our highest redshift samples. Assuming a flat universe, our preliminary detection of the ISW effect provides independent physical evidence for the existence of dark energy.

We consider the possibility that massive primordial black holes are the dominant form of dark matter. Black hole formation generates entropy fluctuations that adds a Poisson noise to the matter power spectrum. We use Lyman-alpha forest observations to constrain this Poisson term in matter power spectrum, then we constrain the mass of black holes to be less than few times 10^4 solar mass. We also find that structures with less than ~ 10^3 primordial black holes evaporate by now.

If there was a weak magnetic field in the early universe, cyclotron emission could play an important role in the thermalization of the CMB. We study this process in the tightly coupled primordial electron-photon plasma and find that if the magnetic field is large enough so that the plasma effects allow emission of cyclotron photons, this process will wipe out deviations from the black body spectrum.

Combining conservation of energy throughout nearly-spherical collapse of galaxy clusters with the virial theorem, we derive the mass-temperature relation for X-ray clusters of galaxies $T=CM^{2/3}$. The normalization factor $C$ and the scatter of the relation are determined from first principles with the additional assumption of initial Gaussian random field. We are also able to reproduce the recently observed break in the M-T relation at $T \sim 3 \keV$, based on the scatter in the underlying density field for a low density $\Lambda$CDM cosmology. Finally, by combining observational data of high redshift clusters with our theoretical formalism, we find a semi-empirical temperature-mass relation which is expected to hold at redshifts up to unity with less than 20% error.

We show that the non-linear evolution of long wavelength perturbations may be important in a wide class of inflationary scenarios. We develop a solution for the evolution of such nonlinear perturbations which is exact to first order in a gradient expansion. As a first application, we demonstrate that in single field models of inflation there can be no parametric amplification of super-Hubble modes during reheating. We consider the implications of the solution for recent discussions of the back-reaction effect of long wavelength perturbations on the background geometry, give a new derivation of the equation of motion of stochastic inflation, and demonstrate that if the (generalized) slow-rolling condition is not satisfied, then inevitably long wavelength vector modes for gravitational fluctuations will be generated.

We write a non-relativistic Lagrangian for a hierarchical universe. The equations of motion are solved numerically and the evolution of the fractal dimension is obtained for different initial conditions. We show that our model is homogeneous at the time of the last scattering, but evolves into a self-similar universe with a remarkably constant fractal dimension. We also show that the Hubble law is implied by this model and make an estimate for the age of the universe.