Source author record

Robert Nikutta

Robert Nikutta appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

From Data to Software to Science with the Rubin Observatory LSST

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science.

preprint2022arXiv

Simulating the Legacy Survey of Space and Time stellar content with TRILEGAL

We describe a large simulation of the stars to be observed by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The simulation is based on the TRILEGAL code, which resorts to large databases of stellar evolutionary tracks, synthetic spectra, and pulsation models, added to simple prescriptions for the stellar density and star formation histories of the main structures of the Galaxy, to generate mock stellar samples through a population synthesis approach. The main bodies of the Magellanic Clouds are also included. A complete simulation is provided for single stars, down to the $r=27.5$ mag depth of the co-added wide-fast-deep survey images. A second simulation is provided for a fraction of the binaries, including the interacting ones, as derived with the BinaPSE module of TRILEGAL. We illustrate the main properties and numbers derived from these simulations, including: comparisons with real star counts; the expected numbers of Cepheids, long-period variables and eclipsing binaries; the crowding limits as a function of seeing and filter; the star-to-galaxy ratios, etc. Complete catalogs are accessible through the NOIRLab Astro Data Lab, while the stellar density maps are incorporated in the LSST metrics analysis framework (MAF).

preprint2022arXiv

The Black Hole-Galaxy Connection: Interplay between Feedback, Obscuration, and Host Galaxy Substructure

There is growing evidence for physical influence between supermassive black holes and their host galaxies. We present a case study of nearby galaxy NGC 7582, for which we find evidence that galactic substructure plays an important role in affecting the collimation of ionized outflows as well as contributing to the heavy active galactic nucleus (AGN) obscuration. This result contrasts with a simple, small-scale AGN torus model, according to which AGN wind collimation may take place inside the torus itself, at subparsec scale. Using 3D spectroscopy with the MUSE instrument, we probe the kinematics of the stellar and ionized gas components as well as the ionization state of the gas from a combination of emission line ratios. We report for the first time a kinematically distinct core (KDC) in NGC 7582, on a scale of ~600pc. This KDC coincides spatially with dust lanes and starbursting complexes previously observed. We interpret it as a circumnuclear ring of stars and dusty, gas-rich material. We obtain a clear view of the outflowing cones over kpc scales, and demonstrate that they are predominantly photoionized by the central engine. We detect the back cone (behind the galaxy), and confirm previous results of a large nuclear obscuration of both the stellar continuum and HII regions. While we tentatively associate the presence of the KDC to a large-scale bar and/or a minor galaxy merger, we stress the importance of gaining a better understanding of the role of galaxy substructure in controlling the fueling, feedback and obscuration of AGN.

preprint2021arXiv

Conditional Neural Process for non-parametric modeling of AGN light curve

The consequences of complex disturbed environments in the vicinity of a supermassive black hole are not well represented by standard statistical models of optical variability in active galactic nuclei (AGN). Thus, developing new methodologies for investigating and modeling AGN light curves is crucial. Conditional Neural Processes (CNPs) are nonlinear function models that forecast stochastic time-series based on a finite amount of known data without the use of any additional parameters or prior knowledge (kernels). We provide a CNP algorithm that is specifically designed for simulating AGN light curves. It was trained using data from the All-Sky Automated Survey for Supernovae, which included 153 AGN. We present CNP modeling performance for a subsample of five AGNs with distinctive difficult-to-model properties. The performance of CNP in predicting temporal flux fluctuation was assessed using a minimizing loss function, and the results demonstrated the algorithm's usefulness. Our preliminary parallelization experiments show that CNP can efficiently handle large amounts of data. These results imply that CNP can be more effective than standard tools in modeling large volumes of AGN data (as anticipated from time-domain surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time).

preprint2020arXiv

ALMA polarimetry measures magnetically aligned dust grains in the torus of NGC 1068

The obscuring structure surrounding active galactic nuclei (AGN) can be explained as a dust and gas flow cycle that fundamentally connects the AGN with their host galaxies. This structure is believed to be associated with dusty winds driven by radiation pressure. However, the role of magnetic fields, which are invoked in almost all models for accretion onto a supermassive black hole and outflows, is not thoroughly studied. Here we report the first detection of polarized thermal emission by means of magnetically aligned dust grains in the dusty torus of NGC 1068 using ALMA Cycle 4 polarimetric dust continuum observations ($0.07"$, $4.2$ pc; 348.5 GHz, $860$ $μ$m). The polarized torus has an asymmetric variation across the equatorial axis with a peak polarization of $3.7\pm0.5$\% and position angle of $109\pm2^{\circ}$ (B-vector) at $\sim8$ pc east from the core. We compute synthetic polarimetric observations of magnetically aligned dust grains assuming a toroidal magnetic field and homogeneous grain alignment. We conclude that the measured 860 $μ$m continuum polarization arises from magnetically aligned dust grains in an optically thin region of the torus. The asymmetric polarization across the equatorial axis of the torus arises from 1) an inhomogeneous optical depth, and 2) a variation of the velocity dispersion, i.e. variation of the magnetic field turbulence at sub-pc scales, from the eastern to the western region of the torus. These observations and modeling constrain the torus properties beyond spectral energy distribution results. This study strongly supports that magnetic fields up to a few pc contribute to the accretion flow onto the active nuclei.

preprint2020arXiv

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 (AGN) 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 ($\leq 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 $\sim 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 a typically factor of $\sim 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.

preprint2020arXiv

The variable and non-variable X-ray absorbers in Compton-thin type-II Active Galactic Nuclei

We have conducted an extensive X-ray spectral variability study of a sample of 20 Compton-thin type II galaxies using broad band spectra from XMM-Newton, Chandra, and Suzaku. The aim is to study the variability of the neutral intrinsic X-ray obscuration along the line of sight and investigate the properties and location of the dominant component of the X-ray-obscuring gas. The observations are sensitive to absorption columns of $N_{\rm H} \sim 10^{20.5-24} {\rm cm^{-2}}$ of fully- and partially-covering neutral and/or lowly-ionized gas on timescales spanning days to well over a decade. We detected variability in the column density of the full-covering absorber in 7/20 sources, on timescales of months-years, indicating a component of compact-scale X-ray-obscuring gas lying along the line of sight of each of these objects. Our results imply that torus models incorporating clouds or overdense regions should account for line of sight column densities as low as $\sim$ a few $\times 10^{21}$ cm$^{-2}$. However, 13/20 sources yielded no detection of significant variability in the full-covering obscurer, with upper limits to $ΔN_{\rm H}$ spanning $10^{21-23}$ cm$^{-2}$. The dominant absorbing media in these systems could be distant, such as kpc-scale dusty structures associated with the host galaxy, or a homogeneous medium along the line of sight. Thus, we find that overall, strong variability in full-covering obscurers is not highly prevalent in Compton-thin type IIs, at least for our sample, in contrast to previous results in the literature. Finally, 11/20 sources required a partial-covering, obscuring component in all or some of their observations, consistent with clumpy near-Compton-thick compact-scale gas.

preprint2016arXiv

A Tale of Three Galaxies: A "Clumpy" View of the Spectroscopically Anomalous Galaxies IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948

We investigate the dust properties in three spectroscopically anomalous galaxies (IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948). Their Spitzer/IRS spectra are characterized by a steep ~5-8 micron emission continuum, strong emission bands from polycyclic aromatic hydrocarbon (PAH) molecules, and prominent 10 micron silicate emission. The steep ~5-8 micron continuum and strong PAH emission features suggest the presence of starbursts, while the silicate emission is indicative of significant heating from AGNs. The simultaneous detection of these two observational properties has rarely been reported on galactic scale. We employ the PAHFIT software to estimate their starlight contributions, and the CLUMPY model for the components contributed by the AGN tori. We find that the CLUMPY model is generally successful in explaining the overall dust infrared emission, although it appears to emit too flat at the ~5-8 micron continuum to be consistent with that observed in IRAS F10398+1455 and IRAS F21013-0739. The flat ~5-8 micron continuum calculated from the CLUMPY model could arise from the adopted specific silicate opacity of Ossenkopf et al. (1992) which exceeds that of the Draine & Lee (1984) "astronomical silicate" by a factor up to 2 in the ~5-8 micron wavelength range. Future models with a variety of dust species incorporated in the CLUMPY radiation transfer regime are needed for a thorough understanding of the dust properties of these spectroscopically anomalous galaxies.

preprint2016arXiv

High-Velocity Bipolar Molecular Emission from an AGN Torus

We have detected in ALMA observations CO J = 6 - 5 emission from the nucleus of the Seyfert galaxy NGC 1068. The low-velocity (up to +/- 70 km/s relative to systemic) CO emission resolves into a 12x7 pc structure, roughly aligned with the nuclear radio source. Higher-velocity emission (up to +/- 400 km/s) is consistent with a bipolar outflow in a direction nearly perpendicular (roughly 80 degrees) to the nuclear disk. The position-velocity diagram shows that in addition to the outflow, the velocity field may also contain rotation about the disk axis. These observations provide compelling evidence in support of the disk-wind scenario for the AGN obscuring torus.

preprint2016arXiv

Infrared spectral energy distribution decomposition of WISE-selected, hyperluminous hot dust-obscured galaxies

We utilize a Bayesian approach to fit the observed mid-IR-to-submm/mm spectral energy distributions (SEDs) of 22 WISE-selected and submm-detected, hyperluminous hot dust-obscured galaxies. By adopting the Torus+GB model, we decompose the observed IR SEDs of Hot DOGs into torus and cold dust components. The main results are: 1) Hot DOGs in our submm-detected sample are hyperluminous, with torus emission dominating the IR energy output. However, cold dust emission is non-negligible, averagely contributing ~24% of total IR luminosity. 2) Compared to QSO and starburst SED templates, the median SED of Hot DOGs shows the highest luminosity ratio between mid-IR and submm at rest-frame, while it is very similar to that of QSOs at 10-50um suggesting that the heating sources of Hot DOGs should be buried AGNs. 3) Hot DOGs have both high dust temperatures ~73K and IR luminosity of cold dust. The T-L relation of Hot DOGs suggests that the increase in IR luminosity for Hot DOGs is mostly due to the increase of the dust temperature, rather than dust mass. Hot DOGs have lower dust masses than those of submillimeter galaxies (SMGs) and QSOs within the similar redshift range. Both high IR luminosity of cold dust and relatively low dust mass in Hot DOGs can be expected by their relatively high dust temperatures. 4) Hot DOGs have high dust covering factors, which deviate the previously proposed trend of the dust covering factor decreasing with increasing bolometric luminosity. Finally, we can reproduce the observed properties in Hot DOGs by employing a physical model of galaxy evolution. The result suggests that Hot DOGs may lie at or close to peaks of both star formation and black hole growth histories, and represent a transit phase during the evolution of massive galaxies, transforming from the dusty starburst dominated phase to the optically bright QSO phase. (abridged)

preprint2015arXiv

A Tale of Three Galaxies: Deciphering the Infrared Emission of the Spectroscopically Anomalous Galaxies IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948

The \textit{Spitzer}/Infrared Spectrograph spectra of three spectroscopically anomalous galaxies (IRAS~F10398+1455, IRAS~F21013-0739 and SDSS~J0808+3948) are modeled in terms of a mixture of warm and cold silicate dust, and warm and cold carbon dust. Their unique infrared (IR) emission spectra are characterized by a steep $\simali$5--8$\mum$ emission continuum, strong emission bands from polycyclic aromatic hydrocarbon (PAH) molecules, and prominent silicate emission. The steep $\simali$5--8$\mum$ emission continuum and strong PAH emission features suggest the dominance of starbursts, while the silicate emission is indicative of significant heating from active galactic nuclei (AGNs). With warm and cold silicate dust of various compositions ("astronomical silicate," amorphous olivine, or amorphous pyroxene) combined with warm and cold carbon dust (amorphous carbon, or graphite), we are able to closely reproduce the observed IR emission of these %spectroscopically anomalous galaxies. We find that the dust temperature is the primary cause in regulating the steep $\sim$5--8$\mum$ continuum and silicate emission, insensitive to the exact silicate or carbon dust mineralogy and grain size $a$ as long as $a\simlt1\mum$. More specifically, the temperature of the $\simali$5--8$\mum$ continuum emitter (which is essentially carbon dust) of these galaxies is $\sim$250--400$\K$, much lower than that of typical quasars which is $\sim$640$\K$. Moreover, it appears that larger dust grains are preferred in quasars. The lower dust temperature and smaller grain sizes inferred for these three galaxies compared with that of quasars could be due to the fact that they may harbor a young/weak AGN which is not maturely developed yet.

preprint2015arXiv

Results from the Wide-field Infrared Survey Explorer (WISE) Future Uses Session at the WISE at 5 Meeting

During the "WISE at 5: Legacy and Prospects" conference in Pasadena, CA -- which ran from February 10 - 12, 2015 -- attendees were invited to engage in an interactive session exploring the future uses of the Wide-field Infrared Survey Explorer (WISE) data. The 65 participants -- many of whom are extensive users of the data -- brainstormed the top questions still to be answered by the mission, as well as the complementary current and future datasets and additional processing of WISE/NEOWISE data that would aid in addressing these most important scientific questions. The results were mainly bifurcated between topics related to extragalactic studies (e.g. AGN, QSOs) and substellar mass objects. In summary, participants found that complementing WISE/NEOWISE data with cross-correlated multiwavelength surveys (e.g. SDSS, Pan-STARRS, LSST, Gaia, Euclid, etc.) would be highly beneficial for all future mission goals. Moreover, developing or implementing machine-learning tools to comb through and understand cross-correlated data was often mentioned for future uses. Finally, attendees agreed that additional processing of the data such as co-adding WISE and NEOWISE and extracting a multi-epoch photometric database and parallax and proper motion catalog would greatly improve the scientific results of the most important projects identified. In that respect, a project such as MaxWISE which would execute the most important additional processing and extraction as well as make the data and catalogs easily accessible via a public portal was deemed extremely important.

preprint2014arXiv

First X-ray-Based Statistical Tests for Clumpy-Torus Models: Eclipse Events from 230 Years of Monitoring of Seyfert AGN

We present an analysis of multi-timescale variability in line-of-sight X-ray absorbing gas as a function of optical classification in a large sample of Seyfert AGN to derive the first X-ray statistical constraints for clumpy-torus models. We systematically search for discrete absorption events in the vast archive of RXTE monitoring of dozens of nearby type I and Compton-thin type II AGN. We are sensitive to discrete absorption events due to clouds of full-covering, neutral or mildly ionized gas with columns >~ 10^(22-25) cm^-2 transiting the line of sight. We detect 12 eclipse events in 8 objects, roughly tripling the number previously published from this archive. Peak column densities span ~ 4-26 x 10^22 cm^-2. Event durations span hours to months. The column density profile for an eclipsing cloud in NGC 3783 is doubly spiked, possibly indicating a cloud that is being tidally sheared. We infer the clouds' distances from the black hole to span ~0.3 -140 x 10^4 R_g. In seven objects, the clouds' distances are commensurate with the outer portions of Broad Line Regions (BLR), or the inner regions of infrared-emitting dusty tori. We discuss implications for cloud distributions in the context of clumpy-torus models. The probability of observing a source undergoing an absorption event, independent of constant absorption due to non-clumpy material, is 0.006(+0.160,-0.003) for type Is and 0.110(+0.461,-0.071) for type IIs.

preprint2014arXiv

Near- to Mid- Infrared Imaging and Spectroscopy of Two Buried AGNs of the Nearby Merging Galaxy NGC 6240 with Subaru/IRCS+AO and GTC/CanariCam

We report near-infrared K', L', and M' band imaging observations of the nearby merging galaxy NGC 6240 with the Infrared Camera and Spectrograph on the Subaru telescope. The observations were performed with the assistance of the Subaru Adaptive Optics System, and the achieved spatial resolutions were around 0.10--0.20$^{\prime\prime}$. We also obtained new mid-infrared imaging in the Si-2 filter band (8.7$μ$m) and N-band (7.5--13$μ$m) spectroscopy of this galaxy with the CanariCam on the Gran Telescopio Canarias with a spatial resolution of 0.4--0.5$^{\prime\prime}$. In the K' band image the two nuclei of the galaxy each show a double-peak suggesting the complex geometry of the source, while the L', M', and Si-2 band images show single compact structures in each of the two nuclei. Assuming that the center core observed at wavelengths longer than 3.8$μ$m is associated with dust heated by the buried AGN, we re-evaluated the spectral energy distributions (SEDs) of the southern nucleus from 2 to 30$μ$m with the additional literature values, and performed the SED+spectroscopy fitting using the clumpy torus models of Nenkova et al. (2008) and a Bayesian fitting approach. The model fit suggests that the high covering factor torus emission in the southern nucleus is also obscured by foreground dust in the host galaxy. The estimated AGN bolometric luminosity of the southern nucleus, $L_{\rm bol}({\rm AGN})\sim1\times10^{45}$ [erg$\cdot$s$^{-1}$], accounts for approximately 40% of the whole luminosity of the system.

preprint2014arXiv

Subaru Spectroscopy and Spectral Modeling of Cygnus A

We present high angular resolution ($\sim$0.5$^\prime$$^\prime$) MIR spectra of the powerful radio galaxy, Cygnus A, obtained with the Subaru telescope. The overall shape of the spectra agree with previous high angular resolution MIR observations, as well as previous Spitzer spectra. Our spectra, both on and off nucleus, show a deep silicate absorption feature. The absorption feature can be modeled with a blackbody obscured by cold dust or a clumpy torus. The deep silicate feature is best fit by a simple model of a screened blackbody, suggesting foreground absorption plays a significant, if not dominant role, in shaping the spectrum of Cygnus A. This foreground absorption prevents a clear view of the central engine and surrounding torus, making it difficult to quantify the extent the torus attributes to the obscuration of the central engine, but does not eliminate the need for a torus in Cygnus A.

preprint2014arXiv

The meaning of WISE colours - I. The Galaxy and its satellites

Through matches with the Sloan Digital Sky Survey (SDSS) catalogue we identify the location of various families of astronomical objects in WISE colour space. We identify reliable indicators that separate Galactic/local from extragalactic sources and concentrate here on the objects in our Galaxy and its closest satellites. We develop colour and magnitude criteria that are based only on WISE data to select asymptotic giant branch (AGB) stars with circumstellar dust shells, and separate them into O-rich and C-rich classes. With these criteria we produce an all-sky map for the count ratio of the two populations. The map reveals differences between the Galactic disc, the Magellanic Clouds and the Sgr Dwarf Spheroidal galaxy, as well as a radial gradient in the Large Magellanic Cloud (LMC) disc. We find that the C:O number ratio for dusty AGB stars increases with distance from the LMC centre about twice as fast as measured for near-IR selected samples of early AGB stars. Detailed radiative transfer models show that WISE colours are well explained by the emission of centrally heated dusty shells where the dust has standard properties of interstellar medium (ISM) grains. The segregation of different classes of objects in WISE colour space arises from differences in properties of the dust shells: those around young stellar objects have uniform density distributions while in evolved stars they have steep radial profiles.

preprint2013arXiv

AGN torus properties with WISE

The Wide-field Infrared Survey Explorer (WISE) has scanned the entire sky with unprecedented sensitivity in four infrared bands, at 3.4, 4.6, 12, and 22 micron. The WISE Point Source Catalog contains more than 560 million objects, among them hundreds of thousands of galaxies with Active Nuclei (AGN). While type 1 AGN, owing to their bright and unobscured nature, are easy to detect and constitute a rather complete and unbiased sample, their type 2 counterparts, postulated by AGN unification, are not as straightforward to identify. Matching the WISE catalog with known QSOs in the Sloan Digital Sky Survey we confirm previous identification of the type 1 locus in the WISE color space. Using a very large database of the popular CLUMPY torus models, we find the colors of the putative type 2 counterparts, and also, for the first time, predict their number vs. flux relation that can be expected to be observed in any given WISE color range. This will allow us to put statistically very significant constraints on the torus parameters. Our results are a successful test of the AGN unification scheme.

preprint2011arXiv

Dusty Tori of Luminous Type 1 Quasars at z \sim 2

We present Spitzer infrared spectra and ultra-violet to mid-infrared spectral energy distributions (SEDs) of 25 luminous type 1 quasars at z \sim 2. In general, the spectra show a bump peaking around 3 μm, and the 10 μm silicate emission feature. The 3 μm emission is identified with hot dust emission at its sublimation temperature. We explore two approaches to modeling the SED: (i) using the Clumpy model SED from Nenkova et al. (2008a), and (ii) the Clumpy model SED, and an additional blackbody component to represent the 3 μm emission. In the first case, a parameter search of \sim 1.25 million Clumpy models shows: (i) if we ignore the UV-to-near-IR SED, models fit the 2-8 μm region well, but not the 10 μm feature; (ii) if we include the UV-to-near-IR SED in the fit, models do not fit the 2-8 μm region. The observed 10 μm features are broader and shallower than those in the best-fit models in the first approach. In the second case, the shape of the 10 μm feature is better reproduced by the Clumpy models. The additional blackbody contribution in the 2-8 μm range allows Clumpy models dominated by cooler temperatures (T < 800K) to better fit the 8-12μm SED. A centrally concentrated distribution of a small number of torus clouds is required in the first case, while in the second case the clouds are more spread out radially. The temperature of the blackbody component is ~ 1200 K as expected for graphite grains.

preprint2011arXiv

Emission from Hot Dust in the Infrared Spectra of Gamma-ray Bright Blazars

A possible source of $γ$-ray photons observed from the jets of blazars is inverse Compton scattering by relativistic electrons of infrared seed photons from a hot, dusty torus in the nucleus. We use observations from the Spitzer Space Telescope to search for signatures of such dust in the infrared spectra of four $γ$-ray bright blazars, the quasars 4C 21.35, CTA102, and PKS 1510$-$089, and the BL Lacertae object ON231. The spectral energy distribution (SED) of 4C 21.35 contains a prominent infrared excess indicative of dust emission. After subtracting a non-thermal component with a power-law spectrum, we fit a dust model to the residual SED. The model consists of a blackbody with temperature $\sim1200$ K, plus a much weaker optically thin component at $\sim660$ K. The total luminosity of the thermal dust emission is $7.9\pm0.2 \times 10^{45}$ erg s$^{-1}$. If the dust lies in an equatorial torus, the density of IR photons from the torus is sufficient to explain the $γ$-ray flux from 4C 21.35 as long as the scattering occurs within a few parsecs of the central engine. We also report a tentative detection of dust in the quasar CTA102, in which the luminosity of the infrared excess is $7 \pm 2 \times 10^{45}$ erg s$^{-1}$. However, in CTA102 the far-IR spectra are too noisy to detect the $10 μ$m silicate feature. Upper limits to the luminosity from thermal emission from dust in PKS 1510-089, and ON231, are, $2.3\times10^{45}$, and $6.6\times10^{43}$ erg s$^{-1}$, respectively. These upper limits do not rule out the possibility of inverse Compton up-scattering of IR photons to $γ$-ray energies in these two sources. The estimated covering factor of the hot dust in 4C 21.35, 22%, is similar to that of non-blazar quasars; however, 4C 21.35 is deficient in cooler dust.

preprint2009arXiv

On the 10-micron silicate feature in Active Galactic Nuclei

The 10-micron silicate feature observed with Spitzer in active galactic nuclei (AGN) reveals some puzzling behavior. It (1) has been detected in emission in type 2 sources, (2) shows broad, flat-topped emission peaks shifted toward long wavelengths in several type 1 sources, and (3) is not seen in deep absorption in any source observed so far. We solve all three puzzles with our clumpy dust radiative transfer formalism. (1) We present the spectral energy distribution (SED) of SST1721+6012, the first type 2 quasar observed to show a clear 10-mic silicate feature in emission. We constructed a large database of clumpy torus models and performed extensive fitting of the observed SED, constraining several of the torus parameters. We find that the source bolometric luminosity is ~3*10^12 L_sun. Our modeling suggests that <35% of objects with tori sharing characteristics and geometry similar to the best fit would have their central engines obscured. This relatively low obscuration probability can explain the clear appearance of the 10-mic emission feature in SST1721+6012 together with its rarity among other QSO2. (2) We also fitted the SED of PG1211+143, one of the first type 1 QSOs with a 10-mic silicate feature in emission. Among similar sources, this QSO appears to display an unusually broadened feature whose peak is shifted toward longer wavelengths. Although this led to suggestions of non-standard dust chemistry in these sources, our analysis fits such SEDs with standard galactic dust; the apparent peak shifts arise from radiative transfer effects. (3) We find that the distribution of silicate feature strengths among clumpy torus models closely resembles the observed distribution, and the feature never occurs deeply absorbed. (abridged)

preprint2008arXiv

AGN Dusty Tori: II. Observational Implications of Clumpiness

From extensive radiative transfer calculations we find that clumpy torus models with \No \about 5--15 dusty clouds along radial equatorial rays successfully explain AGN infrared observations. The dust has standard Galactic composition, with individual cloud optical depth \tV \about 30--100 at visual. The models naturally explain the observed behavior of the 10\mic silicate feature, in particular the lack of deep absorption features in AGN of any type. The weak 10\mic emission feature tentatively detected in type 2 QSO can be reproduced if in these sources \No drops to \about 2 or \tV exceeds \about 100. The clouds angular distribution must have a soft-edge, e.g., Gaussian profile, the radial distribution should decrease as $1/r$ or $1/r^2$. Compact tori can explain all observations, in agreement with the recent interferometric evidence that the ratio of the torus outer to inner radius is perhaps as small as \about 5--10. Clumpy torus models can produce nearly isotropic IR emission together with highly anisotropic obscuration, as required by observations. In contrast with strict variants of unification schemes where the viewing-angle uniquely determines the classification of an AGN into type 1 or 2, clumpiness implies that it is only a probabilistic effect; a source can display type 1 properties even from directions close to the equatorial plane. The fraction of obscured sources depends not only on the torus angular thickness but also on the cloud number \No. The observed decrease of this fraction at increasing luminosity can be explained with a decrease of either torus angular thickness or cloud number, but only the latter option explains also the possible emergence of a 10\mic emission feature in QSO2.