We first show that the effective non-relativistic theory of gravitationally interacting, massive integer-spin fields (spin-$0$, $1$, and $2$ in particular) is described by a $2s+1$ component Schr{o}dinger-Poisson action, where $s$ is the spin of the
field. We then construct $s+1$ distinct, gravitationally supported solitons in this non-relativistic theory from identically polarized plane waves. Such solitons are extremally polarized, with macroscopically large spin, but no orbital angular momentum. These $s+1$ solitons form a basis set, out of which partially polarized solitons can be constructed. All such solitons are ground states, have a spherically symmetric energy density but not field configurations. We discuss how solitons in higher-spin fields can be distinguished from scalar solitons, and potential gravitational and non-gravitational probes of them.
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 nonl
inear 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.
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 pho
ton 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.
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 strin
g 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 axions mass, is introduced if the network collapses before today.
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 in
to 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.
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 r
adiation 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.
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) exceptiona
lly 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.
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.
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.
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.
