Active Galactic Nuclei Winds Research Articles

Q wind code release: a non-hydrodynamical approach to modelling line-driven winds in active galactic nuclei

ABSTRACT Ultraviolet (UV) line-driven winds may be an important part of the active galactic nucleus (AGN) feedback process, but understanding their impact is hindered by the complex nature of the radiation hydrodynamics. Instead, we have taken the approach pioneered by Risaliti & Elvis, calculating only ballistic trajectories from radiation forces and gravity but neglecting gas pressure. We have completely rewritten their Qwind code using more robust algorithms and can now quickly model the acceleration phase of these winds for any AGN spectral energy distribution spanning UV and X-ray wavebands. We demonstrate the code using an AGN with black hole mass $10^8\, \text{M}_\odot$ emitting at half the Eddington rate and show that this can effectively eject a wind with velocities ${\simeq}(0.1-0.2)\, c$. The mass loss rates can be up to ≃0.3M⊙ per year, consistent with more computationally expensive hydrodynamical simulations, though we highlight the importance of future improvements in radiation transfer along the multiple different lines of sight illuminating the wind. The code is fully public and can be used to quickly explore the conditions under which AGN feedback can be dominated by accretion disc winds.

Clumpy AGN Outflows due to Thermal Instability

Abstract One of the main mechanisms that could drive mass outflows on parsec scales in active galactic nuclei (AGN) is thermal driving. The same X-rays that ionize and heat the plasma are also expected to make it thermally unstable. Indeed, it has been proposed that the observed clumpiness in AGN winds is caused by thermal instability (TI). While many studies employing time-dependent numerical simulations of AGN outflows have included the necessary physics for TI, none have so far managed to produce clumpiness. Here we present the first such clumpy wind simulations in 1D and 2D, obtained by simulating parsec-scale outflows irradiated by an AGN. By combining an analysis of our extensive parameter survey with physical arguments, we show that the lack of clumps in previous numerical models can be attributed to the following three effects: (i) insufficient radiative heating or other physical processes that prevent the outflowing gas from entering the TI zone; (ii) the stabilizing effect of stretching (due to rapid radial acceleration) in cases where the gas enters the TI zone; and (iii) a flow speed effect: in circumstances where stretching is inefficient, the flow can still be so fast that it passes through the TI zone too quickly for perturbations to grow. In addition to these considerations, we also find that a necessary condition to trigger TI in an outflow is for the pressure ionization parameter to decrease along a streamline once gas enters a TI zone.

Metal abundances in the MACER simulations of the hot interstellar medium

AbstractHot plasma is the dominant phase of the interstellar medium of early‐type galaxies. It can originate from stellar mass losses, residual gas from the formation epoch, and accretion from outside the galaxies. Its evolution is linked to the dynamical structure of the host galaxy, to the supernova and active galactic nuclei (AGN) feedback, and to (late‐epoch) star formation in a way that has yet to be fully understood. Important clues about the origin and evolution of the hot gas come from the abundances of heavy metals that have been studied with increasing detail through XMM‐Newton and Chandra. We present recent high‐resolution hydrodynamical simulations of the hot gas evolution that include the above processes, where several chemical species, originating in asymptotic giant branch stars and supernovae of types Ia and II, have also been considered. The high resolution, of few parsecs in the central galactic region, allows us to track the metal enrichment, transportation, and dilution throughout the galaxy. The comparison of model results with observed abundances demonstrates good agreement for the region enriched by the AGN wind but also shows discrepancies for the diffuse hot gas; the latter indicates the need for a revision of standard assumptions and/or the importance of neglected effects due to the dust and/or residual uncertainties in deriving abundances from the X‐ray spectra.

A Dissection of Spatially Resolved AGN Feedback across the Electromagnetic Spectrum

Abstract We present optical SuperNova Integral Field Spectrograph integral field spectroscopy, Hubble Space Telescope optical imaging, Chandra X-ray imaging, and Very Large Array radio interferometry of the merging galaxy 2MASX J04234080+0408017, which hosts a Seyfert 2 active galactic nucleus (AGN) at z = 0.046. Our observations reveal that radiatively driven, ionized gas outflows are successful to distances >10 kpc due to the low mass of the host system, encompassing the entirety of the observed optical emission. We also find that at large radii, where observed velocities cannot be reproduced by radiative driving models, high-velocity kinematics are likely due to mechanical driving from AGN winds impacting high-density host material. This impacting deposits sufficient energy to shock the host material, producing thermal X-ray emission and cosmic rays, which in turn promote the formation of in situ radio structure in a pseudo-jet morphology along the high-density lanes.

BALQSO spectra explained by shock disruption of galactic clouds

ABSTRACT Blue-shifted broad absorption lines (BALs) detected in quasar’s spectra are indicative of active galactic nuclei (AGN) outflows. We show, using 2D hydrodynamical simulations, that disruption of interstellar clouds by a fast AGN wind can lead to formation of cold, dense high-speed blobs that give rise to broad absorption features in the transmission spectrum of the AGN continuum source. For a wind velocity of 0.1c and sufficiently high cloud density (nc ≳ 104 cm−3; depending on size), absorption troughs with velocities up to about 3000 km s−1can be produced. For slower winds and/or lower cloud density, the anticipated velocity of the absorbing clouds should be smaller.

ALMA images the many faces of the NGC 1068 torus and its surroundings

Aims. We investigate the fueling and the feedback of nuclear activity in the nearby (D = 14 Mpc) Seyfert 2 barred galaxy NGC 1068 by studying the distribution and kinematics of molecular gas in the torus and its connections to the host galaxy disk. Methods. We used the Atacama Large Millimeter Array (ALMA ) to image the emission of a set of molecular gas tracers in the circumnuclear disk (CND) and the torus of the galaxy using the CO(2–1), CO(3–2), and HCO+(4–3) lines and their underlying continuum emission with high spatial resolutions (0.03″ − 0.09″ ≃ 2 − 6 pc). These transitions, which span a wide range of physical conditions of molecular gas (n(H2)⊂103 − 107 cm−3), are instrumental in revealing the density radial stratification and the complex kinematics of the gas in the torus and its surroundings. Results. The ALMA images resolve the CND as an asymmetric ringed disk of D ≃ 400 pc in size and ≃1.4 × 108 M⊙ in mass. The CND shows a marked deficit of molecular gas in its central ≃130 pc region. The inner edge of the ring is associated with the presence of edge-brightened arcs of NIR polarized emission, which are identified with the current working surface of the ionized wind of the active galactic nucleus (AGN). ALMA proves the existence of an elongated molecular disk/torus in NGC 1068 of Mtorusgas ≃ 3 × 105 M⊙, which extends over a large range of spatial scales D ≃ 10 − 30 pc around the central engine. The new observations evidence the density radial stratification of the torus: the HCO+(4–3) torus, with a full size DHCO+(4 − 3) = 11 ± 0.6 pc, is a factor of between two and three smaller than its CO(2–1) and CO(3–2) counterparts, which have full sizes of DCO(3 − 2) = 26 ± 0.6 pc and DCO(2 − 1) = 28 ± 0.6 pc, respectively. This result brings into light the many faces of the molecular torus. The torus is connected to the CND through a network of molecular gas streamers detected inside the CND ring. The kinematics of molecular gas show strong departures from circular motions in the torus, the gas streamers, and the CND ring. These velocity field distortions are interconnected and are part of a 3D outflow that reflects the effects of AGN feedback on the kinematics of molecular gas across a wide range of spatial scales around the central engine. In particular, we estimate through modeling that a significant fraction of the gas inside the torus (≃ 0.4 − 0.6 × Mtorusgas) and a comparable amount of mass along the gas streamers are outflowing. However, the bulk of the mass, momentum, and energy of the molecular outflow of NGC 1068 is contained at larger radii in the CND region, where the AGN wind and the radio jet are currently pushing the gas assembled at the Inner Lindblad Resonance (ILR) ring of the nuclear stellar bar. Conclusions. In our favored scenario a wide-angle AGN wind launched from the accretion disk of NGC1068 is currently impacting a sizable fraction of the gas inside the torus. However, a large gas reservoir (≃1.2 − 1.8 × 105 M⊙), which lies close to the equatorial plane of the torus, remains unaffected by the feedback of the AGN wind and can therefore continue fueling the AGN for at least ≃1 − 4 Myr. Nevertheless, AGN fueling currently seems thwarted on intermediate scales (15 pc ≤r ≤ 50 pc).

Dust Formation in AGN Winds

Abstract Infrared observations of active galactic nuclei (AGNs) reveal emission from the putative dusty circumnuclear “torus” invoked by AGN unification, which is heated up by radiation from the central accreting black hole (BH). The strong 9.7 and 18 μm silicate features observed in the AGN spectra, in both emission and absorption, further indicate the presence of such dusty environments. We present detailed calculations of the chemistry of silicate dust formation in AGN accretion disk winds. The winds considered herein are magnetohydrodynamic winds driven off the entire accretion disk domain that extends from the BH vicinity to the radius of BH influence, of order ∼1–100 pc depending on the mass of the resident BH. Our results indicate that these winds provide conditions conducive to the formation of significant amounts of dust, especially for objects accreting close to their Eddington limit, making AGNs a significant source of dust in the universe, especially for luminous quasars. Our models justify the importance of an r −1 density law in the winds for efficient formation and survival of dust grains. The dust production rate scales linearly with the mass of the central BH and varies as a power law of index between 2 and 2.5 with the dimensionless mass accretion rate. The resultant distribution of the dense dusty gas resembles a toroidal shape, with high column density and optical depths along the equatorial viewing angles, in agreement with the AGN unification picture.

The role of failed accretion disk winds in active galactic nuclei

AbstractBoth observational and theoretical evidence point at outflows originating from accretion disks as fundamental ingredients of active galactic nuclei (AGN). These outflows can have more than one component, for example an unbound supersonic wind and a failed wind (FW). The latter is a prediction of the simulations of radiation-driven disk outflows which show that the former is accompanied by an inner failed component, where the flow struggles to escape from the strong gravitational pull of the supermassive black hole. This FW component could provide a physical framework to interpret various phenomenological components of AGN. Here we briefly discuss a few of them: the broad line region, the X-ray obscurer, and the X-ray corona.

Tracing AGN feedback, from the SMBH horizon up to cluster scales

AbstractObservations performed in the last decades have shown that supermassive black holes (SMBHs) and cosmic structures are not separate elements of the Universe. While galaxies extend on spatial scales about ten orders of magnitude larger than the horizon of SMBHs, black holes would not exist without matter feeding them, and cosmic structures would not be the same without feedback from SMBHs. Powerful winds/jets in active galactic nuclei (AGN) may be the basis of this co-evolution. Synergistic observations in the X-rays and other wavebands have been proven to be fundamental to map AGN winds from the event horizon up to galaxy scales, providing a promising avenue to study the multi-phase SMBH feeding and feedback processes. Moreover, a spatially resolved, spectroscopic analysis of AGN in clusters will allow us to probe the multiphase medium ranging from galactic up to cluster scales. Revolutionary advances are expected in the upcoming decade with new multi-wavelength observatories, ranging from radio to X-rays.

Outflows in the Disks of Active Galaxies

Abstract Recent advances in observations have provided a wealth of measurements of the expansions of outflows in galactic disks out to large radii in a variety of galactic hosts. To provide an updated baseline for the interpretation of such data, and to assess to what extent the present status of the modeling is consistent with the existing observations, we provide a compact two-dimensional description for the expansion of active galactic nucleus (AGN)-driven shocks in realistic galactic disks with exponential gas density profiles in a disk geometry. We derive solutions for the outflow expansion and the mass outflow rates in different directions with respect to the plane of the disk. These are expressed in terms of the global properties of the host galaxy and of the central AGN to allow for an easy and direct comparison with existing observations in a variety of galactic hosts with measured properties, and out to distances of ∼10 kpc from the center. The results are compared with a state-of-the-art compilation of observed outflows in 19 galaxies with different measured gas and dynamical mass, allowing for a detailed, one-by-one comparison with the model predictions. The agreement we obtain for a wide range of host galaxy gas mass ( ) and AGN bolometric luminosity ( ) provides a quantitative systematic test for the modeling of AGN-driven outflows in galactic disks. We also consider a larger sample of 48 objects in galaxies with no reliable measurements of the gas and dynamical mass. In this case, we perform a comparison of the model predictions for different bins of AGN luminosities assuming different reference values for the gas mass and dynamical mass derived from average scaling relations. Finally, we reconsider the AGN wind scaling laws empirically derived by many authors in light of the results from our updated models. The encouraging, quantitative agreement of the model predictions with a wide set of existing observations constitutes a baseline for the interpretation of forthcoming data, and for a more detailed treatment of AGN feedback in galaxy formation models.

The origin of radio emission from radio-quiet active galactic nuclei

The central nuclei of galaxies, where supermassive black holes (SMBHs) are thought to reside, can experience phases of activity when they become active galactic nuclei (AGNs). An AGN can eject winds and jets and produce radiation across the entire electromagnetic spectrum. The fraction of the bolometric emission in the radio spans a factor of approximately 105 across the different classes of AGNs. The weakest radio sources, radio-quiet (RQ) AGNs, are typically 1,000 times fainter than the radio-loud (RL) AGNs, and represent the majority of the AGN population. In RQ AGNs, the absence of luminous jets allows us to probe radio emission from a wide range of possible mechanisms: star formation, AGN-driven wind, free-free emission from photoionized gas, low-power jets and the innermost accretion disk coronal activity. All these mechanisms can now be probed with unprecedented precision and spatial resolution, owing to the current and forthcoming generation of highly sensitive radio arrays. This article reviews radio emission mechanisms in radio-quiet active galactic nuclei (AGNs), from star formation and AGN winds, to free-free emission from photoionized gas and AGN disk coronal activity. These mechanisms can be probed by sensitive radio observatories.

CHEERS Results from NGC 3393. III. Chandra X-Ray Spectroscopy of the Narrow Line Region

Abstract We present spatially resolved Chandra narrow-band imaging and imaging spectroscopy of NGC 3393. This galaxy hosts a Compton-thick Seyfert 2 active galactic nucleus (AGN) with sub-kpc bipolar outflows that are strongly interacting with the circumnuclear gas. We identify narrow-band excess emission associated with the Ne ix 0.905 keV transition (with likely contributions due to intermediate-state iron emission) that points to strong shocks driven by AGN feedback. Imaging spectroscopy resolves outflow-interstellar medium (ISM) interaction sites and the surrounding ISM at ∼100 pc scales, and suggests the presence of a hot AGN wind above the plane at radii beyond the shock sites. The cross-cone shows evidence for reprocessing of photoionization that has passed through gaps in the torus, and also for collisionally excited plasma that may be powered by a shock-confined equatorial outflow. Deep X-ray observations at subarcsecond resolution (such as may be performed very efficiently by Lynx, which would also energetically resolve the complex line emission) are necessary to eliminate model degeneracies and reduce uncertainties in local feedback properties.

Probing the nature and origin of dust in the reddened quasar IC 4329A with global modelling from X-ray to infrared

Cosmic dust is a key tracer of structure formation and evolution in the universe. In active galactic nuclei (AGN) the origin and role of dust are uncertain. Here, we have studied dust in the X-ray bright and reddened type-1 quasar IC 4329A, which exhibits an ionised AGN wind. We incorporated high-resolution X-ray and mid-IR spectroscopy, combined with broad-band continuum modelling, to investigate the properties of dust in this AGN. We used new chandra HETGS observations taken in June 2017, as well as archival data from XMM-Newton, Swift, HST, Spitzer, IRAS, and Herschel for our IR-optical-UV-X-ray modelling. Two distinct components of dust in IC 4329A are found. One component is in the interstellar medium (ISM) of the host galaxy, and the other is a nuclear component in the AGN torus and its associated wind. The emitting dust in the torus is evident in mid-IR emission (9.7 and 18 μm features), while dust in the wind is present through both reddening and X-ray absorption (O, Si, and Fe edge features). The gas depletion factors into dust for O, Si, and Fe are measured. We derive an intrinsic reddening E(B – V) ≈ 1.0, which is most consistent with a grey (flat) extinction law. The AGN wind consists of three ionisation components. From analysis of long-term changes in the wind, we determine limits on the location of the wind components. The two lowest ionisation components are likely carriers of dust from the AGN torus. We find that the dust in the nuclear component of IC 4329A is different from dust in the Milky Way. The dust grains in the AGN torus and wind are likely larger than the standard Galactic dust, and are in a porous composite form (containing amorphous silicate with iron and oxygen). This can be a consequence of grain coagulation in the dense nuclear environment of the AGN.

Multi-wavelength campaign on NGC 7469

We investigate the physical structure of the active galactic nucleus (AGN) wind in the Seyfert-1 galaxy NGC 7469 through high-resolution X-ray spectroscopy with Chandra HETGS and photoionisation modelling. Contemporaneous data from Chandra, HST, and Swift are used to model the optical-UV-X-ray continuum and determine the spectral energy distribution (SED) at two epochs, 13 yr apart. For our investigation we use new observations taken in December 2015–January 2016, and historical ones taken in December 2002. We study the impact of a change in the SED shape, seen between the two epochs, on the photoionisation of the wind. The HETGS spectroscopy shows that the AGN wind in NGC 7469 consists of four ionisation components, with their outflow velocities ranging from − 400 to − 1800 km s-1. From our modelling we find that the change in the ionising continuum shape between the two epochs results in some variation in the ionisation state of the wind components. However, for the main ions detected in X-rays, the sum of their column densities over all components remains in practice unchanged. For two of the four components, which are found to be thermally unstable in both epochs, we obtain 2 < r < 31 pc and 12 < r < 29 pc using the cooling and recombination timescales. For the other two thermally stable components, we obtain r < 31 pc and r < 80 pc from the recombination timescale. The results of our photoionisation modelling and thermal stability analysis suggest that the absorber components in NGC 7469 are consistent with being a thermally driven wind from the AGN torus. Finally, from analysis of the zeroth-order ACIS/HETG data, we discover that the X-ray emission in the range 0.2–1 keV is spatially extended over 1.5–12′′. This diffuse soft X-ray emission is explained by coronal emission from the nuclear starburst ring in NGC 7469. The star formation rate inferred from this diffuse soft X-ray emission is consistent with those found by far-infrared studies of NGC 7469.

Thick turbulent gas disks with magnetocentrifugal winds in active galactic nuclei

Infrared high-resolution imaging and interferometry show that the dust distribution is frequently elongated along the polar direction of an AGN. In addition, interferometric mm line observations have revealed a bipolar outflow in a direction nearly perpendicular to the nuclear disk. To explain these findings, we developed a model scenario for the inner ~30 pc of an AGN. The structure of the gas within this region is entirely determined by the gas inflow from larger scales. We assumed a rotating thick gas disk between about one and ten parsec. External gas accretion adds mass and injects energy via gas compression into this gas disk and drives turbulence. We extended the description of a massive turbulent thick gas disk developed in a recent paper by adding a magnetocentrifugal wind. Our disks are assumed to be strongly magnetized via equipartition between the turbulent gas pressure and the energy density of the magnetic field. In a second step, we built 3D density cubes based on the analytical model, illuminated them with a central source, and made radiative transfer calculations. In a third step, we calculated mid-infrared (MIR) visibility amplitudes and compared them to available interferometric observations. We show that magnetocentrifugal winds starting from a thin and thick gas disk are viable in active galaxy centers. The magnetic field associated with this thick gas disk plays a major role in driving a magnetocentrifugal wind at a distance of ~1 pc from the central black hole. Once the wind is launched, it is responsible for the transport of angular momentum and the gas disk can become thin. A magnetocentrifugal wind is also expected above the thin magnetized gas disk. The structure and outflow rate of this wind is determined by the properties of the thick gas disk. The outflow scenario can account for the elongated dust structures, outer edges of the thin maser disks, and molecular outflows observed in local AGN. The models reproduce the observed terminal wind velocities, the scatter of the MIR – intrinsic X-ray correlation, and point source fractions. An application of the model to the Circinus galaxy and NGC 1068 shows that the infrared spectral energy distribution, available MIR interferometric observations, and optical polarization can be reproduced in a satisfactory way, provided that (i) a puff-up at the inner edge of the thin disk is present and (ii) a local screen with an optical depth of τV ~ 20 in form of a local gas filament and/or a warp of the thick disk hide a significant fraction of both nuclei. Our thick disk, wind, thin disk model is thus a promising scenario for local Seyfert galaxies.

Can Winds Driven by Active Galactic Nuclei Account for the Extragalactic Gamma-Ray and Neutrino Backgrounds?

Abstract Various observations are revealing the widespread occurrence of fast and powerful winds in active galactic nuclei (AGNs) that are distinct from relativistic jets, likely launched from accretion disks and interacting strongly with the gas of their host galaxies. During the interaction, strong shocks are expected to form that can accelerate nonthermal particles to high energies. Such winds have been suggested to be responsible for a large fraction of the observed extragalactic gamma-ray background (EGB) and the diffuse neutrino background, via the decay of neutral and charged pions generated in inelastic pp collisions between protons accelerated by the forward shock and the ambient gas. However, previous studies did not properly account for processes such as adiabatic losses that may reduce the gamma-ray and neutrino fluxes significantly. We evaluate the production of gamma rays and neutrinos by AGN-driven winds in detail by modeling their hydrodynamic and thermal evolution, including the effects of their two-temperature structure. We find that they can only account for less than ∼30% of the EGB flux, as otherwise the model would violate the independent upper limit derived from the diffuse isotropic gamma-ray background. If the neutrino spectral index is steep with Γ ≳ 2.2, a severe tension with the isotropic gamma-ray background would arise as long as the winds contribute more than 20% of the IceCube neutrino flux in the 10–100 TeV range. At energies ≳ 100 TeV, we find that the IceCube neutrino flux may still be accountable by AGN-driven winds if the spectral index is as small as Γ ∼ 2.0–2.1.

The Effects of Ram Pressure on the Cold Clouds in the Centers of Galaxy Clusters

Abstract We discuss the effect of ram pressure on the cold clouds in the centers of cool-core galaxy clusters, and in particular, how it reduces cloud velocity and sometimes causes an offset between the cold gas and young stars. The velocities of the molecular gas in both observations and our simulations fall in the range of 100–400 km s−1, which is much lower than expected if they fall from a few tens of kiloparsecs ballistically. If the intracluster medium (ICM) is at rest, the ram pressure of the ICM only slightly reduces the velocity of the clouds. When we assume that the clouds are actually “fluffier” because they are co-moving with a warm-hot layer, the velocity becomes smaller. If we also consider the active galactic nucleus wind in the cluster center by adding a wind profile measured from the simulation, the clouds are further slowed down at small radii, and the resulting velocities are in general agreement with the observations and simulations. Because ram pressure only affects gas but not stars, it can cause a separation between a filament and young stars that formed in the filament as they move through the ICM together. This separation has been observed in Perseus and also exists in our simulations. We show that the star-filament offset, combined with line-of-sight velocity measurements, can help determine the true motion of the cold gas, and thus distinguish between inflows and outflows.

NuSTAR View of the Black Hole Wind in the Galaxy Merger IRAS F11119+3257

Abstract Galactic winds driven by active galactic nuclei (AGNs) have been invoked to play a fundamental role in the co-evolution between supermassive black holes and their host galaxies. Finding observational evidence of such feedback mechanisms is of crucial importance and it requires a multi-wavelength approach in order to compare winds at different scales and phases. In Tombesi et al., we reported the detection of a powerful ultra-fast outflow (UFO) in the Suzaku X-ray spectrum of the ultra-luminous infrared galaxy IRAS F11119+3257. The comparison with a galaxy-scale OH molecular outflow observed with Herschel in the same source supported the energy-conserving scenario for AGN feedback. The main objective of this work is to perform an independent check of the Suzaku results using the higher sensitivity and wider X-ray continuum coverage of NuSTAR. We clearly detect a highly ionized Fe K UFO in the 100 ks NuSTAR spectrum with parameters N H = (3.2 ± 1.5) × 1024 cm−2, log ξ = erg s−1 cm, and . The launching radius is likely at a distance of r ≥ 16r s from the black hole. The mass outflow rate is in the range of ≃ 0.5–2 M ⊙ yr−1. The UFO momentum rate and power are ≃ 0.5–2 L AGN/c and ≃ 7%–27% L AGN, respectively. The UFO parameters are consistent between the 2013 Suzaku and the 2015 NuSTAR observations. Only the column density is found to be variable, possibly suggesting a clumpy wind. The comparison with the energetics of molecular outflows estimated in infrared and millimeter wavelengths support a connection between the nuclear and galaxy-scale winds in luminous AGNs.

The origin of fast molecular outflows in quasars: molecule formation in AGN-driven galactic winds

We explore the origin of fast molecular outflows that have been observed in Active Galactic Nuclei (AGN). Previous numerical studies have shown that it is difficult to create such an outflow by accelerating existing molecular clouds in the host galaxy, as the clouds will be destroyed before they can reach the high velocities that are observed. In this work, we consider an alternative scenario where molecules form in-situ within the AGN outflow. We present a series of hydro-chemical simulations of an isotropic AGN wind interacting with a uniform medium. We follow the time-dependent chemistry of 157 species, including 20 molecules, to determine whether molecules can form rapidly enough to produce the observed molecular outflows. We find H$_2$ outflow rates up to 140 M$_\odot$ yr$^{-1}$, which is sensitive to density, AGN luminosity, and metallicity. We compute emission and absorption lines of CO, OH and warm (a few hundred K) H$_2$ from the simulations in post-processing. The CO-derived outflow rates and OH absorption strengths at solar metallicity agree with observations, although the maximum line of sight velocities from the model CO spectra are a factor $\approx$2 lower than is observed. We derive a CO (1-0) to H$_2$ conversion factor of $\alpha_{\rm{CO} (1-0)}$ = 0.13 M$_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, 6 times lower than is commonly assumed in observations of such systems. We find strong emission from the mid-infrared lines of H$_2$. The mass of H$_2$ traced by this infrared emission is within a few per cent of the total H$_2$ mass. This H$_2$ emission may be observable by JWST.

The reverberation signatures of rotating disc winds in active galactic nuclei

The broad emission lines (BELs) in active galactic nuclei (AGN) respond to ionizing continuum variations. The time and velocity dependence of their response depends on the structure of the broad-line region: its geometry, kinematics and ionization state. Here, we predict the reverberation signatures of BELs formed in rotating accretion disc winds. We use a Monte Carlo radiative transfer and ionization code to predict velocity-delay maps for representative high- (C$~IV$) and low-ionization (H$\alpha$) emission lines in both high- and moderate-luminosity AGN. Self-shielding, multiple scattering and the ionization structure of the outflows are all self-consistently taken into account, while small-scale structure in the outflow is modelled in the micro-clumping approximation. Our main findings are: (1) The velocity-delay maps of smooth/micro-clumped outflows often contain significant negative responses. (2)~The reverberation signatures of disc wind models tend to be rotation dominated and can even resemble the classic "red-leads-blue" inflow signature. (3) Traditional "blue-leads-red" outflow signatures can usually only be observed in the long-delay limit. (4) Our models predict lag-luminosity relationships similar to those inferred from observations, but systematically underpredict the observed centroid delays. (5) The ratio between "virial product" and black hole mass predicted by our models depends on viewing angle. Our results imply that considerable care needs to be taken in interpreting data obtained by observational reverberation mapping campaigns. In particular, basic signatures such as "red-leads-blue", "blue-leads-red" and "blue and red vary jointly" are not always reliable indicators of inflow, outflow or rotation. This may help to explain the perplexing diversity of such signatures seen in observational campaigns to date.