Search

Evidence has emerged for a stochastic signal correlated among 67 pulsars within the 15-year pulsar-timing data set compiled by the NANOGrav collaboration. Similar signals have been found in data from the European, Indian, Parkes, and Chinese PTAs. This signal has been interpreted as indicative of the presence of a nanohertz stochastic gravitational wave background. To explore the internal consistency of this result we investigate how the recovered signal strength changes as we remove the pulsars one by one from the data set. We calculate the signal strength using the (noise-marginalized) optimal statistic, a frequentist metric designed to measure correlated excess power in the residuals of the arrival times of the radio pulses. We identify several features emerging from this analysis that were initially unexpected. The significance of these features, however, can only be assessed by comparing the real data to synthetic data sets. After conducting identical analyses on simulated data sets, we do not find anything inconsistent with the presence of a stochastic gravitational wave background in the NANOGrav 15-year data. The methodologies developed here can offer additional tools for application to future, more sensitive data sets. While this analysis provides an internal consistency check of the NANOGrav results, it does not eliminate the necessity for additional investigations that could identify potential systematics or uncover unmodeled physical phenomena in the data., Comment: 21 pages, 11 figures, 2 tables

Pulsar timing array observations have found evidence for an isotropic gravitational wave background with the Hellings-Downs angular correlations, expected from general relativity. This interpretation hinges on the measured shape of the angular correlations, which is predominately quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlations into a sum of Legendre polynomials and use a Bayesian analysis to constrain their coefficients with the 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). When including Legendre polynomials with multipoles $\ell \geq 2$, we only find a significant signal in the quadrupole with an amplitude consistent with general relativity and non-zero at the $\sim 95\%$ confidence level and a Bayes factor of 200. When we include multipoles $\ell \leq 1$, the Bayes factor evidence for quadrupole correlations decreases by more than an order of magnitude due to evidence for a monopolar signal at approximately 4 nHz which has also been noted in previous analyses of the NANOGrav 15-year data. Further work needs to be done in order to better characterize the properties of this monopolar signal and its effect on the evidence for quadrupolar angular correlations., Comment: 15 pages, 6 figures

The detection of a stochastic gravitational wave background by pulsar timing arrays suggests the presence of a supermassive black hole binary population. Although the observed spectrum generally aligns with predictions from orbital evolution driven by gravitational wave emission in circular orbits, there is a discernible preference for a turnover at the lowest observed frequencies. This turnover could indicate a significant hardening phase, transitioning from early environmental influences to later stages predominantly influenced by gravitational wave emission. In the vicinity of these binaries, the ejection of stars or dark matter particles through gravitational three-body slingshots efficiently extracts orbital energy, leading to a low-frequency turnover in the spectrum. By analyzing the NANOGrav 15-year data, we assess how the gravitational wave spectrum depends on the initial inner galactic profile prior to disruption by binary ejections, accounting for a range of initial binary eccentricities. Our findings suggest a parsec-scale galactic center density around $10^6\,M_\odot/\textrm{pc}^3$ across most of the parameter space, offering insights into the environmental effects on black hole evolution and combined matter density near galaxy centers., Comment: 15 pages, 5 figures

The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index \gamma. In this paper, we consider the next logical step beyond this minimal spectral model, allowing for a running (i.e., logarithmic frequency dependence) of the spectral index, \gamma_run(f) = \gamma + \beta \ln(f/f_ref). We fit this running-power-law (RPL) model to the NANOGrav 15-year data and perform a Bayesian model comparison with the minimal constant-power-law (CPL) model, which results in a 95% credible interval for the parameter \beta consistent with no running, \beta \in [-0.80,2.96], and an inconclusive Bayes factor, B(RPL vs. CPL) = 0.69 +- 0.01. We thus conclude that, at present, the minimal CPL model still suffices to adequately describe the NANOGrav signal; however, future data sets may well lead to a measurement of nonzero \beta. Finally, we interpret the RPL model as a description of primordial GWs generated during cosmic inflation, which allows us to combine our results with upper limits from big-bang nucleosynthesis, the cosmic microwave background, and LIGO-Virgo-KAGRA., Comment: 17 pages, 4 figures, 2 tables

Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Those analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residuals across pulsars. These assumptions are idealizations that may not be realized in actuality. We test them in the NANOGrav 15 yr data set using Bayesian posterior predictive checks: after fitting our fiducial model to real data, we generate a population of simulated data-set replications, and use them to assess whether the optimal-statistic significance, inter-pulsar correlations, and spectral coefficients assume extreme values for the real data when compared to the replications. We confirm that the NANOGrav 15 yr data set is consistent with power-law and Hellings--Downs assumptions. We also evaluate the evidence for the stochastic background using posterior-predictive versions of the frequentist optimal statistic and of Bayesian model comparison, and find comparable significance (3.2\ $\sigma$ and 3\ $\sigma$ respectively) to what was previously reported for the standard statistics. We conclude with novel visualizations of the reconstructed gravitational waveforms that enter the residuals for each pulsar. Our analysis strengthens confidence in the identification and characterization of the gravitational-wave background as reported by NANOGrav., Comment: 20 pages, 14 Figures

The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 \sigma - 1.9 \sigma$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 \sigma - 2.1 \sigma$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26~\mathrm{nHz}$ break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe., Comment: 13 pages, 8 figures, 2 appendices, in press at ApJ

Transient gravitational waves (aka gravitational wave bursts) within the nanohertz frequency band could be generated by a variety of astrophysical phenomena such as the encounter of supermassive black holes, the kinks or cusps in cosmic strings, or other as-yet-unknown physical processes. Radio-pulses emitted from millisecond pulsars could be perturbed by passing gravitational waves, hence the correlation of the perturbations in a pulsar timing array can be used to detect and characterize burst signals with a duration of $\mathcal{O}(1\text{-}10)$ years. We propose a fully Bayesian framework for the analysis of the pulsar timing array data, where the burst waveform is generically modeled by piecewise straight lines, and the waveform parameters in the likelihood can be integrated out analytically. As a result, with merely three parameters (in addition to those describing the pulsars' intrinsic and background noise), one is able to efficiently search for the existence and the sky location of {a burst signal}. If a signal is present, the posterior of the waveform can be found without further Bayesian inference. We demonstrate this model by analyzing simulated data sets containing a stochastic gravitational wave background {and a burst signal generated by the parabolic encounter of two supermassive black holes., Comment: 13 pages, 10 figures

The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations., Comment: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

The 15-year 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., Comment: 74 pages, 31 figures, 4 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year 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 $10^{14}$, 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 these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3\sigma$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4\sigma$). 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.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). 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., Comment: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The population of black holes inferred from the detection of gravitational waves by the LIGO-Virgo-KAGRA Collaboration has revealed interesting features in the properties of black holes in the Universe. We analyze the GWTC-3 dataset assuming the detected black holes in each event had an either astrophysical or primordial origin. In particular, we consider astrophysical black holes described by the fiducial \textsc{Power Law + Peak} distribution and primordial black holes whose mass function obeys a broken power law. These primordial black holes can be generated by vacuum bubbles that nucleate during inflation. We find that astrophysical black holes dominate the events with mass less than $\sim 30M_\odot$, whereas primordial black holes are responsible for the massive end, and also for the peak at $\sim 30M_\odot$ in the mass distribution. More than half of the observed events could come from primordial black hole mergers. We also discuss the implications on the primordial black hole formation mechanism and the underlying inflationary model., Comment: The ABH model has been updated

We study the dynamics of a cosmic string loop captured by a rotating black hole, ignoring string reconnections. A loop is numerically evolved in Kerr spacetime, with the result that it turns into one or more growing or contracting double-lines rotating around the black hole in the equatorial plane. This is in good agreement with the approximate analytical treatment of the problem investigated by Xing et al., who studied the evolution of the auxiliary curve associated with the string loop. We confirm that the auxiliary curve deformation can indeed describe the string motion in realistic physical scenarios to a reasonable accuracy, and can thus be used to further study other phenomena such as superradiance and reconnections of the captured loop., Comment: 28 pages, 10 figures

The Peters formula, which tells how the coalescence time of a binary system emitting gravitational radiation is determined by the initial size and shape of the elliptic orbit, is often used in estimating the merger rate of primordial black holes and the gravitational wave background from the mergers. Valid as it is in some interesting scenarios, such as the analysis of the LIGO-Virgo events, the Peters formula fails to describe the coalescence time if the orbital period of the binary exceeds the value given by the formula. This could underestimate the event rate of mergers that occur before the cosmic time $t\sim 10^{13}\ \text{s}$. As a result, the energy density spectrum of the gravitational wave background could develop a peak, which is from mergers occurring at either $t\sim 10^{13}\ \text{s}$ (for black holes with mass $M\gtrsim 10^8 M_\odot$) or $t\sim 10^{26}(M/M_\odot)^{-5/3}\ \text{s}$ (for $10^5 M_\odot \lesssim M\lesssim 10^8 M_\odot$). This can be used to constrain the fraction of dark matter in primordial black holes (denoted by $f$) if potential probes (such as SKA and U-DECIGO) do not discover such a background, with the result $f\lesssim 10^{-6}\text{-}10^{-4}$ for the mass range $10\text{-} 10^9M_\odot$. We then consider the effect of mass accretion onto primordial black holes at redshift $z\sim 10$, and find that the merger rate could drop significantly at low redshifts. The spectrum of the gravitational wave background thus gets suppressed at the high-frequency end. This feature might be captured by future detectors such as ET and CE for initial mass $M= \mathcal{O}(10\text{-}100) M_\odot$ with $f\gtrsim 10^{-4}$., Comment: 22 pages and 4 figures in REVTeX 4

In a variety of mechanisms generating primordial black holes, each black hole is expected to form along with a surrounding underdense region that roughly compensates the black hole mass. This region will propagate outwards and expand as a shell at the speed of sound in the homogeneous background. Dissipation of the shell due to Silk damping could lead to detectable $\mu$-distortion in the CMB spectrum. While the current bound on the average $\mu$-distortion is $\left| \bar{\mu}\right|\lesssim10^{-4}$, the standard $\Lambda$CDM model predicts $\left| \bar{\mu}\right|\sim10^{-8}$, which could possibly be detected in future missions. It is shown in this work that the non-observation of $\bar{\mu}$ beyond $\Lambda$CDM can place a new upper bound on the density of supermassive primordial black holes within the mass range $10^{6}M_{\odot}\lesssim M\lesssim10^{15}M_{\odot}$. Furthermore, black holes with initial mass $M\gtrsim10^{12}M_{\odot}$ could leave a pointlike distortion with $\mu\gtrsim10^{-8}$ at an angular scale $\sim 1^{\circ}$ in CMB, and its non-observation would impose an even more stringent bound on the population of these stupendously large primordial black holes., Comment: 21 pages in REVTeX (V.4), 3 figures

The LIGO-Virgo Collaboration has so far detected around 90 black holes, some of which have masses larger than what were expected from the collapse of stars. The mass distribution of LIGO-Virgo black holes appears to have a peak at $\sim30M_{\odot}$ and two tails on the ends. By assuming that they all have a primordial origin, we analyze the GWTC-1 (O1\&O2) and GWTC-2 (O3a) datasets by performing maximum likelihood estimation on a broken power law mass function $f(m)$, with the result $f\propto m^{1.2}$ for $m<35M_{\odot}$ and $f\propto m^{-4}$ for $m>35M_{\odot}$. This appears to behave better than the popular log-normal mass function. Surprisingly, such a simple and unique distribution can be realized in our previously proposed mechanism of PBH formation, where the black holes are formed by vacuum bubbles that nucleate during inflation via quantum tunneling. Moreover, this mass distribution can also provide an explanation to supermassive black holes formed at high redshifts., Comment: 24 pages in revtex4, 5 figures

The discoveries of LIGO/Virgo black holes in recent years have revitalized the study of primordial black holes. In this work, we investigate a mechanism where primordial black holes are formed by vacuum bubbles that randomly nucleate during inflation through quantum tunneling. After inflation, these bubbles typically run into the ambient radiation fluid with a large Lorentz factor. In our previous work, we assumed the bubble fields are strongly coupled to the standard model particles so that the bubble wall is impermeable. Here we complete this picture by considering bubbles interacting with the fluid only through gravity. By studying the scenario in several limits, we found that black holes could form in either subcritical or supercritical regime. Depending on the model parameters, the resulting mass spectrum of the black holes could be wide or narrow, and may develop two peaks separated by a large mass range. With different spectra, these black holes may account for the LIGO/Virgo black holes, supermassive black holes, and may play an important role in dark matter., Comment: 28 pages, 5 figures

Primordial bubbles that possibly nucleate through quantum tunneling during inflation in a multidimensional potential might have left some relic detectable at the present time. These bubbles turn into black holes during the radiation era, which may account for the LIGO black holes, supermassive black holes, and may play an important role in dark matter. Typically, these black holes are surrounded by an energy deficit in the form of a spherical sound wave packet propagating outwards. In this work we study how this perturbation of the cosmic plasma dissipates before the time of recombination, leading to spectral distortions in CMB. We find that there may exist some rare regions on the last scattering surface containing huge black holes, which have produced potentially detectable point-like signals of $\mu$-type distortions., Comment: 21 pages, 3 figures

The much-discussed swampland conjectures suggest significant constraints on the properties of string theory landscape and on the nature of the multiverse that this landscape can support. The conjectures are especially constraining for models of inflation; in particular, they exclude the existence of de Sitter (dS) vacua. If the conjectures are false and dS vacua do exist, it still appears that their construction in string theory requires a fair amount of fine-tuning, so they may be vastly outnumbered by AdS vacua. Here we explore the multiverse structure suggested by these considerations. We consider two scenarios: (i) a landscape where dS vacua are rare and (ii) a landscape where dS vacua do not exist and the dS potential maxima and saddle points are not flat enough to allow for the usual hilltop inflation, even though slow-roll inflation is possible on the slopes of the potential. We argue that in both scenarios inflation is eternal and all parts of the landscape that can support inflation get represented in the multiverse. The spacetime structure of the multiverse in such models is nontrivial and is rather different from the standard picture., Comment: 36 pages, 11 figures

We use analytic estimates and numerical simulations to explore the stochastic approach to vacuum decay. According to this approach, the time derivative of a scalar field, which is in a local vacuum state, develops a large fluctuation and the field \flies over" a potential barrier to another vacuum. The probability distribution for the initial fluctuation is found quantum mechanically, while the subsequent nonlinear evolution is determined by classical dynamics. We find in a variety of cases that the rate of such flyover transitions has the same parametric form as that of tunneling transitions calculated using the instanton method, differing only by a numerical factor O(1) in the exponent. An important exception is an "upward" transition from a de Sitter vacuum to a higher-energy de Sitter vacuum state. The rate of flyover transitions in this case is parametrically different and can be many orders of magnitude higher than tunneling. This result is in conflict with the conventional picture of quantum de Sitter space as a thermal state. Our numerical simulations indicate that the dynamics of bubble nucleation in flyover transitions is rather different from the standard picture. The difference is especially strong for thin-wall bubbles in flat space, where the transition region oscillates between true and false vacuum until a true vacuum shell is formed which expands both inwards and outwards, and for upward de Sitter transitions, where the inflating new vacuum region is contained inside of a black hole., Comment: 28 pages, 12 figures

Vacuum bubbles may nucleate and expand during the cosmic inflation. When inflation ends, the bubbles run into the ambient plasma, producing strong shocks followed by underdensity waves, which propagate outwards. The bubbles themselves eventually form black holes with a wide distribution of masses. It has been recently suggested that such black holes may account for LIGO observations and may provide seeds for supermassive black holes observed at galactic centers. They may also provide a significant part or even the whole of the dark matter. We estimate the spectral $\mu$-distortion of the CMB induced by expanding shocks and underdensities. The predicted distortions averaged over the sky are well below the current bounds, but localized peaks due to the largest black holes impose constraints on the model parameters., Comment: 27 pages, 6 figures; v2: Added a brief discussion of PBH effect on structure formation

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 $\rho_c$ varies as a function of the core radius $R_c$ as $\rho_c\propto1/R_c^\beta$ with $\beta\approx1$. 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 $\rho_c$ and $R_c$ for ground state solutions is always in one of the following regimes: (i) $\beta\gg1$, or (ii) $\beta\ll1$, 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., Comment: 10 pages in double column format, 3 figures. V2: Updated to resemble version published in PRD

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., Comment: 39 pages, 12 figures

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.