Active Supermassive Black Holes Research Articles

Observational Evidence for a Correlation between the Magnetic Field of Jets and Star Formation Rate in Host Galaxies

Abstract Accretion supermassive black holes in the center of active galaxies usually produce “jet”-collimated bipolar outflows of relativistic particles. Magnetic fields near the black hole event horizon may play a crucial role in the formation of jets/outflows. Both theory and observation indicate that jets/outflows driven by centrally active supermassive black holes have a feedback effect on the overall properties of the host galaxies. Therefore, the magnetic field is a key ingredient for the formation and evolution of galaxies. Here, we report a clear correlation between the magnetic field of jets and star formation rate for a large sample of 96 galaxies hosting supermassive black holes, which suggests that the star formation of active galactic nuclei host galaxies may be powered by the jets.

VODKA-JWST: A 3.8 kpc Dual Quasar at Cosmic Noon in a Powerful Starburst Galaxy with JWST/MIRI Integral Field Unit

Dual quasars—two active supermassive black holes at galactic scales—represent crucial objects for studying the impact of galaxy mergers and quasar activity on the star formation rate (SFR) within their host galaxies, particularly at cosmic noon when SFR peaks. We present JWST/MIRI mid-infrared integral field spectroscopy of J074922.96+225511.7, a dual quasar with a projected separation of 3.8 kpc at a redshift z = 2.17. We detect spatially extended [Fe ii] 5.34 μm and polycyclic aromatic hydrocarbon (PAH) 3.3 μm emissions from the star formation activity in its host galaxy. We derive the SFR of 103.0±0.2 M ⊙ yr−1 using PAH 3.3 μm, which is 5 times higher than that derived from the knee of the infrared luminosity function for galaxies at z ∼ 2. While the SFR of J0749+2255 agrees with that of star-forming galaxies of comparable stellar mass at the same redshifts, its molecular gas content falls short of expectations based on the molecular Kennicutt–Schmidt law. This discrepancy may result from molecular gas depletion due to the longer elevated stage of star formation, even after the molecular gas reservoir is depleted. We do not observe any quasar-driven outflow that impacts PAH and [Fe ii] in the host galaxy based on the spatially resolved maps. From the expected flux in PAH-based star formation, the [Fe ii] line likely originates from the star-forming regions in the host galaxy. Our study highlights the extreme stardust nature of J0749+2255, indicating a potential connection between the dual quasar phase and intense star formation activities.

The high energy X-ray probe (HEX-P): the future of hard X-ray dual AGN science

A fundamental goal of modern-day astrophysics is to understand the connection between supermassive black hole (SMBH) growth and galaxy evolution. Merging galaxies offer one of the most dramatic channels for galaxy evolution known, capable of driving inflows of gas into galactic nuclei, potentially fueling both star formation and central SMBH activity. Dual active galactic nuclei (dual AGNs) in late-stage mergers with nuclear pair separations <10 kpc are thus ideal candidates to study SMBH growth along the merger sequence since they coincide with the most transformative period for galaxies. However, dual AGNs can be extremely difficult to confirm and study. Hard X-ray (>10 keV) studies offer a relatively contamination-free tool for probing the dense obscuring environments predicted to surround the majority of dual AGN in late-stage mergers. To date, only a handful of the brightest and closest systems have been studied at these energies due to the demanding instrumental requirements involved. We demonstrate the unique capabilities of HEX-P to spatially resolve the soft and - for the first time - hard X-ray counterparts of closely-separated (∼2″−5″) dual AGNs in the local Universe. By incorporating physically-motivated obscuration models, we reproduce realistic broadband X-ray spectra expected for deeply embedded accreting SMBHs. Hard X-ray spatially resolved observations of dual AGNs—accessible only to HEX-P—will hence transform our understanding of dual AGN in the nearby Universe.

A JVLA, LOFAR, e-Merlin, VLBA, and EVN study of RBS 797: can binary supermassive black holes explain the outburst history of the central radio galaxy?

The multifaceted central radio galaxy of the cluster RBS 797 shows several episodes of jet activity in multiple directions. We wish to understand the causes behind these dramatic misalignments and measure the timescales of the successive outbursts. We present a multifrequency (144 MHz -- 9 GHz) and multiscale (5 pc -- $50$ kpc) investigation of the central radio galaxy in RBS 797, by means of JVLA, LOFAR (with international stations), e-Merlin, VLBA, and EVN data. We investigate the morphological and spectral properties of the radio lobes, the jets, and the active core. We confirm the co-spatiality of the radio lobes with the four perpendicular X-ray cavities previously discovered. The radiative ages of the east-west lobes ($31.4 Myr) and of the north-south lobes ($32.1 Myr) support a coeval origin of the perpendicular outbursts, which also have similar active phase duration (sim 12 Myr). Based on the analysis of the inner north-south jets (on scales of $ kpc), we (a) confirm the S-shaped jet morphology; (b) show the presence of two hotspots per jet with a similar spectral index; and (c) estimate the age of the twisting north-south jets to be less than $ Myr. Based on these results, we determine that jet precession, with a period sim 9 Myr, half-opening angle sim circ $, and jet advance speed sim 0.01$c$, can explain the properties of the north-south jets. We also find that the synchrotron injection index has steepened from the large, older outbursts ($ to the younger S-shaped jets ($ possibly due to a transition from FR I-like to FR II-like activity. The e-Merlin, VLBA, and EVN data reveal a single, compact core at the heart of RBS 797, surrounded by extended radio emission whose orientation depends on the spatial scale sampled by the data. We explore several engine-based scenarios to explain these results. Piecing together the available evidence, we argue that RBS 797 likely hosts (or hosted) binary active supermassive black holes (SMBHs). The detection of a single component in the VLBA and EVN data is still consistent with this interpretation, since the predicted separation of the binary SMBHs (leq 0.6 pc) is an order of magnitude smaller than the resolution of the available radio data (5 pc).

Joint Analysis of the Iron Emission in the Optical and Near-Infrared Spectrum of I Zw 1

Constraining the physical conditions of the ionized media in the vicinity of an active supermassive black hole (SMBH) is crucial to understanding how these complex systems operate. Metal emission lines such as iron (Fe) are useful probes to trace the gaseous media’s abundance, activity, and evolution in these accreting systems. Among these, the Feii emission has been the focus of many prior studies to investigate the energetics, kinematics, and composition of the broad-emission line region (BELR) from where these emission lines are produced. In this paper, we present the first simultaneous Feii modeling in the optical and near-infrared (NIR) regions. We use cloudy photoionization code to simulate both spectral regions in the wavelength interval 4000–12,000 Å. We compare our model predictions with the observed line flux ratios for I Zw (Zwicky) 1—a prototypical strong Feii-emitting active galactic nuclei (AGNs). This allows setting constraints on the BLR cloud density and metal content that is optimal for the production of the Feii emission, which can be extended to I Zw 1-like sources by examining a broad parameter space. We demonstrate the salient and distinct features of the Feii pseudo-continuum in the optical and NIR, giving special attention to the effect of micro-turbulence on the intensity of the Feii emission.

A Massive Protocluster Anchored by a Luminous Quasar at z = 6.63

Protoclusters, the progenitors of galaxy clusters, trace large scale structures in the early Universe and are important to our understanding of structure formation and galaxy evolution. To date, only a handful of protoclusters have been identified in the Epoch of Reionization. As one of the rarest populations in the early Universe, distant quasars that host active supermassive black holes are thought to reside in the most massive dark matter halos at that cosmic epoch and could thus potentially pinpoint some of the earliest protoclusters. In this Letter, we report the discovery of a massive protocluster around a luminous quasar at z = 6.63. This protocluster is anchored by the quasar and includes three [C ii] emitters at z ∼ 6.63, 12 spectroscopically confirmed Lyα emitters (LAEs) at 6.54 < z ≤ 6.64, and a large number of narrow-band-imaging selected LAE candidates at the same redshift. This structure has an overall overdensity of δ=3.3−0.9+1.1 within ∼35 × 74 cMpc2 on the sky and an extreme overdensity of δ > 30 in its central region (i.e., R ≲ 2 cMpc). We estimate that this protocluster will collapse into a galaxy cluster with a mass of 6.9−1.4+1.2×1015M⊙ at the current epoch, more massive than the most massive clusters known in the local Universe such as Coma. In the quasar vicinity, we discover a double-peaked LAE, which implies that the quasar has a UV lifetime greater than 0.8 Myrs and has already ionized its surrounding intergalactic medium.

The active CGCG 077-102 NED02 galaxy within the Abell 2063 galaxy cluster

Context. Within the framework of investigating the link between the central super massive black holes in the cores of galaxies and the galaxies themselves, we detected a variable X-ray source in the center of CGCG 077-102 NED02, which is a member of the CGCG 077-102 galaxy pair within the Abell 2063 cluster of galaxies. Aims. Our goal is to combine X-ray and optical data to demonstrate that this object harbors an active super massive black hole in its core, and to relate this to the dynamical status of the galaxy pair within the Abell 2063 cluster. Methods. We used Chandra and XMM-Newton archival data to derive the X-ray spectral shape and variability. We also obtained optical spectroscopy to detect the expected emission lines that are typically found in active galactic nuclei. Finally, we used public ZTF imaging data to investigate the optical variability. Results. There is no evidence of multiple X-ray sources or extended components within CGCG 077-102 NED02. Single X-ray spectral models fit the source well. We detect significant, nonrandom inter-observation 0.5–10 keV X-ray flux variabilities, for observations separated by ∼4 days for short-term variations and by up to ∼700 days for long-term variations. Optical spectroscopy points toward a passive galaxy for CGCG 077-102 NED01 and a Seyfert for CGCG 077-102 NED02. The classification of CGCG 077-102 NED02 is also consistent with its X-ray luminosity of over 1042 erg s−1. We do not detect short-term variability in the optical ZTF light curves. However, we find a significant long-term stochastic variability in the g-band that can be well described by the damped random walk model with a best-fit characteristic damping timescale of τDRW = 30−12+28 days. Finally, the CGCG 077-102 galaxy pair is deeply embedded within the Abell 2063 potential, with a long enough history within this massive structure to have been affected by the influence of this cluster for a long time. Conclusions. Our observations point toward a moderately massive black hole in the center of CGCG 077-102 NED02 of ∼106 M⊙. As compared to another similar pair in the literature, CGCG 077-102 NED02 is not heavily obscured, perhaps because of the surrounding intracluster medium ram-pressure stripping.

Blueshifted Fe ii and Balmer Broad Emission lines in a Reddened Quasar

We present detailed analysis of the optical and near-infrared spectrum of the reddened radio-loud quasar F2M J154043.74+492323.9 with strong optical Fe ii emission of . It is found that the peaks of broad emission lines (BELs), including Hα, Hβ, and Hγ, and the peaks of individual Fe ii BELs, are blueshifted as large as v ≈ 1500 km s−1 with respect to the systemic redshift determined by the low-ionization [O ii] narrow emission line. BEL components at the systemic redshift are also significantly detected, though much weaker than the corresponding blueshifted components. Interestingly, the Mg ii BEL also shows a similar double-peaked profile, yet with a weaker blueshifted component relative to the unshifted one. We interpret these exceptional BELs as originating from massive outflows, possibly with both densities and column densities much higher than those of the normal BEL region (BELR). Alternative scenarios, however, cannot be completely ruled out based on the data available, including a gravitationally bound binary of active supermassive black holes (SMBHs), each with an unattached BELR, or a recoiling SMBH caused by anisotropic gravitational radiation during the coalescence of the progenitor binary. Follow-up observations will be able to disentangle these scenarios, and to also provide valuable clues to understand the growth of SMBHs and the assembly of their host galaxies.

Varstrometry for Off-nucleus and Dual Sub-kiloparsec Active Galactic Nuclei (VODKA): Very Long Baseline Array Searches for Dual or Off-nucleus Quasars and Small-scale Jets

Dual and off-nucleus active supermassive black holes are expected to be common in the hierarchical structure formation paradigm, but their identification at parsec scales remains a challenge due to strict angular resolution requirements. We conducted a systematic study using the Very Long Baseline Array (VLBA) to examine 23 radio-bright candidate dual and off-nucleus quasars. The targets are selected by a novel astrometric technique (varstrometry) from Gaia, aiming to identify dual or off-nucleus quasars at (sub)kiloparsec scales. Among these quasars, eight exhibit either multiple radio components or significant (>3σ) positional offsets between the VLBA and Gaia positions. The radio emission from the three candidates, which exhibit multiple radio components, is likely to originate from small-scale jets based on their morphology. Among the remaining five candidates with significant VLBA-Gaia offsets, three are identified as potential dual quasars at parsec scales, one is likely attributed to small-scale jets, and the origin of the last candidate remains unclear. We explore alternative explanations for the observed VLBA-Gaia offsets. We find no evidence for optical jets at kiloparsec scales, nor any contamination to Gaia astrometric noise from the host galaxy; misaligned coordinate systems are unlikely to account for our offsets. Our study highlights the promise of the varstrometry technique in discovering candidate dual or off-nucleus quasars and emphasizes the need for further confirmation and investigation to validate and understand these intriguing candidates.

Tracing obscured galaxy build-up at high redshift using deep radio surveys

Context. A fundamental question of extra-galactic astronomy that is yet to be fully understood, concerns the evolution of the star formation rate (SFR) and supermassive black hole (SMBH) activity with cosmic time, as well as their interplay and how it impacts galaxy evolution. A primary focus that could shed more light on these questions is the study of merging systems, comprising highly star-forming galaxies (SFGs) and active galactic nuclei (AGN) at the earliest stages of galactic formation. However, considering the challenges associated with identifying these objects, it is essential to explore complementary selection methods across multiple wavelengths. Aims. The primary objective of this study is to conduct a comprehensive analysis of a sample of high-redshift (z &gt; 3) far-infrared (far-IR) and radio-emitting galaxies in the highest possible spatial resolution. The aim is to study the properties of this population, such as their morphological characteristics, and to explore the interplay of SFR and SMBH activity at this epoch. Methods. In order to select the galactic population of our interest, we employed two selection criteria that have frequently been used as separate methods in the literature. In more detail, we selected galaxies that present relatively compact radio morphologies at 1.4 GHz (i.e., an angular size smaller than 10 arcsec) as well as a far-IR spectrum that peaks in flux at λ ≥ 350 μm (i.e., flux350 μm &gt; flux250 μm). For these selection criteria, we used the COSMOS and ECDF-S fields, two of the most extensively observed astronomical fields currently available, which provide high spectral and spatial resolution at a multi-wavelength scale. By accepting only galaxies that satisfied these selection criteria, we derived a sample of eight galaxies that were identified either photometrically or spectroscopically at z &gt; 3 from literature studies and by our team. Results. A thorough investigation of available optical, near-IR, and millimetre (mm) imaging reveals a possible merging scenario in five out of eight cases in our sample. Additionally, available multi-wavelength photometry strongly suggests active star formation at the 103 M⊙ yr−1 level in massive systems (stellar masses of M⋆ ∼ 1011 M⊙) co-hosting an active SMBH. Conclusions. Comparison of these results with previous studies, suggests that our selection method preferentially identifies galaxies hosting an active SMBH, as well as a strong SFG component, resulting in high SFR and IR luminosity. An additional examination of the efficacy of the radio and far-IR selection criteria provides further support for their combined application in selecting co-evolving AGN and star formation activity at high redshift. In this regard, future use of these selection criteria on radio and far-IR/mm observations of statistically larger galaxy samples is of high interest.

AMUSE-Antlia. I. Nuclear X-Ray Properties of Early-type Galaxies in a Dynamically Young Galaxy Cluster

To understand the formation and growth of supermassive black holes (SMBHs) and their coevolution with host galaxies, it is essential to know the impact of environment on the activity of active galactic nuclei (AGNs). We present new Chandra X-ray observations of nuclear emission from member galaxies in the Antlia cluster, the nearest non-cool core and the nearest merging galaxy cluster, residing at D = 35.2 Mpc. Its inner region, centered on two dominant galaxies NGC 3268 and NGC 3258, has been mapped with three deep Chandra ACIS-I pointings. Nuclear X-ray sources are detected in 7/84 (8.3%) early-type galaxies (ETG) and 2/8 (25%) late-type galaxies with a median detection limit of 8 × 1038 erg s−1. All nuclear X-ray sources but one have a corresponding radio continuum source detected by MeerKAT at the L band. Nuclear X-ray sources detected in early-type galaxies are considered the genuine X-ray counterpart of low-luminosity AGN. When restricted to a detection limit of and a stellar mass of , six of 11 ETGs are found to contain an X-ray AGN in Antlia, exceeding the AGN occupation fraction of 7/39 (18.0%) and 2/12 (16.7%) in the more relaxed, cool core clusters, Virgo and Fornax, respectively, and rivaling that of the AMUSE-Field ETG of 27/49 (55.1%). Furthermore, more than half of the X-ray AGN in Antlia is hosted by its younger subcluster, centered on NGC 3258. We believe that this is because SMBH activity is enhanced in a dynamically young cluster compared to relatively relaxed clusters.

AT2018dyk revisited: a tidal disruption event candidate with prominent infrared echo and delayed X-ray emission in a LINER galaxy

ABSTRACT The multiwavelength data of nuclear transient AT2018dyk, initially discovered as a changing-look low-ionization nuclear emission-line region (LINER) galaxy, has been revisited by us and found being in agreement with a tidal disruption event (TDE) scenario. The optical light curve of AT2018dyk declines as a power-law form approximately with index –5/3 yet its X-ray emission lags behind the optical peak by ∼140 d, both of which are typical characteristics for TDEs. The X-ray spectra are softer than normal active galactic nuclei (AGNs) although they show a slight trend of hardening. Interestingly, its rising time-scale belongs to the longest among TDEs while it is nicely consistent with the theoretical prediction from its relatively large supermassive black hole (SMBH) mass (∼107.38M⊙). Moreover, a prominent infrared echo with peak luminosity ∼7.4 × 1042 erg s−1 has been also detected in AT2018dyk, implying an unusually dusty subparsec nuclear environment in contrast with other TDEs. In our sample, LINERs share similar covering factors with AGNs, which indicates the existence of the dusty torus in these objects. Our work suggests that the nature of nuclear transients in LINERs needs to be carefully identified and their infrared echoes offer us a unique opportunity for exploring the environment of SMBHs at low accretion rate, which has been so far poorly explored but is crucial for understanding the SMBH activity.

AT 2021loi: A Bowen Fluorescence Flare with a Rebrightening Episode Occurring in a Previously Known AGN

The optical-ultraviolet transient AT 2021loi is located at the center of its host galaxy. Its spectral features identify it as a member of the Bowen fluorescence flare (BFF) class. The first member of this class was considered to be related to a tidal disruption event, but enhanced accretion onto an already active supermassive black hole was suggested as an alternative explanation. Having occurred in a previously known unobscured active galactic nucleus, AT 2021loi strengthens the latter interpretation. Its light curve is similar to those of previous BFFs, showing a rebrightening approximately 1 yr after the main peak (which was not explicitly identified but might be the case in all previous BFFs). An emission feature around 4680 Å, seen in the preflare spectrum, strengthens by a factor of ∼2 around the optical peak of the flare and is clearly seen as a double-peaked feature then, suggesting a blend of N iii λ4640 with He ii λ4686 as its origin. The appearance of O iii λ3133 and possible N iii λλ4097, 4103 (blended with Hδ) during the flare further support a Bowen fluorescence classification. Here we present ZTF, ATLAS, Keck, Las Cumbres Observatory, NEOWISE-R, Swift AMI, and Very Large Array observations of AT 2021loi, making it one of the best-observed BFFs to date. It thus provides some clarity on the nature of BFFs but also further demonstrates the diversity of nuclear transients.

An increase in black hole activity in galaxies with kinematically misaligned gas

External accretion events such as a galaxy merger or the accretion of gas from the immediate environment of a galaxy can create a large misalignment between the gas and the stellar kinematics. Numerical simulations have suggested that misaligned structures may promote the inflow of gas to the nucleus of the galaxy and the accretion of gas by the central supermassive black hole. We show for the first time that galaxies with a strong misalignment between the ionized gas and stellar kinematic angles have a higher observed fraction of active black holes than galaxies with aligned rotation of gas and stars. The increase in black hole activity suggests that the process of formation and/or the presence of misaligned structures are connected with the fuelling of active supermassive black holes.

Discovery of 24 radio-bright quasars at 4.9 ≤z≤ 6.6 using low-frequency radio observations

High-redshift quasars (z &gt; 5) that also shine brightly at radio wavelengths are unique signposts of supermassive black hole activity in the early universe. However, bright radio sources atz ≥ 5 are extremely rare and therefore we have started a campaign to search for new high-zquasars by combining an optical dropout selection driven by theg,r, andzbands from the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys with low-frequency radio observations from the LOFAR Two-metre Sky Survey (LoTSS). Currently, LoTSS covers a large fraction of the northern sky (∼5720 deg2) to such a depth (median noise level ∼83 μJy beam−1) that about 30% of the general quasar population is detected − which is a factor of 5–10 more than previous large sky radio surveys such as NVSS and FIRST, respectively. In this paper, we present the discovery of 20 new quasars (and the independent confirmation of four) between 4.9 ≤ z ≤ 6.6. Out of the 24 quasars, 21 satisfy the traditional radio-loudness criterion ofR = f5 GHz/f4400 Å &gt; 10, with the full sample spanningR ∼ 6–1000, thereby more than doubling the sample of known radio-loud quasars atz ≥ 5. Our radio detection requirement strongly decreases the contamination of stellar sources and allows one to select these quasars in a broad redshift range. Despite selecting our quasar candidates using fewer and less conservative colour restrictions, both the optical and near-infrared colours, Lyαemission line properties, and dust reddening,E(B − V), measurements of our quasar sample do not deviate from the known radio-quiet quasar population, suggesting similar optical quasar properties of the radio-loud and radio-quiet quasar population at high-z. Our campaign demonstrates the potential for discovering new high-zquasar populations through next generation radio continuum surveys.

Multiwavelength Observations of the Obscuring Wind in the Radio-quiet Quasar MR 2251-178

Obscuring winds driven away from active supermassive black holes are rarely seen due to their transient nature. They have been observed with multiwavelength observations in a few Seyfert 1 galaxies and one broad absorption line radio-quiet quasar so far. An X-ray obscuration event in MR 2251-178 was caught in late-2020, which triggered multiwavelength (near-IR (NIR) to X-ray) observations targeting this radio-quiet quasar. In the X-ray band, the obscurer leads to a flux drop in the soft X-ray band from late-2020 to early-2021. X-ray obscuration events might have a quasi-period of two decades considering earlier events in 1980 and 1996. In the UV band, a forest of weak blueshifted absorption features emerged in the blue wing of Lyα λ1216 in late-2020. Our XMM-Newton, NuSTAR, and Hubble Space Telescope/COS observations are obtained simultaneously; hence, the transient X-ray obscuration event is expected to account for the UV outflow, although they are not necessarily caused by the same part of the wind. Both blueshifted and redshifted absorption features were found for He i λ10830, but no previous NIR spectra are available for comparison. The X-ray observational features of MR 2251-178 shared similarities with some other type 1 active galactic nuclei (AGNs) with obscuring wind. However, observational features in the UV to NIR bands are distinctly different from those seen in other AGN with obscuring winds. A general understanding of the observational variety and the nature of obscuring wind is still lacking.

Blazar nature of high-z radio-loud quasars

We report on the Swift/XRT observation and classification of eleven blazar candidates at z &gt; 4. These sources were selected as part of a sample of extremely radio-loud quasars, with a focus on quasars with jets oriented roughly close to our line of sight. Deriving their viewing angles and their jets’ bulk Lorentz factors was crucial for a strict blazar classification, which was made possible only thanks to X-ray observations. Out of eleven sources, five show strong and hard X-ray fluxes that set the foundation for their blazar classification, while two are uncertain and three host relativistic jets that we observe just outside their beaming cone (i.e. are not strictly blazars), while one went undetected by Swift/XRT. Following this approach, we were able to trace the &gt; ​109 M⊙ active super-massive black hole (SMBH) population hosted in jetted active galactic nuclei (AGN). At z ≥ 4, the massive jetted sources are likely predominant in the overall quasar population: this calls for a deep review of our understanding of the first SMBH formations and evolution. Jets are indeed key actors in fast accretion and must be searched for across the whole high-redshift quasar population. A note of caution must be added: radio-loudness and, in general, radio features at high redshifts seem do not seem to perfectly reflect high-energy properties. A strong effect attributed to the interaction with cosmic microwave background (CMB) radiation is surely at work, which quenches the radio emission with respect to the X-rays; however, in addition, more frequent occasions for the jet to be bent seem to play a relevant role in this regard. Thus, classifications and population studies must be carefully performed, so as to avoid interference resulting from these inconsistencies.

Variable structures in M87* from space, time and frequency resolved interferometry

The immediate vicinity of an active supermassive black hole—with its event horizon, photon ring, accretion disk and relativistic jets—is an appropriate place to study physics under extreme conditions, particularly general relativity and magnetohydrodynamics. Observing the dynamics of such compact astrophysical objects provides insights into their inner workings, and the recent observations of M87* by the Event Horizon Telescope1–6 using very-long-baseline interferometry techniques allows us to investigate the dynamical processes of M87* on timescales of days. Compared with most radio interferometers, very-long-baseline interferometry networks typically have fewer antennas and low signal-to-noise ratios. Furthermore, the source is variable, prohibiting integration over time to improve signal-to-noise ratio. Here, we present an imaging algorithm7,8 that copes with the data scarcity and temporal evolution, while providing an uncertainty quantification. Our algorithm views the imaging task as a Bayesian inference problem of a time-varying brightness, exploits the correlation structure in time and reconstructs (2 + 1 + 1)-dimensional time-variable and spectrally resolved images. We apply this method to the Event Horizon Telescope observations of M87*9 and validate our approach on synthetic data. The time- and frequency-resolved reconstruction of M87* confirms variable structures on the emission ring and indicates extended and time-variable emission structures outside the ring itself.

Deep learning Bayesian inference for low-luminosity active galactic nuclei spectra

ABSTRACT Most active supermassive black holes in present-day galaxies are underfed and consist of low-luminosity active galactic nuclei (LLAGN). LLAGNs display complex multiwavelength broadband spectral energy distributions (SED), dominated by non-thermal processes which are explained to first order by a radiatively inefficient accretion flow (RIAF) and a relativistic jet. Due to the computational cost of generating such SEDs, it has not been hitherto possible to perform statistical fits to observed broadband SEDs, since such procedures require generating many thousands of models on-the-fly. Here, we have used a deep learning (DL) method to interpolate a large grid consisting of dozens of thousands of model SEDs for RIAFs and jets covering the parameter space appropriate for LLAGNs. Not only the DL method computes accurate models, it does so hundreds of thousands of times faster than solving the underlying dynamical and radiative transfer equations. This brings RIAF and jet models to the realm of Bayesian inference. We demonstrate that the combination of a DL interpolator and a Markov chain Monte Carlo ensemble sampler can recover the ground truth parameters of Mock LLAGN data. We apply our model to existing radio-to-X-rays observations of three LLAGNs: M87, NGC 315, and NGC 4261. We demonstrate that our model can estimate the relevant parameters of these accreting black holes such as the mass accretion and outflow rate at a small fraction of the computational cost of previous approaches.

X-ray bubbles in the circumgalactic medium of TNG50 Milky Way- and M31-like galaxies: signposts of supermassive black hole activity

ABSTRACT The TNG50 cosmological simulation produces X-ray emitting bubbles, shells, and cavities in the circumgalactic gas above and below the stellar discs of Milky Way- and Andromeda-like galaxies with morphological features reminiscent of the eROSITA and Fermi bubbles in the Galaxy. Two-thirds of the 198 MW/M31 analogues inspected in TNG50 at z = 0 show one or more large-scale, coherent features of overpressurized gas that impinge into the gaseous halo. Some of the galaxies include a succession of bubbles or shells of increasing size, ranging from a few to many tens of kpc. These are prominent in gas pressure, X-ray emission, and gas temperature, and often exhibit sharp boundaries with typical shock Mach numbers of 2–4. The gas in the bubbles outflows with maximum (95th pctl) radial velocities of ∼100–1500 km s−1. TNG50 bubbles expand with speeds as high as 1000–2000 km s−1 (about 1–2 kpc Myr−1), but with a great diversity and with larger bubbles expanding at slower speeds. The bubble gas is at 106.4−7.2 K temperatures and is enriched to metallicities of $0.5-2~\rm Z_{\odot }$. In TNG50, the bubbles are a manifestation of episodic, kinetic, wind-like energy injections from the supermassive black holes at the galaxy centres that accrete at low Eddington ratios. According to TNG50, X-ray, and possibly γ-ray, bubbles similar to those observed in the Milky Way should be a frequent feature of disc-like galaxies prior to, or on the verge of, being quenched. They should be within the grasp of eROSITA in the local Universe.