lsst_inaf_agile.mbh

Tools for assigning MBH.

Attributes

GRAHAM2023_ROWS
X
Y
Z
SPLINE
FUN_LOG_MBH_CONTINUITY_NEW2
REFERENCES
log_mstar

Functions

get_log_mbh_shankar2019(log_mstar)	Return equation 5 from Shankar+2019.
get_log_mbh_sg2016(log_mstar)	Return equation 3 from Shankar+2019.
get_log_mbh_tanaka2024(log_mstar)	Return equation 10 from Tanaka+2024, see also Figure 12.
get_log_mbh_pacucci2023(log_mstar)	Return the first equation from Pacucci+2023 abstract.
get_log_mbh_suh2020(log_mstar[, alpha, beta, dalpha, ...])	Return the first equation of Suh+2020 Sec. 5.
get_log_mbh_reines_volonteri2015(log_mstar[, alpha, beta])	Return the first equation of Reines & Volonteri 2015 Sec. 4.1.
get_log_mbh_kormendy_ho2013(log_mstar[, alpha, beta, ...])	Return Kormendy & Ho 2013.
get_log_mbh_haring_rix2004(log_mstar[, a, b, da, db, ...])	Return Häring & Rix from Angela.
get_log_mbh_const(log_mstar[, A])	Return logMBH according to Haring+Rix.
get_log_mbh_zou2024(log_mstar[, z, is_tng300])	Return Zou+2024 MBH-Mstar(z).
get_log_mbh_graham2023(log_mstar[, which])	Return a MBH relation from Graham+2023 Table 1: https://arxiv.org/pdf/2305.03242
get_delta_log_mbh_shankar2019(log_mstar[, low, high])	Return equation 6 from Shankar+2019.
get_spline(Mst, Mbh)	Take the Mstar and Mbh values at z=0 and creates a cubic spline.
get_log_mbh_continuity_new2(log_mstar, z[, filename])	Return logMBH according to the continuity equation with Zou+24 p(lambda).
get_log_mbh(log_mstar[, reference, z])	Return logMBH according to the scaling relation from the reference.
get_occupation_fraction(logMstar)	Return the black hole occupation fraction.

Module Contents

get_log_mbh_shankar2019(log_mstar)

Return equation 5 from Shankar+2019.

get_log_mbh_sg2016(log_mstar)

Return equation 3 from Shankar+2019.

get_log_mbh_tanaka2024(log_mstar)

Return equation 10 from Tanaka+2024, see also Figure 12.

get_log_mbh_pacucci2023(log_mstar)

Return the first equation from Pacucci+2023 abstract.

get_log_mbh_suh2020(log_mstar, alpha=1.64, beta=10.29, dalpha=0.07, dbeta=0.04, add_error=False)

Return the first equation of Suh+2020 Sec. 5.

get_log_mbh_reines_volonteri2015(log_mstar, alpha=7.45, beta=1.05)

Return the first equation of Reines & Volonteri 2015 Sec. 4.1.

get_log_mbh_kormendy_ho2013(log_mstar, alpha=8.56, beta=1.58, from_rv15=False)

Return Kormendy & Ho 2013.

If from rv15 = True, then return the “scaled” K&H13 relation from their Fig. 8

get_log_mbh_haring_rix2004(log_mstar, a=8.2, b=1.12, da=0.1, db=0.06, add_error=False)

Return Häring & Rix from Angela.

get_log_mbh_const(log_mstar, A=500)

Return logMBH according to Haring+Rix.

get_log_mbh_zou2024(log_mstar, z=0.0, is_tng300=False)

Return Zou+2024 MBH-Mstar(z).

The MBH-Mstar is based on a hybrid approach using the observed BHAR from earlier Zou work combined with mergers from IllustrisTNG. The Mbh-Mstar(z) is then built by following seeded black hole masses starting from z=4.0.

Refer to Fig. 4 and Table 1 of Zou+2024: https://ui.adsabs.harvard.edu/abs/2024ApJ…976….6Z

GRAHAM2023_ROWS = ['E/Es,e', 'E/ES,e/S0(dust=Y)', 'BCG∗', 'major_mergers', 'S', 'S_(w/o_Circinus)', 'S0/Es,b_(dust=Y)']

get_log_mbh_graham2023(log_mstar, which='E/Es,e')

Return a MBH relation from Graham+2023 Table 1: https://arxiv.org/pdf/2305.03242

get_delta_log_mbh_shankar2019(log_mstar, low=0.0, high=np.inf)

Return equation 6 from Shankar+2019.

The reference for the relation is: Shankar F. et al., 2016, MNRAS, 460, 3119

get_spline(Mst, Mbh)

Take the Mstar and Mbh values at z=0 and creates a cubic spline.

The cubic spline extrapolates with dy/dx of the final knot and d2y/dx2 = 0.

X = None

Y = None

Z = None

SPLINE: dict[float, scipy.interpolate.CubicSpline]

FUN_LOG_MBH_CONTINUITY_NEW2

get_log_mbh_continuity_new2(log_mstar, z, filename='data/AGILE_Mstar_Mbh/DECODEmeans_Mstar_Mbh_Zou2024cc.dat')

Return logMBH according to the continuity equation with Zou+24 p(lambda).

Steps to generate the required input file: 1) run opt/AGILE_Mstar_Mbh/Modules/ERDFs/genERDFdata.py 2) run opt/AGILE_Mstar_Mbh/DECODEmeans.py

Zou2024c refers to the case of assuming Zou+2024 from the parameter maps and extrapolating above z>4 and logM<9.5.

The extrapolation scheme assumed is simple boundary extrapolation where z and logM are clipped to their minimum and maximum values.

REFERENCES = ('H&R04', 'K&H13', 'R&V15', 'S&G16', 'Shankar+16', 'Suh+19', 'Pacucci+23', 'Tanaka+24',...

get_log_mbh(log_mstar, reference='Shankar+16', z=0.0)

Return logMBH according to the scaling relation from the reference.

get_occupation_fraction(logMstar)

Return the black hole occupation fraction.

The occupation fraction is the percentage of host galaxies that are expected to host SMBH, which could be below unity for dwarf galaxies (Miller+15, Burke+25, Zou+). The updated AGN duty cycle is then the product occupation fraction and the previous AGN duty cycle.

Parameters:

logMstar (float or array_like) – Log10 of host galaxy stellar mass in Msun.

log_mstar