A number of astronomical standard stars were observed to obtain a valid absolute flux calibration. Stars with a range of fluxes were observed at a number of positions across the array many times throughout the mission, to monitor any changes that may have occurred. Spectral types and accurate absolutely calibrated fluxes in the IRAC channels were determined for the calibration stars. Their flux densities are given in Table 4.1. All of the calibration data taken with these stars are public, and are available in the Spitzer Heritage Archive.
The observational program to measure the flux calibration of IRAC and to ensure its constancy had two components. The primary calibrator stars in the continuous viewing zone (CVZ) were observed in every campaign, and used to derive the overall flux calibration. The secondary calibrator stars were used to monitor short-term variations in the absolute calibration. To avoid slew overheads, they were observed close to downlinks, and, therefore, had to be located near the ecliptic plane, within a tightly constrained window of about 20 degrees in width. Because of the motion of Earth about the Sun, this window constantly moved, and so any one secondary calibrator was visible and near the downlink pointing for only a campaign or two per year. Using the secondary standard stars, the IRAC calibration was measured to be stable over the mission.
The absolute calibration was derived by means of aperture photometry, using a 10 native pixel radius (≈ 12 arcseconds) aperture. The background was measured using a robust average in a 12 - 20 native pixel annulus around the centroid of the star. The steps to remove known systematics are described in detail by Reach et al. (2005), Hora et al. (2008) and Carey et al. (2012). To obtain photometry at the highest possible accuracy, photometric corrections for the location of the peak of the source within a single pixel, and the location of the source within the array, have to be made. If this methodology is not applied, then point source photometry from the Level 1 products (BCDs) can be in error by up to 10%.
Analysis of the flux calibrator data indicated that absolute flux calibration is accurate to 3%. Repeatability of the measurements of the individual stars is better than 1.5% (dispersion), and can be as good as 0.01% with very careful observation design (e.g., Charbonneau et al. 2005). Unfortunately, ground-based infrared calibrators were too bright to use as calibrators for IRAC. Therefore, we had to use models to predict the actual flux for each channel as a function of star spectral type. Uncertainty in those models, and, therefore, in the absolute flux densities of the calibration stars, are 2% – 3% (Cohen et al. 2003). Please note that if you measure aperture photometry on IRAC images, then in our nomenclature you will have measured F*K*(where F* is the flux density of the calibration source at the nominal wavelength, and K*is the color correction factor for the calibration source spectrum, see Section 4.4). If a comparison of values at the nominal IRAC channel wavelengths (listed in
Table 4.3) to measurements of the star at other wavelengths is desired (e.g., the spectral energy distribution), then the F*K* values should be divided by the corresponding K* values in Table 4.1.
Table 4.1: IRAC primary calibrator quoted flux densities and color corrections.
Calibration Source Name(s)
Stellar Type
Ks mag
Quoted Flux Densities F* (mJy) and Color Corrections K*
3.6 µm
F* K*
4.5 µm
F* K*
5.8 µm
F* K*
8.0 µm
F* K*
NPM_67.0533
SAO17718
2MASS J17585466+6747368
K2III
6.4
804.4
1.0037
454.8
1.0614
318.8
1.0023
178.6
1.0554
HD 165459
2MASS J18023073+5837381
A1V
6.6
645.0
1.0055
422.6
0.9992
267.4
1.0114
145.4
1.0266
NPM1+68.0422
BD+68 1022
2MASS J18471980+6816448
K2III
6.8
554.6
1.0031
305.4
1.0711
218.5
1.0005
121.0
1.0684
KF09T1
GSC 04212-01074
2MASS J17592304+6602561
K0III
8.1
165.8
1.0050
97.65
1.0459
66.42
1.0054
37.26
1.043
NPM1+64.0581
HD 180609
2MASS J19124720+6410373
A0V
9.1
62.86
1.0060
41.09
>1.0011
25.97
1.0142
14.16
1.0252
NPM1+60.0581
BD+60 1753
2MASS J17245227+6025508
A1V
9.6
38.10
1.0059
24.82
0.9994
15.67
1.0114
8.506
1.0304
KF06T1
2MASS J17575849+6652293
K1.5III
11.0
13.24
1.0040
7.567
1.0573
5.261
1.0022
2.955
1.051
KF08T3
2MASS J17551622+6610116
K0.5III
11.1
11.85
1.0047
6.922
1.0492
4.738
1.0047
2.662
1.0447
KF06T2
2MASS J17583798+6646522
K1.5III
11.3
10.22
1.0040
5.840
1.0573
4.061
1.0031
2.280
1.051
2MASS J18120957+6329423
A3V
11.3
8.663
1.0055
5.685
0.9986
3.611
1.0084
1.959
1.0252
Table 4.2 lists the calibration factors that were used in the final processing of all IRAC data.
Table 4.2: The photometric calibration and zero magnitude flux densities for IRAC. Warm mission corrections are given as the second set of values in channels 1 and 2. Note that at the beginning of the warm mission the calibration was changing frequently. Different FLUXCONV values for these first months can be found on IRSA’s IRAC web pages (under “Warm IRAC Characteristics”).
λ (μm)
FLUXCONV (MJy/sr)/(DN/second)
Cryogenic/Warm
Fν0 (Jy)
3.6
0.1069±0.0026
0.1257±0.0030
280.9±4.1
4.5
0.1382±0.0028
0.1447±0.0029
179.7±2.6
5.8
0.5858±0.0100
115.0±1.7
8.0
0.2026±0.0043
64.9±0.9
Table 4.2 lists the calibration factors derived from the primary calibration program that were used in the final processing of all IRAC data. These calibration factors can be used to convert native IRAC data units (DN/second) into units of surface brightness (MJy/sr). Note that the zero magnitude flux densities listed below are evaluated at exactly 3.6, 4.5, 5.8, and 8.0 µm and not at the nominal channel wavelengths given in Table 4.3. To estimate the zero magnitude flux densities at the nominal channel wavelengths, scale the values in Table 4.2 by the square of the ratio of the wavelengths. For example, in channel 1 the zero-magnitude flux density at 3.544 μm would be (3.544/3.6)2 × 280.9 Jy = 272.2 Jy. The difference between the two is a few percent.
Note that IRAC was not an absolute background photometer, so the total brightness in IRAC images should be used with great caution. There was a cold shutter in the cold assembly, but it was not operated in flight, in order to minimize the mission risk. Therefore, in the cryogenic mission the offset level in IRAC images was referenced to laboratory measurements before the launch, where the offset level was observed to change very significantly from one laboratory experiment to another. In laboratory tests, the absolute offset of IRAC images was found to vary at levels that are comparable to the minimum celestial background in channels 1 and 2. Furthermore, the offset level changed depending on whether the detector was recently annealed (annealing means heating the detectors to 23 - 27 K for a short time with a small current flowing through the detectors. It restores the distribution of charges that are trapped to the same distribution each time, but does not restore most hot pixels; annealing was done only during the cryogenic mission). Thus, for diffuse surface brightness measurements, we recommend making differential measurements among at least two sky positions, preferably from the same campaign.