NICO RAMOS

Research · Astrophysics research

Estimating Black Hole Masses for Type 1 AGNs

A Python pipeline that estimates black hole masses for about 5,000 active galactic nuclei from SDSS spectra.

My mass estimates against those from Shen et al. 2011
My mass estimates against those from Shen et al. 2011

Background

Active Galactic Nuclei (AGN) are among the most massive and energetic objects in the universe, which makes them key to understanding galactic morphology. One of the fundamental values when characterizing an AGN is its mass.

Methods for getting an AGN's mass are based on the virial theorem and need two values:

AGN structure (image from Ricci 2011, PhD thesis)
AGN structure (image from Ricci 2011, PhD thesis)

The most rigorous of these methods, reverberation mapping, can take years to get the first of those two values for just one object. That's why quicker and more automatable (if less rigorous) methods are desirable.

One such method, the one I use, uses a correlation between the monochromatic luminosity at 5100 Å and Rorbit to estimate Rorbit.

What I did

Starting in June 2023, I wrote a Python script to estimate black hole masses from SDSS spectral data, using a Python package called PyQSOFit (Guo, Shen & Wang 2018). With PyQSOFit I obtained the FWHM of the Hβ emission line and the monochromatic luminosity at 5100 Å for about 5,000 AGN from the catalog described in Souchay et al. 2019.

Example of a Type 1 AGN Hβ emission fit from PyQSOFit
Example of a Type 1 AGN Hβ emission fit from PyQSOFit

Then, using the empirical relationships in Vestergaard & Peterson (2006), I estimated the masses of about 5,000 black holes.

Results

The plot below compares several hundred of these masses with those estimated by Shen et al. The "clean results" are the objects with an Hβ FWHM above 1,000 km/s (Type 1 AGN) and a signal-to-noise ratio above 10.

My results vs those in Shen et al. 2011
My results vs those in Shen et al. 2011