Search

Massive scalar fields provide excellent dark matter candidates, whose dynamics are often explored analytically and numerically using nonrelativistic Schr\"{o}dinger-Poisson (SP) equations in a cosmological context. In this paper, starting from the nonlinear and fully relativistic Klein-Gordon-Einstein (KGE) equations in an expanding universe, we provide a systematic framework for deriving the SP equations, as well as relativistic corrections to them, by integrating out `fast modes' and including nonlinear metric and matter contributions. We provide explicit equations for the leading-order relativistic corrections, which provide insight into deviations from the SP equations as the system approaches the relativistic regime. Upon including the leading-order corrections, our equations are applicable beyond the domain of validity of the SP system, and are simpler to use than the full KGE case in some contexts. As a concrete application, we calculate the mass-radius relationship of solitons in scalar dark matter and accurately capture the deviations of this relationship from the SP system towards the KGE one., Comment: 20 pages, 5 figures, plus appendices

We investigate the production of photons from coherently oscillating, spatially localized clumps of axionic fields (oscillons and axion stars) in the presence of external electromagnetic fields. We delineate different qualitative behaviour of the photon luminosity in terms of an effective dimensionless coupling parameter constructed out of the axion-photon coupling, and field amplitude, oscillation frequency and radius of the axion star. For small values of this dimensionless coupling, we provide a general analytic formula for the dipole radiation field and the photon luminosity per solid angle, including a strong dependence on the radius of the configuration. For moderate to large coupling, we report on a non-monotonic behavior of the luminosity with the coupling strength in the presence of external magnetic fields. After an initial rise in luminosity with the coupling strength, we see a suppression (by an order of magnitude or more compared to the dipole radiation approximation) at moderately large coupling. At sufficiently large coupling, we find a transition to a regime of exponential growth of the luminosity due to parametric resonance. We carry out 3+1 dimensional lattice simulations of axion electrodynamics, at small and large coupling, including non-perturbative effects of parametric resonance as well as backreaction effects when necessary. We also discuss medium (plasma) effects that lead to resonant axion to photon conversion, relevance of the coherence of the soliton, and implications of our results in astrophysical and cosmological settings., Comment: 12 figures, 28 pages + 1 Appendix

The polarization of Cosmic Microwave Background (CMB) photons is rotated as they pass through (ultralight-) axion string loops. Studying this birefringence can reveal valuable information about the axion-photon coupling and the structure of the string network. We develop an approximate analytic formalism and identify a kernel function that can be used to calculate the two-point correlation function for CMB birefringence induced by an arbitrary axion string network. Using this formalism, we evaluate the birefringence signal for some simple loop distributions (including scaling and network collapse). We find that the angular correlation function has a characteristic angular scale set by $\theta_\mathrm{min}$, which corresponds to the angular extent of the loops at the time of recombination. This results in a peak in the birefringence power spectrum around $\ell_p \sim 1/\theta_\mathrm{min}$. An additional scale, controlled by the axion's mass, is introduced if the network collapses before today., Comment: 31 pages + appendices; 9 figures. v2: updated towards the version published in JCAP (added limits from ACTPol n SPTPol in power spectrum figures; added more references)

We investigate the bursts of electromagnetic and scalar radiation resulting from the collision, and merger of oscillons made from axion-like particles using 3+1 dimensional lattice simulations of the coupled axion-gauge field system. The radiation into photons is suppressed before the merger. However, it becomes the dominant source of energy loss after the merger if a resonance condition is satisfied. Conversely, the radiation in scalar waves is dominant during initial merger phase but suppressed after the merger. The backreaction of scalar and electromagnetic radiation is included in our simulations. We evolve the system long enough to see that the resonant photon production extracts a significant fraction of the initial axion energy, and again falls out of the resonance condition. We provide a parametric understanding of the time, and energy scales involved in the process and discuss observational prospects of detecting the electromagnetic signal., Comment: 21 pages, 6 figures + 4 pages of Appendix (v2, updated references and minor changes in text)

We study the problem of constructing a general hybrid quantum-classical bracket from a partial classical limit of a full quantum bracket. Introducing a hybrid composition product, we show that such a bracket is the commutator of that product. From this we see that the hybrid bracket will obey the Jacobi identity and the Leibniz rule provided the composition product is associative. This suggests that the set of hybrid variables belonging to an associative subalgebra with the composition product will have consistent quantum-classical dynamics. This restricts the class of allowed quantum-classical interaction Hamiltonians. Furthermore, we show that pure quantum or classical variables can interact in a consistent framework, unaffected by no-go theorems in the literature or the restrictions for hybrid variables. In the proposed scheme, quantum backreaction appears as quantum-dependent terms in the classical equations of motion., Comment: Published in Physical Review A, minor revisions to v3, 11 pages

CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semi-analytic projection tool, targeted explicitly towards optimizing constraints on the tensor-to-scalar ratio, $r$, in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2--3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments given a desired scientific goal. To form a closed-loop process, we couple this semi-analytic tool with map-based validation studies, which allow for the injection of additional complexity and verification of our forecasts with several independent analysis methods. We document multiple rounds of forecasts for CMB-S4 using this process and the resulting establishment of the current reference design of the primordial gravitational-wave component of the Stage-4 experiment, optimized to achieve our science goals of detecting primordial gravitational waves for $r > 0.003$ at greater than $5\sigma$, or, in the absence of a detection, of reaching an upper limit of $r < 0.001$ at $95\%$ CL., Comment: 24 pages, 8 figures, 9 tables, submitted to ApJ. arXiv admin note: text overlap with arXiv:1907.04473

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., Comment: 67 pages, 18 figures. v2: Discussion and references added. Accepted for publication in The Open Journal of Astrophysics

Oscillons are extremely long-lived, spatially-localized field configurations in real-valued scalar field theories that slowly lose energy via radiation of scalar waves. Before their eventual demise, oscillons can pass through (one or more) exceptionally stable field configurations where their decay rate is highly suppressed. We provide an improved calculation of the non-trivial behavior of the decay rates, and lifetimes of oscillons. In particular, our calculation correctly captures the existence (or absence) of the exceptionally long-lived states for large amplitude oscillons in a broad class of potentials, including non-polynomial potentials that flatten at large field values. The key underlying reason for the improved (by many orders of magnitude in some cases) calculation is the systematic inclusion of a spacetime-dependent effective mass term in the equation describing the radiation emitted by oscillons (in addition to a source term). Our results for the exceptionally stable configurations, decay rates, and lifetime of large amplitude oscillons (in some cases $\gtrsim 10^8$ oscillations) in such flattened potentials might be relevant for cosmological applications., Comment: v2: 27 pages, 7 figures. References added, added short discussion on cosine potentials

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., Comment: 68 pages, 23 figures

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., Comment: 35 pages, 13 figures; minor changes, accepted for publication in JCAP

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., Comment: 19 pages, 12 figures, submitted to MNRAS

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., Comment: 6 pages, 2 figures, accepted to Phys. Rev. Lett

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., Comment: 26 pages, 22 figures, comments welcome

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., Comment: Project White Paper submitted to the 2020 Decadal Survey, 10 pages plus references. arXiv admin note: substantial text overlap with arXiv:1907.04473

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., Comment: 10 pages + references, 2 figures, 1 table; APC white paper submitted to the Astro2020 decadal survey. arXiv admin note: text overlap with arXiv:1810.09572

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., Comment: 17 pages, 2 figures. Comments welcome. v2: Added references and minor corrections

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., Comment: 5 pages, 1 figure, Astro2020 Decadal Survey science white paper

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., Comment: 5 pages + references; 1 figure; science white paper submitted to the Astro2020 decadal survey

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., Comment: 5 pages + references; Submitted to the Astro2020 call for science white papers

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., Comment: White Paper submitted on March 8, 2019 to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics)

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., Comment: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

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., Comment: 80 pages, 44 figures

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., Comment: v3: 17 pages, 8 Figures. References added, extended discussion in the appendices, minor changes in the text. Matches version accepted for publication in PRD

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., Comment: 18 pages, 8 figures

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$)., Comment: 12 pages, 4 figures

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., Comment: v2: 6 pages, 2 figures; matches PRD version

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., Comment: 8 Pages, 8 figures, movies : https://www.youtube.com/playlist?list=PLSkfizpQDrcahgvc5TvBk5qtXAzkSyHPe

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., Comment: v2: 8 pages + 14 pages of appendices, 10 figures. Clarified and expanded discussion, more numerical simulations and evidence, references added. Accepted by Phys Rev D

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$., Comment: 25 pages, 16 figures

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$., Comment: 41 pages, 7 figures, 4 appendices. Corrected typos, added references

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., Comment: 11 pages, 3 figures, accepted by MNRAS

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., Comment: v3: 6 pages, 5 figures; clarifying the comment on asymmetric potentials

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., Comment: 48 pages, 6 figures, 2 appendices, v3: references added, minor changes to text, to appear in JCAP

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., Comment: v3: minor changes and references added; matches published version in JCAP

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., Comment: v1: 44 pages, 9 figures. Invited review for the International Journal of Modern Physics D. For v2: Minor changes in text; 19 references added. For v3: Minor changes in text; 16 references added; similar to version to be published by IJMPD

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., Comment: 12 pages, 3 figures, version 2, added one reference and matching the version to appear on PRD

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., Comment: Accepted for publication in ApJ Letters

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., Comment: 19 pages, 9 figures. v3: references added and minor changes in text. Published in Phys. Rev. D