Black Hole Spectra¶
Black hole spectra can be generated by combining a BlackHoles (for particle, BlackHole for parametric) object with an EmissionModel, translating the physical properties of the blackhole(s) (e.g. mass, accretion_rate, etc.) to a spectral energy distribution.
These models are described in detail in the emission model docs. Here, we’ll use an instance of a UnifiedAGN model for demonstration purposes.
The following sections demonstrate the generation of combined spectra (which is the same for both parametric and particle BlackHoles) and per-particle spectra.
[1]:
import numpy as np
from unyt import K, Mpc, Msun, cm, deg, yr
from synthesizer import Grid
from synthesizer.emission_models import (
Greybody,
UnifiedAGN,
)
from synthesizer.emission_models.attenuation import PowerLaw
from synthesizer.parametric import BlackHole
# Get the NLR and BLR grids
nlr_grid = Grid("test_grid_agn-nlr")
blr_grid = Grid("test_grid_agn-blr")
# Initialise the BlackHole object setting most of the key attributes
blackhole = BlackHole(
mass=1e8 * Msun,
inclination=60 * deg,
accretion_rate_eddington=0.1,
covering_fraction_nlr=0.1,
covering_fraction_blr=0.1,
metallicity=0.01,
theta_torus=20 * deg,
)
# Initialise the UnifiedAGN model
uniagn = UnifiedAGN(
nlr_grid,
blr_grid,
ionisation_parameter_nlr=0.01,
hydrogen_density_nlr=1e4 * cm**-3,
ionisation_parameter_blr=0.1,
hydrogen_density_blr=1e10 * cm**-3,
torus_emission_model=Greybody(1000 * K, 1.5),
)
Integrated spectra¶
To generate integrated spectra we simply call the component’s get_spectra method. This method will populate the component’s spectra attribute with a dictionary containing Sed objects for each spectra in the EmissionModel. It will also return the spectra at the root of the EmissionModel.
[2]:
# Get the spectra using a unified agn model (instantiated elsewhere)
spectra = blackhole.get_spectra(uniagn)
fig, ax = blackhole.plot_spectra(
show=True,
ylimits=(10**27.5, 10**34.0),
figsize=(10, 8),
)
print(blackhole.model_param_cache)
{'full_reprocessed_blr': {'inclination': unyt_array(60., 'degree'), 'theta_torus': unyt_array(20., 'degree'), 'torus_edgeon_cond': unyt_quantity(80., 'degree'), 'mass': unyt_quantity(1.e+08, 'Msun'), 'log10mass': array(8.), 'accretion_rate_eddington': np.float64(3.827e+32), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.1), 'hydrogen_density': unyt_quantity(1.e+09, 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'nebular', 'emitter': 'blackhole', 'masks': 'torus_edgeon_cond < 90 degree'}, 'full_reprocessed_nlr': {'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.01), 'hydrogen_density': unyt_quantity(10000., 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'nebular', 'emitter': 'blackhole'}, 'disc_incident_isotropic': {'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.1), 'hydrogen_density': unyt_quantity(1.e+09, 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'incident', 'emitter': 'blackhole'}, 'disc_transmitted_blr_full': {'inclination': unyt_quantity(60., 'degree'), 'theta_torus': unyt_array(20., 'degree'), 'torus_edgeon_cond': unyt_quantity(80., 'degree'), 'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5000000000000001), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.1), 'hydrogen_density': unyt_quantity(1.e+09, 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'transmitted', 'emitter': 'blackhole', 'masks': 'torus_edgeon_cond < 90 degree'}, 'disc_incident_masked': {'inclination': unyt_quantity(60., 'degree'), 'theta_torus': unyt_array(20., 'degree'), 'torus_edgeon_cond': unyt_quantity(80., 'degree'), 'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5000000000000001), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.1), 'hydrogen_density': unyt_quantity(1.e+09, 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'incident', 'emitter': 'blackhole', 'masks': 'torus_edgeon_cond < 90 degree'}, 'disc_transmitted_nlr_full': {'inclination': unyt_quantity(60., 'degree'), 'theta_torus': unyt_array(20., 'degree'), 'torus_edgeon_cond': unyt_quantity(80., 'degree'), 'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5000000000000001), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.01), 'hydrogen_density': unyt_quantity(10000., 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'transmitted', 'emitter': 'blackhole', 'masks': 'torus_edgeon_cond < 90 degree'}, 'disc_transmitted_nlr_isotropic_full': {'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.01), 'hydrogen_density': unyt_quantity(10000., 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'transmitted', 'emitter': 'blackhole'}, 'disc_transmitted_blr_isotropic_full': {'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.1), 'hydrogen_density': unyt_quantity(1.e+09, 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'transmitted', 'emitter': 'blackhole'}, 'disc_incident': {'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5000000000000001), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.1), 'hydrogen_density': unyt_quantity(1.e+09, 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'incident', 'emitter': 'blackhole'}, 'full_continuum_nlr': {'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.01), 'hydrogen_density': unyt_quantity(10000., 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'nebular_continuum', 'emitter': 'blackhole'}, 'full_continuum_blr': {'inclination': unyt_quantity(60., 'degree'), 'theta_torus': unyt_array(20., 'degree'), 'torus_edgeon_cond': unyt_quantity(80., 'degree'), 'log10mass': array(8.), 'log10accretion_rate_eddington': np.float64(32.5828584622245), 'cosine_inclination': np.float64(0.5), 'metallicity': np.float64(0.01), 'metallicities': np.float64(0.01), 'ionisation_parameter': np.float64(0.1), 'hydrogen_density': unyt_quantity(1.e+09, 'cm**(-3)'), 'bolometric_luminosity': unyt_quantity(1.2567868e+45, 'erg/s'), 'bolometric_luminosities': unyt_quantity(1.2567868e+45, 'erg/s'), 'extract': 'nebular_continuum', 'emitter': 'blackhole', 'masks': 'torus_edgeon_cond < 90 degree'}, 'blr': {'covering_fraction_blr': np.float64(0.1), 'apply_to': 'full_reprocessed_blr', 'transformer': "CoveringFraction(covering_attrs=('covering_fraction_blr',))", 'emitter': 'blackhole'}, 'nlr': {'covering_fraction_nlr': np.float64(0.1), 'apply_to': 'full_reprocessed_nlr', 'transformer': "CoveringFraction(covering_attrs=('covering_fraction_nlr',))", 'emitter': 'blackhole'}, 'torus': {'temperature': unyt_array(1000., 'K'), 'emissivity': np.float64(1.5), 'generator': 'Greybody(scaler=disc_incident_isotropic, temperature=1000.0 K, emissivity=1.5, optically_thin=True, lam_0=100.0 μm)', 'emitter': 'blackhole'}, 'disc_transmitted_blr': {'transmission_fraction_blr': np.float64(0.0), 'apply_to': 'disc_transmitted_blr_full', 'transformer': "CoveringFraction(covering_attrs=('transmission_fraction_blr',))", 'emitter': 'blackhole'}, 'disc_escaped': {'transmission_fraction_escape': np.float64(1.0), 'apply_to': 'disc_incident_masked', 'transformer': "CoveringFraction(covering_attrs=('transmission_fraction_escape',))", 'emitter': 'blackhole'}, 'disc_transmitted_nlr': {'transmission_fraction_nlr': np.float64(0.0), 'apply_to': 'disc_transmitted_nlr_full', 'transformer': "CoveringFraction(covering_attrs=('transmission_fraction_nlr',))", 'emitter': 'blackhole'}, 'disc_transmitted': {'combine': ['disc_escaped', 'disc_transmitted_nlr', 'disc_transmitted_blr'], 'emitter': 'blackhole'}, 'disc': {'combine': ['disc_transmitted'], 'emitter': 'blackhole'}, 'disc_escaped_isotropic': {'covering_fraction_blr': np.float64(0.1), 'covering_fraction_nlr': np.float64(0.1), 'apply_to': 'disc_incident_isotropic', 'transformer': "EscapingFraction(covering_attrs=('covering_fraction_blr', 'covering_fraction_nlr'))", 'emitter': 'blackhole'}, 'disc_transmitted_nlr_isotropic': {'covering_fraction_nlr': np.float64(0.1), 'apply_to': 'disc_transmitted_nlr_isotropic_full', 'transformer': "CoveringFraction(covering_attrs=('covering_fraction_nlr',))", 'emitter': 'blackhole'}, 'disc_transmitted_blr_isotropic': {'covering_fraction_blr': np.float64(0.1), 'apply_to': 'disc_transmitted_blr_isotropic_full', 'transformer': "CoveringFraction(covering_attrs=('covering_fraction_blr',))", 'emitter': 'blackhole'}, 'disc_averaged_without_torus': {'combine': ['disc_escaped_isotropic', 'disc_transmitted_nlr_isotropic', 'disc_transmitted_blr_isotropic'], 'emitter': 'blackhole'}, 'disc_averaged': {'torus_fraction': array(0.22222222), 'apply_to': 'disc_averaged_without_torus', 'transformer': "EscapingFraction(covering_attrs=('torus_fraction',))", 'emitter': 'blackhole'}, 'disc_escaped_weighted': {'escape_fraction': np.float64(0.8), 'apply_to': 'disc_incident_masked', 'transformer': "CoveringFraction(covering_attrs=('escape_fraction',))", 'emitter': 'blackhole'}, 'disc_transmitted_blr_weighted': {'covering_fraction_blr': np.float64(0.1), 'apply_to': 'disc_transmitted_blr_full', 'transformer': "CoveringFraction(covering_attrs=('covering_fraction_blr',))", 'emitter': 'blackhole'}, 'disc_transmitted_nlr_weighted': {'covering_fraction_nlr': np.float64(0.1), 'apply_to': 'disc_transmitted_nlr_full', 'transformer': "CoveringFraction(covering_attrs=('covering_fraction_nlr',))", 'emitter': 'blackhole'}, 'disc_transmitted_weighted_combination': {'combine': ['disc_escaped_weighted', 'disc_transmitted_nlr_weighted', 'disc_transmitted_blr_weighted'], 'emitter': 'blackhole'}, 'line_regions': {'combine': ['nlr', 'blr'], 'emitter': 'blackhole'}, 'nlr_continuum': {'covering_fraction_nlr': np.float64(0.1), 'apply_to': 'full_continuum_nlr', 'transformer': "CoveringFraction(covering_attrs=('covering_fraction_nlr',))", 'emitter': 'blackhole'}, 'blr_continuum': {'covering_fraction_blr': np.float64(0.1), 'apply_to': 'full_continuum_blr', 'transformer': "CoveringFraction(covering_attrs=('covering_fraction_blr',))", 'emitter': 'blackhole'}, 'intrinsic': {'combine': ['disc', 'nlr', 'blr', 'torus'], 'emitter': 'blackhole'}}
Including dust attenuation¶
We can also generate spectra including attenuation and emission from diffuse dust along the line of sight to the black hole. This is now possible directly with the UnifiedAGN by passing a dust curve. The optical depth (tau_v) must be available on the emitter (the blackhole) or set by the emission model.
[3]:
tau_v = 0.5
# Initialise the UnifiedAGN model
uniagn_attenuated = UnifiedAGN(
nlr_grid,
blr_grid,
ionisation_parameter_nlr=0.01,
hydrogen_density_nlr=1e4 * cm**-3,
ionisation_parameter_blr=0.1,
hydrogen_density_blr=1e10 * cm**-3,
torus_emission_model=Greybody(1000 * K, 1.5),
diffuse_dust_curve=PowerLaw(slope=-1.0),
tau_v=tau_v,
)
We then follow the same process of calling get_spectra with the new model. The plot here shows luminosity rather than spectral energy density.
[4]:
spectra = blackhole.get_spectra(uniagn_attenuated)
fig, ax = blackhole.plot_spectra(
quantity_to_plot="luminosity",
figsize=(6, 4),
spectra_to_plot=["intrinsic", "attenuated"],
)
The spectra returned by get_spectra is the “dust_emission” spectra at the root of the emission model.
[5]:
print(spectra)
+----------------------------------------------------------------------------------------------------+
| SED |
+---------------------------+------------------------------------------------------------------------+
| Attribute | Value |
+---------------------------+------------------------------------------------------------------------+
| redshift | 0 |
+---------------------------+------------------------------------------------------------------------+
| ndim | 1 |
+---------------------------+------------------------------------------------------------------------+
| nlam | 9244 |
+---------------------------+------------------------------------------------------------------------+
| shape | (9244,) |
+---------------------------+------------------------------------------------------------------------+
| lam (9244,) | 1.30e-04 Å -> 2.99e+11 Å (Mean: 9.73e+09 Å) |
+---------------------------+------------------------------------------------------------------------+
| nu (9244,) | 1.00e+07 Hz -> 2.31e+22 Hz (Mean: 8.51e+19 Hz) |
+---------------------------+------------------------------------------------------------------------+
| lnu (9244,) | 0.00e+00 erg/(Hz*s) -> 1.45e+31 erg/(Hz*s) (Mean: 1.12e+29 erg/(Hz*s)) |
+---------------------------+------------------------------------------------------------------------+
| bolometric_luminosity | 3.851157354809097e+44 erg/s |
+---------------------------+------------------------------------------------------------------------+
| energy (9244,) | 4.14e-08 eV -> 9.56e+07 eV (Mean: 3.52e+05 eV) |
+---------------------------+------------------------------------------------------------------------+
| frequency (9244,) | 1.00e+07 Hz -> 2.31e+22 Hz (Mean: 8.51e+19 Hz) |
+---------------------------+------------------------------------------------------------------------+
| llam (9244,) | 0.00e+00 erg/(s*Å) -> 3.59e+41 erg/(s*Å) (Mean: 1.97e+39 erg/(s*Å)) |
+---------------------------+------------------------------------------------------------------------+
| luminosity (9244,) | 0.00e+00 erg/s -> 1.80e+45 erg/s (Mean: 1.25e+43 erg/s) |
+---------------------------+------------------------------------------------------------------------+
| luminosity_lambda (9244,) | 0.00e+00 erg/(s*Å) -> 3.59e+41 erg/(s*Å) (Mean: 1.97e+39 erg/(s*Å)) |
+---------------------------+------------------------------------------------------------------------+
| luminosity_nu (9244,) | 0.00e+00 erg/(Hz*s) -> 1.45e+31 erg/(Hz*s) (Mean: 1.12e+29 erg/(Hz*s)) |
+---------------------------+------------------------------------------------------------------------+
| wavelength (9244,) | 1.30e-04 Å -> 2.99e+11 Å (Mean: 9.73e+09 Å) |
+---------------------------+------------------------------------------------------------------------+
However, all the spectra are stored within a dictionary under the spectra attribute.
[6]:
print(blackhole.spectra)
{'full_reprocessed_blr': <synthesizer.emissions.sed.Sed object at 0x7f77eee27880>, 'full_reprocessed_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77eee27d00>, 'disc_incident_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f783ca7fe20>, 'disc_transmitted_blr_full': <synthesizer.emissions.sed.Sed object at 0x7f783ca7fb50>, 'disc_incident_masked': <synthesizer.emissions.sed.Sed object at 0x7f77eee626e0>, 'disc_transmitted_nlr_full': <synthesizer.emissions.sed.Sed object at 0x7f77ef17e530>, 'disc_transmitted_nlr_isotropic_full': <synthesizer.emissions.sed.Sed object at 0x7f783ca7de70>, 'disc_transmitted_blr_isotropic_full': <synthesizer.emissions.sed.Sed object at 0x7f783ca7f4c0>, 'disc_incident': <synthesizer.emissions.sed.Sed object at 0x7f783ca7d840>, 'full_continuum_nlr': <synthesizer.emissions.sed.Sed object at 0x7f783ca7fe50>, 'full_continuum_blr': <synthesizer.emissions.sed.Sed object at 0x7f77f0aff310>, 'blr': <synthesizer.emissions.sed.Sed object at 0x7f77eee61a20>, 'nlr': <synthesizer.emissions.sed.Sed object at 0x7f77eee62110>, 'torus': <synthesizer.emissions.sed.Sed object at 0x7f77eee271c0>, 'disc_transmitted_blr': <synthesizer.emissions.sed.Sed object at 0x7f77eee619f0>, 'disc_escaped': <synthesizer.emissions.sed.Sed object at 0x7f77eee60250>, 'disc_transmitted_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77eee60400>, 'disc_transmitted': <synthesizer.emissions.sed.Sed object at 0x7f77f0affb80>, 'disc': <synthesizer.emissions.sed.Sed object at 0x7f77f0aff0a0>, 'disc_escaped_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77eee61690>, 'disc_transmitted_nlr_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77eee62da0>, 'disc_transmitted_blr_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77eee61510>, 'disc_averaged_without_torus': <synthesizer.emissions.sed.Sed object at 0x7f77eee60fd0>, 'disc_averaged': <synthesizer.emissions.sed.Sed object at 0x7f77eee620b0>, 'disc_escaped_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77eee61420>, 'disc_transmitted_blr_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77eee622c0>, 'disc_transmitted_nlr_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77eee617b0>, 'disc_transmitted_weighted_combination': <synthesizer.emissions.sed.Sed object at 0x7f77eee639a0>, 'line_regions': <synthesizer.emissions.sed.Sed object at 0x7f77eee60b50>, 'nlr_continuum': <synthesizer.emissions.sed.Sed object at 0x7f77eee614e0>, 'blr_continuum': <synthesizer.emissions.sed.Sed object at 0x7f77eee613f0>, 'intrinsic': <synthesizer.emissions.sed.Sed object at 0x7f77eee61de0>, 'attenuated': <synthesizer.emissions.sed.Sed object at 0x7f77eee61480>}
Particle spectra¶
To demonstrate the particle spectra functionality we first generate some mock particle black hole data, and initialise a BlackHoles object.
[7]:
from synthesizer.particle import BlackHoles
# Make fake properties
n = 4
masses = 10 ** np.random.uniform(low=7, high=9, size=n) * Msun
coordinates = np.random.normal(0, 1.5, (n, 3)) * Mpc
accretion_rates = 10 ** np.random.uniform(low=-2, high=1, size=n) * Msun / yr
metallicities = np.full(n, 0.01)
# And get the black holes object
blackholes = BlackHoles(
masses=masses,
coordinates=coordinates,
accretion_rates=accretion_rates,
metallicities=metallicities,
ionisation_parameter_nlr=0.01,
hydrogen_density_nlr=1e4 * cm**-3,
ionisation_parameter_blr=0.1,
hydrogen_density_blr=1e10 * cm**-3,
)
To generate a spectra for each black hole (per particle) we use the same emission model, but we need to tell the model to produce a spectrum for each particle. This is done by setting the per_particle flag to True on the model.
[8]:
uniagn_attenuated.set_per_particle(True)
With that done we just call the same get_spectra method on the component, and the particle spectra will be stored in the particle_spectra attribute of the component.
[9]:
spectra = blackholes.get_spectra(uniagn_attenuated, verbose=True)
Again, the returned spectra is the “dust_emission” spectra from the root of the model.
[10]:
print(spectra)
+------------------------------------------------------------------------------------------------------+
| SED |
+-----------------------------+------------------------------------------------------------------------+
| Attribute | Value |
+-----------------------------+------------------------------------------------------------------------+
| redshift | 0 |
+-----------------------------+------------------------------------------------------------------------+
| ndim | 2 |
+-----------------------------+------------------------------------------------------------------------+
| nlam | 9244 |
+-----------------------------+------------------------------------------------------------------------+
| shape | (4, 9244) |
+-----------------------------+------------------------------------------------------------------------+
| lam (9244,) | 1.30e-04 Å -> 2.99e+11 Å (Mean: 9.73e+09 Å) |
+-----------------------------+------------------------------------------------------------------------+
| nu (9244,) | 1.00e+07 Hz -> 2.31e+22 Hz (Mean: 8.51e+19 Hz) |
+-----------------------------+------------------------------------------------------------------------+
| lnu (4, 9244) | 0.00e+00 erg/(Hz*s) -> 1.96e+31 erg/(Hz*s) (Mean: 7.56e+28 erg/(Hz*s)) |
+-----------------------------+------------------------------------------------------------------------+
| bolometric_luminosity (4,) | 1.62e+44 erg/s -> 6.04e+44 erg/s (Mean: 3.13e+44 erg/s) |
+-----------------------------+------------------------------------------------------------------------+
| energy (9244,) | 4.14e-08 eV -> 9.56e+07 eV (Mean: 3.52e+05 eV) |
+-----------------------------+------------------------------------------------------------------------+
| frequency (9244,) | 1.00e+07 Hz -> 2.31e+22 Hz (Mean: 8.51e+19 Hz) |
+-----------------------------+------------------------------------------------------------------------+
| llam (4, 9244) | 0.00e+00 erg/(s*Å) -> 4.48e+41 erg/(s*Å) (Mean: 2.18e+39 erg/(s*Å)) |
+-----------------------------+------------------------------------------------------------------------+
| luminosity (4, 9244) | 0.00e+00 erg/s -> 2.25e+45 erg/s (Mean: 1.02e+43 erg/s) |
+-----------------------------+------------------------------------------------------------------------+
| luminosity_lambda (4, 9244) | 0.00e+00 erg/(s*Å) -> 4.48e+41 erg/(s*Å) (Mean: 2.18e+39 erg/(s*Å)) |
+-----------------------------+------------------------------------------------------------------------+
| luminosity_nu (4, 9244) | 0.00e+00 erg/(Hz*s) -> 1.96e+31 erg/(Hz*s) (Mean: 7.56e+28 erg/(Hz*s)) |
+-----------------------------+------------------------------------------------------------------------+
| wavelength (9244,) | 1.30e-04 Å -> 2.99e+11 Å (Mean: 9.73e+09 Å) |
+-----------------------------+------------------------------------------------------------------------+
While the spectra produced by get_particle_spectra are stored in a dictionary under the particle_spectra attribute.
[11]:
print(blackholes.particle_spectra)
{'disc_incident_masked': <synthesizer.emissions.sed.Sed object at 0x7f77eee90f40>, 'disc_transmitted_nlr_full': <synthesizer.emissions.sed.Sed object at 0x7f77eee90670>, 'disc_transmitted_blr_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7a920>, 'disc_incident_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b3d0>, 'full_reprocessed_blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b220>, 'full_reprocessed_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79cc0>, 'disc_transmitted_blr_isotropic_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b5b0>, 'disc_transmitted_nlr_isotropic_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b010>, 'full_continuum_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7af50>, 'disc_incident': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7bd60>, 'full_continuum_blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc78df0>, 'disc_escaped': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b400>, 'disc_transmitted_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7a230>, 'disc_transmitted_blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7af20>, 'disc_transmitted': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79a20>, 'disc': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79f30>, 'torus': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79900>, 'blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecd62a70>, 'nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecd62800>, 'line_regions': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7e440>, 'nlr_continuum': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7bb50>, 'disc_escaped_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7ef20>, 'disc_transmitted_blr_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7f9d0>, 'disc_transmitted_nlr_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7e3b0>, 'disc_averaged_without_torus': <synthesizer.emissions.sed.Sed object at 0x7f77ecd63850>, 'disc_averaged': <synthesizer.emissions.sed.Sed object at 0x7f77ecd63550>, 'disc_transmitted_nlr_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7f730>, 'disc_transmitted_blr_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7faf0>, 'disc_escaped_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7d060>, 'disc_transmitted_weighted_combination': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7cca0>, 'blr_continuum': <synthesizer.emissions.sed.Sed object at 0x7f77ecd62590>, 'intrinsic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd63ca0>, 'attenuated': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7fc70>}
Integrating spectra¶
The integrated spectra are automatically produced alongside per particle spectra. However, if we wanted to explictly get the integrated spectra from the particle spectra we just generated (for instance if we had made some modification after generation), we can call the integrate_particle_spectra method. This method will sum the individual spectra and populate the spectra dictionary.
Note, we can also integrate individual spectra using the Sed.sum() method.
[12]:
print(blackholes.spectra)
blackholes.integrate_particle_spectra()
print(blackholes.spectra)
fig, ax = blackholes.plot_spectra(
show=True, ylimits=(10**28.5, 10**34.0), figsize=(9, 7.2)
)
{'disc_incident_masked': <synthesizer.emissions.sed.Sed object at 0x7f77eee90850>, 'disc_transmitted_nlr_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79f60>, 'disc_transmitted_blr_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79d20>, 'disc_incident_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7ab00>, 'full_reprocessed_blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc78ac0>, 'full_reprocessed_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7be20>, 'disc_transmitted_blr_isotropic_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79d50>, 'disc_transmitted_nlr_isotropic_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b040>, 'full_continuum_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b430>, 'disc_incident': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7bb80>, 'full_continuum_blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc78a30>, 'disc_escaped': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79c00>, 'disc_transmitted_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b460>, 'disc_transmitted_blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b190>, 'disc_transmitted': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7a080>, 'disc': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b340>, 'torus': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7a320>, 'blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b730>, 'nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7bf40>, 'line_regions': <synthesizer.emissions.sed.Sed object at 0x7f77ecd602e0>, 'nlr_continuum': <synthesizer.emissions.sed.Sed object at 0x7f77ecd61330>, 'disc_escaped_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7ed70>, 'disc_transmitted_blr_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd63d00>, 'disc_transmitted_nlr_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd63670>, 'disc_averaged_without_torus': <synthesizer.emissions.sed.Sed object at 0x7f77ecd63400>, 'disc_averaged': <synthesizer.emissions.sed.Sed object at 0x7f77ecd62290>, 'disc_transmitted_nlr_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77ecd635e0>, 'disc_transmitted_blr_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77ecd61e10>, 'disc_escaped_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b580>, 'disc_transmitted_weighted_combination': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7c8b0>, 'blr_continuum': <synthesizer.emissions.sed.Sed object at 0x7f77ecc78490>, 'intrinsic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7cc70>, 'attenuated': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7ef80>}
{'disc_incident_masked': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7fb50>, 'disc_transmitted_nlr_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecd633d0>, 'disc_transmitted_blr_full': <synthesizer.emissions.sed.Sed object at 0x7f77eee91060>, 'disc_incident_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77eee90b50>, 'full_reprocessed_blr': <synthesizer.emissions.sed.Sed object at 0x7f77eee90850>, 'full_reprocessed_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77eee90ca0>, 'disc_transmitted_blr_isotropic_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecb17490>, 'disc_transmitted_nlr_isotropic_full': <synthesizer.emissions.sed.Sed object at 0x7f77ecb17e20>, 'full_continuum_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecb17e50>, 'disc_incident': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79d50>, 'full_continuum_blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b430>, 'disc_escaped': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7bb80>, 'disc_transmitted_nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc78a30>, 'disc_transmitted_blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc79c00>, 'disc_transmitted': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b460>, 'disc': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b190>, 'torus': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7a080>, 'blr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b340>, 'nlr': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7a320>, 'line_regions': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b730>, 'nlr_continuum': <synthesizer.emissions.sed.Sed object at 0x7f77ecd602e0>, 'disc_escaped_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd61330>, 'disc_transmitted_blr_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7ed70>, 'disc_transmitted_nlr_isotropic': <synthesizer.emissions.sed.Sed object at 0x7f77ecd63d00>, 'disc_averaged_without_torus': <synthesizer.emissions.sed.Sed object at 0x7f77ecd63670>, 'disc_averaged': <synthesizer.emissions.sed.Sed object at 0x7f77ecd63400>, 'disc_transmitted_nlr_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77ecd62290>, 'disc_transmitted_blr_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77ecd635e0>, 'disc_escaped_weighted': <synthesizer.emissions.sed.Sed object at 0x7f77ecd61e10>, 'disc_transmitted_weighted_combination': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7bf40>, 'blr_continuum': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7c8b0>, 'intrinsic': <synthesizer.emissions.sed.Sed object at 0x7f77ecc7b580>, 'attenuated': <synthesizer.emissions.sed.Sed object at 0x7f77ecd7cc70>}
Printing Used Parameters¶
During spectra generation, emission models cache the parameters they extract and use from the emitter. These cached parameters can be printed in a nicely formatted table to inspect which values were actually used by each model.
[13]:
# Print the cached parameters used by the models
blackholes.print_used_parameters()
+----------------------------------------------------------------------------------------------+
| MODEL: DISC_INCIDENT_MASKED |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 1.00 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.10 |
+------------------------------------+---------------------------------------------------------+
| extract | 'incident' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| masks | 'torus_edgeon_cond < 90 degree' |
+------------------------------------+---------------------------------------------------------+
| inclination | 0.0 degree |
+------------------------------------+---------------------------------------------------------+
| theta_torus | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| torus_edgeon_cond | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| mass (4,) | 5.26e+07 Msun -> 1.55e+08 Msun (Mean: 9.79e+07 Msun) |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| accretion_rate_eddington (4,) | 2.02e+32 -> 8.87e+32 (Mean: 5.18e+32) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000000000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_NLR_FULL |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 1.00 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.01 |
+------------------------------------+---------------------------------------------------------+
| extract | 'transmitted' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| masks | 'torus_edgeon_cond < 90 degree' |
+------------------------------------+---------------------------------------------------------+
| inclination | 0.0 degree |
+------------------------------------+---------------------------------------------------------+
| theta_torus | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| torus_edgeon_cond | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_BLR_FULL |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 1.00 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.10 |
+------------------------------------+---------------------------------------------------------+
| extract | 'transmitted' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| masks | 'torus_edgeon_cond < 90 degree' |
+------------------------------------+---------------------------------------------------------+
| inclination | 0.0 degree |
+------------------------------------+---------------------------------------------------------+
| theta_torus | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| torus_edgeon_cond | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000000000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: DISC_INCIDENT_ISOTROPIC |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 0.50 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.10 |
+------------------------------------+---------------------------------------------------------+
| extract | 'incident' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000000000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: FULL_REPROCESSED_BLR |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 0.50 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.10 |
+------------------------------------+---------------------------------------------------------+
| extract | 'nebular' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| masks | 'torus_edgeon_cond < 90 degree' |
+------------------------------------+---------------------------------------------------------+
| inclination | 0.0 degree |
+------------------------------------+---------------------------------------------------------+
| theta_torus | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| torus_edgeon_cond | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000000000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: FULL_REPROCESSED_NLR |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 0.50 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.01 |
+------------------------------------+---------------------------------------------------------+
| extract | 'nebular' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_BLR_ISOTROPIC_FULL |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 0.50 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.10 |
+------------------------------------+---------------------------------------------------------+
| extract | 'transmitted' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000000000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_NLR_ISOTROPIC_FULL |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 0.50 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.01 |
+------------------------------------+---------------------------------------------------------+
| extract | 'transmitted' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: FULL_CONTINUUM_NLR |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 0.50 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.01 |
+------------------------------------+---------------------------------------------------------+
| extract | 'nebular_continuum' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: DISC_INCIDENT |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 1.00 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.10 |
+------------------------------------+---------------------------------------------------------+
| extract | 'incident' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000000000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+----------------------------------------------------------------------------------------------+
| MODEL: FULL_CONTINUUM_BLR |
+------------------------------------+---------------------------------------------------------+
| Attribute | Value |
+------------------------------------+---------------------------------------------------------+
| cosine_inclination | 0.50 |
+------------------------------------+---------------------------------------------------------+
| ionisation_parameter | 0.10 |
+------------------------------------+---------------------------------------------------------+
| extract | 'nebular_continuum' |
+------------------------------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------------+---------------------------------------------------------+
| masks | 'torus_edgeon_cond < 90 degree' |
+------------------------------------+---------------------------------------------------------+
| inclination | 0.0 degree |
+------------------------------------+---------------------------------------------------------+
| theta_torus | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| torus_edgeon_cond | 10.0 degree |
+------------------------------------+---------------------------------------------------------+
| log10mass (4,) | 7.72e+00 -> 8.19e+00 (Mean: 7.95e+00) |
+------------------------------------+---------------------------------------------------------+
| log10accretion_rate_eddington (4,) | 3.23e+01 -> 3.29e+01 (Mean: 3.27e+01) |
+------------------------------------+---------------------------------------------------------+
| metallicities (4,) | 1.00e-02 -> 1.00e-02 (Mean: 1.00e-02) |
+------------------------------------+---------------------------------------------------------+
| hydrogen_density | 10000000000.0 cm**(-3) |
+------------------------------------+---------------------------------------------------------+
| bolometric_luminosities (4,) | 7.94e+44 erg/s -> 2.96e+45 erg/s (Mean: 1.53e+45 erg/s) |
+------------------------------------+---------------------------------------------------------+
+-----------------------------------------------------------------------------------------------------+
| MODEL: DISC_ESCAPED |
+------------------------------+----------------------------------------------------------------------+
| Attribute | Value |
+------------------------------+----------------------------------------------------------------------+
| transmission_fraction_escape | 1.00 |
+------------------------------+----------------------------------------------------------------------+
| apply_to | 'disc_incident_masked' |
+------------------------------+----------------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('transmission_fraction_escape',))" |
+------------------------------+----------------------------------------------------------------------+
| emitter | 'blackhole' |
+------------------------------+----------------------------------------------------------------------+
+-----------------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_NLR |
+---------------------------+-------------------------------------------------------------------+
| Attribute | Value |
+---------------------------+-------------------------------------------------------------------+
| transmission_fraction_nlr | 0.00e+00 |
+---------------------------+-------------------------------------------------------------------+
| apply_to | 'disc_transmitted_nlr_full' |
+---------------------------+-------------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('transmission_fraction_nlr',))" |
+---------------------------+-------------------------------------------------------------------+
| emitter | 'blackhole' |
+---------------------------+-------------------------------------------------------------------+
+-----------------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_BLR |
+---------------------------+-------------------------------------------------------------------+
| Attribute | Value |
+---------------------------+-------------------------------------------------------------------+
| transmission_fraction_blr | 0.00e+00 |
+---------------------------+-------------------------------------------------------------------+
| apply_to | 'disc_transmitted_blr_full' |
+---------------------------+-------------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('transmission_fraction_blr',))" |
+---------------------------+-------------------------------------------------------------------+
| emitter | 'blackhole' |
+---------------------------+-------------------------------------------------------------------+
+------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED |
+----------+-------------------------------------------------------------+
| Attribute | Value |
+----------+-------------------------------------------------------------+
| emitter | 'blackhole' |
+----------+-------------------------------------------------------------+
| combine | [disc_escaped, disc_transmitted_nlr, disc_transmitted_blr,] |
+----------+-------------------------------------------------------------+
+---------------------------------+
| MODEL: DISC |
+----------+----------------------+
| Attribute | Value |
+----------+----------------------+
| emitter | 'blackhole' |
+----------+----------------------+
| combine | [disc_transmitted, ] |
+----------+----------------------+
+-------------------------------------------------------------------------------------------------------------------------------------+
| MODEL: TORUS |
+-------------+-----------------------------------------------------------------------------------------------------------------------+
| Attribute | Value |
+-------------+-----------------------------------------------------------------------------------------------------------------------+
| emissivity | 1.50 |
+-------------+-----------------------------------------------------------------------------------------------------------------------+
| generator | 'Greybody(scaler=disc_incident_isotropic, temperature=1000.0 K, emissivity=1.5, optically_thin=True, lam_0=100.0 μm)' |
+-------------+-----------------------------------------------------------------------------------------------------------------------+
| emitter | 'blackhole' |
+-------------+-----------------------------------------------------------------------------------------------------------------------+
| temperature | 1000.0 K |
+-------------+-----------------------------------------------------------------------------------------------------------------------+
+---------------------------------------------------------------------------------------+
| MODEL: BLR |
+-----------------------+---------------------------------------------------------------+
| Attribute | Value |
+-----------------------+---------------------------------------------------------------+
| covering_fraction_blr | 0.10 |
+-----------------------+---------------------------------------------------------------+
| apply_to | 'full_reprocessed_blr' |
+-----------------------+---------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('covering_fraction_blr',))" |
+-----------------------+---------------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------------+---------------------------------------------------------------+
+---------------------------------------------------------------------------------------+
| MODEL: NLR |
+-----------------------+---------------------------------------------------------------+
| Attribute | Value |
+-----------------------+---------------------------------------------------------------+
| covering_fraction_nlr | 0.10 |
+-----------------------+---------------------------------------------------------------+
| apply_to | 'full_reprocessed_nlr' |
+-----------------------+---------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('covering_fraction_nlr',))" |
+-----------------------+---------------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------------+---------------------------------------------------------------+
+------------------------+
| MODEL: LINE_REGIONS |
+----------+-------------+
| Attribute | Value |
+----------+-------------+
| emitter | 'blackhole' |
+----------+-------------+
| combine | [nlr, blr] |
+----------+-------------+
+---------------------------------------------------------------------------------------+
| MODEL: NLR_CONTINUUM |
+-----------------------+---------------------------------------------------------------+
| Attribute | Value |
+-----------------------+---------------------------------------------------------------+
| covering_fraction_nlr | 0.10 |
+-----------------------+---------------------------------------------------------------+
| apply_to | 'full_continuum_nlr' |
+-----------------------+---------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('covering_fraction_nlr',))" |
+-----------------------+---------------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------------+---------------------------------------------------------------+
+---------------------------------------------------------------------------------------------------------------+
| MODEL: DISC_ESCAPED_ISOTROPIC |
+-----------------------+---------------------------------------------------------------------------------------+
| Attribute | Value |
+-----------------------+---------------------------------------------------------------------------------------+
| covering_fraction_blr | 0.10 |
+-----------------------+---------------------------------------------------------------------------------------+
| covering_fraction_nlr | 0.10 |
+-----------------------+---------------------------------------------------------------------------------------+
| apply_to | 'disc_incident_isotropic' |
+-----------------------+---------------------------------------------------------------------------------------+
| transformer | "EscapingFraction(covering_attrs=('covering_fraction_blr', 'covering_fraction_nlr'))" |
+-----------------------+---------------------------------------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------------+---------------------------------------------------------------------------------------+
+---------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_BLR_ISOTROPIC |
+-----------------------+---------------------------------------------------------------+
| Attribute | Value |
+-----------------------+---------------------------------------------------------------+
| covering_fraction_blr | 0.10 |
+-----------------------+---------------------------------------------------------------+
| apply_to | 'disc_transmitted_blr_isotropic_full' |
+-----------------------+---------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('covering_fraction_blr',))" |
+-----------------------+---------------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------------+---------------------------------------------------------------+
+---------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_NLR_ISOTROPIC |
+-----------------------+---------------------------------------------------------------+
| Attribute | Value |
+-----------------------+---------------------------------------------------------------+
| covering_fraction_nlr | 0.10 |
+-----------------------+---------------------------------------------------------------+
| apply_to | 'disc_transmitted_nlr_isotropic_full' |
+-----------------------+---------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('covering_fraction_nlr',))" |
+-----------------------+---------------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------------+---------------------------------------------------------------+
+---------------------------------------------------------------------+
| MODEL: DISC_AVERAGED_WITHOUT_TORUS |
+----------+----------------------------------------------------------+
| Attribute | Value |
+----------+----------------------------------------------------------+
| emitter | 'blackhole' |
+----------+----------------------------------------------------------+
| combine | [disc_escaped_isotropic, disc_transmitted_nlr_isotropic, |
| | disc_transmitted_blr_isotropic, ] |
+----------+----------------------------------------------------------+
+-------------------------------------------------------------------------+
| MODEL: DISC_AVERAGED |
+----------------+--------------------------------------------------------+
| Attribute | Value |
+----------------+--------------------------------------------------------+
| apply_to | 'disc_averaged_without_torus' |
+----------------+--------------------------------------------------------+
| transformer | "EscapingFraction(covering_attrs=('torus_fraction',))" |
+----------------+--------------------------------------------------------+
| emitter | 'blackhole' |
+----------------+--------------------------------------------------------+
| torus_fraction | 0.1111111111111111 |
+----------------+--------------------------------------------------------+
+---------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_NLR_WEIGHTED |
+-----------------------+---------------------------------------------------------------+
| Attribute | Value |
+-----------------------+---------------------------------------------------------------+
| covering_fraction_nlr | 0.10 |
+-----------------------+---------------------------------------------------------------+
| apply_to | 'disc_transmitted_nlr_full' |
+-----------------------+---------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('covering_fraction_nlr',))" |
+-----------------------+---------------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------------+---------------------------------------------------------------+
+---------------------------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_BLR_WEIGHTED |
+-----------------------+---------------------------------------------------------------+
| Attribute | Value |
+-----------------------+---------------------------------------------------------------+
| covering_fraction_blr | 0.10 |
+-----------------------+---------------------------------------------------------------+
| apply_to | 'disc_transmitted_blr_full' |
+-----------------------+---------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('covering_fraction_blr',))" |
+-----------------------+---------------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------------+---------------------------------------------------------------+
+---------------------------------------------------------------------------+
| MODEL: DISC_ESCAPED_WEIGHTED |
+-----------------+---------------------------------------------------------+
| Attribute | Value |
+-----------------+---------------------------------------------------------+
| escape_fraction | 0.80 |
+-----------------+---------------------------------------------------------+
| apply_to | 'disc_incident_masked' |
+-----------------+---------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('escape_fraction',))" |
+-----------------+---------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------+---------------------------------------------------------+
+-------------------------------------------------------------------+
| MODEL: DISC_TRANSMITTED_WEIGHTED_COMBINATION |
+----------+--------------------------------------------------------+
| Attribute | Value |
+----------+--------------------------------------------------------+
| emitter | 'blackhole' |
+----------+--------------------------------------------------------+
| combine | [disc_escaped_weighted, disc_transmitted_nlr_weighted, |
| | disc_transmitted_blr_weighted, ] |
+----------+--------------------------------------------------------+
+---------------------------------------------------------------------------------------+
| MODEL: BLR_CONTINUUM |
+-----------------------+---------------------------------------------------------------+
| Attribute | Value |
+-----------------------+---------------------------------------------------------------+
| covering_fraction_blr | 0.10 |
+-----------------------+---------------------------------------------------------------+
| apply_to | 'full_continuum_blr' |
+-----------------------+---------------------------------------------------------------+
| transformer | "CoveringFraction(covering_attrs=('covering_fraction_blr',))" |
+-----------------------+---------------------------------------------------------------+
| emitter | 'blackhole' |
+-----------------------+---------------------------------------------------------------+
+------------------------------------+
| MODEL: INTRINSIC |
+----------+-------------------------+
| Attribute | Value |
+----------+-------------------------+
| emitter | 'blackhole' |
+----------+-------------------------+
| combine | [disc, nlr, blr, torus] |
+----------+-------------------------+
+--------------------------------------+
| MODEL: ATTENUATED |
+-------------+------------------------+
| Attribute | Value |
+-------------+------------------------+
| tau_v | 0.50 |
+-------------+------------------------+
| apply_to | 'intrinsic' |
+-------------+------------------------+
| transformer | 'PowerLaw(slope=-1.0)' |
+-------------+------------------------+
| emitter | 'blackhole' |
+-------------+------------------------+