Search

The occupation fraction of massive black holes (MBHs) in low mass galaxies offers interesting insights into initial black hole seeding mechanisms and their mass assembly history, though disentangling these two effects remains challenging. Using the {\sc Romulus} cosmological simulations we examine the impact of environment on the occupation fraction of MBHs in low mass galaxies. Unlike most modern cosmological simulations, {\sc Romulus} seeds MBHs based on local gas properties, selecting dense (n>3 cm^-3, 15 times the threshold for star formation), pristine (Z3), a difference arises as late-forming dwarfs -- which do not exist in the cluster environment -- begin to dominate in the field and pull the MBH occupation fraction down for low mass galaxies. Additionally, prior to in-fall some cluster dwarfs are similar to progenitors of massive, isolated galaxies, indicating that they might have grown to higher masses had they not been impeded by the cluster environment. While the population of MBHs in dwarf galaxies is already widely understood to be important for understanding MBH formation, this work demonstrates that environmental dependence is important to consider as future observations search for low mass black holes in dwarf galaxies.

The discovery of Little Red Dots (LRDs)—a population of compact, high-redshift, dust-reddened galaxies—is one of the most surprising results from JWST. However, the nature of LRDs is still debated: does the near-infrared emission originate from accreting supermassive black holes (SMBHs), or intense star formation? In this work, we utilize ultra-deep Chandra observations and study LRDs residing behind the lensing galaxy cluster, A2744. We probe the X-ray emission from individual galaxies but find that they remain undetected and provide SMBH mass upper limits of ≲(1.5–16) × 10 ^6 M _⊙ assuming Eddington limited accretion. To increase the signal-to-noise ratios, we conduct a stacking analysis of the full sample with a total lensed exposure time of ≈87 Ms. We also bin the galaxies based on their stellar mass, lensing magnification, and detected broad-line H α emission. For the LRDs exhibiting broad-line H α emission, there is a hint of a stacked signal (∼2.6 σ ), corresponding to an SMBH mass of ∼3.2 × 10 ^6 M _⊙ . Assuming unobscured, Eddington-limited accretion, this black hole (BH) mass is at least 1.5 orders of magnitude lower than that inferred from virial mass estimates using JWST spectra. Given galaxy-dominated stellar mass estimates, our results imply that LRDs do not host overmassive SMBHs and/or accrete at a few percent of their Eddington limit. However, alternative stellar mass estimates may still support that LRDs host overmassive BHs. The significant discrepancy between the JWST and Chandra data hints that the scaling relations used to infer the SMBH mass from the H α line and virial relations may not be applicable for high-redshift LRDs.

The first-generation metal-free stars, referred to as Population III (Pop III) stars, are believed to be the first objects to form out of the pristine gas in the very early Universe. Pop III stars have different structures from the current generation of stars and are important for generating heavy elements and shaping subsequent star formation. However, it is very challenging to directly detect Pop III stars given their high redshifts and short lifetimes. In this Letter, we propose a novel method for detecting Pop III stars through their tidal disruption events (TDEs) by massive black holes. We model the emission properties and calculate the expected rates for these unique TDEs in the early Universe at z ∼ 10. We find that Pop III star TDEs have much higher mass fallback rates and longer evolution timescales compared to solar-type star TDEs in the local Universe, which enhances the feasibility of their detection, although a good survey strategy will be needed for categorizing these sources as transients. We further demonstrate that a large fraction of the flare emissions are redshifted to infrared wavelengths, which can be detected by the JWST and the Nancy Grace Roman Space Telescope (Roman). Last but not least, we find a promising Pop III star TDE detection rate of up to a few tens per year using Roman, based on our current understanding of the black hole mass function in the early Universe.

We report for the first time a sample of 12 supermassive black holes (SMBHs) hosted by low-mass galaxies at cosmic noon, i.e., in a redshift range consistent with the peak of star formation history: z ∼ 1–3. These black holes are 2 orders of magnitude too massive for the stellar content of their hosts when compared with the local relation for active galaxies. These overmassive systems at cosmic noon share similar properties with the high- z sources found ubiquitously in recent James Webb Space Telescope (JWST) surveys (same range of black-hole-to-stellar-mass ratio, bolometric luminosity, and Eddington ratio). We argue that black hole feedback processes, for which there is possible evidence in five of the sources, and the differing environments in galactic nuclei at these respective epochs play a key role in these overmassive systems. These findings contribute to our understanding of the growth and coevolution of SMBHs and their host galaxies across cosmic time, offering a link between the early Universe ( z > 4) observed by JWST and observations of the present-day Universe ( z ≲ 1).

While supermassive black holes (SMBHs) are widely observed in the nearby and distant Universe, their origin remains debated with two viable formation scenarios with light and heavy seeds. In the light seeding model, the seed of the first SMBHs form from the collapse of massive stars with masses of 10–100 M _⊙ , while the heavy seeding model posits the formation of 10 ^4–5 M _⊙ seeds from direct collapse. The detection of SMBHs at redshifts z ≳ 10, edging closer to their formation epoch, provides critical observational discrimination between these scenarios. Here, we focus on the JWST-detected galaxy, GHZ 9, at z ≈ 10 that is lensed by the foreground cluster, A2744. Based on 2.1 Ms deep Chandra observations, we detect a candidate X-ray active galactic nucleus (AGN), which is spatially coincident with the high-redshift galaxy, GHZ 9. The SMBH candidate is inferred to have a bolometric luminosity of $({1.0}_{-0.4}^{+0.5})\times {10}^{46}\ \mathrm{erg}\,{{\rm{s}}}^{-1}$ , which corresponds to a black hole (BH) mass of $({8.0}_{-3.2}^{+3.7})\times {10}^{7}\ {M}_{\odot }$ assuming Eddington-limited accretion. This extreme mass at such an early cosmic epoch suggests the heavy seed origin for this BH candidate. Based on the Chandra and JWST discoveries of extremely high-redshift quasars, we have constructed the first simple AGN luminosity function extending to z ≈ 10. Comparison of this luminosity function with theoretical models indicates an overabundant z ≈ 10 SMBH population, consistent with a higher-than-expected seed formation efficiency.

We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy ( ${f}_{\mathrm{lim},\mathrm{line},5\sigma }\approx 1.2\,\times \,$ 10 ^−18 erg s ^−1 cm ^−2 ) to measure metallicities and star formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.5 < z < 4). In parallel, NGDEEP targets the HUDF-Par2 parallel field with NIRCam ( ${m}_{\mathrm{lim},5\sigma }=30.6-30.9$ ) to discover galaxies to z > 12, constraining the slope of the faint end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST Advanced Camera for Surveys optical imaging in the sky, F435W in the HUDF ( ${m}_{\mathrm{lim},{\rm{F}}435{\rm{W}}}=29.6$ ) and F814W in HUDF-Par2 ( ${m}_{\mathrm{lim},{\rm{F}}814{\rm{W}}}=30$ ), making this a premier HST+JWST deep field. As a treasury survey, NGDEEP data are public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep-field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products.

Recently we found compelling evidence for a gravitational-wave background with Hellings and Downs (HD) correlations in our 15 yr data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes, which produce different interpulsar correlations. In this work, we search the NANOGrav 15 yr data set for evidence of a gravitational-wave background with quadrupolar HD and scalar-transverse (ST) correlations. We find that HD correlations are the best fit to the data and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors ∼2 when comparing HD to ST correlations, and ∼1 for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis.

Supermassive black holes (SMBHs) can be ejected from their galactic centers due to gravitational wave recoil or the slingshot mechanism following a galaxy merger. If an ejected SMBH retains its inner accretion disk, it may be visible as an off-nuclear active galactic nucleus (AGN). At present, only a handful of offset AGNs that are recoil or slingshot candidates have been found, and none have been robustly confirmed. Compiling a large sample of runaway SMBHs would enable us to constrain the mass and spin evolution of binary SMBHs and study feedback effects of displaced AGNs. We adapt the method of varstrometry—which was developed for Gaia observations to identify off-center, dual, and lensed AGNs—in order to quickly identify off-nuclear AGNs in optical survey data by looking for an excess of blue versus red astrometric jitter. We apply this to the Pan-STARRS1 3 π Survey and report on five new runaway AGN candidates. We focus on ZTF18aajyzfv: a luminous quasar offset by 6.7 ± 0.2 kpc from an adjacent galaxy at z = 0.224, and conclude after Keck LRIS spectroscopy and comparison to ASTRID simulation analogs that it is likely a dual AGN. This selection method can be easily adapted to work with data from the soon-to-be commissioned Vera C. Rubin Telescope Legacy Survey of Space and Time (LSST). LSST will have a higher cadence and deeper magnitude limit than Pan-STARRS1, and should permit detection of many more runaway SMBH candidates.

The Local Volume Complete Cluster Survey is an ongoing program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts 0.03 < z < 0.12 and X-ray luminosities in the 0.1–2.4 keV band L _X500c > 10 ^44 erg s ^−1 ) with the Dark Energy Camera, capturing data in the u , g , r , i , z bands with a 5 σ point source depth of approximately 25th–26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19° within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32° within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales.

Several large JWST blank field observing programs have not yet discovered the first galaxies expected to form at 15 ≤ z ≤ 20. This has motivated the search for more effective survey strategies that will be able to effectively probe this redshift range. Here, we explore the use of gravitationally lensed cluster fields, which have historically been the most effective discovery tool with the Hubble Space Telescope. In this paper, we analyze the effectiveness of the most massive galaxy clusters that provide the highest median magnification factor within a single JWST NIRCam module in uncovering this population. The results of exploiting these lensing clusters to break the z > 15 barrier are compared against the results from large-area, blank-field surveys such as JADES and CEERS in order to determine the most effective survey strategy for JWST. We report that the fields containing massive foreground galaxy clusters specifically chosen to occupy the largest fraction of a single NIRCam module with high magnification factors in the source plane while containing all multiple images in the image plane within a single module provide the highest probability of both probing the 15 ≤ z ≤ 20 regime as well as discovering the highest-redshift galaxy possible with JWST. We also find that using multiple massive clusters in exchange for shallower survey depths is a more time-efficient method of probing the z > 15 regime.

We study the black hole mass–host galaxy stellar mass relation, M _BH – M _* , of a sample of z < 4 optically variable active galactic nuclei (AGNs) in the COSMOS field. The parent sample of 491 COSMOS AGNs were identified by optical variability from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) program. Using publicly available catalogs and spectra, we consolidate their spectroscopic redshifts and estimate virial black hole masses using broad-line widths and luminosities. We show that variability searches with deep, high-precision photometry like the HSC-SSP can identity AGNs in low-mass galaxies up to z ∼ 1. However, their black holes are more massive given their host galaxy stellar masses than predicted by the local relation for active galaxies. We report that z ∼ 0.5–4 variability-selected AGNs are meanwhile more consistent with the M _BH – M _* relation for local inactive early-type galaxies. This result is in agreement with most previous studies of the M _BH – M _* relation at similar redshifts and indicates that AGNs selected from variability are not intrinsically different from the broad-line Type 1 AGN population at similar luminosities. Our results demonstrate the need for robust black hole and stellar mass estimates for intermediate-mass black hole candidates in low-mass galaxies at similar redshifts to anchor this scaling relation. Assuming that these results do not reflect a selection bias, they appear to be consistent with self-regulated feedback models wherein the central black hole and stars in galaxies grow in tandem.

Strong gravitational lensing offers a powerful probe of the detailed distribution of matter in lenses, while magnifying and bringing faint background sources into view. Observed strong lensing by massive galaxy clusters, which are often in complex dynamical states, has also been used to map their dark matter (DM) substructures on smaller scales. Deep high-resolution imaging has revealed the presence of strong lensing events associated with these substructures, namely galaxy-scale sub-halos. However, an inventory of these observed galaxy–galaxy strong lensing (GGSL) events is noted to be discrepant with state-of-the-art ΛCDM simulations. Cluster sub-halos appear to be over-concentrated compared to their simulated counterparts yielding an order-of-magnitude higher value of GGSL. In this paper, we explore the possibility of resolving this observed discrepancy by redistributing the mass within observed cluster sub-halos in ways that are consistent within the ΛCDM paradigm of structure formation. Lensing mass reconstructions from data provide constraints on the mass enclosed within apertures and are agnostic to the detailed mass profile within them. Therefore, as the detailed density profile within cluster sub-halos currently remains unconstrained by data, we are afforded the freedom to redistribute the enclosed mass. We investigate if rearranging the mass to a more centrally concentrated density profile helps alleviate the GGSL discrepancy. We report that refitting cluster sub-halos to the ubiquitous ΛCDM-motivated Navarro–Frenk–White profile, and further modifying them to include significant baryonic components, does not resolve this tension. A resolution to this persisting GGSL discrepancy may require more careful exploration of alternative DM models.

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph. We present early results obtained from the first set of observations (Epoch 1, 50% of the allocated orbits) of this program (NGDEEP-NISS1). Using a set of independently developed calibration files designed to deal with a complex combination of overlapping spectra, multiple position angles, and multiple cross filters and grisms, in conjunction with a robust and proven algorithm for quantifying contamination from overlapping dispersed spectra, NGDEEP-NISS1 has achieved a 3 σ sensitivity limit of 2 × 10 ^−18 erg s ^−1 cm ^−2 . We demonstrate the power of deep wide field slitless spectroscopy (WFSS) to characterize the star formation rates, and metallicity ([O iii ]/H β ), and dust content, of galaxies at 1 < z < 3.5. The latter showing intriguing initial results on the applicability and assumptions made regarding the use of Case B recombination. Further, we identify the presence of active galactic nuclei and infer the mass of their supermassive black holes using broadened restframe Mg ii and H β emission lines. The spectroscopic results are then compared with the physical properties of galaxies extrapolated from fitting spectral energy distribution models to photometry alone. The results clearly demonstrate the unique power and efficiency of WFSS at near-infrared wavelengths over other methods to determine the properties of galaxies across a broad range of redshifts.

The recent Chandra-JWST discovery of a quasar in the z ≈ 10.1 galaxy UHZ1 reveals that accreting supermassive black holes were already in place 470 million years after the Big Bang. The Chandra X-ray source detected in UHZ1 is a Compton-thick quasar with a bolometric luminosity of L _bol ∼ 5 × 10 ^45 erg s ^−1 , which corresponds to an estimated black hole (BH) mass of ∼4 × 10 ^7 M _⊙ , assuming accretion at the Eddington rate. JWST NIRCAM and NIRSpec data yield a stellar mass estimate for UHZ1 comparable to its BH mass. These characteristics are in excellent agreement with prior theoretical predictions for a unique class of transient, high-redshift objects, overmassive black hole galaxies (OBGs) by Natarajan et al., that harbor a heavy initial black hole seed that likely formed from the direct collapse of the gas. Given the excellent agreement between the observed multiwavelength properties of UHZ1 and theoretical model template predictions, we suggest that UHZ1 is the first detected OBG candidate. Our assertion rests on multiple lines of concordant evidence between model predictions and the following observed properties of UHZ1: its X-ray detection and the estimated ratio of the X-ray flux to the IR flux, which is consistent with theoretical expectations for a heavy initial BH seed; its high measured redshift of z ≈ 10.1, as predicted for the transient OBG stage (9 < z < 12); the amplitude and shape of the detected JWST spectral energy distribution (SED) between 1 and 5 μ m, which is in very good agreement with simulated template SEDs for OBGs; and the extended JWST morphology of UHZ1, which is suggestive of a recent merge and is also expected for the formation of transient OBGs. As the first OBG candidate, UHZ1 provides compelling evidence for the formation of heavy initial seeds from direct collapse in the early Universe.

Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius of R _B ≈ 2 × 10 ^5 GM _• / c ^2 , where M _• is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ∼ r ^−1 rather than r ^−3/2 and the accretion rate $\dot{M}$ is consequently suppressed by over 2 orders of magnitude below the Bondi rate ${\dot{M}}_{{\rm{B}}}$ . We find continuous energy feedback from the accretion flow to the external medium at a level of $\sim {10}^{-2}\dot{M}{c}^{2}\sim 5\,\times \,{10}^{-5}{\dot{M}}_{{\rm{B}}}{c}^{2}$ . Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.

The James Webb Space Telescope is now detecting early black holes (BHs) as they transition from “seeds” to supermassive BHs. Recently, Bogdan et al. reported the detection of an X-ray luminous supermassive BH, UHZ-1, with a photometric redshift at z > 10. Such an extreme source at this very high redshift provides new insights on seeding and growth models for BHs given the short time available for formation and growth. Harnessing the exquisite sensitivity of JWST/NIRSpec, here we report the spectroscopic confirmation of UHZ-1 at z = 10.073 ± 0.002. We find that the NIRSpec/Prism spectrum is typical of recently discovered z ≈ 10 galaxies, characterized primarily by star formation features. We see no clear evidence of the powerful X-ray source in the rest-frame UV/optical spectrum, which may suggest heavy obscuration of the central BH, in line with the Compton-thick column density measured in the X-rays. We perform a stellar population fit simultaneously to the new NIRSpec spectroscopy and previously available photometry. The fit yields a stellar-mass estimate for the host galaxy that is significantly better constrained than prior photometric estimates ( ${M}_{\star }\sim {1.4}_{-0.4}^{+0.3}\times {10}^{8}$ M _⊙ ). Given the predicted BH mass ( M _BH ∼ 10 ^7 –10 ^8 M _⊙ ), the resulting ratio of M _BH / M _⋆ remains 2 to 3 orders of magnitude higher than local values, thus lending support to the heavy seeding channel for the formation of supermassive BHs within the first billion years of cosmic evolution.

We present a robust sample of very high redshift galaxy candidates from the first epoch of JWST/NIRCam imaging from the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey. The NGDEEP NIRCam imaging, spanning 9.7 arcmin ^2 in the Hubble Ultra Deep Field Parallel Field 2, reaches m = 30.4 (5 σ , point-source, 2″ diameter apertures corrected to total) in F277W, making it the deepest public JWST GO imaging data set to date. We describe our detailed data reduction process of the six-filter broadband JWST/NIRCam imaging, incorporating custom corrections for systematic effects to produce high-quality calibrated images. Using robust photometric redshift selection criteria, we identify a sample of 38 z ≳ 9 galaxy candidates. These objects span a redshift range of z = 8.5–15.8 and apparent magnitudes of m _F277W = 27–30.5 AB mag, reaching ∼1.5 mag deeper than previous public JWST imaging surveys. We calculate the rest-frame ultraviolet luminosity function at z ∼ 9 and 11 and present a new measurement of the luminosity function faint-end slope at z ∼ 11. We find a faint-end slope of α = −2.5 ± 0.4 and −2.2 ± 0.2 at z ∼ 9 and 11, respectively. This is consistent with no significant evolution in the faint-end slope and number density from z = 9 to 11. Comparing our results with theoretical predictions, we find that some models produce better agreement at the faint end than the bright end. These results will help to constrain how stellar feedback impacts star formation at these early epochs.

We look for simulated star-forming linear features such as the one recently discovered by van Dokkum et al. in the cosmological hydrodynamical simulation ASTRID . Among the runaway black holes in ASTRID , none are able to produce clear star-forming wakes. Meanwhile, flyby encounters, typically involving a compact galaxy (with a central black hole) and a star-forming galaxy (with a duo of black holes), reproduce remarkably well many of the key properties (length and linearity, recent star formation, etc.) of the observed star-forming linear feature. We predict that the feature will persist for approximately 100 Myr in such a system and hence constitute a rare event. The feature contains a partly stripped galaxy (with M _gal = 10 ^9 –10 ^10 M _⊙ ) and a dual black hole system ( M _BH = 10 ^5 –10 ^7 M _⊙ ) in its brightest knot. The X-ray emission from AGN in the knot should be detectable in such systems. After 100–200 Myr from the first flyby, the galaxies merge, leaving behind a triple black hole system in a (still) actively star-forming early-type remnant of mass ∼5 × 10 ^10 M _⊙ . Follow-up JWST observations may be key for revealing the nature of these linear features by potentially detecting the older stellar populations constituting the bright knot. Confirmation of such detections may therefore help discriminate a flyby encounter from a massive black hole wake to reveal the origin of such features.

The NANOGrav 15 yr 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 yr 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 subparsec 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.

Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ∼10–100 Hz, over a characteristic jet activity timescale ∼10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ∼10 events yr ^−1 . The detection rate drops to ∼1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events.

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 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 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}\times {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.

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.

The interaction of a runaway supermassive black hole (SMBH) with the circumgalactic medium (CGM) can lead to the formation of a wake of shocked gas and young stars behind it. Here we report the serendipitous discovery of an extremely narrow linear feature in Hubble Space Telescope (HST) Advanced Camera for Surveys images that may be an example of such a wake. The feature extends 62 kpc from the nucleus of a compact star-forming galaxy at z = 0.964. Keck Low-resolution Imaging Spectrometer spectra show that the [O iii ]/H β ratio varies from ∼1 to ∼10 along the feature, indicating a mixture of star formation and fast shocks. The feature terminates in a bright [O iii ] knot with a luminosity of ≈1.9 × 10 ^41 erg s ^−1 . The stellar continuum colors vary along the feature and are well fit by a simple model that has a monotonically increasing age with the distance from the tip. The line ratios, colors, and overall morphology are consistent with an ejected SMBH moving through the CGM at a high speed while triggering star formation. The best-fit time since ejection is ∼39 Myr, and the implied velocity is v _BH ∼ 1600 km s ^−1 . The feature is not perfectly straight in the HST images, and we show that the amplitude of the observed spatial variations is consistent with the runaway SMBH interpretation. Opposite the primary wake is a fainter and shorter feature, marginally detected only in [O iii ] and the rest-frame far-ultraviolet. This feature may be shocked gas behind a binary SMBH that was ejected at the same time as the SMBH that produced the primary wake.

The Event Horizon Telescope (EHT) has led to the first images of a supermassive black hole, revealing the central compact objects in the elliptical galaxy M87 and the Milky Way. Proposed upgrades to this array through the next-generation EHT (ngEHT) program would sharply improve the angular resolution, dynamic range, and temporal coverage of the existing EHT observations. These improvements will uniquely enable a wealth of transformative new discoveries related to black hole science, extending from event-horizon-scale studies of strong gravity to studies of explosive transients to the cosmological growth and influence of supermassive black holes. Here, we present the key science goals for the ngEHT and their associated instrument requirements, both of which have been formulated through a multi-year international effort involving hundreds of scientists worldwide.

This white paper outlines the plans of the History Philosophy Culture Working Group of the Next Generation Event Horizon Telescope Collaboration.

Supermassive black hole binary systems (SMBHBs) should be the most powerful sources of gravitational waves (GWs) in the universe. Once pulsar timing arrays (PTAs) detect the stochastic GW background from their cosmic merger history, searching for individually resolvable binaries will take on new importance. Since these individual SMBHBs are expected to be rare, here we explore how strong gravitational lensing can act as a tool for increasing their detection prospects by magnifying fainter sources and bringing them into view. Unlike for electromagnetic waves, when the geometric optics limit is nearly always valid, for GWs the wave-diffraction-interference effects can become important when the wavelength of the GWs is larger than the Schwarzchild radius of the lens, i.e., ${M}_{\mathrm{lens}}\sim {10}^{8}{\left(\tfrac{f}{\mathrm{mHz}}\right)}^{-1}\,{M}_{\odot }$ . For the GW frequency range explored in this work, the geometric optics limit holds. We investigate GW signals from SMBHBs that might be detectable with current and future PTAs under the assumption that quasars serve as bright beacons that signal a recent merger. Using the black hole mass function derived from quasars and a physically motivated magnification distribution, we expect to detect a few strongly lensed binary systems out to z ≈ 2. Additionally, for a range of fixed magnifications 2 ≤ μ ≤ 100, strong lensing adds up to ∼30 more detectable binaries for PTAs. Finally, we investigate the possibility of observing both time-delayed electromagnetic signals and GW signals from these strongly lensed binary systems—that will provide us with unprecedented multi-messenger insights into their orbital evolution.

We use the Galaxy Morphology Posterior Estimation Network (GaMPEN) to estimate morphological parameters and associated uncertainties for ∼8 million galaxies in the Hyper Suprime-Cam Wide survey with z ≤ 0.75 and m ≤ 23. GaMPEN is a machine-learning framework that estimates Bayesian posteriors for a galaxy’s bulge-to-total light ratio ( L _B / L _T ), effective radius ( R _e ), and flux ( F ). By first training on simulations of galaxies and then applying transfer learning using real data, we trained GaMPEN with

We present QUOTAS, a novel research platform for the data-driven investigation of supermassive black hole (SMBH) populations. While SMBH data—observations and simulations—have grown in complexity and abundance, our computational environments and tools have not matured commensurately to exhaust opportunities for discovery. To explore the BH, host galaxy, and parent dark matter halo connection—in this pilot version—we assemble and colocate the high-redshift, z > 3 quasar population alongside simulated data at the same cosmic epochs. As a first demonstration of the utility of QUOTAS, we investigate correlations between observed Sloan Digital Sky Survey (SDSS) quasars and their hosts with those derived from simulations. Leveraging machine-learning algorithms (ML), to expand simulation volumes, we show that halo properties extracted from smaller dark-matter-only simulation boxes successfully replicate halo populations in larger boxes. Next, using the Illustris-TNG300 simulation that includes baryonic physics as the training set, we populate the larger LEGACY Expanse dark-matter-only box with quasars, and show that observed SDSS quasar occupation statistics are accurately replicated. First science results from QUOTAS comparing colocated observational and ML-trained simulated data at z 3 are presented. QUOTAS demonstrates the power of ML, in analyzing and exploring large data sets, while also offering a unique opportunity to interrogate theoretical assumptions that underpin accretion and feedback models. QUOTAS and all related materials are publicly available at the Google Kaggle platform. (The full data set—observational data and simulation data—are available at: https://www.kaggle.com/ and the codes are available at: https://www.kaggle.com/datasets/quotasplatform/quotas )

Exploiting the fundamentally achromatic nature of gravitational lensing, we present a lens model for the massive galaxy cluster SMACS J0723.3−7323 (SMACS J0723; z = 0.388) that significantly improves upon earlier work. Building on strong-lensing constraints identified in prior Hubble Space Telescope (HST) observations, the mass model utilizes 21 multiple-image systems, 17 of which were newly discovered in Early Release Observation data from the JWST. The resulting lens model maps the cluster mass distribution to an rms spatial precision of 0.″32, and is publicly available. Consistent with previous analyses, our study shows SMACS J0723.3 to be well described by a single large-scale component centered on the location of the brightest cluster galaxy. However, satisfying all lensing constraints provided by the JWST data, the model points to the need for the inclusion of an additional, diffuse component west of the cluster. A comparison of the galaxy, mass, and gas distributions in the core of SMACS J0723 based on HST, JWST, and Chandra data reveals a concentrated regular elliptical profile along with tell-tale signs of a recent merger, possibly proceeding almost along our line of sight. The exquisite sensitivity of JWST’s NIRCam reveals in spectacular fashion both the extended intracluster light distribution and numerous star-forming clumps in magnified background galaxies. The high-precision lens model derived here for SMACS J0723 demonstrates the unprecedented power of combining HST and JWST data for studies of structure formation and evolution in the distant universe.

We present new diagnostic metrics to probe the dynamical state of galaxy clusters. These novel metrics rely on the computation of the power spectra of the matter and gas distributions and their cross-correlation derived from cluster observations. This analysis permits us to cross-correlate the fluctuations in the matter distribution, inferred from high-resolution lensing mass maps derived from Hubble Space Telescope (HST) data, with those derived from the emitted X-ray surface brightness distribution of the hot intracluster medium from the Chandra X-ray Observatory. These methodological tools allow us to quantify with unprecedented resolution the coherence with which the gas traces the mass and interrogate the assumption that the gas is in hydrostatic equilibrium with the underlying gravitational potential. We characterize departures from equilibrium as a function of scale with a new gas-mass coherence parameter. The efficacy of these metrics is demonstrated by applying them to the analysis of two representative clusters known to be in different dynamical states: the massive merging cluster A2744, from the HST Frontier Fields sample, and the dynamically relaxed cluster A383, from the Cluster Lensing and Supernova Survey with the Hubble sample. Using lensing mass maps in combination with archival Chandra data, and simulated cluster analogs available from the OMEGA500 suite, we quantify the fluctuations in the mass and X-ray surface brightness and show that new insights into the dynamical state of the clusters can be obtained from our gas-mass coherence analysis.

Dwarf galaxies are promising test beds for constraining models of supermassive and intermediate-mass black holes (MBHs) via their BH occupation fraction (BHOF). Disentangling seeding from the confounding effects of mass assembly over a Hubble time is a challenging problem that we tackle in this study with a suite of semianalytical models (SAMs). We show how the measured BHOF depends on the lowest BH mass or active galactic nucleus (AGN) luminosity achieved by a survey. To tell seeding models apart, we need to detect or model all AGNs brighter than 10 ^37 erg s ^−1 in galaxies of M _* ∼ 10 ^8−10 M _⊙ . Shallower surveys, like eRASS, cannot distinguish between seed models even with the compensation of a much larger survey volume. We show that the AMUSE survey, with its inference of the MBH population underlying the observed AGNs, strongly favors heavy seed models, growing with either a power-law Eddington ratio distribution function or one in which BH accretion is tied to the star formation rate (i.e., the AGN-main sequence, AGN-MS, model). These two growth channels can then be distinguished by the AGN luminosity function at >10 ^40 erg s ^−1 , with the AGN-MS model requiring more accretion than observed at z ∼ 0. Thus, current X-ray observations favor heavy seeds whose Eddington ratios follow a power-law distribution. The different models also predict different radio scaling relations, which we quantify using the fundamental plane of BH activity. We close with recommendations for the design of upcoming multiwavelength campaigns that can optimally detect MBHs in dwarf galaxies.

We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models, including a constant electron to magnetic pressure (constant βe model) with a population of non-thermal electrons as well as an R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations are enhanced with every pair added. On the contrary, we show that, at higher frequencies, a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of a large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function.

The next generation Event Horizon Telescope (ngEHT) will observe multiple supermassive black hole (SMBH) candidates down to a few tens of mJy, and profoundly transform our understanding of the local SMBH population. Given the impossibility of large-area high-resolution millimeter surveys, multi-frequency spectral energy densities (SEDs), and models are required to both identify source samples tailored to specific science goals, and to predict the feasibility of detection of individual interesting sources. Here, we present the Event Horizon and Environs (ETHER) source and SED model database whose primary use is to enable the selection and optimization of targets for EHT and ngEHT science. The living ETHER database currently consolidates 1.6 million black hole mass estimates, ∼15,500 milliarcsec-scale radio fluxes, ∼14,000 hard X-ray fluxes (expected to grow by factor ≳40 with the eROSITA data release) and SED information as obtained from catalogs and database queries, the literature, and our own new observations. Jet and accretion flow models are fit to individual SEDs in an automated way in order to predict the ngEHT observable fluxes from the jet base and accretion inflow. The database can be filtered by parameters or cross matched to a user source list, with the automated SED fitting models optionally fine tuned by the user. We have identified an initial ngEHT ‘gold sample’ for jet base studies and potentially black hole shadows; this sample will grow significantly in the coming years. While the ngEHT requires and will best exploit the ETHER database, six (eleven) ETHER sources have already been observed (scheduled) with the EHT in 2022 (2023), and the database has wide ranging applications in galaxy and black hole mass evolution studies.

We propose the tracing of the motion of a shearing hot spot near the Sgr A* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its boundary during its orbit, resulting in a distinct signature. We subdivide the motion into two different phases: the first phase refers to the appearance of the hot spot modeled as a bright blob, followed by a subsequent shearing phase. We employ different observational array configurations, including EHT (2017, 2022) and the next-generation Event Horizon Telescope (ngEHTp1, ngEHT) arrays, with several new sites added, and make dynamical image reconstructions for each of them. Subsequently, we infer the hot spot angular image location in the first phase, followed by the axes ratio and the ellipse area in the second phase. We focus on the direct observability of the orbiting hot spot in the sub-mm wavelength. Our analysis demonstrates that for this particular simulation, the newly added dishes are better able to trace the first phase as well as part of the second phase before the flux is reduced substantially, compared to the EHT arrays. The algorithm used in this work can be easily extended to other types of dynamics, as well as different shearing timescales. More simulations are required to prove whether the current set of newly proposed sites are sufficient to resolve any motions near variable sources, such as Sgr A*.

We present estimates for the number of supermassive black holes (SMBHs) for which the next-generation Event Horizon Telescope (ngEHT) can identify the black hole “shadow”, along with estimates for how many black hole masses and spins the ngEHT can expect to constrain using measurements of horizon-resolved emission structure. Building on prior theoretical studies of SMBH accretion flows and analyses carried out by the Event Horizon Telescope (EHT) collaboration, we construct a simple geometric model for the polarized emission structure around a black hole, and we associate parameters of this model with the three physical quantities of interest. We generate a large number of realistic synthetic ngEHT datasets across different assumed source sizes and flux densities, and we estimate the precision with which our defined proxies for physical parameters could be measured from these datasets. Under April weather conditions and using an observing frequency of 230 GHz, we predict that a “Phase 1” ngEHT can potentially measure ∼50 black hole masses, ∼30 black hole spins, and ∼7 black hole shadows across the entire sky.

While supermassive black-hole masses have been cataloged across cosmic time, only a few dozen of them have robust spin measurements. By extending and improving the existing Event Horizon Telescope (EHT) array, the next-generation Event Horizon Telescope (ngEHT) will enable multifrequency, polarimetric movies on event-horizon scales, which will place new constraints on the space-time and accretion flow. By combining this information, it is anticipated that the ngEHT may be able to measure tens of supermassive black-hole masses and spins. In this white paper, we discuss existing spin measurements and many proposed techniques with which the ngEHT could potentially measure spins of target supermassive black holes. Spins measured by the ngEHT would represent a completely new sample of sources that, unlike pre-existing samples, would not be biased towards objects with high accretion rates. Such a sample would provide new insights into the accretion, feedback, and cosmic assembly of supermassive black holes.

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.