Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GALAXY NUCLEI”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Long-term follow-up observations of extreme coronal line emitting galaxies

We present new spectroscopic and photometric follow-up observations of the known sample of extreme coronal line-emitting galaxies (ECLEs) identified in the Sloan Digital Sky Survey (SDSS). With these new data, observations of the ECLE sample now span a period of two decades following their initial SDSS detections. We confirm the non-recurrence of the iron coronal line signatures in five of the seven objects, further supporting their identification as the transient light echoes of tidal disruption events (TDEs). Photometric observations of these objects in optical bands show little overall evolution. In contrast, mid-infrared (MIR) observations show ongoing long-term declines consistent with power-law decay. The remaining two objects had been classified as active galactic nuclei (AGNs) with unusually strong coronal lines rather than being TDE related, given the persistence of the coronal lines in earlier follow-up spectra. We confirm this classification, with our spectra continuing to show the presence of strong, unchanged coronal line features and AGN-like MIR colors and behavior. We have constructed spectral templates of both subtypes of ECLE to aid in distinguishing the likely origin of newly discovered ECLEs. We highlight the need for higher cadence, and more rapid, follow-up observations of such objects to better constrain their properties and evolution. We also discuss the relationships between ECLEs, TDEs, and other identified transients having significant MIR variability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Spectroscopic Classification of Astronomical Transients (SCAT) Survey: Overview, Pipeline Description, Initial Results, and Future Plans

Here, we present the Spectroscopic Classification of Astronomical Transients (SCAT) survey, which is dedicated to spectrophotometric observations of transient objects such as supernovae and tidal disruption events. SCAT uses the SuperNova Integral-Field Spectrograph (SNIFS) on the University of Hawai’i 2.2 m (UH2.2m) telescope. SNIFS was designed specifically for accurate transient spectrophotometry, including absolute flux calibration and host-galaxy removal. We describe the data reduction and calibration pipeline including spectral extraction, telluric correction, atmospheric characterization, nightly photometricity, and spectrophotometric precision. We achieve ≲5% spectrophotometry across the full optical wavelength range (3500–9000 Å) under photometric conditions. The inclusion of photometry from the SNIFS multi-filter mosaic imager allows for decent spectrophotometric calibration (10%–20%) even under unfavorable weather/atmospheric conditions. SCAT obtained ≈640 spectra of transients over the first 3 yr of operations, including supernovae of all types, active galactic nuclei, cataclysmic variables, and rare transients such as superluminous supernovae and tidal disruption events. These observations will provide the community with benchmark spectrophotometry to constrain the next generation of hydrodynamic and radiative transfer models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The LSST AGN Data Challenge: Selection Methods

Abstract Development of the Rubin Observatory Legacy Survey of Space and Time (LSST) includes a series of Data Challenges (DCs) arranged by various LSST Scientific Collaborations that are taking place during the project's preoperational phase. The AGN Science Collaboration Data Challenge (AGNSC-DC) is a partial prototype of the expected LSST data on active galactic nuclei (AGNs), aimed at validating machine learning approaches for AGN selection and characterization in large surveys like LSST. The AGNSC-DC took place in 2021, focusing on accuracy, robustness, and scalability. The training and the blinded data sets were constructed to mimic the future LSST release catalogs using the data from the Sloan Digital Sky Survey Stripe 82 region and the XMM-Newton Large Scale Structure Survey region. Data features were divided into astrometry, photometry, color, morphology, redshift, and class label with the addition of variability features and images. We present the results of four submitted solutions to DCs using both classical and machine learning methods. We systematically test the performance of supervised models (support vector machine, random forest, extreme gradient boosting, artificial neural network, convolutional neural network) and unsupervised ones (deep embedding clustering) when applied to the problem of classifying/clustering sources as stars, galaxies, or AGNs. We obtained classification accuracy of 97.5% for supervised models and clustering accuracy of 96.0% for unsupervised ones and 95.0% with a classic approach for a blinded data set. We find that variability features significantly improve the accuracy of the trained models, and correlation analysis among different bands enables a fast and inexpensive first-order selection of quasar candidates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Sensitivity of the Cherenkov Telescope Array to TeV photon emission from the Large Magellanic Cloud

ABSTRACT A deep survey of the Large Magellanic Cloud at ∼0.1–100 TeV photon energies with the Cherenkov Telescope Array is planned. We assess the detection prospects based on a model for the emission of the galaxy, comprising the four known TeV emitters, mock populations of sources, and interstellar emission on galactic scales. We also assess the detectability of 30 Doradus and SN 1987A, and the constraints that can be derived on the nature of dark matter. The survey will allow for fine spectral studies of N 157B, N 132D, LMC P3, and 30 Doradus C, and half a dozen other sources should be revealed, mainly pulsar-powered objects. The remnant from SN 1987A could be detected if it produces cosmic-ray nuclei with a flat power-law spectrum at high energies, or with a steeper index 2.3–2.4 pending a flux increase by a factor of >3–4 over ∼2015–2035. Large-scale interstellar emission remains mostly out of reach of the survey if its >10 GeV spectrum has a soft photon index ∼2.7, but degree-scale 0.1–10 TeV pion-decay emission could be detected if the cosmic-ray spectrum hardens above >100 GeV. The 30 Doradus star-forming region is detectable if acceleration efficiency is on the order of 1−10 per cent of the mechanical luminosity and diffusion is suppressed by two orders of magnitude within <100 pc. Finally, the survey could probe the canonical velocity-averaged cross-section for self-annihilation of weakly interacting massive particles for cuspy Navarro–Frenk–White profiles.

Astronomy & Astrophysics↗
DOE OSTI2023

A multifrequency characterization of the extragalactic hard X-ray sky

Context: Nowadays, we know that the origin of the cosmic X-ray background (CXB) is mainly due to the integrated emission of active galactic nuclei. Therefore, in order to obtain a precise estimate of the contribution of different source classes to the CXB, it is crucial to achieve full characterization of the hard-X ray sky. Aims: We present a multifrequency analysis of all sources listed in the third release of the Palermo Swift-BAT hard X-ray catalog (3PBC) with the goal of (i) identifying and classifying the largest number of sources adopting multifrequency criteria, with particular emphasis on extragalactic populations and (ii) extracting sources belonging to the class of Seyfert galaxies to present here the release of the second version of the Turin-SyCAT. Methods: We outline a classification scheme based on radio, infrared (IR), and optical criteria that allows us to distinguish between unidentified and unclassified hard X-ray sources, as well as to classify those sources belonging to the Galactic and the extragalactic populations. Results: Our revised version of the 3PBC lists 1176 classified, 820 extragalactic, and 356 Galactic sources, as well as 199 unclassified and 218 unidentified sources. According to our analysis, the hard X-ray sky is mainly populated by Seyfert galaxies and blazars. For the blazar population, we report trends between the hard X-ray and the gamma-ray emissions based on the fact that a large fraction of them also have a counterpart detected by the Fermi satellite. These trends are all in agreement with the expectations of inverse Compton models which are widely adopted to explain the blazar broadband emission. For the Seyfert galaxies, we present the second version of the Turin-SyCAT, including a total of 633 Seyfert galaxies, with 282 new sources corresponding to an increase of ~80% with respect to the previous release. Comparing the hard X-ray and the infrared emissions of Seyfert galaxies, we confirm that there is no clear difference between the flux distribution of the infrared-to-hard X-ray flux ratio of Seyfert galaxies Type 1 and Type 2. However, there is a significant trend between the mid-IR flux and hard X-ray flux, confirming previous statistical results in the literature. Conclusions: We provide two catalog tables. The first is the revised version of the 3PBC catalog based on our multifrequency analyses. The second catalog table is a release of the second version of the Turin-SyCAT catalog. Finally, we highlight that extensive soft X-ray data are already available in the form of the Swift archive which can be used to search for potential counterparts of unidentified hard X-ray sources. All these datasets will be reduced and analyzed in a forthcoming analysis to determine the precise position of low-energy counterparts in the 0.5–10 keV energy range for 3PBC sources that can be targets of future optical spectroscopic campaigns; this is necessary to obtain their precise classification.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The radio loudness of SDSS quasars from the LOFAR Two-metre Sky Survey: ubiquitous jet activity and constraints on star formation

ABSTRACT We examine the distribution of radio emission from ∼42 000 quasars from the Sloan Digital Sky Survey, as measured in the LOFAR Two-metre Sky Survey (LoTSS). We present a model of the radio luminosity distribution of the quasars that assumes that every quasar displays a superposition of two sources of radio emission: active galactic nuclei (jets) and star formation. Our two-component model provides an excellent match to the observed radio flux density distributions across a wide range of redshifts and quasar optical luminosities; this suggests that the jet-launching mechanism operates in all quasars but with different powering efficiency. The wide distribution of jet powers allows for a smooth transition between the ‘radio-quiet’ and ‘radio-loud’ quasar regimes, without need for any explicit bimodality. The best-fitting model parameters indicate that the star formation rate of quasar host galaxies correlates strongly with quasar luminosity and also increases with redshift at least out to z ∼ 2. For a model where star formation rate scales as $L_{\rm bol}^{\alpha } (1+z)^{\beta }$, we find α = 0.47 ± 0.01 and β = 1.61 ± 0.05, in agreement with far-infrared studies. Quasars contribute ≈0.15 per cent of the cosmic star formation rate density at z = 0.5, rising to 0.4 per cent by z ∼ 2. The typical radio jet power is seen to increase with both increasing optical luminosity and black hole mass independently, but does not vary with redshift, suggesting intrinsic properties govern the production of the radio jets. We discuss the implications of these results for the triggering of quasar activity and the launching of jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dwarf AGNs from Optical Variability for the Origins of Seeds (DAVOS): insights from the dark energy survey deep fields

ABSTRACT We present a sample of 706, z < 1.5 active galactic nuclei (AGNs) selected from optical photometric variability in three of the Dark Energy Survey (DES) deep fields (E2, C3, and X3) over an area of 4.64 deg2. We construct light curves using difference imaging aperture photometry for resolved sources and non-difference imaging PSF photometry for unresolved sources, respectively, and characterize the variability significance. Our DES light curves have a mean cadence of 7 d, a 6-yr baseline, and a single-epoch imaging depth of up to g ∼ 24.5. Using spectral energy distribution (SED) fitting, we find 26 out of total 706 variable galaxies are consistent with dwarf galaxies with a reliable stellar mass estimate ($M_{\ast }\lt 10^{9.5}\, {\rm M}_\odot$; median photometric redshift of 0.9). We were able to constrain rapid characteristic variability time-scales (∼ weeks) using the DES light curves in 15 dwarf AGN candidates (a subset of our variable AGN candidates) at a median photometric redshift of 0.4. This rapid variability is consistent with their low black hole (BH) masses. We confirm the low-mass AGN nature of one source with a high S/N optical spectrum. We publish our catalogue, optical light curves, and supplementary data, such as X-ray properties and optical spectra, when available. We measure a variable AGN fraction versus stellar mass and compare to results from a forward model. This work demonstrates the feasibility of optical variability to identify AGNs with lower BH masses in deep fields, which may be more ‘pristine’ analogues of supermassive BH seeds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Computational general relativistic force-free electrodynamics: I. Multi-coordinate implementation and testing

General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole (BH) shadows from the Galactic Center and the M87 galaxy provide a first direct glimpse into the physics of accretion flows in the most extreme environments of the universe. The efficient extraction of energy in the form of collimated outflows or jets from a rotating BH is directly linked to the topology of the surrounding magnetic field. We aim at providing a tool to numerically model the dynamics of such fields in magnetospheres around compact objects, such as BHs and neutron stars. To do so, we probe their role in the formation of high energy phenomena such as magnetar flares and the highly variable teraelectronvolt emission of some active galactic nuclei. In this work, we present numerical strategies capable of modeling fully dynamical force-free magnetospheres of compact astrophysical objects. Here, we provide implementation details and extensive testing of our implementation of general relativistic force-free electrodynamics in Cartesian and spherical coordinates using the infrastructure of the E INSTEIN T OOLKIT . The employed hyperbolic/parabolic cleaning of numerical errors with full general relativistic compatibility allows for fast advection of numerical errors in dynamical spacetimes. Such fast advection of divergence errors significantly improves the stability of the general relativistic force-free electrodynamics modeling of BH magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Optical and Infrared Are Connected

Galaxies are often modeled as composites of separable components with distinct spectral signatures, implying that different wavelength ranges are only weakly correlated. They are not. We present a data-driven model that exploits subtle correlations between physical processes to accurately predict infrared (IR) Wide-field Infrared Survey Explorer (WISE) photometry from a neural summary of optical Sloan Digital Sky Survey spectra. The model achieves accuracies of $χ^{2}_{N} ≈ 1$ for all photometric bands in WISE, as well as good colors. We are able to tightly constrain typically IR-derived properties, e.g., the bolometric luminosities of active galactic nuclei (AGN) and dust parameters such as q PAH . We also test whether current spectral energy distribution (SED) fitting methods reproduce such panchromatic relations, but find their predictions biased and overconfident, likely due to model misspecification, with correlated biases in star-formation rates (SFRs) and AGN luminosities being most evident. To help improve SED models, we determine which features of the optical spectrum are responsible for our improved predictions, and identify several lines (Ca II , Sr II , Fe I , [O II ], and Hα), which point to the complex chronology of star formation and chemical enrichment being incorrectly modeled.

Jespersen, Christian Kragh [Princeton Univ., NJ (U↗
DOE OSTI2022

Multiple giant eruptions and X-ray emission in the recoiling AGN/LBV candidate SDSS1133

We present a comprehensive analysis of 20 yr worth of multicolour photometric light curves, multiepoch optical spectra, and X-ray data of an off-nuclear variable object SDSS1133 in Mrk 177 at z = 0.0079. The UV-optical light curves reveal that SDSS1133 experienced four outbursts in 2001, 2014, 2019, and 2021. The persistent UV-optical luminosity in the non-outbursting state is ~1041 erg s -1 with small-scale flux variations, and peak luminosities during the outbursts reach ~10 42 erg s -1 . The optical spectra exhibit enduring broad hydrogen Balmer P-Cygni profiles with the absorption minimum at ~-2000 km s -1 , indicating the presence of fast-moving ejecta. Chandra detected weak X-ray emission at a 0.3-10-keV luminosity of L X = 4 × 10 38 erg s -1 after the 2019 outburst. These lines of evidence suggests that SDSS1133 is an extreme luminous blue variable (LBV) star experiencing multiple giant eruptions with interactions of the ejected shell with different shells and/or circumstellar medium (CSM), and disfavours the recoiling active galactic nuclei scenario suggested in the literature. We suggest that pulsational pair-instability may provide a viable explanation for the multiple energetic eruptions in SDSS1133. If the current activity of SDSS1133 is a precursor of a supernova explosion, we may be able to observe a few additional giant eruptions and then the terminal supernova explosion or collapse to a massive black hole in future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Revisiting AGN as the Source of IceCube's Diffuse Neutrino Flux

The origin of the astrophysical neutrino flux reported by the IceCube Collaboration remains an open question. In this study, we use three years of publicly available IceCube data to search for evidence of neutrino emission from the blazars and non-blazar Active Galactic Nuclei (AGN) contained the Fermi 4LAC catalog. We find no evidence that these sources produce high-energy neutrinos, and conclude that blazars can produce no more than 15% of IceCube's observed flux. The constraint we derive on the contribution from non-blazar AGN, which are less luminous and more numerous than blazars, is significantly less restrictive, and it remains possible that this class of sources could produce the entirety of the diffuse neutrino flux observed by IceCube. We anticipate that it will become possible to definitively test such scenarios as IceCube accumulates and releases more data, and as gamma-ray catalogs of AGN become increasingly complete. We also comment on starburst and other starforming galaxies, and conclude that these sources could contribute substantially to the signal observed by IceCube, in particular at the lowest detected energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Neutrino and pair creation in reconnection-powered coronae of accreting black holes

A ubiquitous feature of accreting black hole systems is their hard X-ray emission which is thought to be produced through Comptonization of soft photons by electrons and positrons in the vicinity of the black hole, in a region with optical depth of order unity. The origin and composition of this Comptonizing region, known as the corona, is a matter open for debate. In this paper we investigate the role of relativistic protons accelerated in black-hole magnetospheric current sheets for the pair enrichment and neutrino emission of AGN coronae. Our model has two free parameters, namely the proton plasma magnetization σ p , which controls the peak energy of the neutrino spectrum, and the Eddington ratio λ X,Edd (defined as the ratio between X-ray luminosity L X and Eddington luminosity L Edd ), which controls the amount of energy transferred to secondary particles. For sources with λ X,Edd ≳ λ Edd,crit (where λ Edd,crit ∼ 10 -1 for σ p = 10 5 or ∼ 10 -2 for σ p = 10 7 ), proton-photon interactions and γγ annihilation produce enough secondary pairs to achieve Thomson optical depths τ T ∼ 0.1-10. In the opposite case of λ X,Edd ≲ λ Edd,crit , the coronal pairs cannot originate only from hadronic interactions. Additionally, we find that the neutrino luminosity scales as L 2 X /L Edd for λ X,Edd ≲ λ Edd,crit , while it is proportional to L X for higher λ X,Edd values. We apply our model to four Seyfert galaxies, including NGC 1068, and discuss our results in light of recent IceCube observations.

active galactic nuclei↗
DOE OSTI2024

3D hybrid fluid-particle jet simulations and the importance of synchrotron radiative losses

Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Aims. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. Methods. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. Results. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Chemical abundances of Seyfert 2 AGNs – I. Comparing oxygen abundances from distinct methods using SDSS

ABSTRACT We compare the oxygen abundance (O/H) of the narrow-line regions (NLRs) of Seyfert 2 AGNs obtained through strong-line methods and from direct measurements of the electron temperature (Te-method). The aim of this study is to explore the effects of the use of distinct methods on the range of metallicity and on the mass–metallicity relation of active galactic nuclei (AGNs) at low redshifts (z ≲ 0.4). We used the Sloan Digital Sky Survey (SDSS) and NASA/IPAC Extragalactic Database (NED) to selected optical (3000 < λ(Å) < 7000) emission line intensities of 463 confirmed Seyfert 2 AGNs. The oxygen abundances of the NLRs were estimated using the theoretical Storchi-Bergmann et al. calibrations, the semi-empirical N2O2 calibration, the Bayesian H ii-Chi-mistry code and the Te-method. We found that the oxygen abundance estimations via the strong-line methods differ from each other up to ∼0.8 dex, with the largest discrepancies in the low-metallicity regime ($\rm 12+\log (O/H) \: \lesssim \: 8.5$). We confirmed that the Te-method underestimates the oxygen abundance in NLRs, producing unreal subsolar values. We did not find any correlation between the stellar mass of the host galaxies and the metallicity of their AGNs. This result is independent of the method used to estimate Z.

Dors, O. L.↗
DOE OSTI2021

A New Sample of Gamma-Ray Emitting Jetted Active Galactic Nuclei—Preliminary Results

We are compiling a new list of gamma-ray jetted active galactic nuclei (AGN), starting from the fourth catalog of point sources of the Fermi Large Area Telescope (LAT). Our aim is to prepare a list of jetted AGN with known redshifts and classifications to be used to calibrate jet power. We searched in the available literature for all the published optical spectra and multiwavelength studies useful to characterize the sources. We found new, missed, or even forgotten information leading to a substantial change in the redshift values and classification of many sources. We present here the preliminary results of this analysis and some statistics based on the gamma-ray sources with right ascension within the interval 0 h --12 h (J2000). Although flat-spectrum radio quasars and BL Lac objects are still the dominant populations, there is a significant increase in the number of other objects, such as misaligned AGN, narrow-line Seyfert 1 galaxies, and Seyfert galaxies. We also introduced two new classes of objects: changing-look AGN and ambiguous sources. About one third of the sources remain unclassified.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Cosmic-ray confinement in radio bubbles by micromirrors

Radio bubbles, ubiquitous features of the intracluster medium around active galactic nuclei, are known to rise buoyantly for multiple scale heights through the intracluster medium (ICM). It is an open question how the bubbles can retain their high-energy cosmic-ray content over such distances. We propose that the enhanced scattering of cosmic rays due to micromirrors generated in the ICM is a viable mechanism for confining the cosmic rays within bubbles and can qualitatively reproduce their morphology. We discuss the observational implications of such a model of cosmic-ray confinement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The contribution of magnetized galactic outflows to extragalactic Faraday rotation

ABSTRACT Galactic outflows driven by star formation and active galactic nuclei blow bubbles into their local environments, causing galactic magnetic fields to be carried into intergalactic space. We explore the redshift-dependent effect of these magnetized bubbles on the Faraday rotation measure (RM) of extragalactic radio sources. Using the IllustrisTNG cosmological simulations, we separate the contribution from magnetic bubbles from that of the volume-filling magnetic component expected to be due to the seed field originating in the early universe. We use this separation to extract the redshift dependence of each component and to compare IllustrisTNG model predictions with observation measurements of the NRAO VLA sky survey (NVSS). We find that magnetized bubbles provide a sizeable contribution to the extragalactic RM, with redshift-independent 〈|RM|〉 ≃ 13 rad/m2 for sources at redshifts z ≥ 2. This is close to the mean residual RM of 16 rad/m2 found from NVSS data in this redshift range. Using the IllustrisTNG simulations, we also evaluate a simple model for the contribution to residual RM from individual host galaxies and show that this contribution is negligible at high-redshift. While the contribution from magnetic bubbles in the IllustrisTNG model is currently compatible with observational measurements of residual RM, the next-generation RM sky surveys, which will be free from the wrapping uncertainty, have larger statistics and better sensitivity should be able to observe predicted flat contribution from magnetic bubbles at large redshifts. This should allow to experimentally probe magnetic bubbles and check models of galaxy feedback in cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Jets with a Twist: The Emergence of FR0 Jets in a 3D GRMHD Simulation of Zero-angular-momentum Black Hole Accretion

Abstract Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius, R B = 10 3 R g . The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outside R B and suppress BH accretion to 1.5% of the Bondi rate, M ̇ B . Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% of M ̇ B reach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2% M ̇ B . Because the BAD and RAD states tangle up the jets and destroy them well inside R B , they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies.

Astronomy & Astrophysics↗