We study the robustness of single-field inflation against inhomogeneities. We derive a simple analytic criterion on the shape of the potential for successful inflation in the presence of inhomogeneities, and demonstrate its validity using full 3+1 dimensional numerical relativity simulations on several classes of popular models of single-field inflation. We find that models with convex potentials are more robust to inhomogeneities than those with concave potentials, and that concave potentials that vary on super-Planckian scales are significantly more robust than those that vary on sub-Planckian scales.

Setting aside anthropic arguments, there is no reason why the universe should initially favour a net expanding phase rather than one experiencing a net contraction. However, a collapsing universe containing "normal" matter will end at a singularity in a finite time. We point out that there is a mechanism, derived from non-perturbative effects in Quantum Field Theory in a finite volume, which may provide a bias towards expansion when the spacetime volume shrinks, by dynamically violating the null energy condition, without the need for modified gravity or exotic matter. We describe a scalar field component subjected to this non-perturbative effect in a cosmological background and consider its impact on a contracting phase. We discuss how this could dynamically generate the necessary initial conditions for inflation to get started, or form part of the mechanism for a non-singular cosmological bounce.

Through numerical simulations of a minimally coupled massive Klein-Gordon scalar field, we show that it is possible to grow hair on a Schwarzschild black hole if one assumes an initial periodically time-varying but spatially homogeneous scalar background. By "hair", we mean a non-trivial profile in the scalar field. We find that this profile emerges on a timescale related to the mass of the black hole, with features related to the mass of the scalar particle. We undertake simulations with and without backreaction on the metric and see that the essential, qualitative features remain consistent. We also contrast the results from higher mass scalars to the case of a low mass with a large Compton wavelength. The results are particularly relevant for scalar-tensor theories of gravity and dark matter models consisting of a massive scalar, e.g. axions.

The QCD axion is a hypothetical particle motivated by the Strong CP problem of particle physics. One of the primary ways in which its existence can be inferred is via its function as an additional cooling channel in stars, with some of the strongest constraints coming from the supernova observation SN1987A. Multimessenger observations of binary neutron star mergers (such as those of GW170817, AT2017gfo, and GRB170817A) may provide another scenario in which such constraints could be obtained. In particular, the axion could potentially alter the lifetime, the ejection of material, and the emitted gravitational wave signal of the postmerger remnant. In this article, we perform numerical relativity simulations of a binary neutron star merger, including a phenomenological description of the nucleon-nucleon-axion bremsstrahlung to quantify the effects of such a cooling channel on the dynamical evolution. While our simulations show a difference in the temperature profile of the merger remnant, the imprint of the axion via nucleon-nucleon-axion bremsstrahlung on the emitted gravitational wave signal and the ejecta mass is too small to improve constraints on the axion mass with current or future planned detectors. Whilst we consider a limited number of cases, and a simplified cooling model, these broadly represent the "best case" scenario, thus, a more thorough investigation is unlikely to change the conclusions, at least for this particular interaction channel.

Axions are a potential dark matter candidate, which may condense and form self gravitating compact objects, called axion stars (ASs). In this work, we study for the first time head-on collisions of relativistic ASs with black holes (BHs) and neutron stars (NSs). In the case of BH-AS mergers we find that, in general, the largest scalar clouds are produced by mergers of low compactness ASs and spinning BHs. Although in most of the cases which we study the majority of the mass is absorbed by the BH within a short time after the merger, in favourable cases the remaining cloud surrounding the final BH remnant can be as large as 30% of the initial axion star mass, with a bosonic cloud mass of order ~ 0.1 M_BH and peak energy density comparable to that obtained in a superradiant build up. This provides a dynamical mechanism for the formation of long lived scalar hair, which could lead to observable signals in cases where the axion interacts with baryonic matter around the BH, or where it forms the seed of a future superradiant build up in highly spinning cases. Considering NS-AS collisions we find two possible final states (i) a BH surrounded by a (small) scalar cloud, or (ii) a stable NS enveloped in an axion cloud of roughly the same mass as the initial AS. Whilst for low mass ASs the NS is only mildly perturbed by the collision, a larger mass AS gives rise to a massive ejection of baryonic mass from the system, purely due to gravitational effects. Therefore, even in the absence of a direct axion coupling to baryonic matter, NS-AS collisions could give rise to electromagnetic observables in addition to their gravitational wave signatures.

Axions are increasingly favoured as a candidate particle for the dark matter in galaxies, since they satisfy the observational requirements for cold dark matter and are theoretically well motivated. Fluctuations in the axion field give rise to stable localised overdensities known as axion stars, which, for the most massive, compact cases, are potential neutron star mimickers. In principle, there are no fundamental arguments against the multi-messenger observations of GW170817/GRB170817A/AT2017gfo arising from the merger of a neutron star with a neutron star mimicker, rather than from a binary neutron star. To constrain this possibility and better understand the astrophysical signatures of a neutron star--axion star (NSAS) merger, we present in this work a detailed example case of a NSAS merger based on full 3D numerical relativity simulations, and give an overview of the many potential observables - ranging from gravitational waves, to optical and near-infrared electromagnetic signals, radio flares, fast radio bursts, gamma ray bursts, and neutrino emission. We discuss the individual channels and estimate to which distances current and future observatories might be able to detect such a NSAS merger. Such signals could constrain the unknown axion mass and its couplings to standard baryonic matter,thus enhancing our understanding of the dark matter sector of the Universe.

With the first direct detections of gravitational waves (GWs) from the coalescence of compact binaries observed by the advanced LIGO and VIRGO interferometers, the era of GW astronomy has begun. Whilst there is strong evidence that the observed GWs are connected to the merger of two black holes (BH) or two neutron stars (NS), future detections may present a less consistent picture. Indeed, the possibility that the observed GW signal was created by a merger of exotic compact objects (ECOs) such as boson stars (BS) or axion stars (AS) has not yet been fully excluded. For a detailed understanding of the late stages of the coalescence full 3D numerical relativity simulations are essential. In this paper, we extend the infrastructure of the numerical relativity code BAM, to permit the simultaneous simulation of baryonic matter with bosonic scalar fields, thus enabling the study of BS-BS, BS-NS, and BS-BH mergers. We present a large number of single star evolutions to test the newly implemented routines, and to quantify the numerical challenges of such simulations, which we find to partially differ from the default NS case. We also compare head-on BS-BS simulations with independent numerical relativity codes, namely the SpEC and the GRChombo codes, and find good general agreement. Finally, we present what are, to the best of our knowledge, the first full NR simulations of BS-NS mergers, a first step towards identifying the hallmarks of BS-NS interactions in the strong gravity regime, as well as possible GW and electromagnetic observables.

Extending our previous work on the robustness of inflation to perturbations in the scalar field, we investigate the effects of perturbations in the transverse traceless part of the extrinsic curvature on the evolution of an inhomogeneous inflaton field. Focusing on small field models, we show that these additional metric inhomogeneities initially reduce the total number of e-folds as the amplitude increases, but that the reduction saturates and even reverses above a certain amplitude. We present an argument that this is due to the presence of a large initial Hubble friction when metric perturbations are large.

We consider the effects of inhomogeneous initial conditions in both the scalar field profile and the extrinsic curvature on different inflationary models. In particular, we compare the robustness of small field inflation to that of large field inflation, using numerical simulations with Einstein gravity in 3+1 dimensions. We find that small field inflation can fail in the presence of subdominant gradient energies, suggesting that it is much less robust to inhomogeneities than large field inflation, which withstands dominant gradient energies. However, we also show that small field inflation can be successful even if some regions of spacetime start out in the region of the potential that does not support inflation. In the large field case, we confirm previous results that inflation is robust if the inflaton occupies the inflationary part of the potential. Furthermore, we show that increasing initial scalar gradients will not form sufficiently massive inflation-ending black holes if the initial hypersurface is approximately flat. Finally, we consider the large field case with a varying extrinsic curvature K, such that some regions are initially collapsing. We find that this may again lead to local black holes, but overall the spacetime remains inflationary if the spacetime is open, which confirms previous theoretical studies.

The classical equations of motion for an axion with potential $V(\phi)=m_a^2f_a^2 [1-\cos (\phi/f_a)]$ possess quasi-stable, localized, oscillating solutions, which we refer to as "axion stars". We study, for the first time, collapse of axion stars numerically using the full non-linear Einstein equations of general relativity and the full non-perturbative cosine potential. We map regions on an "axion star stability diagram", parameterized by the initial ADM mass, $M_{\rm ADM}$, and axion decay constant, $f_a$. We identify three regions of the parameter space: i) long-lived oscillating axion star solutions, with a base frequency, $m_a$, modulated by self-interactions, ii) collapse to a BH and iii) complete dispersal due to gravitational cooling and interactions. We locate the boundaries of these three regions and an approximate "triple point" $(M_{\rm TP},f_{\rm TP})\sim (2.4 M_{pl}^2/m_a,0.3 M_{pl})$. For $f_a$ below the triple point BH formation proceeds during winding (in the complex $U(1)$ picture) of the axion field near the dispersal phase. This could prevent astrophysical BH formation from axion stars with $f_a\ll M_{pl}$. For larger $f_a\gtrsim f_{\rm TP}$, BH formation occurs through the stable branch and we estimate the mass ratio of the BH to the stable state at the phase boundary to be $\mathcal{O}(1)$ within numerical uncertainty. We discuss the observational relevance of our findings for axion stars as BH seeds, which are supermassive in the case of ultralight axions. For the QCD axion, the typical BH mass formed from axion star collapse is $M_{\rm BH}\sim 3.4 (f_a/0.6 M_{pl})^{1.2} M_\odot$.