Mid-IR Variability from the Spitzer Deep Wide-Field Survey

Mid-Infrared Variability from the Spitzer Deep Wide-Field Survey

Szymon Kozłowski11affiliation: Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA; e-mail: simkoz@astronomy.ohio-state.edu , Christopher S. Kochanek11affiliation: Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA; e-mail: simkoz@astronomy.ohio-state.edu 22affiliation: The Center for Cosmology and Astroparticle Physics, The Ohio State University, Columbus, OH 43210, USA , Daniel Stern33affiliation: Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Drive, Pasadena, CA 91109, USA , Matthew L. N. Ashby44affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA , Roberto J. Assef11affiliation: Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA; e-mail: simkoz@astronomy.ohio-state.edu , J. J. Bock55affiliation: California Institute of Technology, Pasadena, CA 91125, USA , C. Borys55affiliation: California Institute of Technology, Pasadena, CA 91125, USA , K. Brand66affiliation: Space Telescope Science Institute, Baltimore, MD 21218, USA , M. Brodwin44affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 2121affiliation: W. M. Keck Postdoctoral Fellow at the Harvard-Smithsonian Center for Astrophysics , M. J. I. Brown77affiliation: School of Physics, Monash University, Clayton 3800, Victoria, Australia , R. Cool88affiliation: Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA , A. Cooray99affiliation: University of California, Irvine, CA 92697,USA , S. Croft1010affiliation: University of California, Berkeley, CA 94720,USA , Arjun Dey1111affiliation: National Optical Astronomical Observatory, 950 North Cherry Avenue, Tucson, AZ 85719, USA , P. R. Eisenhardt33affiliation: Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Drive, Pasadena, CA 91109, USA , A. Gonzalez1212affiliation: Department of Astronomy, University of Florida, Gainesville, FL 32611, USA , V. Gorjian33affiliation: Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Drive, Pasadena, CA 91109, USA , R. Griffith33affiliation: Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Drive, Pasadena, CA 91109, USA , N. Grogin66affiliation: Space Telescope Science Institute, Baltimore, MD 21218, USA , R. Ivison1313affiliation: UK Astronomy Technology Centre, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK 1414affiliation: Institute of Astronomy, University of Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK , J. Jacob33affiliation: Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Drive, Pasadena, CA 91109, USA , B. T. Jannuzi1111affiliation: National Optical Astronomical Observatory, 950 North Cherry Avenue, Tucson, AZ 85719, USA , A. Mainzer33affiliation: Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Drive, Pasadena, CA 91109, USA , L. Moustakas33affiliation: Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Drive, Pasadena, CA 91109, USA , H. Röttgering1515affiliation: Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands , N. Seymour1616affiliation: Mullard Space Science Laboratory, University College London, UK , H. A. Smith44affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA , S. A. Stanford1717affiliation: University of California, Davis, CA 95616, USA , J. R. Stauffer1818affiliation: Spitzer Science Center, California Institute of Technology, Pasadena, CA 91125, USA , I. S. Sullivan66affiliation: Space Telescope Science Institute, Baltimore, MD 21218, USA , W. van Breugel1919affiliation: University of California, Merced, CA 95344, USA , S. P. Willner44affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA , and E. L. Wright2020affiliation: University of California, Los Angeles, CA 90095-1562, USA

We use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to investigate the variability of objects in 8.1 deg of the NOAO Deep Wide Field Survey Boötes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 µm bands. Out of analyzed sources, 1.1% meet our standard variability selection criteria, that the two light curves are strongly correlated () and that their joint variance () exceeds that for all sources with the same magnitude by . We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBoötes survey, radio catalogs, 24 µm  selected active galactic nucleus (AGN) candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects, although most of the stellar, galaxy, and unclassified sources are false positives. For our standard selection criteria, 11%–12% of the mid-IR counterparts to X-ray sources, 24 µm AGN candidates, and spectroscopically identified AGNs show variability. The exact fractions depend on both the search depth and the selection criteria. For example, 12% of the 1131 known AGNs in the field and 14%–17% of the known AGNs with well-measured fluxes in all four IRAC bands meet our standard selection criteria. The mid-IR AGN variability can be well described by a single power-law structure function with an index of at both 3.6 and 4.5 µm, and an amplitude of mag on rest-frame timescales of 2 years. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities.

cosmology: observations — galaxies: active — infrared: galaxies — quasars: general

1 Introduction

A substantial fraction of objects in the universe change in brightness with time. This apparent variability can be due to physical periodic changes in the objects (e.g., Cepheids, RR Lyrae), motion (e.g., eclipsing binaries, microlensing, rotating stars), explosive events (e.g., supernovae, novae), and accretion (e.g., cataclysmic variables, active galactic nuclei (AGNs)). In a high-latitude extragalactic survey, the most common variable sources are active galactic nuclei (e.g., Sesar et al. 2007). Optical variability is one method for identifying quasars111We will use words “quasar” and “AGN” interchangeably throughout this paper. (e.g., Eyer 2002; Geha et al. 2003; Rengstorf et al. 2004), although it is only recently that a fully quantitative approach to variability selection has been developed (Kozłowski et al. 2010).

The emission from quasars generally has three components. The UV to near-IR radiation is dominated by a hot accretion disk extending from an inner edge of a few gravitational radii from the black hole outward with, in simple thin disk theory, a temperature profile (Shakura & Sunyaev 1973). Near the inner edge, there is a corona of hotter gas that produces the non-thermal X-ray emission (Haardt et al. 1994). On scales where the temperature is below the dust sublimation temperature ( K), dust absorbs radiation from the disk and reradiates the energy in the mid-IR and far-IR (Barvainis 1987). The overall spectrum typically shows a minimum near 1 µm, with the emission from the disk rising toward the UV and the emission reradiated by dust rising toward the far-IR (e.g., Sanders et al. 1989). There is an increasing evidence that many physical relations between AGNs, galaxies and their large scale clustering have to be taken into account if their formation and evolution is to be understood (e.g., Tasse et al. 2008).

We know a great deal about the optical variability of quasars both from large studies of the variability seen in ensembles of sparsely monitored quasars and from detailed studies of individual quasars. Ensemble studies (e.g., Vanden Berk et al. 2004; de Vries et al. 2005) have shown that variability increases with decreasing optical wavelength, decreasing luminosity, and potentially decreasing black-hole mass. The structure function of the ensemble variability is a power law with smaller variability amplitudes on short timescales, with some evidence for saturation on timescales of order a few decades. Until recently, there were few studies of individual quasars (e.g., Cutri et al. 1985; Clavel et al. 1989; Hook et al. 1994), but this has changed dramatically in the last year. The light curves of individual quasars are well modeled by a damped random walk, a stochastic process described by the amplitude of the random walk and a damping timescale for returning to the mean luminosity (Kelly et al. 2009; Kozłowski et al. 2010; MacLeod et al. 2010). While preliminary indications from Kelly et al. (2009), based on quasars, suggest that these two process parameters are related to the quasar luminosity and black hole mass, MacLeod et al. (2010) used Sloan Digital Sky Survey (SDSS) Stripe 82 quasars to find a number of clear trends. For example, the asymptotic variability on long timescales decreases with increasing luminosity and rest-frame wavelength, and is correlated with black hole mass. The timescale for returning to the mean luminosity increases with wavelength and also with increasing black hole mass, but remains constant with redshift and luminosity.

Far less is known about the near-IR and mid-IR variability of quasars. Since AGN luminosities do vary in time at optical wavelengths, we expect to find some variability in the IR, but it could be heavily smoothed by averaging the response over the large scales of either the cooler parts of the disk or the absorbing dust. Glass (2004) found that the majority (39/41) of the surveyed low-luminosity nearby Seyferts do vary in the near-IR and mid-IR, and that the variability is most apparent at longer wavelengths due to diminishing flux from the host galaxy. The most extensive work to date on near-IR and mid-IR variability of quasars is that of Neugebauer & Matthews (1999). These authors studied 25 low-redshift quasars (), which were observed in five bands (1.3 – 10.6 µm) for years. The amplitude of variations for radio loud (quiet) quasars is mag (0.1 mag) on rest-frame time lags (time differences between any two epochs) of 2 years at 10.6 µm. For one quasar, PG 1226+023, they measured structure functions in five bands, all showing mag variations for a rest-frame time lag of 2 years. Another example is the work of Suganuma et al. (2006), who found a delayed response between the - and -band variability in four Seyfert 1 galaxies. These authors showed that -band light curves have a fairly tight time lag–luminosity relation, , while broad line region (BLR) time lags are more scattered. The interpretation is that the -band emission is mainly coming from dusty clouds bounding the BLR, and that the size scale of these dust clouds is very dependent on the AGN luminosity, since they have to be colder than the dust sublimation temperature.

In this paper, for the first time, we investigate the variability of a truly large number of mid-IR sources, using the four-epoch Spitzer Deep Wide Field Survey (SDWFS; Ashby et al. 2009) of the Boötes field of the NOAO Deep Wide Field Survey (NDWFS; Jannuzi & Dey 1999). The two goals are to characterize the mid-IR variability of a large sample of AGNs, and to use the variability to identify additional AGNs. While mid-IR quasar selection (Stern et al. 2005) works very well when the AGN dominates the luminosity (even for quasars behind dense stellar fields such as LMC/SMC; Kozłowski & Kochanek 2009), the method begins to fail as the luminosity becomes comparable to that of the host (Gorjian et al. 2008; Assef et al. 2010b). With mid-IR variability, we should be able to identify such sources and extend the sample of mid-IR selected quasars to lower luminosities. We also make use of the extensive wavelength coverage of the field, including X-ray (XBoötes; Murray et al. 2005), the earlier Infrared Array Camera (IRAC) Shallow Survey that became the first epoch of SDWFS (Eisenhardt et al., 2004), the 24 µm survey of the field (see Houck et al. 2005), and the FIRST (Becker et al., 1995) and WSRT (de Vries et al., 2002) radio surveys. Based on these data, the AGN and Galaxy Evolution Survey (AGES; C. S. Kochanek et al. 2010, in preparation) selected galaxies and AGNs at all these wavelengths for spectroscopic follow up using the Hectospec instrument on the MMT (Fabricant et al. 2005), providing roughly 23,500 redshifts in the field. In Section 2 we describe the data analysis, and in Section 3 we describe our variability selection criteria. In Section 4, we present the results of this study. The mid-IR structure function analysis is presented in Section 5 followed by a summary in Section 7. Throughout this paper, we use a standard CDM model with and . Where needed, we use the galaxy and AGN templates of Assef et al. (2010b) to describe the mid-IR colors of various populations and to make any K-corrections or estimates of absolute luminosities.

2 Data Analysis

The SDWFS is described in detail by Ashby et al. (2009). The images were taken with the IRAC instrument (Fazio et al. 2004) on the Spitzer Space Telescope, covering wavelengths of 3.6, 4.5, 5.8, and 8.0 m, hereafter referred to as the [3.6], [4.5], [5.8], and [8.0] bands, respectively. A single IRAC frame is pixels. We work with full-field SDWFS mosaics constructed from the 20,000 IRAC images/epoch. Each approximately pixel mosaic covers 10.2 deg, centered at (R.A., decl.14:31:27, +34:10:43), with a pixel scale of 0.84 arcsec. Counting the original IRAC Shallow Survey (Eisenhardt et al. 2004) as epoch 1, the SDWFS data span a period of four years between 2004 and 2008, where the intervals between the epochs are approximately 3.5 years, 0.5 years, and 1 month (observations were done on 2004 January 10-14, 2007 August 8-13, 2008 February 2-6, and 2008 March 6-10). For each individual observation of a patch of the sky ( pixels), there is only a s offset between the channel 1/3 and 2/4 observations. For a given SDWFS epoch, each patch of sky was observed times over a 1–2 day period. Each SDWFS epoch of the entire 10 deg area took days, and we treat it as a single image with a single Julian Date. This does not affect our final conclusions, as we show later that there is little mid-IR AGN variability on these timescales.

To study the variability of objects in the SDWFS field, we used the Difference Image Analysis (DIA) method (Alard & Lupton 1998; Woźniak 2000). While the variability itself is the subject of this paper, the DIA method also proved to be a useful diagnostic tool for constructing artifact-free mosaics (see Ashby et al. 2009). The basic idea behind the DIA method is to match one image, called the template or reference image, both astrometrically and photometrically to another image of the same area of the sky. Once they are astrometrically aligned, the template image (which is chosen to be of better quality) is convolved with a kernel function to match the point spread function (PSF) of the analyzed image. The two images are then subtracted from each other (Figure 1) to leave only variable objects, noise, and systematic problems (cosmic rays, artifacts, etc.). After some experimentation, we concluded that the [5.8] and [8.0] maps were too shallow to perform DIA and obtain reliable light curves, so our variability analysis is restricted to the [3.6] and [4.5] bands.

Spitzer PSFs have a “triangular” shape with extended diffraction-like patterns (more precisely they are “multiplexer bleed” or “muxbleed,” “column pull-down,” and “banding”222For details see the IRAC Data Handbook at http://ssc.spitzer.caltech.edu/irac/dh/.), and the PSF rotates between the epochs. DIA works best on simple quasi-Gaussian PSFs, so we pre-filtered the images to have a simpler PSF structure (Figure 1). For each region, we scaled and median combined the brightest stars to build an empirical model of the PSF, . We then linearly filtered the image to create a filtered image defined by the combination of Fourier transforms , that transforms the PSF from to , where is a Gaussian broad enough to avoid noise amplification (i.e., faster than for small scales/large wave numbers). The filtering process conserves fluxes and was carried out on a grid of pixel regions that overlap by 20 pixels along their edges. To save some computation time, each of these regions was then cut into smaller pixel sub-images (again each with 20 pixels of margin) for the DIA analysis. We combined all four epochs to build the template images and these templates were scaled and subtracted from each epoch using DIA (Woźniak 2000). The brightest stars in the field produce artifacts due to subtraction residuals, particularly from the “diffraction” spikes in the PSF. We therefore masked a region around all stars with mag in the 2MASS catalog of point sources (Cutri et al. 2003). The masking radius was arcsec, which is large enough to include the spikes. We used the Two Micron All Sky Survey (2MASS) catalogs because these bright stars were saturated and not included in the SDWFS catalogs. The masked regions constitute % of the area.

The [3.6], [4.5], [5.8], and [8.0] Vega magnitudes, and the coordinates (J2000.0) are taken from the SDWFS catalogs (Ashby et al. 2009). Throughout this paper, we use 4 arcsec diameter aperture magnitudes with aperture corrections to 24 arcsec. The original 3.6 µm full SDWFS catalog contains 677,522 objects (catalog DR1.1), but we cropped it to 515,055 objects (hereafter called the “cropped sample”) as variations in coverage between epochs reduced the usable area for variability analysis by about 20%, to 8.1 deg. Using DIA’s PSF fitting routines, we could construct light curves (at least three epochs) for 492,420 objects simultaneously in both [3.6] and [4.5]. The final number of objects used in this analysis after masking regions around bright stars is 474,179. While we include the results for all objects from the cropped sample in Table 1, in the analysis that follows we consider only the sources with at least three epochs. Table 2 provides the DIA light curves for all SDWFS sources from the cropped sample.

In order to understand the nature of the variable sources, we matched the SDWFS sources to several of the available data sets for the Boötes field, focusing on information related to AGN activity. Thus, we matched to the X-ray and radio catalogs, to the subset of 24 µm sources with high probabilities of being AGNs, and to the AGES spectroscopic data. We used the X-ray catalogs from the XBoötes Chandra survey covering 9.3 deg of the Boötes field in the 0.5-7 keV band (Murray et al. 2005; Kenter et al. 2005; Brand et al. 2006; Gorjian et al. 2008). We matched a version of the Murray et al. (2005) catalogs extended downward in flux from sources with 4 or more counts to sources with 2 or more counts. While using 2 count (about erg cm s) X-ray sources sounds peculiar, the backgrounds in XBoötes are so low that almost all such sources with mid-IR identifications will be real, particularly if the match is to one of the relatively rare mid-IR variable sources. Using the Bayesian matching method of Brand et al. (2006), we identified X-ray sources with at least two epochs of SDWFS data. We used the same procedures to match to the 20 cm FIRST (Becker et al., 1995) and 1.4 GHz WSRT (de Vries et al. 2002) radio catalogs. Restricting the matching to likely radio point sources in the FIRST and WSRT surveys, we found and matches, respectively, to SDWFS sources with at least two epochs of data. We know from AGES that most 24 µm, optical point sources satisfying the criterion that are quasars (Brown et al. 2006; Eckart et al. 2010; Assef et al. 2010a), where the optical criterion is used to eliminate stars. Finally, we matched to the AGES redshift survey, finding (, , ) SDWFS sources (with redshifts ). The modified SDSS pipeline used for the AGES redshifts (HSRED; Cool et al. 2008) provides an estimate of the best fit spectral template used for the fits that can be used to classify objects. More generally, we have all the galaxy and AGN target selection codes used by AGES to select and prioritize spectroscopic targets.

3 Variability

For each light curve (from Table 2) with at least three epochs in both the [3.6] and [4.5] bands, we computed the standard deviation for that light curve as a function of (Vega) apparent magnitude , as shown in Figure 2. We also calculated the variability covariance and Pearson’s correlation coefficients between the two channels, where


and is the average magnitude for the th light curve in channel . The correlation coefficient is limited to the range , where () means that the [3.6] and [4.5] light curves are perfectly correlated (anti-correlated). The distribution of is dominated by contributions from noise (and any systematic errors) rather than variability, so we need to define a threshold for selecting candidate variables from these distributions. Let be the median of the as a function of magnitude and be the standard deviation of the values around this median. Sources are more likely to be true variables as their distance, , from the median increases, and if the [3.6] and [4.5] variabilities are strongly correlated. A reasonable metric for the joint significance of the variances in the light curves is


which quantifies the degree to which the object deviates from the median variances in both bands. Objects with positive excess variances (large and ) and strong correlations between the two bands are the best candidates for true variables.

Figure 3 shows the distribution of correlation strengths for sources with exceeding a range of thresholds. We see that there is a peak for highly correlated sources whose strength increases as the variability significance increases. This demonstrates that the variability in the two bands is tightly correlated for real variables and that we are likely to produce a relatively clean sample of variables (the false positive rate is 6%–7% for AGNs in the Photometry Group 3 defined below) if we restrict ourselves to and . Variability selection is always a tradeoff between completeness and contamination, so we will consider sources in three different variability levels, all restricted to have correlations between the bands. To avoid confusion, the Variability Levels correspond directly to the values of . Variability Levels 2, 3 and 4 correspond to sources with excess variances and significances exceeding , and , respectively, with Variability Level 2 being the most complete and the most contaminated, and Variability Level 4 being the least complete and least contaminated.

We estimate the fraction of false positives by assuming that the density of uncorrelated sources () in Figure 3 represents the density of false positives amidst the highly correlated () sources. If the number of identified variables is , that is sources with (or 3, 4) and , and the number of sources with (or 3, 4) and is , then we estimate that the number of false positives is , where is the ratio of the intervals and . The fraction of the identified variables that are false positives is then .

Our ability to characterize sources also depends on the magnitudes of the object, because this determines how well we can measure the mid-IR colors. We divide sources into three photometry groups based on their mid-IR magnitudes.

  • Photometry Group 1 consists of all objects that after masking have three or four epochs of data, with no constraints on either their magnitudes or magnitude uncertainties. There are 474,179 Photometry Group 1 objects after masking, of which 5107, 2208, and 1071 meet the Variability Level 2, 3, and 4 criteria.

  • Photometry Group 2 is the subset of the Photometry Group 1 sources for which we can measure the [3.6][4.5] color. These sources must satisfy mag and mag with uncertainties in these bands smaller than 0.1 mag. There are 213,594 Photometry Group 2 objects, of which 1557, 668, and 359 meet the Variability Level 2, 3, and 4 criteria.

  • Photometry Group 3 is the subset of Photometry Group 2 objects for which we can also measure the [5.8][8.0] color. These sources are required to have uncertainties smaller than mag for the 5.8 and 8.0 µm bands. There are 39,522 Photometry Group 3 objects, of which 775, 428, and 260 meet the Variability Level 2, 3, and 4 criteria.

Table 1 provides these significance statistics for the SDWFS sources. Note, however, that our catalogs are limited to the objects detected on the template image – any transient which does not produce a detectable source in the template image, which is just the average of four epochs, will not be identified. The most variable, highly correlated () sources in the field are AGNs, peaking at .

We will make extensive use of the [3.6][4.5] versus [5.8][8.0] mid-IR color-color distribution to characterize the sources. Figure 4 shows the distribution of all Photometry Group 3 SDWFS sources in this space. Eisenhardt et al. (2004) suggested that for bright sources the red plume in [3.6][4.5] color was likely to be dominated by AGNs. This was confirmed by Stern et al. (2005) based on the first year of spectroscopic observations by the AGES survey. Stern et al. (2005) defined a selection region, the “AGN wedge”, for mid-IR AGN selection that is also shown in Figure 4. For these Vega colors, stars lie at color 0 unless they have dusty winds or disks to produce a mid-IR excess. Using the spectral energy distribution (SED) models of Assef et al. (2010b), which are derived by fitting the SEDs of Boötes sources, we can illustrate the locations of various extragalactic populations as a function of redshift. Early-type galaxies at low redshifts have colors similar to stars and then become redder in [3.6][4.5] for at roughly fixed [5.8][8.0], tracking along the left edge of the AGN wedge. They become red in [3.6][4.5] due to the flux “peak” near rest-frame wavelength of 1.6 µm, which passes through the 4.5 µm band at a redshift of 1.5–2.1. Late-type galaxies at low redshift have red [5.8][8.0] colors at low redshift due to the polycyclic aromatic hydrocarbon (PAH) emission in the 8 µm band, then move horizontally to bluer [5.8][8.0] along the bottom of the AGN wedge as the feature is redshifted out of the 8 µm band, and finally follow the early-type galaxies up along the left edge of the selection region. Pure AGNs remain inside the AGN wedge except near where they make an excursion out of the AGN wedge when the H line lies in the 3.6 µm band. If we use an equal combination of an AGN and a host galaxy, defined by the total emission from 0.1 to 30 µm, the combined SED is somewhat bluer. With further increases in the host luminosity relative to the AGN, the source drops below the bottom of the AGN wedge.

In fact, the distribution of the IRAC colors for the XBoötes X-ray sources does extend below the Stern et al. (2005) AGN wedge (Gorjian et al., 2008). Based on Gorjian et al. (2008), we define a “modified AGN wedge” that encompasses most X-ray sources. This modified AGN wedge is defined by mag (left side), mag (right side), and mag (bottom). Because we are identifying AGN candidates based on their mid-IR variability, we should be able to identify them in this blueward extension of the Stern et al. (2005) selection region despite any host contamination. The caveat, however, is that this blueward extension is dominated by false positives (galaxies) if used to select AGNs based solely on their mid-IR colors without additional information such as variability or X-ray emission. We also define a “stellar” box which should contain most normal stars. The edges of this box are defined by mag mag and mag mag. The coldest stars, reaching into the brown dwarf regime, become red in due to methane absorption in the 3.6 µm band, but have colors near zero and thus stay to the left of the AGN wedge (Stern et al. 2007; Eisenhardt et al. 2010).

Beyond these mid-IR color diagnostics, we also use the extensive multi-wavelength and spectroscopic data for the field to characterize the sources. We know that the variable sources must be dominated by AGNs and a small fraction of stars, so the basic check on our selection criteria is that our selection limits on and preferentially select sources with AGN signatures, without introducing large numbers of contaminating sources.

4 Results

We will now characterize the variable sources using three approaches. First, we simply examine the distribution of objects in mid-IR color. Second, we examine the variability of AGN candidates identified at other wavelengths. Finally, we consider the spectral properties of the objects with known redshifts.

Figure 5 and Table 3 characterize the distribution of objects in mid-IR color. The Photometry Group 3, Variability Level 2 candidates are concentrated in either the modified AGN wedge, with a distribution very similar to that of X-ray sources (Gorjian et al. 2008; see Figure 4), or in the stellar box. While only 0.5% of the Photometry Group 3 objects outside these two regions are Variability Level 2 sources, 7.4%, 4.6%, and 1.1% of the sources in the old AGN, modified AGN, and stellar regions are Variability Level 2 sources. The estimated false positive rates are 44%, 6%, 7%, and 58% for these four regions (outside, old AGN wedge, modified AGN wedge, stellar box) in the color-color space, respectively. This simply means that, for example, of the 7.4% variable sources in the old AGN wedge 6% are false positives, so in fact 7.4%% of all objects are real variables and 7.4%% of all objects are falsely identified as variable. Increasing the variability significance level reduces the variable fractions roughly in the same proportion, suggesting that even the Variability Level 2 candidates are dominated by true variables. The distribution also validates our assumption that the variability is dominated by AGN – 75% of the candidates lie in the modified AGN region for Variability Level 4. The selection remains pure even for Variability Level 2 because the requirement that the [3.6] and [4.5] light curves are tightly correlated () eliminates most of the contamination even when we only require (see Figure 3). Most variable sources in the modified AGN wedge are known to be AGNs based on other evidence. Only 16%, 10%, and 7% of the Photometry Group 3, Variability Level 2, 3, and 4 sources lack photometric or spectroscopic signs of AGN activity. Most of these, however, are likely to be false positives given our estimated false positive rates of 7%, 4%, and 3%. For the Photometry Group 2 objects, where we lack accurate [5.8][8.0] color, we can only investigate the distribution in [3.6][4.5] color, and the side panel to Figure 5 shows the fraction of variable sources as a function of this color. Not surprisingly, the fraction is high (%) when the color is red, corresponding to AGNs, and then drops dramatically for colors consistent with normal galaxies. The fraction of variables in the other two photometry groups are dramatically lower than for Photometry Group 3 for two reasons. First, a source must be intrinsically more variable in order to meet any variability criterion because the intrinsic measurement errors are larger. Second, with less well-measured colors, objects are scattered from their true location in the color-color space, and the differential scattering of galaxies into the AGN regions and AGNs out of these regions dilutes the selection criterion.

The 86 (Photometry Group 3, Variability Level 2) variable objects inside the stellar box are red, with optical colors of (i.e., ), so they must be dominated by red giants or main-sequence stars. The red and asymptotic giant branch stars include the OGLE Small Amplitude Red Giants (OSARGs; e.g., Soszyński et al. 2004) with -band amplitudes of mag and periods between 10 and 1000 days, the Long Period Variables (LPVs; e.g., Soszyński et al. 2009) with amplitudes of up to mag and periods 10-1000 days, and the Long Secondary Period variables (LSPs; e.g., Soszyński 2007) with periods 30–1000 days. However, given their absolute magnitudes of mag, they would have to be at a distance of kpc to match the observed , , and IRAC magnitudes. Since this scenario seems unlikely, the variable objects inside the stellar box could be the main-sequence flaring M dwarfs (e.g., Kowalski et al. 2009) with absolute magnitudes of mag, and distances of up to 1 kpc. Based on Kowalski et al. (2009), we expect to find M dwarf flares in our data. The and colors of the variable objects inside the stellar box are inconsistent with the galaxy templates from Assef et al. (2010b). We have not pursued these sources in any detail because they are almost certainly quotidian variable stars and false positives. The estimated false positive rate of 58% is very high.

Next, we examine the variability of sources whose X-ray, MIPS 24 µm or radio fluxes indicate that the source is an AGN. The color distributions of these sources are shown in Figure 6 and the variability statistics are summarized in Table 4. Here we focus on Photometry Group 1 since we are not dependent on having mid-IR colors to interpret the results. First, we see that a significant fraction of the X-ray sources show mid-IR variability, with 12% ( accounting for false positives) being variable even with no limits on the mid-IR photometric uncertainties. The higher fractions in Photometry Groups 2 and 3 are largely due to our improved ability to detect a given level of variability as we require smaller photometric uncertainties. The MIPS quasar candidates, which are found to be either AGNs () or optically unresolved star forming galaxies upon spectroscopic examination, also show high variability rates, with 11% () of the Photometry Group 1 objects showing variability at Variability Level 2. The variable fraction again increases as we go to the brighter photometric groups. The two sets of radio sources show significantly lower levels of variability. This is consistent with the results of Hickox et al. (2009) and Griffith & Stern (2010) who show that, as compared to X-ray and mid-IR selected AGNs, radio-selected AGNs have lower Eddington ratios and are much less likely to be AGN dominated. Thus, variability in radio-selected AGNs will be more challenging to discern relative to the photometric noise due to the host galaxy.

Finally, we compare the variable sources with the template types used to estimate the source redshifts in the AGES survey, as shown in Figure 7 and summarized in Table 5. This is somewhat redundant with the previous two tests because the objects targeted as AGNs are chosen to be in the IRAC AGN wedge, to show X-ray or radio emission, or to be optical point sources with non-stellar 24 µm fluxes. The success rate for the spectroscopic identification of candidate AGNs in AGES is quite high. The sample of stars with spectra is very heterogeneous. Roughly half of the stars are F and G stars observed by AGES as spectrophotometric calibrators. Most of the remainder are K and M stars which were flagged as candidate AGNs either because of their X-ray fluxes or the presence of a mid-IR excess. We again see that the fraction of variable sources is highest for the AGNs, with roughly 11% meeting the Level 2 variability criterion.

In summary, small relatively shallow mid-IR variability surveys like SDWFS are a workable but not a highly efficient means of identifying AGNs. Roughly speaking, the efficiencies are of order 10% for variability selection criterion providing relatively low false positive rates. It does, however, overcome the weaknesses of mid-IR color selection to identify AGNs outside the Stern et al. (2005) selection region. With difference imaging, the host galaxy flux that moves the source outside this color region ceases to be a problem other than through its effects on the photometric noise.

5 Mid-IR Structure Functions

The structure function is a measure of the expected changes in an object’s flux as a function of the time between the measurements. The time lag is defined as the time difference between any two data points in a light curve, where the observed time lag is and the rest-frame time lag is . For each time lag, one wants to find all pairs of points in a light curve and calculate the variance of the magnitude differences, . Such a variance can also be calculated by averaging over many objects of the same class (i.e., AGNs) for each time lag between any two epochs, as we do in this paper.

Given the depth of the AGES redshift survey, all objects are spectroscopically confirmed AGNs (see Figure 2 in Stern et al. 2005). We measured the mid-IR structure function of Photometry Group 3 AGNs to fixed magnitude limits of , , and mag. Only of these objects were actually detected as significantly variable sources (see Table 5), so most of the observed variance in the magnitudes is at the noise level. This is a rather different regime from most optical studies, where the variance is dominated by variability, and makes it critically important to correctly characterize the contribution to the variance from the measurement errors. We do this by matching each AGN with four unresolved galaxies of nearly identical observed flux. We chose these galaxies to be Photometry Group 2 and 3 objects with colors mag and 0.5 mag mag, thus avoiding AGNs. We then compute the structure function as


In addition to the simple calculation over the , , and AGNs with [3.6], , and  mag, we also calculated the structure functions for 1000 bootstrap resamplings of both the AGN and galaxy lists. These give consistent results, and we generally report the median result of the bootstrap samples and the region encompassing 68.3% of the samples about the median. The small values at short lags are very sensitive to any systematic problems in the noise estimate, and while this empirical approach to correctly characterizing the noise appears to work well, our estimates of the amplitude of the structure function on the longer timescales are probably more robust than our estimate of the slope of the structure function.

Figure 8 shows the resulting observed and rest-frame structure functions for the [3.6] and [4.5] bands as a function of magnitude limit. We generally see more variability on longer timescales. We fit our structure functions with a power-law


using  years ( years) for the observed (rest-frame) estimates. Table 6 presents the resulting fits. The [3.6] and [4.5] structure functions are consistent with each other, but also significantly steeper than the -band structure function from an ensemble of SDSS quasars (Vanden Berk et al. 2004), with , and higher amplitude, with  mag for  years.

The steeper mid-IR slope could be explained by the emission being dominated by larger physical scales than the optical emission. At these redshifts, the observed frame mid-IR corresponds to the rest-frame near-IR, where the emission has contributions from both the disk and hot dust. If we focus on the more compact disk emission, the characteristic scale of an accretion disk for emission at wavelength , defined by the point where the disk temperature matches the photon energy, is


where is the Eddington ratio and defines the radiative efficiency . Taking , , and , the optical emission for comes from a region only a light day ( light years) in size, while the mid-IR scale is approximately 10 light days ( light years). Thus, our shortest lags correspond to the fastest possible response time given the size of the emission region since the typical AGES quasar has (Kollmeier et al. 2006). The orbital, thermal and viscous timescales are even longer, and Kelly et al. (2009, also see Kozłowski et al. 2010) find that the optical variability can be characterized by an amplitude and a timescale that might be associated with thermal timescales. In any case, the mid-IR emission is coming from scales so large that it is hard to respond to any driving process on the timescale of our shortest observed lags. This should suppress variability on these short timescales compared to the optical emission, which can respond more easily, leading to a steeper mid-IR structure function.

Next we subdivided the sample into two redshift bins, and , and then divided the objects in these bins at the median value of their absolute magnitudes. Figure 9 shows the structure functions of these subdivided samples as compared to the best fits to the full samples. At least for the well-determined values at longer lags, there is a pattern that shorter rest-frame wavelengths vary more than longer rest-frame wavelengths. Also for AGNs, the lower luminosity systems vary more than the higher luminosity systems. This is consistent with earlier findings in the optical part of the spectrum (e.g., Hook et al. 1994; Hawkins 2002; Vanden Berk et al. 2004; de Vries et al. 2005). In general, the trend of more variability with longer lags also holds, although there are some exceptions in the [4.5] results.

We also subdivided the data into eight smaller redshift–absolute magnitude bins and performed a global fit that includes a dependence on wavelength in the following form (see also Ivezic et al. 2004)


where , and are parameters of the fit, is the average absolute -band magnitude in a bin (the rest-frame magnitudes are taken from Assef et al. (2010a)), is the rest-frame time lag, and is the average rest-frame wavelength in each bin. We either fit (1) all the data or (2) exclude the short time-lags ( year) as the structure function contains little or no signal in this regime. For case (1) we found , , , , while for case (2) we found , , , . We find no detectable dependence of the structure function on wavelength, even though the data cover the redshift range of , the absolute magnitude range of to mag, and the rest-frame wavelengths µm. As compared to optical studies (e.g., MacLeod et al. 2010), the mid-IR structure function appears to have stronger dependence on the luminosity (parameter ), and no dependence on wavelength (parameter ). Note that the large uncertainty in () can potentially hide, within the 2–, the wavelength dependence of observed in optical studies. Another possibility is that the mid-IR wavelength dependence is much weaker than that observed in optical.

Table 7 shows the variability amplitudes at the observed time lag of 4 years for various classes of objects from Figure 4 assuming a fixed slope of . Objects with are AGNs, and show uniform and consistent variability amplitudes of mag as a function of magnitude in both channels for fixed . The X-ray selected sources inside the modified AGN wedge are expected to be AGNs, and not surprisingly, they also show similar variability amplitudes of mag to the objects. The mid-IR counterparts of the 24 µm (MIPS) selected sources inside the modified AGN wedge, again, have almost identical variability amplitudes. The radio selected objects, however, have a drastically different distribution in their mid-IR colors (top right panel in Figure 4) from that of AGNs, suggesting that a large fraction of them may not be AGNs. The structure functions based on the radio-selected objects inside the modified AGN wedge may be driven by non-AGN objects, and we find somewhat smaller variability amplitudes of mag and mag in [3.6] and [4.5], respectively.

6 Increasing Completeness with Future Surveys

Given the structure function, we can now estimate the expected statistical properties of the AGNs in the current survey or any extension by adding new epochs. The structure function is related to the correlation function by , where is the variance as . This makes it straightforward to compute the mean of an AGN since , where is the measurement uncertainty. Because we fit for the mean, , the exact value of is irrelevant. For our observed frame power-law structure function, the result for the SDWFS data and cadence is , where , mag, and is the typical photometric uncertainty. The constant is simply the number of degrees of freedom after fitting the mean, so a non-variable source would have . This expression makes it clear why our AGN completeness is so low and drops rapidly with magnitude (increasing , see Figure 2).

We can, in fact, simulate the completeness. We can decompose the correlation matrix into eigenvectors and eigenvalues . If is a Gaussian deviate of dispersion , then random realizations of light curves consistent with the structure function are


(see Rybicki & Press 1992). If we assume the variability in both bands is identical but the noise is independent, then we can compute the fraction of such AGNs we would select as a function of the noise using Monte Carlo simulations. Figure 10 shows the approximate completeness as a function of magnitude for sources with and either or . We use the observed structure function and the normalized measurement uncertainties to track the median values for [3.6] in Figure 2. For simplicity, the [4.5] band is treated as being identical to the [3.6] band. As might be expected from the structure of Figure 2, the completeness is high when the measurement errors are small compared to the amplitude of the structure function and then drops rapidly for fainter sources.

In general, the statistics of these simulated light curves agree with observations, indicating that our estimates of the structure function are correct. However, they also show that we are “pushing the envelope” of variability, because for most AGNs the variability amplitudes are merging into the measurement uncertainties, as can be seen in Figure 2.

7 Summary

In this paper, we performed DIA photometry on the four epochs of SDWFS data. The common area of the Spitzer mosaics covers 8.1 deg of the NDWFS Boötes field, and contains mostly extragalactic objects. We concentrate on the deeper 3.6 and 4.5 µm channels, and provide variability catalogs and light curves. These variability data are cross-correlated with the AGES redshift survey, X-ray and radio data, and 24 µm AGN candidates from the same area of the sky to study the variability of various classes of objects.

We study the mid-IR variability of objects based on the significance level of their joint variability () in both the [3.6] and [4.5] channels focusing on the sources that also show a high correlation () between the two channels. Our standard selection criteria with identify objects as variables, which constitute 1.1% of all detected sources in the field. We find that the majority (76%) of the mid-IR variable sources in the NDWFS Boötes field are AGNs, where we define AGNs to be objects that satisfy any of the following criteria: (1) lie inside the modified AGN wedge, (2) are X-ray, 24 µm or radio AGN candidates, (3) have , or (4) have an AGN spectroscopic template in any region.

Amongst all the extragalactic classes of objects, AGNs are the primary sources able to produce significant changes in their mid-IR luminosity over a several year timescale. Thus, mid-IR variability may be used as a new tool for selecting AGNs, especially if combined with other methods. However, in our relatively shallow few-epoch survey, only 15% of AGNs are sufficiently variable to be detected; therefore, variability-selected samples of AGNs are highly incomplete. This incompleteness is caused by the typical variability amplitude being comparable to the survey photometric errors, and is not due to the small number of available epochs. For our standard criteria, 14% of X-ray sources, 17% of 24 µm selected AGN candidates, and 15% of spectroscopically confirmed AGNs are detected by their variability.

The variable AGNs primarily occupy the mid-IR AGN selection wedge of Stern et al. (2005), with an extension to bluer colors where the host galaxy dominates the mid-IR colors of low-luminosity AGNs (see Gorjian et al. 2008; Assef et al. 2010b). If we redefine the original Stern et al. (2005) selection wedge to encompass X-ray or variability-selected AGNs, we find that variable objects bluer than mag are likely to be AGNs. This opens a new window for the Spitzer Warm Mission, where the mag color cut combined with the variability criterion will be a solid indication of AGN activity.

The structure function of AGNs is well described by a power law with a logarithmic slope of and amplitude mag in [3.6] and , mag in [4.5] for a rest-frame time lag of years. These structure functions are steeper than those observed in optical bands, and we argue that this is due to the inability of either disk or dust emission to respond on the short timescales seen in optical studies. The mid-IR structure functions go approximately as , which corresponds to the short timescale () behavior in the damped random walk model used by Kelly et al. (2009), Kozłowski et al. (2010), and MacLeod et al. (2010). If we examine the variability of various subsamples, we detect more variability in low luminosity objects observed at shorter rest-frame wavelengths. The mid-IR sources corresponding to X-ray and 24 µm AGN have structure functions identical to that of AGNs. However, this is not the case for the counterparts of radio sources which have mag in both [3.6] and [4.5] bands, as compared to for either objects, X-ray, or MIPS sources. One explanation may be that the radio sample has more contamination by non-AGN (e.g., due to confusion between radio cores and lobes), another is that the radio-selected AGNs are accreting in a different mode (lower accretion rates or efficiencies; e.g., Tasse et al. 2008) and that variability in such a mode is inherently lower.

Expanding the SDWFS survey with new epochs during the Spitzer Warm Mission would represent a significant improvement. Given our structure functions, we can estimate the results from expanding the SDWFS survey with new epochs. We considered two cases of (1) doing one new epoch in 2011 March and (2) doing new epochs in both 2011 and 2012 March. The improvements in the mean are significant, with and , respectively, for AGNs, compared to or 5 for non-variable objects. These represent a significant improvement over the SDWFS baseline, where for AGN and 3 for non-variable objects. Figure 10 shows the increases in completeness for these five and six epoch scenarios. Most of the gain for detecting variable sources comes from adding one additional epoch. The greatest astrophysical gain from two additional epochs comes from being able to better map out and subdivide the structure function by luminosity, redshift, and source type to better characterize the physics of AGN variability at these wavelengths. In particular, it would confirm the different slope of the mid-IR structure function from that observed in the optical.

This work is based on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of the Technology under contract with the National Aeronautics and Space Administration (NASA). Support for this work was provided by NASA through award numbers 1310744 (C.S.K. and S.K.), 1314516 (M.L.N.A.) and 1317692 (H.A.S.) issued by JPL/CalTech. CSK & SK are also supported by NSF grant AST-0708082. The NDWFS and the research of A.D. and B.T.J. are supported by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation. Support for M.B. was provided by the W. M. Keck Foundation.


  • Alard & Lupton (1998) Alard, C., & Lupton R. H., 1998, ApJ, 503, 325
  • Ashby et al. (2009) Ashby, M. L. N., et al. 2009, ApJ, 701, 428
  • Assef et al. (2010a) Assef, R. J., et al. 2010a, arXiv:1001.4529
  • Assef et al. (2010b) Assef, R. J., et al. 2010b, ApJ, 713, 970
  • Barvainis (1987) Barvainis, R. 1987, ApJ, 320, 537
  • Becker et al. (1995) Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559
  • Brand et al. (2006) Brand, K. et al., 2006, ApJ, 641, 140
  • Brown et al. (2006) Brown, M. J. I., et al. 2006, ApJ, 638, 88
  • Clavel et al. (1989) Clavel, J., Wamsteker, W., Glass, I. S., 1989, ApJ, 337, 236
  • Cool et al. (2008) Cool, R. J., et al., 2008, ApJ, 682, 919
  • Cutri et al. (1985) Cutri, R. M., Wisniewski, W. Z., Rieke, G. H., Lebofsky, M. J., 1985, ApJ, 296, 423
  • Cutri et al. (2003) Cutri, R. M., et al. 2003, The IRSA 2MASS All-Sky Point Source Catalog, NASA/IPAC Infrared Science Archive. http://irsa.ipac.caltech.edu/applications/Gator/,
  • de Vries et al. (2002) de Vries, W. H., et al., 2002, AJ, 123, 1784
  • de Vries et al. (2005) de Vries, W. H., Becker, R. H., White, R. L., Loomis, C. 2005, AJ, 129, 615
  • Eckart et al. (2010) Eckart, M. E., McGreer, I. D., Stern, D., Harrison, F. A., & Helfand, D. J. 2010, ApJ, 708, 584
  • Eisenhardt et al. (2004) Eisenhardt, P. R., et al., 2004, ApJS, 154, 48
  • Eisenhardt et al. (2010) Eisenhardt, P. R. M., et al. 2010, AJ, 139, 2455
  • Eyer (2002) Eyer, L. 2002, Acta Astron., 52, 241
  • Fabricant et al. (2005) Fabricant, D., et al. 2005, PASP, 117, 1411
  • Fazio et al. (2004) Fazio, G. G., et al., 2004, ApJS, 154, 10
  • Geha et al. (2003) Geha, M., et al. 2003, AJ, 125, 1
  • Glass (2004) Glass, I. S. 2004, MNRAS, 350, 1049
  • Gorjian et al. (2008) Gorjian, V., et al., 2008, ApJ, 679, 1040
  • Griffith & Stern (2010) Griffith, R. L., & Stern, D. 2010, AJ, submitted
  • Haardt et al. (1994) Haardt, F., Maraschi, L., & Ghisellini, G. 1994, ApJ, 432, L95
  • Hawkins (2002) Hawkins, M. R. S. 2002, MNRAS, 329, 76
  • Hickox et al. (2009) Hickox, R. C., et al. 2009, ApJ, 696, 891
  • Hook et al. (1994) Hook, I. M., McMahon, R. G., Boyle, B. J., Irvin, M. J., 1994, MNRAS, 268, 305
  • Houck et al. (2005) Houck, J. R., et al. 2005, ApJ, 622, L105
  • Ivezic et al. (2004) Ivezic, Ž., et al. 2004, The Interplay Among Black Holes, Stars and ISM in Galactic Nuclei, 222, 525
  • Jannuzi & Dey (1999) Jannuzi, B. T. & Dey, A., 1999, ASPC, 191, 111
  • Kelly et al. (2009) Kelly, B. C., Bechtold, J., & Siemiginowska, A. 2009, ApJ, 698, 895
  • Kenter et al. (2005) Kenter, A., et al. 2005, ApJS, 161, 9
  • Kollmeier et al. (2006) Kollmeier, J. A., et al. 2006, ApJ, 648, 128
  • Kowalski et al. (2009) Kowalski, A. F., Hawley, S. L., Hilton, E. J., Becker, A. C., West, A. A., Bochanski, J. J., & Sesar, B. 2009, AJ, 138, 633
  • Kozłowski & Kochanek (2009) Kozłowski, S., & Kochanek, C. S. 2009, ApJ, 701, 508
  • Kozłowski et al. (2010) Kozłowski, S., et al. 2010, ApJ, 708, 927
  • MacLeod et al. (2010) MacLeod, C. L., et al. 2010, arXiv:1004.0276
  • Murray et al. (2005) Murray, S., et al., 2005, ApJS, 161, 1
  • Neugebauer & Matthews (1999) Neugebauer, G. & Matthews, K., 1999, AJ, 118, 35
  • Rengstorf et al. (2004) Rengstorf, A. W., et al. 2004a, ApJ, 606, 741
  • Rybicki & Press (1992) Rybicki, G. B., & Press, W. H. 1992, ApJ, 398, 169
  • Sanders et al. (1989) Sanders, D. B., Phinney, E. S., Neugebauer, G., Soifer, B. T., Matthews, K., 1989, ApJ, 247, 29
  • Sesar et al. (2007) Sesar, B., et al. 2007, AJ, 134, 2236
  • Shakura & Sunyaev (1973) Shakura, N. I. & Sunyaev, R. A., 1973, A&A, 24, 337
  • Soszyński et al. (2004) Soszyński, I., Udalski, A., Kubiak, M., Szymański, M., Pietrzyński, G., Żebruń, K., Szewczyk, O., & Wyrzykowski, Ł. 2004, Acta Astronomica, 54, 129
  • Soszyński (2007) Soszyński, I. 2007, ApJ, 660, 1486
  • Soszyński et al. (2009) Soszyński, I., et al. 2009, Acta Astronomica, 59, 239
  • Stern et al. (2005) Stern, D., et al., 2005, ApJ, 631, 163
  • Stern et al. (2007) Stern, D., et al. 2007, ApJ, 663, 677
  • Suganuma et al. (2006) Suganuma, M., et al. 2006, ApJ, 639, 46
  • Tasse et al. (2008) Tasse, C., Best, P. N., Röttgering, H., & Le Borgne, D. 2008, A&A, 490, 893
  • Vanden Berk et al. (2004) Vanden Berk, D. E., et al., 2004, ApJ, 601, 692
  • Woźniak (2000) Woźniak, P. R., 2000, Acta Astronomica, 50, 421
Figure 1: Examples of data filtering and image subtraction. All six stamps cover the same area (1.21 arcmin) and are centered on one of the most variable AGNs in the Boötes field, SDWFS J143032.42+332158.8, with , , and (see §3). Top row: stamp of the original 3.6 µm image (a). Before we do the image subtraction, we filter the images in order to convert the Spitzer IRAC’s triangle-shaped PSFs into more Gaussian-like ones (b). Middle row: subtracted images for the 3.6 µm channel. The AGN brightened by about 0.5 mag on a baseline of 4 years between epochs 1 (c) and 3 (d). Black (white) color here reflects the source emitting less (more) flux than in the template image defined by the average of all four epochs. All non-variable objects are subtracted out. Bottom row: subtracted images for the 4.5 µm channel. The AGN brightened by mag between epochs 1 (e) and 3 (f), consistent with what is seen in the 3.6 µm channel.
Figure 2: Variability in the 3.6 µm (top) and 4.5 µm (bottom) channels as a function of magnitude for all objects with four epochs of data shown as contours. The objects are binned into 0.1 mag and 0.05 log() bins. The contours are drawn for 2, 10, and 100 objects per bin, counting from the outer contour. The lower solid line shows the median calculated in 0.1 mag bins, and the upper three dashed, dash-dotted, and dotted lines show variability significances , 2, and 3 (from bottom to top) with respect to this median. We limit our Photometry Group 2 and 3 samples to objects with mag, and mag (vertical dashed lines). Dots mark Photometry Group 3 AGNs from the modified AGN wedge, selected as either X-ray or 24 µm AGN candidates or having a QSO spectroscopic template code. Note that they tend to have higher average variability than the typical sources, but with amplitudes easily masked by the photometric uncertainties. Open circles mark objects with and .
Figure 3: Correlation between the 3.6 µm and 4.5 µm light curves for sources with measured colors (“Photometry Group 3” sources; see Section 3) with four epochs of data for Variability Level 2, 3, and 4 (, , and ) sources. We only consider sources in the shaded, highly correlated region as true variables, and we use the density of the poorly correlated sources to estimate our false positive rates (Tables 35).
Figure 4: SDWFS mid-IR color-color diagrams for the Boötes field. The green contours show the distribution of all 39,522 Photometry Group 3 sources. The objects are binned in 0.05 mag bins in both axes, and the contours are drawn at levels of 2, 10, 20, 50, 100, and 150 objects per bin, counting from the outer contour. Top-left panel: black dashed and solid line “wedges” show the original Stern et al. (2005) AGN selection region and the modified AGN region we define to encompass the distribution of XBoötes X-ray sources discussed in Gorjian et al. (2008). The black box at color zero is the region occupied by normal stars and also nearby galaxies without PAH emission. The blue solid, dotted, and dashed lines show the evolution of the mid-IR colors of late-type spiral, irregular, and early-type galaxies with redshift based on the SED models of Assef et al. (2010b). The red dotted and solid lines show the evolution of a pure AGN and a composite with equal contributions from the AGN and its late-type host to the bolometric luminosity. The dots on each track mark redshift increments of 0.5 (1.0) from to 3 (to 6) for galaxies (AGNs) that start at the bulls eye. Other panels: Photometry Group 3 mid-IR counterparts of radio sources (top right), X-ray sources (bottom left), and MIPS 24 µm AGN candidates (bottom right) are shown.
Figure 5: Left panel: mid-IR color-color distribution of all Photometry Group 3, Variability Level 2 sources in SDWFS. The contours show the distribution of all Photometry Group 3 sources, as described in Figure 4. The Stern et al. (2005) AGN selection region is shown by the dashed wedge, and the modified AGN wedge is shown with the solid line. The stellar box encompasses sources with (stellar) colors near zero. The majority of mid-IR variable objects (76%) are AGNs. Also 11% of objects in the stellar box are variable. Right panel: fraction of Photometry Group 2 sources (e.g., sources lacking reliable and/or photometry) that are Variability Level 2 variables as a function of color.
Figure 6: Left panel: mid-IR color-color distribution of the Photometry Group 3, Variability Level 2 sources selected as AGNs based on their X-ray (squares), radio (triangles), or 24 µm emission (dots). The contours and selection regions are the same as in Figure 5. Right panel: histograms of the fraction of the Photometry Group 2 X-ray (solid), 24 µm (dashed), and radio (dotted) sources identified as Variability Level 2 variables as a function of color.
Figure 7: Left panel: mid-IR color-color distribution of the Photometry Group 3, Variability Level 2 sources with spectroscopic redshifts from AGES. Stars, circles, and squares correspond to objects with stellar, galaxy, or AGN spectroscopic template matches. The contours show the distribution of galaxies (gray), stars, and AGNs with AGES redshifts. Right panel: fraction of objects with stellar, galaxy, or AGN template matches (dotted, dashed, and solid line, respectively) that are the Photometry Group 2, Variability Level 2 variables as a function of color.
Figure 8: Mid-IR AGN structure functions (Equation (5)) in the observed frame (top) and rest-frame (bottom). They are based on , 690, and 1100 SDWFS objects with mag (dots), mag (triangles), and mag (squares), respectively. We also show the i-band rest-frame structure function from Vanden Berk et al. (2004) in the lower panels.

Figure 9: Rest-frame mid-IR AGN structure functions for the objects with mag. We divide the objects into two redshift bins, (filled symbols) and (open symbols), and also into two brightness levels, objects brighter (circles) and fainter (squares) than the median for the redshift bin. For time lags of half a year and longer, fainter and lower redshift sources show greater variability. The solid lines are the mid-IR structure function fits from the bottom panel of Figure 8, and the dotted lines are the i-band rest-frame structure functions from Vanden Berk et al. (2004).

Figure 10: Detection completeness for variable AGNs as a function of magnitude for sources with high correlation and either or (Variability Level 2 or 4). We show the completeness for the current four epoch SDWFS survey, for the current survey plus one epoch in 2011 March (five epochs), and for the current survey plus two epochs in 2011 and 2012 March (six epochs).
Object Name R.A. Decl. [3.6] [4.5] [5.8] [8.0] Mask
(deg) (deg) (mag) (mag) (mag) (mag) (mag) (mag) (mag)
SDWFS J143040.36+322721.8 217.668148 32.456048 17.666 17.488 16.792 16.070 4 0.051 0.062 0.001 0.178 0.894 1
SDWFS J143034.08+322729.0 217.642013 32.458058 17.937 17.466 17.142 16.928 4 0.067 0.185 0.002 0.127 0.758 1
SDWFS J143033.71+322730.3 217.640453 32.458420 18.725 18.352 17.576 16.869 4 0.181 0.223 0.029 0.718 0.392 1
SDWFS J143047.11+322729.5 217.696285 32.458194 17.121 16.995 16.898 14.735 4 0.036 0.096 0.002 0.511 0.295 1
SDWFS J142951.47+322726.0 217.464455 32.457230 18.361 17.875 17.538 16.447 4 0.135 0.102 0.009 0.627 0.734 1
SDWFS J143015.91+322733.8 217.566284 32.459402 19.732 20.200 18.262 99.999 1 99.999 99.999 99.999 99.999 99.999 1

Note. – This is an abridged version of the full table, presenting the variability statistics for the combined channels [3.6] and [4.5] only. The electronic table includes additional columns showing the variability statistics for the separate channels [3.6] and [4.5], and the photometric uncertainties for the IRAC channels. Here, quantities and , are respectively, the variability standard deviation in channel and the common number of epochs in both channels used in the calculation of the joint variability. With , , and we denote the covariance and the correlation between the channels 1 and 2, and the joint variability significance, respectively. The mask is 0 (1) if there is (is no) bright star with the 2mass -band magnitude within arcsec affecting the measurement. The Vega magnitudes presented in the table are measured in 4 arcsec apertures (see Ashby et al. 2009). The error code for magnitudes, reflecting no measurement, is 99.999. The error code for , , , and is 99.999.
(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.)

Table 1: SDWFS-DIA Variability Catalog.
Channel 1 Channel 2
R.A. Decl. Epoch 1 Epoch 2 Epoch 3 Epoch 4 Epoch 1 Epoch 2 Epoch 3 Epoch 4
(deg) (deg) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag)
217.668148 32.456048 17.703 17.823 17.741 17.771 17.457 17.518 17.547 17.409
217.642013 32.458058 17.894 18.027 17.881 17.952 17.290 17.455 17.566 17.728
217.640453 32.458420 19.037 18.793 18.769 18.598 18.721 18.220 18.288 18.371
217.696285 32.458194 17.107 17.128 17.132 17.191 17.019 16.995 17.195 17.143
217.464455 32.457230 18.385 18.070 18.166 18.148 17.925 17.780 17.696 17.876
217.566284 32.459402 20.055 19.751 20.271 19.649 99.999 99.999 19.692 99.999

Note. – The error code for magnitudes, reflecting no measurement, is 99.999.
(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.)

Table 2: The Light Curves for the SDWFS-DIA Variability Catalog.
Variability Variability Old AGN Wedge Modified AGN Wedge Stellar Box Outside All Abjects
Number (fraction/false positive rate) of variable objects: Photometry Group 1
Total 43468 76005 45939 144263 474179
2 746 (1.7%/28%) 1015 (1.3%/33%) 555 (1.2%/58%) 832 (0.6%/71%) 5107 (1.1%/57%)
3 397 (0.9%/17%) 501 (0.7%/22%) 226 (0.5%/54%) 326 (0.2%/56%) 2208 (0.5%/47%)
4 246 (0.6%/11%) 293 (0.4%/17%) 102 (0.2%/51%) 156 (0.1%/53%) 1071 (0.2%/44%)
Number (fraction/false positive rate) of variable objects: Photometry Group 2
Total 29181 54920 18150 113299 213594
2 598 (2.0%/18%) 774 (1.4%/24%) 133 (0.7%/65%) 512 (0.5%/74%) 1557 (0.7%/50%)
3 340 (1.2%/9%) 405 (0.7%/13%) 58 (0.3%/56%) 166 (0.1%/56%) 668 (0.3%/30%)
4 219 (0.8%/5%) 250 (0.5%/9%) 24 (0.1%/42%) 71 (0.1%/55%) 359 (0.2%/24%)
Number (fraction/false positive rate) of variable objects: Photometry Group 3
Total 6677 12741 7567 19214 39522
2 496 (7.4%/6%) 584 (4.6%/7%) 86 (1.1%/58%) 105 (0.5%/44%) 775 (2.0%/19%)
3 310 (4.6%/3%) 349 (2.7%/4%) 44 (0.6%/48%) 35 (0.2%/44%) 428 (1.1%/12%)
4 204 (3.1%/2%) 224 (1.8%/3%) 21 (0.3%/33%) 15 (0.1%/30%) 260 (0.7%/9%)

Note. – Photometry Groups 1 and 2 are defined in Section 3; only 49% and 86% of objects in these groups have all four IRAC band measurements, respectively. The remaining objects lack either [5.8], [8.0], or both measurements, so we are unable to place them on the mid-IR color-color diagram (i.e., we were unable to tell if they fall into the modified/old AGN wedge or stellar box). By definition, Photometry Groups 3 includes only objects with all four IRAC measurements.
Column “Outside” includes objects with all four IRAC band measurements only that are both outside of the “Modified AGN Wedge” and the “Stellar Box.” Column “All Objects” contains all objects, regardless of the number of IRAC band measurements.
The fraction of variable sources is a ratio of the number of variable objects above a given level of and to the total number of objects in a respective region. The false positive rate is estimated from the number of sources above the same level of but with (see Section 3 and Figure 3).

Table 3: Mid-IR Variable Objects in the SDWFS.
Variability Variability MIPS QSOs X-ray FIRST WSRT
Number (fraction/false positive rate) of variable objects: Photometry Group 1
Total 1964 3231 313 2046
2 230 (11.7%/5%) 392 (12.1%/4%) 12 (3.8%/8%) 48 (2.3%/13%)
3 151 (7.7%/3%) 253 (7.8%/2%) 8 (2.6%/3%) 27 (1.3%/8%)
4 107 (5.4%/2%) 164 (5.1%/1%) 4 (1.3%/1%) 17 (0.8%/4%)
Number (fraction/false positive rate) of variable objects: Photometry Group 2
Total 1942 2950 292 1883
2 229 (11.8%/5%) 391 (13.3%/3%) 12 (4.1%/7%) 48 (2.5%/10%)
3 151 (7.8%/3%) 253 (8.6%/2%) 8 (2.7%/3%) 27 (1.4%/6%)
4 107 (5.5%/2%) 164 (5.6%/1%) 4 (1.4%/1%) 17 (0.9%/3%)
Number (fraction/false positive rate) of variable objects: Photometry Group 3
Total 1637 2217 195 1233
2 225 (13.7%/4%) 374 (16.9%/2%) 11 (5.6%/4%) 43 (3.5%/9%)
3 151 (9.2%/3%) 245 (11.1%/1%) 8 (4.1%/1%) 25 (2.0%/5%)
4 107 (6.5%/1%) 160 (7.2%/1%) 4 (2.1%/1%) 15 (1.2%/4%)

Note. – The fraction of variable sources is a ratio of the number of variable objects above a given level of and to the total number of objects in a respective class. The false positive rate is estimated from the number of sources above the same level of but with (see Section 3 and Figure 3).

Table 4: Mid-IR Variability of Photometrically-Selected AGNs.
Variability Variability QSOs Stars Galaxies
Number (fraction/false positive rate) variable objects: Photometry Group 1
Total 2950 1121 18161
2 334 (11.3%/5%) 20 (1.8%/20%) 209 (1.1%/25%)
3 216 (7.3%/3%) 11 (1.0%/7%) 108 (0.6%/15%)
4 144 (4.9%/2%) 6 (0.5%/5%) 60 (0.3%/10%)
Number (fraction/false positive rate) of variable objects: Photometry Group 2
Total 2820 966 17405
2 333 (11.8%/4%) 18 (1.9%/14%) 200 (1.1%/24%)
3 216 (7.7%/2%) 10 (1.0%/5%) 102 (0.6%/13%)
4 144 (5.1%/2%) 6 (0.6%/4%) 58 (0.3%/8%)
Number (fraction/false positive rate) of variable objects: Photometry Group 3
Total 2192 491 11214
2 319 (14.6%/4%) 14 (2.9%/9%) 175 (1.6%/17%)
3 211 (9.6%/2%) 10 (2.0%/3%) 93 (0.8%/10%)
4 142 (6.5%/1%) 6 (1.2%/3%) 56 (0.5%/5%)

Note. – The fraction of variable sources is a ratio of the number of variable objects above a given level of and to the total number of objects in a respective class. The false positive rate is estimated from the number of sources above the same level of but with (see Section 3 and Figure 3).

Table 5: Mid-IR Variability by AGES Spectroscopic Classification.
Structure Function Parameters
Level Channel
Observed Frame ( years)
mag [3.6]
mag [3.6]
mag [3.6]
mag [4.5]
mag [4.5]
mag [4.5]
Rest-Frame ( years)
mag [3.6]
mag [3.6]
mag [3.6]
mag [4.5]
mag [4.5]
mag [4.5]
Rest-Frame ( years)
bright [3.6]
faint [3.6]
bright [3.6]
faint [3.6]
bright [4.5]
faint [4.5]
bright [4.5]
faint [4.5]
Table 6: Mid-IR AGN Structure Functions
Amplitude Parameters
Level [3.6] (mag) [4.5] (mag)
objects ( years)
X-ray sources inside the modified AGN wedge ( years)
MIPS sources inside the modified AGN wedge ( years)
Radio (WSRTFIRST) sources inside the modified AGN wedge ( years)
Table 7: Mid-IR Variability Amplitudes At Fixed .
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description