skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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. We present new Very Large Telescope Interferometer (VLTI)/GRAVITY near-infrared interferometric measurements of the angular size of the innermost hot dust continuum for 14 type 1 active galactic nuclei (AGNs). The angular sizes are resolved on scales of ∼0.7 mas and the inferred ring radii range from 0.028 to 1.33 pc, comparable to those reported previously and a factor of 10−20 smaller than the mid-infrared sizes in the literature. Combining our new data with previously published values, we compiled a sample of 25 AGNs with bolometric luminosity ranging from 1042to 1047erg s−1, with which we studied the radius-luminosity (R − L) relation for the hot dust structure. Our interferometric measurements of radius are offset by a factor of 2 from the equivalent relation derived through reverberation mapping. Using a simple model to explore the dust structure’s geometry, we conclude that this offset can be explained if the 2 μm emitting surface has a concave shape. Our data show that the slope of the relation is in line with the canonicalR ∝ L0.5when using an appropriately non-linear correction for bolometric luminosity. In contrast, using optical luminosity or applying a constant bolometric correction to it results in a significant deviation in the slope, suggesting a potential luminosity dependence on the spectral energy distribution. Over four orders of magnitude in luminosity, the intrinsic scatter around theR − Lrelation is 0.2 dex, suggesting a tight correlation between the innermost hot dust structure size and the AGN luminosity. 
    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 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