We present the first simultaneous spectral energy distribution (SED) of M87 core at a scale of 0.4 arcsec ( ∼ 32 pc) across the electromagnetic spectrum. Two separate, quiescent, and active states ...are sampled that are characterized by a similar featureless SED of power-law form, and that are thus remarkably different from that of a canonical active galactic nuclei or a radiatively inefficient accretion source. We show that the emission from a jet gives an excellent representation of the core of M87 core covering ten orders of magnitude in frequency for both the active and the quiescent phases. The inferred total jet power is, however, one to two orders of magnitude lower than the jet mechanical power reported in the literature. The maximum luminosity of a thin accretion disc allowed by the data yields an accretion rate of < 6 × 10− 5 M⊙ yr− 1, assuming 10 per cent efficiency. This power suffices to explain M87 radiative luminosity at the jet frame, it is however two to three order of magnitude below that required to account for the jet's kinetic power. The simplest explanation is variability, which requires the core power of M87 to have been two to three orders of magnitude higher in the last 200 yr. Alternatively, an extra source of power may derive from black hole spin. Based on the strict upper limit on the accretion rate, such spin power extraction requires an efficiency an order of magnitude higher than predicted from magnetohydrodynamic simulations, currently in the few hundred per cent range.
Abstract
We present new calculations of the Clumpy AGN Tori in a 3D geometry (CAT3D) clumpy torus models, which now include a more physical dust sublimation model as well as active galactic nucleus ...(AGN) anisotropic emission. These new models allow graphite grains to persist at temperatures higher than the silicate dust sublimation temperature. This produces stronger near-infrared emission and bluer mid-infrared (MIR) spectral slopes. We make a statistical comparison of the CAT3D model MIR predictions with a compilation of sub-arcsecond resolution ground-based MIR spectroscopy of 52 nearby Seyfert galaxies (median distance of 36 Mpc) and 10 quasars. We focus on the AGN MIR spectral index αMIR and the strength of the 9.7 μm silicate feature S
Sil. As with other clumpy torus models, the new CAT3D models do not reproduce the Seyfert galaxies with deep silicate absorption (S
Sil < −1). Excluding those, we conclude that the new CAT3D models are in better agreement with the observed αMIR and S
Sil of Seyfert galaxies and quasars. We find that Seyfert 2 are reproduced with models with low photon escape probabilities, while the quasars and the Seyfert 1–1.5 require generally models with higher photon escape probabilities. Quasars and Seyfert 1–1.5 tend to show steeper radial cloud distributions and fewer clouds along an equatorial line of sight than Seyfert 2. Introducing AGN anisotropic emission besides the more physical dust sublimation models alleviates the problem of requiring inverted radial cloud distributions (i.e. more clouds towards the outer parts of the torus) to explain the MIR spectral indices of type 2 Seyferts.
We study the multiphase feedback processes in the central ∼3 kpc of the barred Seyfert 2 galaxy NGC 5643. We used observations of the cold molecular gas (ALMA CO(2−1) transition) and ionized gas ...(MUSE IFU optical emission lines). We studied different regions along the outflow zone, which extends out to ∼2.3 kpc in the same direction (east-west) as the radio jet, as well as nuclear and circumnuclear regions in the host galaxy disk. The CO(2−1) line profiles of regions in the outflow and spiral arms show two or more different velocity components: one associated with the host galaxy rotation, and the others with out- or inflowing material. In the outflow region, the O
III
λ
5007 Å emission lines have two or more components: the narrow component traces rotation of the gas in the disk, and the others are related to the ionized outflow. The deprojected outflowing velocities of the cold molecular gas (median
V
central
∼ 189 km s
−1
) are generally lower than those of the outflowing ionized gas, which reach deprojected velocities of up to 750 km s
−1
close to the active galactic nucleus (AGN), and their spatial profiles follow those of the ionized phase. This suggests that the outflowing molecular gas in the galaxy disk is being entrained by the AGN wind. We derive molecular and ionized outflow masses of ∼5.2 × 10
7
M
⊙
(
α
CO
Galactic
) and 8.5 × 10
4
M
⊙
and molecular and ionized outflow mass rates of ∼51
M
⊙
yr
−1
(
α
CO
Galactic
) and 0.14
M
⊙
yr
−1
, respectively. This means that the molecular phase dominates the outflow mass and outflow mass rate, while the kinetic power and momentum of the outflow are similar in both phases. However, the wind momentum loads (
Ṗ
out
/
Ṗ
AGN
) for the molecular and ionized outflow phases are ∼27−5 (
α
CO
Galactic
and
α
CO
ULIRGs
) and < 1, which suggests that the molecular phase is not momentum conserving, but the ionized phase most certainly is. The molecular gas content (
M
east
∼ 1.5 × 10
7
M
⊙
;
α
CO
Galactic
) of the eastern spiral arm is approximately 50−70% of the content of the western one. We interpret this as destruction or clearing of the molecular gas produced by the AGN wind impacting in the eastern side of the host galaxy (negative feedback process). The increase in molecular phase momentum implies that part of the kinetic energy from the AGN wind is transmitted to the molecular outflow. This suggests that in Seyfert-like AGN such as NGC 5643, the radiative or quasar and the kinetic or radio AGN feedback modes coexist and may shape the host galaxies even at kiloparsec scales through both positive and (mild) negative feedback.
We present J′ and K′ imaging linear polarimetric adaptive optics observations of NGC 1068 using MMT-Pol on the 6.5-m MMT. These observations allow us to study the torus from a magnetohydrodynamical ...(MHD) framework. In a 0.5 arcsec (30 pc) aperture at K′, we find that polarization arising from the passage of radiation from the inner edge of the torus through magnetically aligned dust grains in the clumps is the dominant polarization mechanism, with an intrinsic polarization of 7.0 ± 2.2 per cent. This result yields a torus magnetic field strength in the range of 4–82 mG through paramagnetic alignment, and 139
$^{+11}_{-20}$
mG through the Chandrasekhar–Fermi method. The measured position angle (P.A.) of polarization at K′ is found to be similar to the P.A. of the obscuring dusty component at few parsec scales using infrared interferometric techniques. We show that the constant component of the magnetic field is responsible for the alignment of the dust grains, and aligned with the torus axis on to the plane of the sky. Adopting this magnetic field configuration and the physical conditions of the clumps in the MHD outflow wind model, we estimate a mass outflow rate ≤0.17 M⊙ yr−1 at 0.4 pc from the central engine for those clumps showing near-infrared dichroism. The models used were able to create the torus in a time-scale of ≥105 yr with a rotational velocity of ≤1228 km s−1 at 0.4 pc. We conclude that the evolution, morphology and kinematics of the torus in NGC 1068 can be explained within a MHD framework.
ABSTRACT
We present a polarization variability analysis of a sample of 26 γ-ray blazars monitored by the Steward Observatory between 2008 and 2018 in the optical band. We investigate the properties ...and long-term variability of their optical polarization, searching for differences between blazar types. We observe that BL Lac objects are typically less polarized and less variable than flat spectrum radio quasars (FSRQs). Moreover, BL Lacs display a distribution of their polarization angle typically oriented in a preferential direction, contrary to the rather random distribution of FSRQs. For the latter blazar type, as well as those sources showing a bright stellar emission, we take into account the depolarizing effect introduced by the broad line region and the host galaxy on the measured polarization degree. In this sample, we also observe that BL Lacs present an uncorrelated evolution of the flux and the polarization. Contrary, FSRQs show a correlation before the depolarization correction, that is lost however after considering this effect. In addition, we study the behaviour of the polarization angle, searching for angle rotations in its long-term evolution. We derive that the FSRQs studied here show rotations more frequently than BL Lac objects by a factor ∼1.5. During these periods, we also observe a systematic decrease of the polarization fraction, as well as a marginal flux increase, not significant however to connect rotations with optical flares. We interpret these results within the extended shock-in-jet scenario, able to explain the overall features observed here for the polarization of the blazar sample.
ABSTRACT
Blazars optical emission is generally dominated by relativistic jets, although the host galaxy, accretion disc, and broad-line region (BLR) may also contribute significantly. Disentangling ...their contributions has been challenging for years due to the dominance of the jet. To quantify the contributions to the spectral variability, we use the statistical technique for dimensionality reduction non-negative matrix factorization on a spectroscopic data set of 26 γ-ray blazars. This technique allows to model large numbers of spectra in terms of a reduced number of components. We use a priori knowledge to obtain components associated with meaningful physical processes. The sources are classified according to their optical spectrum as host-galaxy dominated BL Lac objects (BL Lacs), BL Lacs, or flat spectrum radio quasars (FSRQs). Host-galaxy sources show less variability, as expected, and bluer-when-brighter (BWB) trends, as the other BL Lacs. For FSRQs, more complicated colour-flux behaviours are observed: redder-when-brighter for low states saturating above a certain level and, in some cases, turning to BWB. We are able to reproduce the variability observed during 10 yr using only two to four components, depending on the type. The simplest scenario corresponds to host-galaxy blazars, whose spectra are reconstructed using the stellar population and a power law (PL) for the jet. BL Lac spectra are reproduced using from two to four PLs. Different components can be associated with acceleration/cooling processes taking place in the jet. The reconstruction of FSRQs also incorporates a QSO-like component to account for the BLR, plus a very steep PL, associated with the accretion disc.
Abstract
Deep exposure imaging of early-type galaxies (ETGs) are revealing the second-order complexity of these objects, which have been long considered uniform, dispersion-supported spheroidals. ...‘Fine structure’ features (e.g. ripples, plumes, tidal tails, rings) as well as depleted stellar cores (i.e. central light deficits) characterize a number of massive ETG galaxies, and can be interpreted as the result of galaxy–galaxy interactions. We discuss how the time-scale for the evolution of cores and fine structures are comparable, and hence it is expected that they develop in parallel after the major interaction event which shaped the ETG. Using archival data, we compare the ‘depleted stellar mass’ (i.e. the mass missing from the depleted stellar core) against the prominence of the fine structure features, and observe that they correlate inversely. This result confirms our expectation that, while the supermassive black hole (SMBH) binary (constituted by the SMBHs of the merger progenitors) excavates the core via three-body interactions, the gravitational potential of the newborn galaxy relaxes, and the fine structures fade below detection levels. We expect the inverse correlation to hold at least within the first Gyr from the merger which created the SMBH binary; after then, the fine structure evolves independently.
We measure the 10 and 18 m silicate features in a sample of 67 local (z < 0.1) type 1 active galactic nuclei (AGN) with available Spitzer spectra dominated by nonstellar processes. We find that the ...10 m silicate feature peaks at with a strength (Sip = ln fp(spectrum)/fp(continuum)) of , while the 18 m one peaks at with a strength of . We select from this sample sources with the strongest 10 m silicate strength ( , 10 objects). We carry out a detailed modeling of the infrared spectrometer/Spitzer spectra by comparing several models that assume different geometries and dust composition: a smooth torus model, two clumpy torus models, a two-phase medium torus model, and a disk+outflow clumpy model. We find that the silicate features are well modeled by the clumpy model of Nenkova et al., and among all models, those including outflows and complex dust composition are the best. We note that even in AGN-dominated galaxies, it is usually necessary to add stellar contributions to reproduce the emission at the shortest wavelengths.
ABSTRACT
A complete demographic of active galactic nuclei (AGNs) is essential to understand the evolution of the Universe. Optical surveys estimate the population of AGNs in the local Universe to be ...of ∼ 4 per cent. However, these results could be biased towards bright sources, not affected by the host galaxy attenuation. An alternative method for detecting these objects is through the X-ray emission. In this work, we aim to complement the AGN population of the optical CALIFA survey (941 sources), by using X-ray data from Chandra, which provides the best spatial resolution to date, essential to isolate the nuclear emission from the host galaxy. We study a total of 138 sources with available data. We find 34 new bona fide AGNs and 23 AGN candidates, which could increase the AGN population to 7–10 per cent among the CALIFA survey. X-rays are particularly useful for low-luminosity AGNs since they are excluded by the criterion of large equivalent width of the $\rm {H\alpha }$ emission line when applied to optical selections. Indeed, placing such a restrictive criteria might cause a loss of up to 70 per cent of AGN sources. X-ray detected sources are preferentially located in the right side of the ${\mathrm{O}}\, {\small {\rm III}}$ /Hβ versus ${\mathrm{N}}\, {\small {\rm II}}$ /Hα diagram, suggesting that this diagram might be the most reliable at classifying AGN sources. Our results support the idea that multiwavelength studies are the best way to obtain a complete AGN population.