We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings–Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of 1014, and this same model is favored over an uncorrelated common power-law spectrum model with Bayes factors of 200–1000, depending on spectral modeling choices. We have built a statistical background distribution for the latter Bayes factors using a method that removes interpulsar correlations from our data set, finding p = 10−3 (≈3σ) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of interpulsar correlations yields p = 5 × 10−5 to 1.9 × 10−4 (≈3.5σ–4σ). Assuming a fiducial f −2/3 characteristic strain spectrum, as appropriate for an ensemble of binary supermassive black hole inspirals, the strain amplitude is 2.4 − 0.6 + 0.7 × 10 − 15 (median + 90% credible interval) at a reference frequency of 1 yr−1. The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings–Downs correlations points to the gravitational-wave origin of this signal.

In theories with a broken discrete symmetry, Hubble sized spherical domain walls may spontaneously nucleate during inflation. These objects are subsequently stretched by the inflationary expansion, resulting in a broad distribution of sizes. The fate of the walls after inflation depends on their radius. Walls smaller than a critical radius fall within the cosmological horizon early on and collapse due to their own tension, forming ordinary black holes. But if a wall is large enough, its repulsive gravitational field becomes dominant much before the wall can fall within the cosmological horizon. In this "supercritical" case, a wormhole throat develops, connecting the ambient exterior FRW universe with an interior baby universe, where the exponential growth of the wall radius takes place. The wormhole pinches off in a time-scale comparable to its light-crossing time, and black holes are formed at its two mouths. As discussed in previous work, the resulting black hole population has a wide distribution of masses and can have significant astrophysical effects. The mechanism of black hole formation has been previously studied for a dust-dominated universe. Here we investigate the case of a radiation-dominated universe, which is more relevant cosmologically, by using numerical simulations in order to find the initial mass of a black hole as a function of the wall size at the end of inflation. For large supercritical domain walls, this mass nearly saturates the upper bound according to which the black hole cannot be larger than the cosmological horizon. We also find that the subsequent accretion of radiation satisfies a scaling relation, resulting in a mass increase by about a factor of 2.

The 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons.

Vacuum bubbles may nucleate during the inflationary epoch and expand, reaching relativistic speeds. After inflation ends, the bubbles are quickly slowed down, transferring their momentum to a shock wave that propagates outwards in the radiation background. The ultimate fate of the bubble depends on its size. Bubbles smaller than certain critical size collapse to ordinary black holes, while in the supercritical case the bubble interior inflates, forming a baby universe, which is connected to the exterior region by a wormhole. The wormhole then closes up, turning into two black holes at its two mouths. We use numerical simulations to find the masses of black holes formed in this scenario, both in subcritical and supercritical regime. The resulting mass spectrum is extremely broad, ranging over many orders of magnitude. For some parameter values, these black holes can serve as seeds for supermassive black holes and may account for LIGO observations.

Recently there has been much interest in light dark matter, especially ultra-light axions, as they may provide a solution to the core-cusp problem at the center of galaxies. Since very light bosons can have a de Broglie wavelength that is of astrophysical size, they can smooth out the centers of galaxies to produce a core, as opposed to vanilla dark matter models, and so it has been suggested that this solves the core-cusp problem. In this work, we critically examine this claim. While an ultra-light particle will indeed lead to a core, we examine whether the relationship between the density of the core and its radius matches the data over a range of galaxies. We first review data that shows the core density of a galaxy ρc varies as a function of the core radius Rc as ρc ∝ 1/R β c with β ≈ 1. We then compare this to theoretical models. We examine a large class of light scalar dark matter models, governed by some potential V . For simplicity, we take the scalar to be complex with a global U (1) symmetry in order to readily organize solutions by a conserved particle number. However, we expect our central conclusions to persist even for a real scalar, and furthermore, a complex scalar matches the behavior of a real scalar in the non-relativistic limit which is the standard regime of interest. For any potential V , we find the relationship between ρc and Rc for ground state solutions is always in one of the following regimes: (i) β 1, or (ii) β 1, or (iii) unstable, and so it never matches the data. We also find similar conclusions for virialized dark matter, more general scalar field theories, degenerate fermion dark matter, superfluid dark matter, and general polytropes. We conclude that the solution to the core-cusp problem is more likely due to either complicated baryonic effects or some other type of dark matter interactions.