skip to main content

Title: Variability-selected Intermediate-mass Black Hole Candidates in Dwarf Galaxies from ZTF and WISE
Abstract While it is difficult to observe the first black hole seeds in the early universe, we can study intermediate-mass black holes (IMBHs) in local dwarf galaxies for clues about their origins. In this paper we present a sample of variability-selected active galactic nuclei (AGN) in dwarf galaxies using optical photometry from the Zwicky Transient Facility (ZTF) and forward-modeled mid-IR photometry of time-resolved Wide-field Infrared Survey Explorer (WISE) co-added images. We found that 44 out of 25,714 dwarf galaxies had optically variable AGN candidates and 148 out of 79,879 dwarf galaxies had mid-IR variable AGN candidates, corresponding to active fractions of 0.17% ± 0.03% and 0.19% ± 0.02%, respectively. We found that spectroscopic approaches to AGN identification would have missed 81% of our ZTF IMBH candidates and 69% of our WISE IMBH candidates. Only nine candidates have been detected previously in radio, X-ray, and variability searches for dwarf galaxy AGN. The ZTF and WISE dwarf galaxy AGN with broad Balmer lines have virial masses of 10 5 M ⊙ < M BH < 10 7 M ⊙ , but for the rest of the sample, BH masses predicted from host galaxy mass range between 10 5.2 M ⊙ < M BH < 10 7.25 M ⊙ . We found that only 5 of 152 previously reported variability-selected AGN candidates from the Palomar Transient Factory in common with our parent sample were variable in ZTF. We also determined a nuclear supernova fraction of 0.05% ± 0.01% yr −1 for dwarf galaxies in ZTF. Our ZTF and WISE IMBH candidates show the promise of variability searches for the discovery of otherwise hidden low-mass AGN.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We present a sample of 706, z < 1.5 active galactic nuclei (AGNs) selected from optical photometric variability in three of the Dark Energy Survey (DES) deep fields (E2, C3, and X3) over an area of 4.64 deg2. We construct light curves using difference imaging aperture photometry for resolved sources and non-difference imaging PSF photometry for unresolved sources, respectively, and characterize the variability significance. Our DES light curves have a mean cadence of 7 d, a 6-yr baseline, and a single-epoch imaging depth of up to g ∼ 24.5. Using spectral energy distribution (SED) fitting, we find 26 out of total 706 variable galaxies are consistent with dwarf galaxies with a reliable stellar mass estimate ($M_{\ast }\lt 10^{9.5}\, {\rm M}_\odot$; median photometric redshift of 0.9). We were able to constrain rapid characteristic variability time-scales (∼ weeks) using the DES light curves in 15 dwarf AGN candidates (a subset of our variable AGN candidates) at a median photometric redshift of 0.4. This rapid variability is consistent with their low black hole (BH) masses. We confirm the low-mass AGN nature of one source with a high S/N optical spectrum. We publish our catalogue, optical light curves, and supplementary data, such as X-ray properties and optical spectra, when available. We measure a variable AGN fraction versus stellar mass and compare to results from a forward model. This work demonstrates the feasibility of optical variability to identify AGNs with lower BH masses in deep fields, which may be more ‘pristine’ analogues of supermassive BH seeds. 
    more » « less

    We present a phenomenological forward Monte Carlo model for forecasting the population of active galactic nuclei (AGNs) in dwarf galaxies observable via their optical variability. Our model accounts for expected changes in the spectral energy distribution of AGNs in the intermediate-mass black hole (IMBH) mass range and uses observational constraints on optical variability as a function of black hole (BH) mass to generate mock light curves. Adopting several different models for the BH occupation function, including one for off-nuclear IMBHs, we quantify differences in the predicted local AGN mass and luminosity functions in dwarf galaxies. As a result, we are able to model the fraction of variable AGNs as a function of important galaxy host properties, such as host galaxy stellar mass, in the presence of selection effects. We find that our adopted occupation fractions for the ‘heavy’ and ‘light’ initial BH seeding scenarios can be distinguished with variability at the 2–3σ level for galaxy host stellar masses below ∼108M⊙ with data from the upcoming Vera C. Rubin Observatory. We also demonstrate the prevalence of a selection bias whereby recovered IMBH masses fall, on average, above the predicted value from the local host galaxy–BH mass scaling relation with the strength of this bias dependent on the survey sensitivity. Our methodology can be used more broadly to calibrate AGN demographic studies in synoptic surveys. Finally, we show that a targeted ∼ hourly cadence program over a few nights with the Rubin Observatory can provide strong constraints on IMBH masses given their expected rapid variability time-scales.

    more » « less

    We study the optical variability of a sample of candidate low-mass (dwarf and Seyfert) active galactic nuclei (AGNs) using Zwicky Transient Facility g-band light curves. Our sample is compiled from broad-line AGNs in dwarf galaxies reported in the literature with single-epoch virial black hole (BH) masses in the range MBH ∼ 104–108 M⊙. We measure the characteristic ‘damping’ time-scale of the optical variability τDRW, beyond which the power spectral density flattens, of a final sample of 79 candidate low-mass AGNs with high-quality light curves. Our results provide further confirmation of the MBH–τDRW relation from previous work within 1σ agreement, adding 78 new low-mass AGNs to the relation. The agreement suggests that the virial BH mass estimates for these AGNs are generally reasonable. We expect that the optical light curve of an accreting intermediate-mass black hole (IMBH) to vary with a rest-frame damping time-scale of ∼tens of hours, which could enable detection and direct mass estimation of accreting IMBHs in wide-field time-domain imaging surveys with sufficient cadence like with the Vera C. Rubin Observatory.

    more » « less

    We used Transiting Exoplanet Survey Satellite (TESS) data to identify 29 candidate active galactic nuclei (AGNs) through their optical variability. The high-cadence, high-precision TESS light curves present an opportunity for the identification of AGNs, including those not selected through other methods. Of the candidates, we found that 18 have either previously been identified as AGNs in the literature or could have been selected based on emission-line diagnostics, mid-IR colours, or X-ray luminosity. AGNs in low-mass galaxies offer a unique window into supermassive black hole and galaxy co-evolution and 8 of the 29 candidates have estimated black hole masses ≲ 106 M⊙. The low-mass galaxies NGC 4395 and NGC 4449 are two of our five ‘high-confidence’ candidates. Since our initial sample largely draws from just nine TESS sectors, we expect to identify at least ∼45 more candidates in the TESS primary and extended mission data sets, of which ∼60 per cent will be new AGNs and ∼20 per cent will be in low-mass galaxies.

    more » « less
  5. Abstract Current observational facilities have yet to conclusively detect 10 3 –10 4 M ⊙ intermediate-mass black holes (IMBHs) that fill in the evolutionary gap between seed black holes in the early universe and z ∼ 0 supermassive black holes. Dwarf galaxies present an opportunity to reveal active IMBHs amidst persistent star formation. We introduce photoionization simulations tailored to address key physical uncertainties: coincident versus noncoincident mixing of IMBH and starlight excitation, open versus closed geometries of surrounding gas clouds, and different shapes of the spectral energy distribution of active galactic nuclei (AGN). We examine possible AGN emission line diagnostics in the optical and mid-IR, and find that the diagnostics are often degenerate with respect to the investigated physical uncertainties. In spite of these setbacks, and in contrast to recent work, we are able to show that [O iii ]/H β typically remains bright for dwarf AGN powered by IMBHs down to 10 3 M ⊙ . Dwarf AGN are predicted to have inconsistent star-forming and Seyfert/LINER classifications using the most common optical diagnostics. In the mid-IR, [O iv ] 25.9 μ m and [Ar ii ] 6.98 μ m are less sensitive to physical uncertainties than are optical diagnostics. Based on these emission lines, we provide several diagrams of mid-IR emission line diagnostic diagrams with demarcations for separating starbursts and AGN with varying levels of activity. The diagrams are valid over a wide range of ionization parameters and metallicities out to z ∼ 0.1, so will prove useful for future JWST observations of local dwarf AGN in the search for IMBHs. We make our photoionization simulation suite freely available. 
    more » « less