How Quasar Feedback May Shape the Co-evolutionary Paths

Ishibashi, Wako

doi:10.3389/fspas.2017.00026

PERSPECTIVE article

Front. Astron. Space Sci., 17 October 2017
Sec. Extragalactic Astronomy
Volume 4 - 2017 | https://doi.org/10.3389/fspas.2017.00026

How Quasar Feedback May Shape the Co-evolutionary Paths

Wako Ishibashi^*

Physik-Institut, University of Zurich, Zürich, Switzerland

Observations point toward some form of “co-evolutionary sequence,” from dust-enshrouded starbursts to luminous unobscured quasars. Active galactic nucleus (AGN) feedback is generally invoked to expel the obscuring dusty gas in a blow-out event, eventually revealing the hidden central quasar. However, the physical mechanism driving AGN feedback, either due to winds or radiation, remains uncertain and is still a source of much debate. We consider quasar feedback, based on radiation pressure on dust, which directly acts on the obscuring dusty gas. We show that AGN radiative feedback is capable of efficiently removing the obscuring cocoon, and driving powerful outflows on galactic scales, consistent with recent observations. I will discuss how such quasar feedback may provide a natural physical interpretation of the observed evolutionary path, and the physical implications in the broader context of black hole-host galaxy co-evolution.

1. Introduction

Active galactic nuclei (AGN) and nuclear starbursts are observed to be intimately coupled, both fueled by the rapid infall of matter, possibly triggered by galaxy mergers (e.g., Alexander and Hickox, 2012, and references therein). A “co-evolutionary sequence,” starting from dust-obscured starbursts (e.g., ultraluminous infrared galaxies or ULIRG) leading to unobscured optical quasars, has been widely discussed since Sanders et al. (1988). In most evolutionary scenarios, some form of AGN feedback is required to expel the obscuring gas and dust in a short-lived “blow-out” event.

Numerical simulations try to reproduce such an evolutionary sequence, with AGN feedback usually implemented by coupling a fixed fraction of the accretion luminosity to the surrounding medium (Di Matteo et al., 2005; Hopkins et al., 2005). However, the physical mechanism for AGN feedback is not specified, and in particular no explicit connection is made with the dust component (which is the fundamental parameter determining both the source visibility and the spectral energy distribution). On the observational side, there has been renewed interest in the topic, following the discovery of new populations of dust-reddened quasars at high redshifts (close to the peak epoch of both AGN and starburst activities), likely caught in the short-lived blow-out phase (e.g., Banerji et al., 2012, 2015).

Here we propose a physical mechanism that directly makes use of the dust in order to drive AGN feedback, and we briefly discuss how such AGN radiative feedback may naturally explain the observed co-evolutionary path (Ishibashi and Fabian, 2016b).

2. AGN Radiative Feedback: Radiation Pressure on Dust

We consider AGN feedback driven by radiation pressure on dust. The ambient dusty gas is swept up into an outflowing shell, with the corresponding equation of motion given by:

\begin{array}{l} \frac{d}{d t} [M_{s h} (r) v] = \frac{L}{c} (1 + τ_{I R} - e^{- τ_{U V}}) - \frac{G M (r) M_{s h} (r)}{r^{2}} & (1) \end{array}

where L is the central luminosity, M(r) is the total mass distribution, M_sh(r) is the shell mass, and $τ_{I R, U V} (r) = \frac{κ_{I R, U V} M_{s h} (r)}{4 π r^{2}}$ are the infrared (IR) and ultraviolet (UV) optical depths (Ishibashi and Fabian, 2015; Thompson et al., 2015). We recall that there are 3 distinct physical regimes depending on the optical depth of the medium: optically thick to both IR and UV, optically thick to UV but optically thin to IR (single scattering limit), and optically thin to UV. The associated shell column density is defined as:

\begin{array}{l} N_{s h} (r) = \frac{M_{s h} (r)}{4 π r^{2} m_{p}} & (2) \end{array}

where m_p is the proton mass.

By equating the outward force due to radiation pressure to the inward force due to gravity, we can define a critical luminosity, $L_{E}^{'}$ , which may be considered as a generalized form of the effective Eddington luminosity. We recall that the effective Eddington luminosity for dusty gas is distinct from the standard Eddington luminosity relevant for ionized gas (Fabian et al., 2006). The corresponding effective Eddington ratio is then given by:

\begin{array}{l} Γ = \frac{L}{L_{E}^{'}} = \frac{L r^{2}}{c G M (r) M_{s h} (r)} (1 + τ_{I R} - e^{- τ_{U V}}) & (3) \end{array}

The effective Eddington ratios in the 3 optical depth regimes are respectively given by:

\begin{array}{l} Γ_{I R} = \frac{κ_{I R} L}{4 π G c M (r)}; Γ_{S S} = \frac{L}{4 π G c m_{p} M (r) N_{s h} (r)}; \\ Γ_{U V} = \frac{κ_{U V} L}{4 π G c M (r)} & (4) \end{array}

Figure 1 shows the effective Eddington ratio as a function of radius, and its dependence on the underlying physical parameters. The effective Eddington ratio basically corresponds to the ratio of the radiative force to the gravitational force, and must exceed unity for an outflowing shell (the standard Eddington limit does not apply here). From Figure 1, we see that an increase in the shell mass leads to a lower Eddington ratio in the single scattering regime (cyan curve); while in the IR-optically thick and UV-optically thin regimes, Γ_IR and Γ_UV are independent of the shell mass (cf. Equation 4), hence the cyan and black curves overlap at small and large radii. In contrast, enhanced opacities (e.g., due to a higher dust-to-gas ratio) lead to higher Eddington ratios in the IR-optically thick and UV-optically thin regimes (magenta curve); whereas in the single scattering regime, Γ_SS is independent of the opacity (cf. Equation 4), and thus the magenta and black curves overlap at intermediate radii. We note that there are two interesting trends in the IR-optically thick regime:

1. The effective Eddington ratio is independent of the shell column density (or shell mass).

2. The effective Eddington ratio directly scales with the opacity, and hence dust-to-gas ratio (as κ_IR ∝ f_dg).

FIGURE 1

Figure 1. Effective Eddington ratio as a function of radius: L = 5 × 10⁴⁶erg/s, $M_{s h} = 5 \times 1 0^{8} M_{⊙}$ , f_dg = 1/150 (black, fiducial); L = 5 × 10⁴⁶erg/s, $M_{s h} = 1 \times 1 0^{9} M_{⊙}$ , f_dg = 1/150 (cyan); L = 5 × 10⁴⁶erg/s, $M_{s h} = 5 \times 1 0^{8} M_{⊙}$ , f_dg = 1/50 (magenta).

The first point implies that even dense material can potentially be disrupted; while the second point implies that the more dusty gas is more easily ejected (Ishibashi and Fabian, 2016b).

3. Dust Obscuration and Blow-Out

In physical terms, dusty gas surrounding the central source absorbs UV radiation and re-emits in the IR band. If the reprocessed IR photons remain trapped in the nuclear regions, the system is effectively in the IR-optically thick regime, where the Eddington ratio is independent of the column density (Equation 4). Such conditions are likely reached in the nuclear regions of ULIRG-like systems, characterized by high densities and high dust-to-gas fractions. This leads to both higher IR optical depth and higher IR-Eddington ratio, which combine to facilitate the blow-out of dusty gas.

In Figure 2 (left panel) we plot the temporal evolution of the shell column density for different values of the central luminosity. We see that the obscuration falls off with time, with the decline being more rapid for brighter sources. In fact, an increase in the luminosity leads to an increase in the effective Eddington ratio in all 3 optical depth regimes, and thus efficient acceleration, and resulting high velocity. The right panel in Figure 2 shows the shell column density as a function of velocity for different values of the dust-to-gas ratio. We observe that a given column is accelerated to higher velocities for larger dust-to-gas fractions: at a given time, the column density is lower for higher dust-to-gas ratios. Indeed, an increase in the dust-to-gas ratio leads to enhanced IR opacity, which in turn leads to a higher Eddington ratio in the IR-optically thick regime. Thus the more dusty gas (which provides much of the obscuration) is preferentially ejected by radiative feedback, and therefore the central quasar has a natural tendency to remove its own obscuring cocoon and reveal itself. We have previously discussed how such AGN-driven dusty outflows may propagate on larger scales beyond the host galaxy, and contribute to the enrichment of the circum-galactic medium (Ishibashi and Fabian, 2016a).

FIGURE 2

Figure 2. Left: shell column density vs. time for variations in luminosity ( $N_{s h, 0} = 1 0^{24} {cm}^{- 2}$ , f_dg = 1/150): L = 3 × 10⁴⁶erg/s (blue), L = 5 × 10⁴⁶erg/s (black), L = 1 × 10⁴⁷erg/s (green). Right: shell column density vs. velocity for variations in the dust-to-gas ratio (L = 2 × 10⁴⁶erg/s, $N_{s h, 0} = 1 0^{24} {cm}^{- 2}$ ): f_dg = 1/150 (black), f_dg = 1/50 (blue), f_dg = 1/30 (green). High-luminosity case: L = 10⁴⁷erg/s, f_dg = 1/50 (magenta).

The model results can be compared with observational samples of the recently uncovered populations of dust-reddened quasars and dust-obscured galaxies (Assef et al., 2015; Banerji et al., 2015; Zakamska et al., 2016). These sources are likely observed close to the point of ejecting their obscuring cocoon, and thus may represent the long-sought “sources in transition” (as discussed in Ishibashi and Fabian, 2016b).

4. The Dual Role of Dust: Hiding and Revealing the Central Quasar

The large amounts of gas and dust surrounding the nuclear source provide potential fuel for the central black hole, but that same accreting material is also responsible for significant obscuration. Although AGN feedback is generally invoked to clear the obscuring gas and dust in co-evolutionary scenarios, the actual physical mechanism remains uncertain, with no clear connection to the dust component. In our picture, AGN radiative feedback directly acts on the dusty gas, which forms the obscuring medium. Large amounts of dust imply heavy obscuration, but also powerful feedback, suggesting a causal link between dust obscuration and blow-out. In this scenario, dust plays a dual role: it is responsible for both hiding, but eventually revealing the central quasar.

An important requirement in our model is the presence of dust, since the whole AGN feedback process relies on radiation pressure on dust. The large amount of dust required in order to sustain feedback may be provided by supernovae in the starburst phase. In fact, recent observations indicate that significant quantities of dust can be produced in core-collapse supernovae (Owen and Barlow, 2015; Wesson et al., 2015). Continued star formation in the starburst is beneficial, as it helps keeping the ambient medium dusty, and thus supporting the overall AGN feedback process. Stellar feedback likely disrupts the starburst activity, but the nuclear gas reservoir may also be replenished by re-accretion episodes. Furthermore, we previously suggested that star formation (with subsequent supernova explosions and associated dust release) may be triggered within the AGN feedback-driven outflow itself (Ishibashi and Fabian, 2012). Interestingly, the first direct detection of star formation occurring inside a galactic outflow has now been observationally confirmed (Maiolino et al., 2017), nicely supporting our model predictions. In our framework, AGN feedback and starburst phenomena are intrinsically coupled through the production of dust in supernova explosions, leading to a natural interpretation of the co-evolutionary path.

Author Contributions

The author confirms being the sole contributor of this work and approved it for publication.

Conflict of Interest Statement

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer GR and handling Editor declared their shared affiliation.

References

Alexander, D. M., and Hickox, R. C. (2012). What drives the growth of black holes?. New Astron. Rev. 56, 93–121. doi: 10.1016/j.newar.2011.11.003

PERSPECTIVE article

How Quasar Feedback May Shape the Co-evolutionary Paths

1. Introduction

2. AGN Radiative Feedback: Radiation Pressure on Dust

3. Dust Obscuration and Blow-Out

4. The Dual Role of Dust: Hiding and Revealing the Central Quasar

Author Contributions

Conflict of Interest Statement

References

This article is part of the Research Topic

People also looked at