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

    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
  3. 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
  4. ABSTRACT

    Multiwavelength variability studies of active galactic nuclei can be used to probe their inner regions that are not directly resolvable. Dust reverberation mapping (DRM) estimates the size of the dust emitting region by measuring the delays between the infrared (IR) response to variability in the optical light curves. We measure DRM lags of Zw229-015 between optical ground-based and Kepler light curves and concurrent IR Spitzer 3.6 and 4.5 µm light curves from 2010 to 2015, finding an overall mean rest-frame lag of 18.3 ± 4.5 d. Each combination of optical and IR light curve returns lags that are consistent with each other within 1σ, which implies that the different wavelengths are dominated by the same hot dust emission. The lags measured for Zw229-015 are found to be consistently smaller than predictions using the lag–luminosity relationship. Also, the overall IR response to the optical emission actually depends on the geometry and structure of the dust emitting region as well, so we use Markov chain Monte Carlo modelling to simulate the dust distribution to further estimate these structural and geometrical properties. We find that a large increase in flux between the 2011–2012 observation seasons, which is more dramatic in the IR light curve, is not well simulated by a single dust component. When excluding this increase in flux, the modelling consistently suggests that the dust is distributed in an extended flat disc, and finds a mean inclination angle of 49$^{+3}_{-13}$ deg.

     
    more » « less
  5. null (Ed.)
    Using VLTI/GRAVITY and SINFONI data, we investigate the subparsec gas and dust structure around the nearby type 1 active galactic nucleus (AGN) hosted by NGC 3783. The K -band coverage of GRAVITY uniquely allows simultaneous analysis of the size and kinematics of the broad line region (BLR), the size and structure of the near-infrared(near-IR)-continuum-emitting hot dust, and the size of the coronal line region (CLR). We find the BLR, probed through broad Br γ emission, to be well described by a rotating, thick disc with a radial distribution of clouds peaking in the inner region. In our BLR model, the physical mean radius of 16 light-days is nearly twice the ten-day time-lag that would be measured, which closely matches the ten-day time-lag that has been measured by reverberation mapping. We measure a hot dust full-width at half-maximum (FWHM) size of 0.74 mas (0.14 pc) and further reconstruct an image of the hot dust, which reveals a faint (5% of the total flux) offset cloud that we interpret as an accreting or outflowing cloud heated by the central AGN. Finally, we directly measure the FWHM size of the nuclear CLR as traced by the [Ca  VIII ] and narrow Br γ line. We find a FWHM size of 2.2 mas (0.4 pc), fully in line with the expectation of the CLR located between the BLR and narrow line region. Combining all of these measurements together with larger scale near-IR integral field unit and mid-IR interferometry data, we are able to comprehensively map the structure and dynamics of gas and dust from 0.01 to 100 pc. 
    more » « less