skip to main content


Title: Modeling the Infrared Reverberation Response of the Circumnuclear Dusty Torus in AGNs: An Investigation of Torus Response Functions
The size and structure of the dusty circumnuclear torus in active galactic nuclei (AGNs) can be investigated by analyzing the temporal response of the torus's infrared (IR) dust emission to variations in the AGN ultraviolet/optical luminosity. This method, reverberation mapping, is applicable over a wide redshift range, but the IR response is sensitive to several poorly constrained variables relating to the dust distribution and its illumination, complicating the interpretation of measured reverberation lags. We have used an enhanced version of our torus reverberation mapping code (TORMAC) to conduct a comprehensive exploration of the torus response functions at selected wavelengths, for the standard interstellar medium grain composition. The shapes of the response functions vary widely over the parameter range covered by our models, with the largest variations occurring at shorter wavelengths (≤4.5 μm). The reverberation lag, quantified as the response-weighted delay (RWD), is most affected by the radial depth of the torus, the steepness of the radial cloud distribution, the degree of anisotropy of the AGN radiation field, and the volume filling factor. Nevertheless, we find that the RWD provides a reasonably robust estimate, to within a factor of ~3, of the luminosity-weighted torus radius, confirming the basic assumption underlying reverberation mapping. However, overall, the models predict radii at 2.2 μm that are typically a factor of ~2 larger than those derived from K-band reverberation mapping. This is likely an indication that the innermost region of the torus is populated by clouds dominated by large graphite grains.  more » « less
Award ID(s):
1659740
NSF-PAR ID:
10227733
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
The astrophysical journal
Volume:
891
Issue:
1
ISSN:
2041-8205
Page Range / eLocation ID:
26
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We have recently initiated the first spectroscopic dust reverberation programme on active galactic nuclei in the near-infrared. Spectroscopy enables measurement of dust properties, such as flux, temperature, and covering factor, with higher precision than photometry. In particular, it enables measurement of both luminosity-based dust radii and dust response times. Here we report results from a 1 yr campaign on NGC 5548. The hot dust responds to changes in the irradiating flux with a lag time of ∼70 light-days, similar to what was previously found in photometric reverberation campaigns. The mean and rms spectra are similar, implying that the same dust component dominates both the emission and the variations. The dust lag time is consistent with the luminosity-based dust radius only if we assume a wavelength-independent dust emissivity law, i.e. a blackbody, which is appropriate for grains of large sizes (of a few μm). For such grains the dust temperature is ∼1450 K. Therefore, silicate grains have most likely evaporated and carbon is the main chemical component. But the hot dust is not close to its sublimation temperature, contrary to popular belief. This is further supported by our observation of temperature variations largely consistent with a heating/cooling process. Therefore, the inner dust-free region is enlarged and the dusty torus rather a ‘dusty wall’, whose inner radius is expected to be luminosity-invariant. The dust-destruction mechanism that enlarges the dust-free region seems to also partly affect the dusty region. We observe a cyclical decrease in dust mass with implied dust reformation times of ∼5–6 months. 
    more » « less
  2. Abstract

    The variable continuum emission of an active galactic nucleus (AGN) produces corresponding responses in the broad emission lines, which are modulated by light travel delays, and contain information on the physical properties, structure, and kinematics of the emitting gas region. The reverberation mapping technique, a time series analysis of the driving light curve and response, can recover some of this information, including the size and velocity field of the broad-line region (BLR). Here we introduce a new forward-modeling tool, the Broad Emission Line MApping Code, which simulates the velocity-resolved reverberation response of the BLR to any given input light curve by setting up a 3D ensemble of gas clouds for various specified geometries, velocity fields, and cloud properties. In this work, we present numerical approximations to the transfer function by simulating the velocity-resolved responses to a single continuum pulse for sets of models representing a spherical BLR with a radiatively driven outflow and a disklike BLR with Keplerian rotation. We explore how the structure, velocity field, and other BLR properties affect the transfer function. We calculate the response-weighted time delay (reverberation “lag”), which is considered to be a proxy for the luminosity-weighted radius of the BLR. We investigate the effects of anisotropic cloud emission and matter-bounded (completely ionized) clouds and find the response-weighted delay is only equivalent to the luminosity-weighted radius when clouds emit isotropically and are radiation-bounded (partially ionized). Otherwise, the luminosity-weighted radius can be overestimated by up to a factor of 2.

     
    more » « less
  3. ABSTRACT

    To unravel the dominant cause of the weak emission line in a subset of optically selected radio-quiet ‘weak emission line quasars’ (WLQs), we have investigated the possibility of an underdeveloped broad line region (BLR). For this, we have modelled spectral energy distributions (SED) of 61 WLQs by using their optical and infrared (IR) photometric observations from Sloan Digital Sky Survey (SDSS) and WISE, respectively. SED fit consists of various emission components, including the luminosity from the dusty torus (Ltor). For comparison with the normal quasar, we have used a control sample of 55 Quasi-Stellar Objects (QSOs) for each WLQ matching in emission redshift and SDSS r-band. Based on our measurement of Ltor, we found a decrement of 42 ± 2 per cent in IR luminosity in WLQs with respect to the control sample of normal QSOs. Using Ltor/Lbol as the measure of torus covering factor (CFtor), we found a similar decrement in WLQ covering factor, with their CFtor distribution being significantly different with respect to the normal QSOs with a KS test Pnull of 4.27 × 10−14. As dusty torus and BLR covering factors are expected to be of a similar order in active galactic nucleus (AGN), our results suggest that the BLR in the WLQs is underdeveloped and could be a dominant cause of the weakness of their emission line. As a result, our analysis gives support to the models of WLQs based on the evolution scenario being in an early stage of AGNs.

     
    more » « less
  4. ABSTRACT

    Near IR spectroscopic reverberation of Active Galactic Nuclei (AGN) potentially allows the infrared (IR) broad line region (BLR) to be reverberated alongside the disc and dust continua, while the spectra can also reveal details of dust astro-chemistry. Here, we describe results of a short pilot study (17 near-IR spectra over a 183 d period) for Mrk 509. The spectra give a luminosity-weighted dust radius of 〈Rd,lum〉 = 186 ± 4 light-days for blackbody (large grain dust), consistent with previous (photometric) reverberation campaigns, whereas carbon and silicate dust give much larger radii. We develop a method of calibrating spectral data in objects where the narrow lines are extended beyond the slit width. We demonstrate this by showing our resultant photometric band light curves are consistent with previous results, with a hot dust lag at >40 d in the K band, clearly different from the accretion disc response at <20 d in the z band. We place this limit of 40 d by demonstrating clearly that the modest variability that we do detect in the H and K band does not reverberate on time-scales of less than 40 d. We also extract the Pa β line light curve, and find a lag which is consistent with the optical BLR H β line of ∼70–90 d. This is important as direct imaging of the near-IR BLR is now possible in a few objects, so we need to understand its relation to the better studied optical BLR.

     
    more » « less
  5. ABSTRACT

    We examine the long-term optical/near-infrared (NIR) flux variability of a ‘changing-look’ active galactic nucleus (AGN) Mrk 590 between 1998 and 2007. Multiband multi-epoch optical/NIR photometry data from the SDSS Stripe 82 data base and the Multicolor Active Galactic Nuclei Monitoring (MAGNUM) project reveal that Mrk 590 experienced a sudden luminosity decrease during the period from 2000 to 2001. Detection of dust reverberation lag signals between V- and K-band light curves obtained by the MAGNUM project during the faint state in 2003–2007 suggests that the dust torus innermost radius Rdust of Mrk 590 had become very small [Rdust ≃ 32 light-days (lt-days)] by the year 2004 according to the aforementioned significant decrease in AGN luminosity. The Rdust in the faint state is comparable to the H β broad-line region (BLR) radius of RH β, BLR ≃ 26 lt-days measured by previous reverberation mapping observations during the bright state of Mrk 590 in 1990–1996. These observations indicate that the innermost radius of the dust torus in Mrk 590 decreased rapidly after the AGN ultraviolet-optical luminosity drop, and that the replenishment time-scale of the innermost dust distribution is less than 4 yr, which is much shorter than the free fall time-scale of BLR gas or dust clouds. We suggest that rapid replenishment of the innermost dust distribution can be accomplished either by new dust formation in radiatively cooled BLR gas clouds or by new dust formation in the disc atmosphere and subsequent vertical wind from the dusty disc as a result of radiation pressure.

     
    more » « less