It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

We present a new numerical algorithm and code, ${\sf GFiRe}$, for solving the non-linear evolution of Abelian gauge fields coupled to complex scalar fields in homogeneous and isotropic spacetimes. We adopt a hybrid approach to solving the system: the spatial derivatives are discretized using standard Lattice Gauge Field Theory techniques, whereas the time evolution of the fields and scale factor is implemented with explicit, composite, symplectic integrators. An important property of our compound algorithm is that the discretized Gauss constraint is respected exactly, regardless of the order of the symplectic integrator. This remains true even when the background expansion is computed "self-consistently"; that is, when the expansion history is computed using spatial averaged components of the energy momentum tensor in the simulation volume. Hence, our code can also be used in cases where the fields dominate the energy density of the universe, for example, during reheating after inflation. We test the algorithm in scenarios of reheating where the inflaton is a complex scalar field with a potential $\propto(2|\varphi|^2-v^2)^2$ and is coupled to an Abelian gauge field. Tracing the evolution of the system through complex dynamics (including resonant excitation of fields, backreaction, formation of solitons, and changes in the equation of state) in a self-consistently expanding universe, we find the energy conservation violation ($<10^{-4}$) to be very stable and the Gauss constraint violation ($<10^{-6}$) to be dominated by differencing noise.

We calculate the curvature power spectrum sourced by spectator fields that are excited repeatedly and non-adiabatically during inflation. In the absence of detailed information of the nature of spectator field interactions, we consider an ensemble of models with intervals between the repeated interactions and interaction strengths drawn from simple probabilistic distributions. We show that the curvature power spectra of each member of the ensemble shows rich structure with many features, and there is a large variability between different realizations of the same ensemble. Such features can be probed by the cosmic microwave background (CMB) and large scale structure observations. They can also have implications for primordial black hole formation and CMB spectral distortions. The geometric random walk behavior of the spectator field allows us to calculate the ensemble-averaged power spectrum of curvature perturbations semi-analytically. For sufficiently large stochastic sourcing, the ensemble-averaged power spectrum shows a scale dependence arising from the time spent by modes outside the horizon during the period of particle production, in spite of there being no preferred scale in the underlying model. We find that the magnitude of the ensemble-averaged power spectrum overestimates the typical power spectra in the ensemble because the ensemble distribution of the power spectra is highly non-Gaussian with fat tails.

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec (kpc) scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter we investigate how the first galaxies are assembled in a `fuzzy' dark matter (FDM) cosmology where dark matter is an ultralight $\sim 10^{-22}$~eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical soliton-like cores which are unstable under gravity and collapse into kpc-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter -- a smoking gun signature of FDM.

In this white paper, we discuss the prospects for characterizing and identifying dark matter using gravitational waves, covering a wide range of dark matter candidate types and signals. We argue that present and upcoming gravitational wave probes offer unprecedented opportunities for unraveling the nature of dark matter and we identify the most urgent challenges and open problems with the aim of encouraging a strong community effort at the interface between these two exciting fields of research.

We investigate the formation, gravitational clustering and interactions of solitons in a self-interacting, non-relativistic scalar field in an expanding universe. Rapid formation of large number of solitons is driven by attractive self-interactions of the field, whereas the slower clustering of solitons is driven by gravitational forces. Driven closer together by gravity, we see a rich plethora of dynamics in the soliton "gas" including mergers, scatterings and formation of soliton binaries. The numerical simulations are complemented by analytic calculations and estimates of (i) the relevant instability length and time scales, (ii) individual soliton profiles and their stability, (iii) number density of produced solitons, and (iv) the two point correlation function of soliton positions as evidence for gravitational clustering.

We present a detailed investigation of a sub-dominant oscillating scalar field ('early dark energy', EDE) in the context of resolving the Hubble tension. Consistent with earlier work, but without relying on fluid approximations, we find that a scalar field frozen due to Hubble friction until ${\rm log}_{10}(z_c)\sim3.5$, reaching $\rho_{\rm EDE}(z_c)/\rho_{\rm tot}\sim10$%, and diluting faster than matter afterwards can bring cosmic microwave background (CMB), baryonic acoustic oscillations, supernovae luminosity distances, and the late-time estimate of the Hubble constant from the SH0ES collaboration into agreement. A scalar field potential which scales as $V(\phi) \propto \phi^{2n}$ with $2\lesssim n\lesssim 3.4$ around the minimum is preferred at the 68% confidence level, and the {\em Planck} polarization places additional constraints on the dynamics of perturbations in the scalar field. In particular, the data prefers a potential which flattens at large field displacements. An MCMC analysis of mock data shows that the next-generation CMB observations (i.e., CMB-S4) can unambiguously detect the presence of the EDE at very high significance. This projected sensitivity to the EDE dynamics is mainly driven by improved measurements of the $E$-mode polarization. We also explore new observational signatures of EDE scalar field dynamics: (i) We find that depending on the strength of the tensor-to-scalar ratio, the presence of the EDE might imply the existence of isocurvature perturbations in the CMB. (ii) We show that a strikingly rapid, scale-dependent growth of EDE field perturbations can result from parametric resonance driven by the anharmonic oscillating field for $n\approx 2$. This instability and ensuing potentially nonlinear, spatially inhomogenoues, dynamics may provide unique signatures of this scenario.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

The axion is a hypothetical, well-motivated dark-matter particle whose existence would explain the lack of charge-parity violation in the strong interaction. In addition to this original motivation, an `axiverse' of ultra-light axions (ULAs) with masses $10^{-33}\,{\rm eV}\lesssim m_{\rm a}\lesssim 10^{-10}\,{\rm eV}$ also emerges from string theory. Depending on the mass, such a ULA contributes to the dark-matter density, or alternatively, behaves like dark energy. At these masses, ULAs' classical wave-like properties are astronomically manifested, potentially mitigating observational tensions within the $\Lambda$CDM paradigm on local-group scales. ULAs also provide signatures on small scales such as suppression of structure, interference patterns and solitons to distinguish them from heavier dark matter candidates. Through their gravitational imprint, ULAs in the presently allowed parameter space furnish a host of observational tests to target in the next decade, altering standard predictions for microwave background anisotropies, galaxy clustering, Lyman-$\alpha$ absorption by neutral hydrogen along quasar sightlines, pulsar timing, and the black-hole mass spectrum.

We numerically investigate the gravitational waves generated by the head-on collision of equal-mass, self-gravitating, real scalar field solitons (oscillatons) as a function of their compactness $\mathcal{C}$. We show that there exist three different possible outcomes for such collisions: (1) an excited stable oscillaton for low $\mathcal{C}$, (2) a merger and formation of a black-hole for intermediate $\mathcal{C}$, and (3) a pre-merger collapse of both oscillatons into individual black-holes for large $\mathcal{C}$. For (1), the excited, aspherical oscillaton continues to emit gravitational waves. For (2), the total energy in gravitational waves emitted increases with compactness, and possesses a maximum which is greater than that from the merger of a pair of equivalent mass black-holes. The initial amplitudes of the quasi-normal modes in the post-merger ring-down in this case are larger than that of corresponding mass black-holes -- potentially a key observable to distinguish black-hole mergers with their scalar mimics. For (3), the gravitational wave output is indistinguishable from a similar mass, black-hole--black-hole merger.

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.

The next generation of instruments designed to measure the polarization of the cosmic microwave background (CMB) will provide a historic opportunity to open the gravitational wave window to the primordial Universe. Through high sensitivity searches for primordial gravitational waves, and tighter limits on the energy released in processes like phase transitions, the CMB polarization data of the next decade has the potential to transform our understanding of the laws of physics underlying the formation of the Universe.

There are two big questions cosmologists would like to answer -- How does the Universe work, and what are its origin and destiny? A long wavelength gravitational wave detector -- with million km interferometer arms, achievable only from space -- gives a unique opportunity to address both of these questions. A sensitive, mHz frequency observatory could use the inspiral and merger of massive black hole binaries as standard sirens, extending our ability to characterize the expansion history of the Universe from the onset of dark energy-domination out to a redshift z ~ 10. A low-frequency detector, furthermore, offers the best chance for discovery of exotic gravitational wave sources, including a primordial stochastic background, that could reveal clues to the origin of our Universe.

We explore non-adiabatic particle production in a de Sitter universe for a scalar spectator field, by allowing the effective mass $m^2(t)$ of this field and the cosmic time interval between non-adiabatic events to vary stochastically. Two main scenarios are considered depending on the (non-stochastic) mass $M$ of the spectator field: the conformal case with $M^2=2H^2$, and the case of a massless field. We make use of the transfer matrix formalism to parametrize the evolution of the system in terms of the "occupation number", and two phases associated with the transfer matrix; these are used to construct the evolution of the spectator field. Assuming short-time interactions approximated by Dirac-delta functions, we numerically track the change of these parameters and the field in all regimes: sub- and super-horizon with weak and strong scattering. In all cases a log-normally distributed field amplitude is observed, and the logarithm of the field amplitude approximately satisfies the properties of a Wiener process outside the horizon. We derive a Fokker-Planck equation for the evolution of the transfer matrix parameters, which allows us to calculate analytically non-trivial distributions and moments in the weak-scattering limit.

We study the scalar and tensor perturbations generated by the fragmentation of the inflaton condensate into oscillons or transients after inflation, using nonlinear classical lattice simulations. Without including the backreaction of metric perturbations, we find that the magnitude of scalar metric perturbations never exceeds a few $\times 10^{-3}$, whereas the maximal strength of the gravitational wave signal today is $\mathcal{O}(10^{-9})$ for standard post-inflationary expansion histories. We provide parameter scalings for the $\alpha$-attractor models of inflation, which can be easily applied to other models. We also discuss the likelihood of primordial black hole formation, as well as conditions under which the gravitational wave signal can be at observationally interesting frequencies and amplitudes. Finally, we provide an upper bound on the frequency of the peak of the gravitational wave signal, which applies to all preheating scenarios.

The Higgs potential appears to be fine-tuned, hence very sensitive to values of other scalar fields that couple to the Higgs. We show that this feature can lead to a new epoch in the early universe featuring violent dynamics coupling the Higgs to a scalar modulus. The oscillating modulus drives tachyonic Higgs particle production. We find a simple parametric understanding of when this process can lead to rapid modulus fragmentation, resulting in gravitational wave production. A nontrivial equation-of-state arising from the nonlinear dynamics also affects the time elapsed from inflation to the CMB, influencing fits of inflationary models. Supersymmetric theories automatically contain useful ingredients for this picture.

At the end of inflation, the inflaton field decays into an initially nonthermal population of relativistic particles which eventually thermalize. We consider the production of dark matter from this relativistic plasma, focusing on the prethermal phase. We find that for a production cross section $\sigma(E)\sim E^n$ with $n> 2$, the present dark matter abundance is produced during the prethermal phase of its progenitors. For $n\le 2$, entropy production during reheating makes the nonthermal contribution to the present dark matter abundance subdominant compared to that produced thermally. As specific examples, we verify that the nonthermal contribution is irrelevant for gravitino production in low scale supersymmetric models ($n=0$) and is dominant for gravitino production in high scale supersymmetry models ($n=6$).

It has been recently suggested that oscillons produced in the early universe from certain asymmetric potentials continue to emit gravitational waves for a number of $e$-folds of expansion after their formation, leading to potentially detectable gravitational wave signals. We revisit this claim by conducting a convergence study using graphics processing unit (GPU)-accelerated lattice simulations and show that numerical errors accumulated with time are significant in low-resolution scenarios, or in scenarios where the run-time causes the resolution to drop below the relevant scales in the problem. Our study determines that the dominant, growing high frequency peak of the gravitational wave signals in the fiducial "hill-top model" in [arXiv:1607.01314] is a numerical artifact. This finding prompts the need for a more careful analysis of the numerical validity of other similar results related to gravitational waves from oscillon dynamics.

Homogeneous oscillations of the inflaton after inflation can be unstable to small spatial perturbations even without coupling to other fields. We show that for inflaton potentials $\propto |\phi|^{2n}$ near $|\phi|=0$ and flatter beyond some $|\phi|=M$, the inflaton condensate oscillations can lead to self-resonance, followed by its complete fragmentation. We find that for non-quadratic minima ($n>1$), shortly after backreaction, the equation of state parameter, $w\rightarrow1/3$. If $M\ll m_{pl}$, radiation domination is established within less than an e-fold of expansion after the end of inflation. In this case self-resonance is efficient and the condensate fragments into transient, localised spherical objects which are unstable and decay, leaving behind them a virialized field with mean kinetic and gradient energies much greater than the potential energy. This end-state yields $w=1/3$. When $M\ll m_{pl}$ we observe slow and steady, self-resonace that can last many {\it e}-folds before backreaction eventually shuts it off, followed by fragmentation and $w\rightarrow 1/3$. We provide analytical estimates for the duration to $w\rightarrow 1/3$ after inflation, which can be used as an upper bound (under certain assumptions) on the duration of the transition between the inflationary and the radiation dominated states of expansion. This upper bound can reduce uncertainties in CMB observables such as the spectral tilt $n_{\rm{s}}$, and the tensor-to-scalar ratio $r$. For quadratic minima ($n=1$), $w\rightarrow0$ regardless of the value of $M$.

We explore non-adiabatic particle production for $N_{\rm f}$ coupled scalar fields in a time-dependent background with stochastically varying effective masses, cross-couplings and intervals between interactions. Under the assumption of weak scattering per interaction, we provide a framework for calculating the typical particle production rates after a large number of interactions. After setting up the framework, for analytic tractability, we consider interactions (effective masses and cross couplings) characterized by series of Dirac-delta functions in time with amplitudes and locations drawn from different distributions. Without assuming that the fields are statistically equivalent, we present closed form results (up to quadratures) for the asymptotic particle production rates for the $N_{\rm f}=1$ and $N_{\rm f}=2$ cases. We also present results for the general $N_{\rm f} >2$ case, but with more restrictive assumptions. We find agreement between our analytic results and direct numerical calculations of the total occupation number of the produced particles, with departures that can be explained in terms of violation of our assumptions. We elucidate the precise connection between the maximum entropy ansatz (MEA) used in Amin & Baumann (2015) and the underlying statistical distribution of the self and cross couplings. We provide and justify a simple to use (MEA-inspired) expression for the particle production rate, which agrees with our more detailed treatment when the parameters characterizing the effective mass and cross-couplings between fields are all comparable to each other. However, deviations are seen when some parameters differ significantly from others. We show that such deviations become negligible for a broad range of parameters when $N_{\rm f}\gg 1$.

Anisotropies in distortions to the frequency spectrum of the cosmic microwave background (CMB) can be created through spatially varying heating processes in the early Universe. For instance, the dissipation of small-scale acoustic modes does create distortion anisotropies, in particular for non-Gaussian primordial perturbations. In this work, we derive approximations that allow describing the associated distortion field. We provide a systematic formulation of the problem using Fourier-space window functions, clarifying and generalizing previous approximations. Our expressions highlight the fact that the amplitudes of the spectral-distortion fluctuations induced by non-Gaussianity depend also on the homogeneous value of those distortions. Absolute measurements are thus required to obtain model-independent distortion constraints on primordial non-Gaussianity. We also include a simple description for the evolution of distortions through photon diffusion, showing that these corrections can usually be neglected. Our formulation provides a systematic framework for computing higher order correlation functions of distortions with CMB temperature anisotropies and can be extended to describe correlations with polarization anisotropies.

We calculate the equation of state after inflation and provide an upper bound on the duration before radiation domination by taking the nonlinear dynamics of the fragmented inflaton field into account. A broad class of single-field inflationary models with observationally consistent flattening of the potential at a scale $M$ away from the origin, $V(\phi)\propto |\phi|^{2n}$ near the origin, and where the couplings to other fields are ignored are included in our analysis. We find that the equation of state parameter $w\rightarrow 0$ for $n=1$ and $w\rightarrow 1/3$ (after sufficient time) for $n\gtrsim 1$. We calculate how the number of $e$-folds to radiation domination depends on both $n$ and $M$ when $M\sim m_{\rm pl}$, whereas when $M\ll m_{\rm pl}$, we find that the duration to radiation domination is negligible. Our results are explained in terms of a linear instability analysis in an expanding universe, scaling arguments, and are supported by detailed 3+1 dimensional lattice simulations. We show how our work significantly reduces the uncertainty in inflationary observables, even after including couplings to additional light fields.

We calculate particle production during inflation and in the early stages of reheating after inflation in models with a charged scalar field coupled to Abelian and non-Abelian gauge fields. A detailed analysis of the power spectra of primordial electric fields, magnetic fields and charge fluctuations at the end of inflation and preheating is provided. We carefully account for the Gauss constraints during inflation and preheating, and clarify the role of the longitudinal components of the electric field. We calculate the timescale for the back-reaction of the produced gauge fields on the inflaton condensate, marking the onset of non-linear evolution of the fields. We provide a prescription for initial conditions for lattice simulations necessary to capture the subsequent nonlinear dynamics. On the observational side, we find that the primordial magnetic fields generated are too small to explain the origin of magnetic fields on galactic scales and the charge fluctuations are well within observational bounds for the models considered in this paper.

We provide a statistical framework for characterizing stochastic particle production in the early universe via a precise correspondence to current conduction in wires with impurities. Our approach is particularly useful when the microphysics is uncertain and the dynamics are complex, but only coarse-grained information is of interest. We study scenarios with multiple interacting fields and derive the evolution of the particle occupation numbers from a Fokker-Planck equation. At late times, the typical occupation numbers grow exponentially which is the analog of Anderson localization for disordered wires. Some statistical features of the occupation numbers show hints of universality in the limit of a large number of interactions and/or a large number of fields. For test cases, excellent agreement is found between our analytic results and numerical simulations.

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope routinely detects the highly dust-absorbed, reddened, and MeV-peaked flat spectrum radio quasar PKS 1830-211 (z=2.507). Its apparent isotropic gamma-ray luminosity (E>100 MeV) averaged over $\sim$ 3 years of observations and peaking on 2010 October 14/15 at 2.9 X 10^{50} erg s^{-1}, makes it among the brightest high-redshift Fermi blazars. No published model with a single lens can account for all of the observed characteristics of this complex system. Based on radio observations, one expects time delayed variability to follow about 25 days after a primary flare, with flux about a factor 1.5 less. Two large gamma-ray flares of PKS 1830-211 have been detected by the LAT in the considered period and no substantial evidence for such a delayed activity was found. This allows us to place a lower limit of about 6 on the gamma rays flux ratio between the two lensed images. Swift XRT observations from a dedicated Target of Opportunity program indicate a hard spectrum and with no significant correlation of X-ray flux with the gamma-ray variability. The spectral energy distribution can be modeled with inverse Compton scattering of thermal photons from the dusty torus. The implications of the LAT data in terms of variability, the lack of evident delayed flare events, and different radio and gamma-ray flux ratios are discussed. Microlensing effects, absorption, size and location of the emitting regions, the complex mass distribution of the system, an energy-dependent inner structure of the source, and flux suppression by the lens galaxy for one image path may be considered as hypotheses for understanding our results.

Our understanding of the state of the universe between the end of inflation and big bang nucleosynthesis (BBN) is incomplete. The dynamics at the end of inflation are rich and a potential source of observational signatures. Reheating, the energy transfer between the inflaton and Standard Model fields (possibly through intermediaries) and their subsequent thermalization, can provide clues to how inflation fits in with known high-energy physics. We provide an overview of our current understanding of the nonperturbative, nonlinear dynamics at the end of inflation, some salient features of realistic particle physics models of reheating, and how the universe reaches a thermal state before BBN. In addition, we review the analytical and numerical tools available in the literature to study preheating and reheating and discuss potential observational signatures from this fascinating era.

The dynamics at the end of inflation can generate an asymmetry between particles and anti-particles of the inflaton field. This asymmetry can be transferred to baryons via decays, generating a baryon asymmetry in our Universe. We explore this idea in detail for a complex inflaton governed by an observationally consistent -"flatter than quadratic"- potential with a weakly broken global U(1) symmetry. We find that most of the inflaton asymmetry is locked in non-topological soliton like configurations (oscillons) produced copiously at the end of inflation. These solitons eventually decay into baryons and generate the observed matter-antimatter asymmetry for a range of model parameters. Through a combination of three dimensional lattice simulations and a detailed linearized analysis, we show how the inflaton asymmetry depends on the fragmentation, the magnitude of the symmetry breaking term and initial conditions at the end of inflation. We discuss the final decay into baryons, but leave a detailed analysis of the inhomogeneous annihilation, reheating and thermalization to future work. As part of our work, we pay particular attention to generating multifield initial conditions for the field fluctuations (including metric perturbations) at the end of inflation for lattice simulations.

Spectral distortions of the cosmic microwave background (CMB) may become a powerful probe of primordial perturbations at small scales. Existing studies of spectral distortions focus almost exclusively on primordial scalar metric perturbations. Similarly, vector and tensor perturbations should source CMB spectral distortions. In this paper, we give general expressions for the effective heating rate caused by these types of perturbations, including previously neglected contributions from polarization states and higher multipoles. We then focus our discussion on the dissipation of tensors, showing that for nearly scale invariant tensor power spectra, the overall distortion is some six orders of magnitudes smaller than from the damping of adiabatic scalar modes. We find simple analytic expressions describing the effective heating rate from tensors using a quasi-tight coupling approximation. In contrast to adiabatic modes, tensors cause heating without additional photon diffusion and thus over a wider range of scales, as recently pointed out by Ota et. al 2014. Our results are in broad agreement with their conclusions, but we find that small-scale modes beyond k< 2x10^4 Mpc^{-1} cannot be neglected, leading to a larger distortion, especially for very blue tensor power spectra. At small scales, also the effect of neutrino damping on the tensor amplitude needs to be included.

In this short note we investigate the ultra-relativistic collisions of small amplitude oscillons in 1+1 dimensions. Using the amplitude of the oscillons and the inverse relativistic boost factor $\gamma^{-1}$ as the perturbation variables, we analytically calculate the leading order spatial and temporal phase shifts, and the change in the amplitude of the oscillons after the collisions. At leading order, we find that only the temporal phase shift receives a nonzero contribution, and that the collision is elastic. This work is also the first application of the general kinematic framework for understanding ultra-relativistic collisions (arXiv:1308.0606) to intrinsically time-dependent solitons.

The origin of hot Jupiters -- gas giant exoplanets orbiting very close to their host stars -- is a long-standing puzzle. Planet formation theories suggest that such planets are unlikely to have formed in-situ but instead may have formed at large orbital separations beyond the snow line and migrated inward to their present orbits. Two competing hypotheses suggest that the planets migrated either through interaction with the protoplanetary disk during their formation, or by disk-free mechanisms such as gravitational interactions with a third body. Observations of eccentricities and spin-orbit misalignments of hot Jupiter systems have been unable to differentiate between the two hypotheses. In the present work, we suggest that chemical depletions in hot Jupiter atmospheres might be able to constrain their migration mechanisms. We find that sub-solar carbon and oxygen abundances in Jovian-mass hot Jupiters around Sun-like stars are hard to explain by disk migration. Instead, such abundances are more readily explained by giant planets forming at large orbital separations, either by core accretion or gravitational instability, and migrating to close-in orbits via disk-free mechanisms involving dynamical encounters. Such planets also contain solar or super-solar C/O ratios. On the contrary, hot Jupiters with super-solar O and C abundances can be explained by a variety of formation-migration pathways which, however, lead to solar or sub-solar C/O ratios. Current estimates of low oxygen abundances in hot Jupiter atmospheres may be indicative of disk-free migration mechanisms. We discuss open questions in this area which future studies will need to investigate.

Global, symmetry-breaking phase transitions in the early universe can generate scaling seed networks which lead to metric perturbations. The acoustic waves in the photon-baryon plasma sourced by these metric perturbations, when Silk damped, generate spectral distortions of the cosmic microwave background (CMB). In this work, the chemical potential distortion ($\mu$) due to scaling seed networks is computed and the accompanying Compton $y$-type distortion is estimated. The specific model of choice is the $O(N)$ nonlinear $\sigma$-model for $N\gg 1$, but the results remain the same order of magnitude for other scaling seeds. If CMB anisotropy constraints to the $O(N)$ model are saturated, the resulting chemical potential distortion $\mu \lesssim 2\times 10^{-9}$.

Using data from the Fermi Large Area Telescope (LAT), we report the first clear gamma-ray measurement of a delay between flares from the gravitationally lensed images of a blazar. The delay was detected in B0218+357, a known double-image lensed system, during a period of enhanced gamma-ray activity with peak fluxes consistently observed to reach >20-50 times its previous average flux. An auto-correlation function analysis identified a delay in the gamma-ray data of 11.46 +/- 0.16 days (1 sigma) that is ~1 day greater than previous radio measurements. Considering that it is beyond the capabilities of the LAT to spatially resolve the two images, we nevertheless decomposed individual sequences of superposing gamma-ray flares/delayed emissions. In three such ~8-10 day-long sequences within a ~4-month span, considering confusion due to overlapping flaring emission and flux measurement uncertainties, we found flux ratios consistent with ~1, thus systematically smaller than those from radio observations. During the first, best-defined flare, the delayed emission was detailed with a Fermi pointing, and we observed flux doubling timescales of ~3-6 hrs implying as well extremely compact gamma-ray emitting regions.

We construct a perturbative framework for understanding the collision of solitons (more precisely, solitary waves) in relativistic scalar field theories. Our perturbative framework is based on the suppression of the space-time interaction area proportional to $1/(\gamma v)$, where $v$ is the relative velocity of an incoming solitary wave and $\gamma = 1/\sqrt{1-v^2} \gg 1$. We calculate the leading order results for collisions of (1+1) dimensional kinks in periodic potentials, and provide explicit, closed form expressions for the phase shift and the velocity change after the collisions. We find excellent agreement between our results and detailed numerical simulations. Crucially, our perturbation series is controlled by a kinematic parameter, and hence not restricted to small deviations around integrable cases such as the Sine-Gordon model.

We derive a closed-form expression for the phase shift experienced by 1+1 dimensional kinks colliding at ultra-relativistic velocities (gamma v >> 1), valid for arbitrary periodic potentials. Our closed-form expression is the leading order result of a more general scattering theory of solitary waves described in a companion paper [Amin, Lim and Yang, arXiv:1308.0606]. This theory relies on a small kinematic parameter 1/(gamma v) << 1 rather than a small parameter in the Lagrangian. Our analytic results can be directly extracted from the Lagrangian without solving the equation of motion. Based on our closed-form expression, we prove that kink-kink and kink-antikink collisions have identical phase shifts at leading order.

Oscillons are long-lived, localized, oscillatory scalar field configurations. In this work we derive a condition for the existence of small-amplitude oscillons (and provide solutions) in scalar field theories with non-canonical kinetic terms. While oscillons have been studied extensively in the canonical case, this is the first example of oscillons in scalar field theories with non-canonical kinetic terms. In particular, we demonstrate the existence of oscillons supported solely by the non-canonical kinetic terms, without any need for nonlinear terms in the potential. In the small-amplitude limit, we provide an explicit condition for their stability in d+1 dimensions against long-wavelength perturbations. We show that for d > 2, there exists a long-wavelength instability which can lead to radial collapse of small-amplitude oscillons.

Nonlinear propagation of electrostatic modes in ultrarelativistic dense elelectron-positron gravito-plasma at the polar cap region of pulsar magnetosphere is considered. A nonlinear Schr\"{o}dinger equation is obtained from the reductive perturbation method which predicts the existence of Langmuir dark solitons. Relevance of the propagating dark solitons to the pulsar radio emission is discussed.

We investigate the observational consequences of the quintessence field rolling to and oscillating near a minimum in its potential, "if" it happens close to the present epoch (z<0.2). We show that in a class of models, the oscillations lead to a rapid growth of the field fluctuations and the gravitational potential on subhorizon scales. The growth in the gravitational potential occurs on timescales << H^(1). This effect is present even when the quintessence parameters are chosen to reproduce an expansion history consistent with observations. For linearized fluctuations, we find that although the gravitational potential power spectrum is enhanced in a scale-dependent manner, the shape of the dark matter/galaxy power spectrum is not significantly affected. We find that the best constraints on such a transition in the quintessence field is provided via the integrated Sachs-Wolfe (ISW) effect in the CMB temperature power spectrum. Going beyond the linearized regime, the quintessence field can fragment into large, localized, long lived excitations (oscillons) with sizes comparable to galaxy clusters; this fragmentation could provide additional observational constraints. Two quoted "signatures" of modified gravity are a scale-dependent growth of the gravitational potential and a difference between the matter power spectrum inferred from measurements of lensing and galaxy clustering. Here, both effects are achieved by a minimally coupled scalar field in general relativity with a canonical kinetic term.

Oscillons are massive, long-lived, localized excitations of a scalar field. We show that in a large class of well-motivated single-field models, inflation is followed by self-resonance, leading to copious oscillon generation and a lengthy period of oscillon domination. These models are characterized by an inflaton potential which has a quadratic minimum and is shallower than quadratic away from the minimum. This set includes both string monodromy models and a class of supergravity inspired scenarios, and is in good agreement with the current central values of the concordance cosmology parameters. We assume that the inflaton is weakly coupled to other fields, so as not to quickly drain energy from the oscillons or prevent them from forming. An oscillon-dominated universe has a greatly enhanced primordial power spectrum on very small scales relative to that seen with a quadratic potential, possibly leading to novel gravitational effects in the early universe.

It is shown that a strongly magnetized isothermal pair plasma near the surface of a pulsar supports low-frequency (in comparison to electron cyclotron frequency) toroidal electrostatic plasma modes in the equatorial region. Physically, the thermal pressure coupled with the magnetic pressure creates the low frequency oscillations which may grow for particular case of inhomogeneities of the equilibrium magnetic field and the pair plasma density.

Analytical arguments suggest that a large class of scalar field potentials permit the existence of oscillons -- pseudo-stable, non-topological solitons -- in three spatial dimensions. In this paper we numerically explore oscillon solutions in three dimensions. We confirm the existence of these field configurations as solutions to the Klein-Gorden equation in an expanding background, and verify the predictions of Amin and Shirokoff for the characteristics of individual oscillons for their model. Further, we demonstrate that significant numbers of oscillons can be generated via fragmentation of the inflaton condensate, consistent with the analysis of Amin. These emergent oscillons can easily dominate the post-inflationary universe. Finally, both analytic and numerical results suggest that oscillons are stable on timescales longer than the post-inflationary Hubble time. Consequently, the post-inflationary universe can contain an effective matter-dominated phase, during which it is dominated by localized concentrations of scalar field matter.

We investigate the emergence of large, localized, pseudo-stable configurations (oscillons) from inflaton fragmentation at the end of inflation. We predict the number density of large oscillons, and the conditions necessary for their emergence in a class of inflationary models. Analytic estimates are provided for a 3+1 and 1+1-dimensional universe. We test our predictions with detailed numerical simulations in 1+1-dimensions. We see a zoo of oscillons emerging from the simulations, including the usual small amplitude "sech" oscillons as well as large "flat-topped" oscillons. The emergent oscillons account for approximately 80 per cent of the energy density of the inflaton.

Oscillons are extremely long lived, oscillatory, spatially localized field configurations that arise from generic initial conditions in a large number of non-linear field theories. With an eye towards their cosmological implications, we investigate their properties in an expanding universe. We (1) provide an analytic solution for one dimensional oscillons (for the models under consideration) and discuss their generalization to 3 dimensions, (2) discuss their stability against long wavelength perturbations and (3) estimate the effects of expansion on their shapes and life-times. In particular, we discuss a new, extended class of oscillons with surprisingly flat tops. We show that these flat topped oscillons are more robust against collapse instabilities in (3+1) dimensions than their usual counterparts. Unlike the solutions found in the small amplitude analysis, the width of these configurations is a non-monotonic function of their amplitudes.

The relationship between the metric and nonrelativistic matter distribution depends on the theory of gravity and additional fields, providing a possible way of distinguishing competing theories. With the assumption that the geometry and kinematics of the homogeneous universe have been measured to sufficient accuracy, we present a procedure for understanding and testing the relationship between the cosmological matter distribution and metric perturbations (along with their respective evolution) using the ratio of the physical size of the perturbation to the size of the horizon as our small expansion parameter. We expand around Newtonian gravity on linear, subhorizon scales with coefficient functions in front of the expansion parameter. Our framework relies on an ansatz which ensures that (i) the Poisson equation is recovered on small scales (ii) the metric variables (and any additional fields) are generated and supported by the nonrelativistic matter overdensity. The scales for which our framework is intended are small enough so that cosmic variance does not significantly limit the accuracy of the measurements and large enough to avoid complications from nonlinear effects and baryon cooling. The coefficient functions provide a general framework for contrasting the consequences of Lambda CDM and its alternatives. We calculate the coefficient functions for general relativity with a cosmological constant and dark matter, GR with dark matter and quintessence, scalar-tensor theories, f(R) gravity and braneworld models. We identify a possibly unique signature of braneworld models. Constraining the coefficient functions provides a streamlined approach for testing gravity in a scale dependent manner. We briefly discuss the observations best suited for an application of our framework.

In a large region of the supersymmetry parameter space, the annihilation cross section for neutralino dark matter is strongly dependent on the relative velocity of the incoming particles. We explore the consequences of this velocity dependence in the context of indirect detection of dark matter from the galactic center. We find that the increase in the annihilation cross section at high velocities leads to a flattening of the halo density profile near the galactic center and an enhancement of the annihilation signal.

We present and employ a new kinematical approach to cosmological `dark energy' studies. We construct models in terms of the dimensionless second and third derivatives of the scale factor a(t) with respect to cosmic time t, namely the present-day value of the deceleration parameter q_0 and the cosmic jerk parameter, j(t). An elegant feature of this parameterization is that all LCDM models have j(t)=1 (constant), which facilitates simple tests for departures from the LCDM paradigm. Applying our model to the three best available sets of redshift-independent distance measurements, from type Ia supernovae and X-ray cluster gas mass fraction measurements, we obtain clear statistical evidence for a late time transition from a decelerating to an accelerating phase. For a flat model with constant jerk, j(t)=j, we measure q_0=-0.81+-0.14 and j=2.16+0.81-0.75, results that are consistent with LCDM at about the 1sigma confidence level. A standard `dynamical' analysis of the same data, employing the Friedmann equations and modeling the dark energy as a fluid with an equation of state parameter, w (constant), gives Omega_m=0.306+0.042-0.040 and w=-1.15+0.14-0.18, also consistent with LCDM at about the 1sigma level. In comparison to dynamical analyses, the kinematical approach uses a different model set and employs a minimum of prior information, being independent of any particular gravity theory. The results obtained with this new approach therefore provide important additional information and we argue that both kinematical and dynamical techniques should be employed in future dark energy studies, where possible. Our results provide further interesting support for the concordance LCDM paradigm.

We present a set of new characteristic frequencies associated with accretion disks around compact objects. These frequencies arise from persistent rotating patterns in the disk that are finite in radial extent and driven purely by the gravity of the central body. Their existence depends on general relativistic corrections to orbital motion and, if observed, could be used to probe the strong gravity region around a black hole. We also discuss a possible connection to the puzzle of quasi-periodic oscillations.

Our fundamental lack of understanding of the acceleration of the Universe suggests that we consider a kinematic description. The simplest formalism involves the third derivative of the scale factor through a jerk parameter. A new approach is presented for describing the results of astronomical observations in terms of the contemporary jerk parameter and this is related to the equation of state approach. Simple perturbative expansions about $\Lambda$CDM are given.