On January 14, 2019, the Major Atmospheric Gamma Imaging Cherenkov telescopes detected GRB 190114C above 0.2 TeV, recording the most energetic photons ever observed from a gamma-ray burst. We use this unique observation to probe an energy dependence of the speed of light in vacuo for photons as predicted by several quantum gravity models. Based on a set of assumptions on the possible intrinsic spectral and temporal evolution, we obtain competitive lower limits on the quadratic leading order of speed of light modification.

The disagreement between direct late-time measurements of the Hubble constant from the SH0ES collaboration, and early-universe measurements based on the $\Lambda$CDM model from the Planck collaboration might, at least in principle, be explained by new physics in the early universe. Recently, the application of the Effective Field Theory of Large-Scale Structure to the full shape of the power spectrum of the SDSS/BOSS data has revealed a new, rather powerful, way to measure the Hubble constant and the other cosmological parameters from Large-Scale Structure surveys. In light of this, we analyze two models for early universe physics, Early Dark Energy and Rock 'n' Roll, that were designed to significantly ameliorate the Hubble tension. Upon including the information from the full shape to the Planck, BAO, and Supernovae measurements, we find that the degeneracies in the cosmological parameters that were introduced by these models are well broken by the data, so that these two models do not significantly ameliorate the tension.

Gamma-ray bursts (GRBs) of the long-duration class are the most luminous sources of electromagnetic radiation known in the Universe. They are generated by outflows of plasma ejected at near the speed of light by newly formed neutron stars or black holes of stellar mass at cosmological distances. Prompt flashes of MeV gamma rays are followed by longer-lasting afterglow emission from radio waves to GeV gamma rays, due to synchrotron radiation by energetic electrons in accompanying shock waves. Although emission of gamma rays at even higher, TeV energies by other radiation mechanisms had been theoretically predicted, it had never been detected previously. Here we report the clear detection of GRB 190114C in the TeV band, achieved after many years of dedicated searches for TeV emission from GRBs. Gamma rays in the energy range 0.2--1 TeV are observed from about 1 minute after the burst (at more than 50 standard deviations in the first 20 minutes). This unambiguously reveals a new emission component in the afterglow of a GRB, whose power is comparable to that of the synchrotron component. The observed similarity in the radiated power and temporal behaviour of the TeV and X-ray bands points to processes such as inverse Compton radiation as the mechanism of the TeV emission, while processes such as synchrotron emission by ultrahigh-energy protons are disfavoured due to their low radiative efficiency.

We present a critical assessment of the SN1987A supernova cooling bound on axions and other light particles. Core-collapse simulations used in the literature to substantiate the bound omitted from the calculation the envelope exterior to the proto-neutron star (PNS). As a result, the only source of neutrinos in these simulations was, by construction, a cooling PNS. We show that if the canonical delayed neutrino mechanism failed to explode SN1987A, and if the pre-collapse star was rotating, then an accretion disk would form that could explain the late-time ($t\gtrsim5$ sec) neutrino events. Such accretion disk would be a natural feature if SN1987A was a collapse-induced thermonuclear explosion. Axions do not cool the disk and do not affect its neutrino output, provided the disk is optically-thin to neutrinos, as it naturally is. These considerations cast doubt on the supernova cooling bound.

Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

An accurate theoretical template for the galaxy power spectrum is a key for the success of ongoing and future spectroscopic surveys. We examine to what extent the Effective Field Theory of Large Scale Structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a hundred times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge task was to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within sub-percent statistical errors corresponding to the total simulation volume.

We apply the Effective Field Theory of Large-Scale Structure to analyze the $w$CDM cosmological model. By using the full shape of the power spectrum and the BAO post-reconstruction measurements from BOSS, the Supernovae from Pantheon, and a prior from BBN, we set the competitive CMB-independent limit $w=-1.046_{-0.052}^{+0.055}$ at $68\%$ C.L.. After adding the Planck CMB data, we find $w=-1.023_{-0.030}^{+0.033}$ at $68\%$ C.L.. Our results are obtained using PyBird, a new, fast Python-based code which we make publicly available.

We consider gravitational wave echoes from black hole merger ringdowns which may arise from local Lorentz symmetry violations that modify graviton dispersion relations. If the corrections are sufficiently soft so they do not remove the horizon, the reflection of the infalling waves which trigger the echoes is very weak. As an example, we look at the dispersion relation corrected by roton-like operators depending only on spatial momenta, in Gullstrand-Painlev\'e coordinates. The near-horizon regions of a black hole do become reflective, but only very weakly. The resulting "bounces" of infalling waves can yield repetitive gravity wave emissions but their power is very small. This implies that to see any echoes from black holes we really need an egregious departure from either standard GR or effective field theory, or both.

The precision of the cosmological data allows us to accurately approximate the predictions for cosmological observables by Taylor expanding up to a low order the dependence on the cosmological parameters around a reference cosmology. By applying this observation to the redshift-space one-loop galaxy power spectrum of the Effective Field Theory of Large-Scale Structure, we analyze the BOSS DR12 data by scanning over all the parameters of $\Lambda$CDM cosmology with massive neutrinos. We impose several sets of priors, the widest of which is just a Big Bang Nucleosynthesis prior on the current fractional energy density of baryons, $\Omega_b h^2$, and a bound on the sum of neutrino masses to be less than 0.9 eV. In this case we measure the primordial amplitude of the power spectrum, $A_s$, the abundance of matter, $\Omega_m$, the Hubble parameter, $H_0$, and the tilt of the primordial power spectrum, $n_s$, to about $19\%$, $5.7\%$, $2.2\%$ and $7.3\%$ respectively, obtaining $\ln ( 10^{10} A_s) =2.91\pm 0.19$, $\Omega_m=0.314\pm 0.018$, $H_0=68.7\pm 1.5$ km/(s Mpc) and $n_s=0.979\pm 0.071$ at $68\%$ confidence level. A public code is released with this preprint.

Core-collapse supernovae presumably explode because trapped neutrinos push the material out of the stellar envelope. This process is directly controlled by the weak scale $v$: we argue that supernova explosions happen only if fundamental constants are tuned within a factor of few as $v \sim \Lambda_{\rm QCD}^{3/4} M_{\rm Pl}^{1/4}$, such that neutrinos are trapped in supernovae for a time comparable to the gravitational time-scale. We provide analytic arguments and simulations in spherical approximation, that need to be validated by more comprehensive simulations. The above result can be important for fundamental physics, because core-collapse supernova explosions seem anthropically needed, as they spread intermediate-mass nuclei presumably necessary for `life'. We also study stellar burning, finding that it does not provide anthropic boundaries on $v$.

The Effective Field Theory of Large-Scale Structure (EFTofLSS) is a formalism that allows us to predict the clustering of Cosmological Large-Scale Structure in the mildly non-linear regime in an accurate and reliable way. After validating our technique against several sets of numerical simulations, we perform the analysis for the cosmological parameters of the DR12 BOSS data. We assume $\Lambda$CDM, a fixed value of the baryon/dark-matter ratio, $\Omega_b/\Omega_c$, and of the tilt of the primordial power spectrum, $n_s$, and no significant input from numerical simulations. By using the one-loop power spectrum multipoles, we measure the primordial amplitude of the power spectrum, $A_s$, the abundance of matter, $\Omega_m$, and the Hubble parameter, $H_0$, to about $13\%$, $3.2\%$ and $3.2\%$ respectively, obtaining $\ln(10^{10}As)=2.72\pm 0.13$, $\Omega_m=0.309\pm 0.010$, $H_0=68.5\pm 2.2$ km/(s Mpc) at 68\% confidence level. If we then add a CMB prior on the sound horizon, the error bar on $H_0$ is reduced to $1.6\%$. These results are a substantial qualitative and quantitative improvement with respect to former analyses, and suggest that the EFTofLSS is a powerful instrument to extract cosmological information from Large-Scale Structure.

Neutrinos can gain mass from coupling to an ultralight field in slow roll. When such a field is displaced from its minimum, its vev acts just like the Higgs vev in spontaneous symmetry breaking. Although these masses may eventually vanish, they do it over a very long time. The theory is technically natural, with the ultralight field-dependent part being the right-handed Majorana mass. The mass variation induced by the field correlates with the cosmological evolution. The change of the mass term changes the mixing matrix, and therefore suppresses the fraction of sterile neutrinos at earlier times and increases it at later times. Since the issue of quantum gravity corrections to field theories with large field variations remains open, this framework may give an observational handle on the Weak Gravity Conjecture.

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

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

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Black holes with masses of $\rm 10^6-10^9~M_{\odot}$ dwell in the centers of most galaxies, but their formation mechanisms are not well known. A subdominant dissipative component of dark matter with similar properties to the ordinary baryons, known as mirror dark matter, may collapse to form massive black holes during the epoch of first galaxies formation. In this study, we explore the possibility of massive black hole formation via this alternative scenario. We perform three-dimensional cosmological simulations for four distinct halos and compare their thermal, chemical and dynamical evolution in both the ordinary and the mirror sectors. We find that the collapse of halos is significantly delayed in the mirror sector due to the lack of $\rm H_2$ cooling and only halos with masses above $ \rm \geq 10^7~ M_{\odot}$ are formed. Overall, the mass inflow rates are $\rm \geq 10^{-2}~M_{\odot}/yr$ and there is less fragmentation. This suggests that the conditions for the formation of massive objects, including black holes, are more favorable in the mirror sector.

We study the ratio of neutrino-proton elastic scattering to inverse beta decay event counts, measurable in a scintillation detector like JUNO, as a key observable for identifying the explosion mechanism of a galactic core-collapse supernova. If the supernova is not powered by the core but rather, e.g., by collapse-induced thermonuclear explosion, then a prolonged period of accretion-dominated neutrino luminosity is predicted. Using 1D numerical simulations, we show that the distinct resulting flavour composition of the neutrino burst can be tested in JUNO with high significance, overcoming theoretical uncertainties in the progenitor star profile and equation of state.

We present a family of consistent quantum field theories of monodromy quintessence in strong coupling, which can serve as benchmarks in modeling dark energy different from cosmological constant. These theories have discrete gauge symmetries which can protect them from quantum field theory and quantum gravity corrections, both perturbative and nonperturbative. The strong coupling effects, at scales $\ga {\rm mm}^{-1}$, flatten the potential and activate operators with higher powers of derivatives. The predicted equation of state is close to, but not exactly equal to $-1$, thus being within reach of the (near!) future programs to explore the nature of dark energy.

The observation of an absorption feature in the 21 cm spectrum at redshift $z\approx 17$ implies bounds on Dark Matter annihilations for a broad range of masses, given that significant heating of the intergalactic medium would have erased such feature. The resulting bounds on the DM annihilation cross sections are comparable to the strongest ones from all other observables.

We present a flavor and energy inference analysis for each high-energy neutrino event observed by the IceCube observatory during six years of data taking. Our goal is to obtain, for the first time, an estimate of the posterior probability distribution for the most relevant properties, such as the neutrino energy and flavor, of the neutrino-nucleon interactions producing shower and track events in the IceCube detector. For each event the main observables in the IceCube detector are the deposited energy and the event topology (showers or tracks) produced by the Cherenkov light by the transit through a medium of charged particles created in neutrino interactions. It is crucial to reconstruct from these observables the properties of the neutrino which generated such event. Here we describe how to achieve this goal using Bayesian inference and Markov chain Monte Carlo methods.

We present a simple effective field theory formulation of a general family of single field flux monodromy models for which strong coupling effects at large field values can flatten the potential and activate operators with higher powers of derivatives. These models are radiatively and non-perturbatively stable and can easily sustain $\ga 60$ efolds of inflation. The dynamics combines features of both large field chaotic inflation and $k$-inflation, both of which can suppress the tensor amplitude. Reducing the tensor-scalar ratio below the observational bound $r \lesssim 0.1$ while keeping the scalar spectral index $n_s$ within experimental bounds either yields equilateral nongaussianity $f_{NL}^{eq} \simeq {\cal O}(1)$, close to the current observational bounds, or ultimately gives very small $r$.

We show that a subdominant component of dissipative dark matter resembling the Standard Model can form many intermediate-mass black hole seeds during the first structure formation epoch. We also observe that, in the presence of this matter sector, the black holes will grow at a much faster rate with respect to the ordinary case. These facts can explain the observed abundance of supermassive black holes feeding high-redshift quasars. The scenario will have interesting observational consequences for dark substructures and gravitational wave production.

Some recent studies exposed rather strong statistical evidence of in-vacuo-dispersion-like spectral lags for gamma-ray bursts (GRBs), a linear correlation between time of observation and energy of GRB particles. Those results focused on testing in-vacuo dispersion for the most energetic GRB particles, and in particular only included photons with energy at emission greater than 40 GeV. We here extend the window of the statistical analysis down to 5 GeV and find results that are consistent with what had been previously noticed at higher energies.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

Recently two of the authors proposed a mechanism of vacuum energy sequester as a means of protecting the observable cosmological constant from quantum radiative corrections. The original proposal was based on using global Lagrange multipliers, but later a local formulation was provided. Subsequently other interesting claims of a different non-local approach to the cosmological constant problem were made, based again on global Lagrange multipliers. We examine some of these proposals and find their mutual relationship. We explain that the proposals which do not treat the cosmological constant counterterm as a dynamical variable require fine tunings to have acceptable solutions. Furthermore, the counterterm often needs to be retuned at every order in the loop expansion to cancel the radiative corrections to the cosmological constant, just like in standard GR. These observations are an important reminder of just how the proposal of vacuum energy sequester avoids such problems.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We discuss how to formulate a quantum field theory of dark energy interacting with dark matter. We show that the proposals based on the assumption that dark matter is made up of heavy particles with masses which are very sensitive to the value of dark energy are strongly constrained. Quintessence-generated long range forces and radiative stability of the quintessence potential require that such dark matter and dark energy are completely decoupled. However, if dark energy and a fraction of dark matter are very light axions, they can have significant mixings which are radiatively stable and perfectly consistent with quantum field theory. Such models can naturally occur in multi-axion realizations of monodromies. The mixings yield interesting signatures which are observable and are within current cosmological limits but could be constrained further by future observations.

Over the last 15 years there has been considerable interest in the possibility of quantum-gravity-induced in-vacuo dispersion, the possibility that spacetime itself might behave essentially like a dispersive medium for particle propagation. Two very recent studies have exposed what might be in-vacuo dispersion features for GRB (gamma-ray-burst) neutrinos of energy in the range of 100 TeV and for GRB photons with energy in the range of 10 GeV. We here show that these two features are roughly compatible with a description such that the same effects apply over 4 orders of magnitude in energy. We also characterize quantitatively how rare it would be for such features to arise accidentally, as a result of (still unknown) aspects of the mechanisms producing photons at GRBs or as a result of background neutrinos accidentally fitting the profile of a GRB neutrino affected by in-vacuo dispersion.

Momentum-space curvature, which is expected in some approaches to the quantum-gravity problem, can produce dual redshift, a feature which introduces energy dependence of the travel times of ultrarelativistic particles, and dual lensing, a feature which mainly affects the direction of observation of particles. In our recent arXiv:1605.00496 we explored the possibility that dual redshift might be relevant in the analysis of IceCube neutrinos, obtaining results which are preliminarily encouraging. Here we explore the possibility that also dual lensing might play a role in the analysis of IceCube neutrinos. In doing so we also investigate issues which are of broader interest, such as the possibility of estimating the contribution by background neutrinos and some noteworthy differences between candidate "early neutrinos" and candidate "late neutrinos".

We explore general scalar-tensor models in the presence of a kinetic mixing between matter and the scalar field, which we call Kinetic Matter Mixing. In the frame where gravity is de-mixed from the scalar this is due to disformal couplings of matter species to the gravitational sector, with disformal coefficients that depend on the gradient of the scalar field. In the frame where matter is minimally coupled, it originates from the so-called beyond Horndeski quadratic Lagrangian. We extend the Effective Theory of Interacting Dark Energy by allowing disformal coupling coefficients to depend on the gradient of the scalar field as well. In this very general approach, we derive the conditions to avoid ghost and gradient instabilities and we define Kinetic Matter Mixing independently of the frame metric used to described the action. We study its phenomenological consequences for a $\Lambda$CDM background evolution, first analytically on small scales. Then, we compute the matter power spectrum and the angular spectra of the CMB anisotropies and the CMB lensing potential, on all scales. We employ the public version of COOP, a numerical Einstein-Boltzmann solver that implements very general scalar-tensor modifications of gravity. Rather uniquely, Kinetic Matter Mixing weakens gravity on short scales, predicting a lower $\sigma_8$ with respect to the $\Lambda$CDM case. We propose this as a possible solution to the tension between the CMB best-fit model and low-redshift observables.

Two recent publications have reported intriguing analyses, tentatively suggesting that some aspects of IceCube data might be manifestations of quantum-gravity-modified laws of propagation for neutrinos. We here propose a strategy of data analysis which has the advantage of being applicable to several alternative possibilities for the laws of propagation of neutrinos in a quantum spacetime. In all scenarios here of interest one should find a correlation between the energy of an observed neutrino and the difference between the time of observation of that neutrino and the trigger time of a GRB. We select accordingly some GRB-neutrino candidates among IceCube events, and our data analysis finds a rather strong such correlation. This sort of studies naturally lends itself to the introduction of a "false alarm probability", which for our analysis we estimate conservatively to be of 1%. We therefore argue that our findings should motivate a vigorous program of investigation following the strategy here advocated.

Ultralight axions which couple sufficiently strongly to photons can leave imprints on the sky at diverse frequencies by mixing with cosmic light in the presence of background magnetic fields. We explore such direction dependent grey-body distortions of the CMB spectrum, enhanced by resonant conditions in the IGM plasma. We also find that if such axions are produced in the early universe and represent a subdominant dark radiation component today, they could convert into X-rays in supervoids, and brighten them at X-ray frequencies.

We analyze primordial non-Gaussianity in single field inflationary models when the tensor/scalar ratio is large. Our results show that detectable levels of non-Gaussianity $f_{NL} \sim 50$ are still possible in the simplest class of models described by the effective theory of inflation. However, the shape is very tightly constrained, making a sharp prediction that could be confirmed or falsified by a future detection of non-Gaussianity.

We study the dynamics of parallel brane-brane and brane-antibrane scattering in string theory in flat spacetime, focusing on the pair production of open strings that stretch between the branes. We are particularly interested in the case of scattering at small impact parameter $b < l_s$, where there is a tachyon in the spectrum when a brane and an antibrane approach within a string length. Our conclusion is that despite the tachyon, branes and antibranes can pass through each other with only a very small probability of annihilating, so long as $g_s$ is small and the relative velocity $v$ is neither too small nor too close to 1. Our analysis is relevant also to the case of charged open string production in world-volume electric fields, and we make use of this T-dual scenario in our analysis. We briefly discuss the application of our results to a stringy model of inflation involving moving branes.

We analyze an inflationary model in which part of the power in density perturbations arises due to particle production. The amount of particle production is modulated by an auxiliary field. Given an initial gradient for the auxiliary field, this model produces a hemispherical power asymmetry and a suppression of power at low multipoles similar to those observed by WMAP and Planck in the CMB temperature. It also predicts an additive contribution to $\delta T$ with support only at very small $l$ that is aligned with the direction of the power asymmetry and has a definite sign, as well as small oscillations in the power spectrum at all $l$.

A theory of the quasidilaton is an extension of massive gravity by a scalar field, nonlinearly realizing a certain new global symmetry of the Lagrangian. It has been shown that unlike pure massive gravity, this theory does admit homogeneous and isotropic spatially flat solutions. Among the latter, selfaccelerated solutions attract a special attention. Previous studies of perturbations, performed in the decoupling limit, revealed one healthy scalar mode, while the second relevant scalar was not captured in that limit. Here we study full cosmological perturbations above the simplest selfaccelerated background. We show that the fluctuations of a mixed state of the quasidilaton and the helicity-0 graviton necessarily have a negative kinetic term at short distances, making this background unphysical. In addition, these cosmologies exhibit an order one sensitivity to higher dimensional terms suppressed by an energy scale that is parametrically higher than the strong coupling scale of the quasidilaton effective theory: such terms include Galileons, Goldstone-like selfinteractions and derivatives of the quasidilaton coupled to curvature, none of which introduce extra Ostrogradsky states. As one consequence, cosmology at the Hubble distances for this particular class of solutions depends on an unknown extension of the quasidilaton below its strong coupling distance scale. We note that non-FRW solutions that are similar to those of pure massive gravity should not necessarily suffer from these pathologies.

We study perturbations around some cosmological backgrounds in the dRGT theory of massive gravity. We develop a general formalism to calculate the perturbations around any background. We derive the Lagrangian for fluctuations in the small scale limit, and for the open FRW solution we repeat the analysis around the full background. We find that the perturbations display similar properties: the longitudinal modes of the massive graviton are instantaneous at quadratic level, but they acquire a time-kinetic term at cubic order.

Higher-form flux that extends in all 3+1 dimensions of spacetime is a source of positive vacuum energy that can drive meta-stable eternal inflation. If the flux also threads compact extra dimensions, the spontaneous nucleation of a bubble of brane charged under the flux can trigger a classical cascade that steadily unwinds many units of flux, gradually decreasing the vacuum energy while inflating the bubble, until the cascade ends in the self-annihilation of the brane into radiation. With an initial number of flux quanta Q_{0} \simgeq N, this can result in N efolds of inflationary expansion while producing a scale-invariant spectrum of adiabatic density perturbations with amplitude and tilt consistent with observation. The power spectrum has an oscillatory component that does not decay away during inflation, relatively large tensor power, and interesting non-Gaussianities. Unwinding inflation fits naturally into the string landscape, and our preliminary conclusion is that consistency with observation can be attained without fine-tuning the string parameters. The initial conditions necessary for the unwinding phase are produced automatically by bubble formation, so long as the critical radius of the bubble is smaller than at least one of the compact dimensions threaded by flux.

General Relativity (GR), with or without matter fields, admits a natural extension to a scale invariant theory that requires a dilaton. Here we show that the recently formulated massive GR, minimally coupled to matter, possesses a new global symmetry related to scaling of the reference coordinates w.r.t. the physical ones. The field enforcing this symmetry, dubbed here quasi-dilaton, coincides with an ordinary dilaton if only pure gravity is considered, but differs from it when the matter Lagrangian is present. We study: (1) Theoretical consistency of massive GR with the quasi-dilaton; (2) Consistency with observations for spherically symmetric sources on (nearly) flat backgrounds; (3) Cosmological implications of this theory. We find that: (I) The theory with the quasi-dilaton is as consistent as massive GR is. (II) The Vainshtein mechanism is generically retained, owing to the fact that in the decoupling limit there is an enhanced symmetry, which turns the quasi-dilaton into a second galileon, consistently coupled to a tensor field. (III) Unlike in massive GR, there exist flat FRW solutions. In particular, we find self-accelerated solutions and discuss their quadratic perturbations. These solutions are testable by virtue of the different effective Newton's constants that govern the Hubble expansion and structure growth.

When electric-type flux threads compact extra dimensions, a quantum nucleation event can break a flux line and initiate a cascade that unwinds many units of flux. Here, we present a novel mechanism for inflation based on this phenomenon. From the 4D point of view, the cascade begins with the formation of a bubble containing an open Friedmann-Robertson-Walker cosmology, but the vacuum energy inside the bubble is initially only slightly reduced, and subsequently decreases gradually throughout the cascade. If the initial flux number Q_0 ~ O(100), during the cascade the universe can undergo N ~ 60 efolds of inflationary expansion with gradually decreasing Hubble constant, producing a nearly scale-invariant spectrum of adiabatic density perturbations with amplitude and tilt consistent with observation, and a potentially observable level of non-Gaussianity and tensor modes. The power spectrum has a small oscillatory component that does not decay away during inflation, with a period set approximately by the light-crossing time of the compact dimension(s). Since the ingredients are fluxes threading compact dimensions, this mechanism fits naturally into the string landscape, but does not appear to suffer from the eta problem or require fine-tuning (beyond the usual anthropic requirement of small vacuum energy after reheating).

We prove that, in a generic single-field model, the consistency relation for the 3-point function in the squeezed limit receives corrections that vanish quadratically in the ratio of the momenta, i.e. as (k_L/k_S)^2. This implies that a detection of a bispectrum signal going as 1/k_L^2 in the squeezed limit, that is suppressed only by one power of k_L compared with the local shape, would rule out all single-field models. The absence of this kind of terms in the bispectrum holds also for multifield models, but only if all the fields have a mass much smaller than H. The detection of any scale dependence of the bias, for scales much larger than the size of the haloes, would disprove all single-field models. We comment on the regime of squeezing that can be probed by realistic surveys.

We explore the cosmological solutions of a recently proposed extension of General Relativity with a Lorentz-invariant mass term. We show that the same constraint that removes the Boulware-Deser ghost in this theory also prohibits the existence of homogeneous and isotropic cosmological solutions. Nevertheless, within domains of the size of inverse graviton mass we find approximately homogeneous and isotropic solutions that can well describe the past and present of the Universe. At energy densities above a certain crossover value, these solutions approximate the standard FRW evolution with great accuracy. As the Universe evolves and density drops below the crossover value the inhomogeneities become more and more pronounced. In the low density regime each domain of the size of the inverse graviton mass has essentially non-FRW cosmology. This scenario imposes an upper bound on the graviton mass, which we roughly estimate to be an order of magnitude below the present-day value of the Hubble parameter. The bound becomes especially restrictive if one utilizes an exact self-accelerated solution that this theory offers. Although the above are robust predictions of massive gravity with an explicit mass term, we point out that if the mass parameter emerges from some additional scalar field condensation, the constraint no longer forbids the homogeneous and isotropic cosmologies. In the latter case, there will exist an extra light scalar field at cosmological scales, which is screened by the Vainshtein mechanism at shorter distances.

We study the nonlinear evolution of density perturbations in cosmologies where the late-time accelerated expansion is driven by a quintessence field with vanishing speed of sound. For these models matter and quintessence perturbations are comoving and it is possible to write a single continuity equation for the total density fluctuations given by a weighted sum of the two components. Including the Euler equation for the common velocity field we solve the evolution equations for the nonlinear, total density power spectrum in the Time-Renormalization Group approach. In fact any cosmological observable is directly related by gravity only to the total density perturbations, with the two components being individually unobservable. We show that the clustering of quintessence perturbations induces small corrections with respect to the nonlinear evolution of power spectrum in smooth quintessence models described by the same equation of state. Such small corrections, however, contrast with the large effect of a vanishing speed of sound on the linear growth function at low redshift. For this reason, models with the same normalization of the linear density power spectrum can present significantly different nonlinear corrections depending on the value of the sound speed. Although such differences vanish in the w \rightarrow -1, {\Lambda}CDM limit, we argue that the relation between linear and nonlinear growth of structures should be properly taken into account in constraining models with inhomogeneous dark energy.

The abundance of collapsed objects in the universe, or halo mass function, is an important theoretical tool in studying the effects of primordially generated non-Gaussianities on the large scale structure. The non-Gaussian mass function has been calculated by several authors in different ways, typically by exploiting the smallness of certain parameters which naturally appear in the calculation, to set up a perturbative expansion. We improve upon the existing results for the mass function by combining path integral methods and saddle point techniques (which have been separately applied in previous approaches). Additionally, we carefully account for the various scale dependent combinations of small parameters which appear. Some of these combinations in fact become of order unity for large mass scales and at high redshifts, and must therefore be treated non-perturbatively. Our approach allows us to do this, and to also account for multi-scale density correlations which appear in the calculation. We thus derive an accurate expression for the mass function which is based on approximations that are valid over a larger range of mass scales and redshifts than those of other authors. By tracking the terms ignored in the analysis, we estimate theoretical errors for our result and also for the results of others. We also discuss the complications introduced by the choice of smoothing filter function, which we take to be a top-hat in real space, and which leads to the dominant errors in our expression. Finally, we present a detailed comparison between the various expressions for the mass functions, exploring the accuracy and range of validity of each.

Primordial non-Gaussianity (NG) affects the large scale structure (LSS) of the universe by leaving an imprint on the distribution of matter at late times. Much attention has been focused on using the distribution of collapsed objects (i.e. dark matter halos and the galaxies and galaxy clusters that reside in them) to probe primordial NG. An equally interesting and complementary probe however is the abundance of extended underdense regions or voids in the LSS. The calculation of the abundance of voids using the excursion set formalism in the presence of primordial NG is subject to the same technical issues as the one for halos, which were discussed e.g. in arXiv:1005.1203. However, unlike the excursion set problem for halos which involved random walks in the presence of one barrier $\delta_c$, the void excursion set problem involves two barriers $\delta_v$ and $\delta_c$. This leads to a new complication introduced by what is called the "void-in-cloud" effect discussed in the literature, which is unique to the case of voids. We explore a path integral approach which allows us to carefully account for all these issues, leading to a rigorous derivation of the effects of primordial NG on void abundances. The void-in-cloud issue in particular makes the calculation conceptually rather different from the one for halos. However, we show that its final effect can be described by a simple yet accurate approximation. Our final void abundance function is valid on larger scales than the expressions of other authors, while being broadly in agreement with those expressions on smaller scales.

We study the consequences of imposing an approximate Galilean symmetry on the Effective Theory of Inflation, the theory of small perturbations around the inflationary background. This approach allows us to study the effect of operators with two derivatives on each field, which can be the leading interactions due to non-renormalization properties of the Galilean Lagrangian. In this case cubic non-Gaussianities are given by three independent operators, containing up to six derivatives, two with a shape close to equilateral and one peaking on flattened isosceles triangles. The four-point function is larger than in models with small speed of sound and potentially observable with the Planck satellite.

We study the spherical collapse model in the presence of quintessence with negligible speed of sound. This case is particularly motivated for w<-1 as it is required by stability. As pressure gradients are negligible, quintessence follows dark matter during the collapse. The spherical overdensity behaves as a separate closed FLRW universe, so that its evolution can be studied exactly. We derive the critical overdensity for collapse and we use the extended Press-Schechter theory to study how the clustering of quintessence affects the dark matter mass function. The effect is dominated by the modification of the linear dark matter growth function. A larger effect occurs on the total mass function, which includes the quintessence overdensities. Indeed, here quintessence constitutes a third component of virialized objects, together with baryons and dark matter, and contributes to the total halo mass by a fraction ~ (1+w) Omega_Q / Omega_m. This gives a distinctive modification of the total mass function at low redshift.

We calculate the Cosmic Microwave Background anisotropy bispectrum on large angular scales in the absence of primordial non-Gaussianities, assuming exact matter dominance and extending at second order the classic Sachs-Wolfe result \delta T/T=\Phi/3. The calculation is done in Poisson gauge. Besides intrinsic contributions calculated at last scattering, one must consider integrated effects. These are associated to lensing, and to the time dependence of the potentials (Rees-Sciama) and of the vector and tensor components of the metric generated at second order. The bispectrum is explicitly computed in the flat-sky approximation. It scales as l^(-4) in the scale invariant limit and the shape dependence of its various contributions is represented in 3d plots. Although all the contributions to the bispectrum are parametrically of the same order, the full bispectrum is dominated by lensing. In the squeezed limit it corresponds to f_NL^local = -1/6 - cos(2 \theta), where \theta is the angle between the short and the long modes; the angle dependent contribution comes from lensing. In the equilateral limit it corresponds to f_NL^equil ~ 3.13.

These lectures are intended to provide a brief pedagogical review of dark matter for the newcomer to the subject. We begin with a discussion of the astrophysical evidence for dark matter. The standard weakly-interacting massive particle (WIMP) scenario--the motivation, particle models, and detection techniques--is then reviewed. We provide a brief sampling of some recent variations to the standard WIMP scenario as well as some alternatives (axions and sterile neutrinos). Exercises are provided for the reader.

We study generic single-field dark energy models, by a parametrization of the most general theory of their perturbations around a given background, including higher derivative terms. In appropriate limits this approach reproduces standard quintessence, k-essence and ghost condensation. We find no general pathology associated to an equation of state w_Q < -1 or in crossing the phantom divide w_Q = -1. Stability requires that the w_Q < -1 side of dark energy behaves, on cosmological scales, as a k-essence fluid with a virtually zero speed of sound. This implies that one should set the speed of sound to zero when comparing with data models with w_Q < -1 or crossing the phantom divide. We summarize the theoretical and stability constraints on the quintessential plane (1+w_Q) vs. speed of sound squared.

In this paper we present a complete computation of the Cosmic Microwave Background (CMB) anisotropies up to third order from gravitational perturbations accounting for scalar, vector and tensor perturbations. We then specify our results to the large scale limit, providing the evolution of the gravitational potentials in a flat universe filled with matter and cosmological constant which characterizes the Integrated Sachs-Wolfe effect. As a byproduct in the large scale approximation we are able to give non-perturbative solutions for the photon geodesic equations. Our results are the first step to provide a complete theoretical prediction for cubic non-linearities which are particularly relevant for characterizing the level of non-Gaussianity in the CMB through the detection of the four-point angular connected correlation function (trispectrum). For this purpose we also allow for generic initial conditions due to primordial non-Gaussianity.

This work investigates the dark matters structures that form on the smallest cosmological scales. We find that the types and abundances of structures which form at approximately Earth-mass scales are very sensitive to the nature of dark matter. We explore various candidates for dark matter and determine the corresponding properties of small-scale structure. In particular, we discuss possibilities for indirect detection of dark matter through small-scale structure, and comment on the potential of these methods for discriminating between dark matter candidates.