Skip to content
Visualisation of the Structure of Active Galactic Nuclei
Python JavaScript HTML
Branch: master
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Permalink
Type Name Latest commit message Commit time
Failed to load latest commit information.
img
.nojekyll
LICENSE
README.rst
alphadisk.py
blrtor.py
combinedstructure.py
index.html
jquery-1.9.1.min.js
luminosities.js
radfountain-angle.py
radfountain.py
radfountain_2d.py
radfountain_v2.py
simcones.py
windyblr.py
windyblr_2d.py

README.rst

Visualisation of the Structure of Active Galactic Nuclei

Model-driven illustrations of centres of galaxies.

See it in action at https://johannesbuchner.github.io/agnviz/

Ingredients and Assumptions

This visualisation allows extrapolation outside the range of validity of the used models. This includes going to very high and low Eddington rates. Most reliable it should be at Eddington rates of 0.01-1.

This visualisation is an idealisation of a steady-state situation. In reality, the host galaxy morphology, chaotic accretion events leading to re- and mis-alignments and variability make the picture more complex.

Click to play with visualisation

The following components are implemented:

  • BH: Black hole

    • Schwarzschild radius (r_g)
  • AD: Accretion disk

  • Corona: X-ray emitting corona

  • BLR: Broad line emission region

    • BLCH model (Baskin & Laor (2017), Czerny & Hryniewicz (2011))
    • describes where accretion disk radiation pressure lifts dusty material up (can be clumpy)
    • The left edge is the dust sublimation radius of the smallest graphites (2000K).
    • This model is known to produce a little too low covering factors
    • Note that at low L/M^(2/3), the BLR goes away (see Elitzur & Netzer (2016))
  • RF TOR: Toroidal Obscuring Region

    • Material is likely clumpy
    • analytic form of the radiative fountain model from Wada (2016)
    • describes where dust is no longer accelerated by radiation pressure. Material can pile up near and to the right of the line up to >>10 pc
    • Note that at low L/M^(2/3), the region is closed (high covering fraction).
    • This region starts near the sublimation radius of most dust grains (~1.3pc sqrt(L46)), but will be stratified by dust grain size (possibly down to the BLR, see above).
  • NLR: Narrow line emission region

    • Appears to be shaped substantially by the host galaxy (see e.g. Prieto et al. 2014) as well as the nuclear obscurer.
  • Jet

    • Width is an empirical intermediate (from Algaba et al. (2017)) between

      • quasi-conical streamline (Blandford & Znajek 1977)
      • force-free steady jet solution (Narayan et al 2007; Tchekhovskoy et al 2008)
    • Base is assumed at r_g, maximum 10^8 pc.

  • Sphere of influence

    • Where gravity from the BH dominates over the host galaxy
    • Computed the host gravity with M_BH-sigma from Kormendy & Ho (2013)
  • Outflows & Inflows

    • Inflow onto BH is given by bolometric luminosity and accretion efficiency (assumed 10%)
    • Outflow relation from Fluetsch et al. (2018), when setting SFR=0.
  • Viewing angles

    • Are a bit difficult to interpret in log-log plots.
    • Note that views to the corona, to the BLR and to the TOR are not necessarily identical.

Usage

Clone this repository, and look at the end of combinedstructure.py how to use the plot_log_agn_postcard function:

from combinedstructure import plot_log_agn_postcard

import astropy.units as u

MBH = 1e8 * u.Msun
eddrate = 0.1

import matplotlib.pyplot as plt
plt.figure(figsize=(10,7))

plot_log_agn_postcard(MBH, eddrate,
        # for Baskin & Laor (2017)
        rad_efficiency = 0.1,
        Z = 5, # metallicity
        # Radiative Fountain system assumptions
        dust_to_gas_ratio = 1/20.,
        kappa = 1e3 * u.cm**2/u.g,
        rho_g = 300 * u.Msun / u.pc**3,
        show_BH = True,
        show_corona = True,
        show_disk = True,
        show_BLR = True,
        show_NLR = True,
        show_jet = True,
        show_TOR = True,
        show_SOI = True,
        show_viewing = True,
        colored_disk = True,
        show_flows = True,
)

plt.show()

If useful for your work and paper, please cite this repository URL. Code is BSD licensed, see LICENSE file.

Conferences and talks about AGN

Conferences:

Talks:

More general:

You can’t perform that action at this time.