Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.stsci.edu/hst/wfc3/documents/ISRs/WFC3-2004-07.pdf Дата изменения: Fri Apr 1 01:15:13 2005 Дата индексирования: Sun Dec 23 01:16:44 2007 Кодировка: Поисковые слова: m 35

Instrument Science Report WFC3 2004-07

Predicted Countrates for the UV WFC3 Calibration Subsystem using Deuterium Lamps
S.Baggett, J. Sullivan, and M. Quijada May 17,2004

ABSTRACT Predicted WFC3 calibration subsystem countrates have been computed for a number of candidate deuterium (D2) lamps. Based on the calibration system component throughput estimates from Ball, the detector curves from the GSFC Detector Characterization Lab (DCL), and the lamp spectra measured at GSFC, use of any of the lamps in the calibration subsystem should provide sufficient flux to satisfy the CEI specification of 16.7 e-/s/pix in all filters blueward of and including F390N. This report summarizes the various lamps available and provides details of the procedure used for deriving the estimated countrates.

Introduction
Deuterium lamps manufactured by IST (Imaging and Sensing Technology) are planned for use in the WFC3 calibration subsystem (calsystem) to provide flatfields from ~2000е to 4000е. Initially, four D2 lamps were delivered to GSFC as possible candidates for the calsystem; since then, two have suffered weld problems and are no longer viable. As a result of the welding issues, IST revised their production procedures in order to fabricate more rugged lamps; a batch of five D2 candidate lamps with improved welds has been delivered to GSFC and one is currently undergoing vibration tests. This report summarizes the various lamps available for the calsystem, presents predicted calsystem countrates for filters

Operated by the Association of Universities for Research in Astronomy, Inc., for the National Aeronautics and Space Administration


Instrument Science Report WFC3 2004-07 out to F390N, and provides details on the procedure used to derive the throughput estimates.

Data
Figure 1 presents spectra of the nine D2 lamps in question. The first four delivered were 31302, 31305, 31306, and 31298; the lamps with improved welds are 31774, 31775, 31776, 31886, and 31887. Of the first four lamps, number 31302 had one of the highest fluxes and was shipped to Ball for possible integration into the calsystem; the other high flux lamp (31305) suffered a broken anode weld during vibration testing and is no longer viable. Lamp 31298 was not considered a candidate for the calsystem, as not only did it have the lowest flux of the initial four lamps, it also exhibited some spectral features of unknown origin; since then, it too has suffered from a broken weld. The intermediate flux lamp 31306 was used in a lifetime test (Boucarut et al., 2003; Baggett & Quijada, 2003) and thus is also no longer a candidate for use in the calsystem. The five new design D2 lamps all have similar spectra but with relatively low output; their spectra cluster around the level of lamp 31298, the lowest flux lamp of the first four delivered. Figure 1: Spectra of nine D2 lamps. The top plot shows the entire wavelength range as measured; the bottom plot highlights the WFC3 calsystem wavelength regime.
2 1.8 1.6

Spectral Scans All Viable D2 Lamps
31302 (Flight Candidate) 31305(Broken) 31306 (After Burn-in) 31298(Broken) 31774 31775 31776 31886 31887

Radiance ( W/cm^2/nm)

1.4 1.2 1 0.8 0.6 0.4 0.2 0 150

200

250

300

350

400

Wavelength (nm)
0.5
31302 (Flight Candidate) 31305(Broken) 31306 (After Burn-in) 31298(Broken) 31774 31775 31776 31886 31887

Radiance ( W/cm^2/nm)

0.4

0.3

0.2

0.1

0 200

250

300

350

400

Wavelength (nm)

2


Instrument Science Report WFC3 2004-07 Initial calsystem throughputs from Ball The Ball Systems Engineering Report OPT-060 (Sullivan, 2001) provided predictions of the calibration system fluxes at the UVIS detector for filters out to F393N using an IST D2 spectral curve from the Cosmic Origins Spectrograph (COS) program. Table 1 in that report listed predicted countrates in each filter; these estimates were computed via an Excel spreadsheet, using NIST spectral radiance estimates, radiometric analysis by D.Ebbets, optical design as revised by M. Dittman, filter data from JPL, and additional updates by J.Sullivan in Jan 2003. In that particular analysis, each calsystem observing mode was represented by a set of terms at the central wavelength of the filter in question (e.g., detector quantum efficiency, optics throughputs, and lamp intensity), which were multiplied by other necessary corrective factors (inclination of the detector to the line of sight, f/# of the light from the D2 lamp lens, and the total number of pixels) to derive the countrate estimates. Transfer of Ball calsystem throughput estimates into SYNPHOT For ease of use at STScI, the Excel spreadsheet computations were transferred into SYNPHOT, a Synthetic Photometry package available as part of STSDAS/IRAF1. SYNPHOT is commonly used at STScI for providing photometric calibration for HST data but it can also be used to provide countrate predictions. The SYNPHOT system consists of a set of throughput tables containing optical element throughput values as a function of wavelength, along with a graph table which defines which types of throughput tables to use for a given observing mode, and a component table, which provides specific names of throughput tables to use. To compute a predicted countrate, SYNPHOT identifies the necessary throughput tables for the specified mode, multiplies the files together, and sums the results across the bandpass. Once a SYNPHOT system has been established, it becomes relatively easy to swap out specific components and re-run countrate predictions. The current WFC3 SYNPHOT system, designed for estimating countrates from external sources, was used as a basis for generating a separate, calsystem-specific SYNPHOT system. Two terms are specified in the WFC3 SYNPHOT external observing mode: channel (IR or UVIS) and filter; these terms map to several tables for computing total throughput, namely, the HST OTA throughput table, WFC3 optics, filter throughput, and detector QE tables. For the calsystem, the HST OTA term is not needed; it was removed from the SYNPHOT throughput path in the calsystem graph table. A new optics table was generated, to include the calsystem-specific components such as the UVIS beam splitter and the UVIS parabola, using the values given in the SER OPT-060-based Excel spreadsheet (Sullivan, STSDAS is the Space Telescope Science Data Analysis System, which is part of IRAF, the Image Reduction and Analysis Facility, a general purpose software for the analysis of astronomical data

1.

3


Instrument Science Report WFC3 2004-07 2003). Those contributions, listed in the SER as a function of filter, were converted for SYNPHOT to a table of throughput as a function of filter central wavelength. Individual filter curves from the current SYNPHOT used for estimating external source countrates were kept as-is, since those contain the actual measured throughputs across each filter as a function of wavelength. For validation purposes, a SYNPHOT detector QE table and a D2 lamp spectrum were populated with the values used in the Excel spreadsheet (Sullivan, 2003), also as a function of filter central wavelength. Once the calsystem SYNPHOT setup could be shown to reproduce the countrates given in OPT-060 and the Excel spreadsheet, the QE and lamp tables could be updated to reflect the most recent data available for the flight build. Since the Excel terms were given as a function of filter central wavelength while SYNPHOT requires the tables to be populated out to the shortest and longest wavelengths in the system, some extrapolations in throughput and wavelength were necessary, from 2222е down to 1800е and from 3933е up to 4500е. A straight extrapolation to shorter wavelengths of the calsystem optics and the initial lamp curve gave rise to very steep curves, resulting in predicted SYNPHOT countrates significantly higher than the rates predicted in SER OPT-060, so the slope of those extrapolations were reduced somewhat; the extensions that were used in the calsystem SYNPHOT tables are shown in Figure 2 below. Figure 2: Initial SYNPHOT curves, based upon SER-060; dashed lines mark areas outside of which extensions were necessary (the QE required extrapolation on the blue side only). Calsystem optics, initial QE, and initial lamp curves are shown in the left, middle, and right plots, respectively.
wfc3_uvis_d2
initial QE .6

initial lamp

2.00E-10

3.00E11
.5

throughput

1.80E-10

throughput

.4

1.60E-10

.3

.2

1.40E-10 2000 2500 3000 3500 4000 wavelength (angstroms) 4500

2000

2500 3000 3500 4000 wavelength (angstroms)

4500

flux

2.00E11

1.00E11

2000

2500

3000 3500 4000 wavelength (angstroms)

4500

Using these extended curves along with the SYNPHOT filter throughput tables, a first set of SYNPHOT countrate predictions were computed. Some discrepancies were found between these initial SYNPHOT countrates and the Excel spreadsheet countrates, which were traced to differences between SYNPHOT filter throughput and/or bandwidth and filter throughput assumed in the Excel spreadsheet. Since the SYNPHOT filter tables are values measured at GSFC (Lupie and Boucarut, 2003), the throughputs and bandwidths

4


Instrument Science Report WFC3 2004-07 from the SYNPHOT tables were used to update the Excel spreadsheet quantities and the Excel countrates were recomputed.

Results
Verification of Migration to SYNPHOT system One observing mode As an example, the countrate estimates for F280N are shown below. The prediction from the Excel spreadsheet (Sullivan, 2003) is 60.7 e-/s/pix, computed using the following terms
filter bandwidth = 42е lamp intensity = 1e11 photons/s/е optics efficiency = 0.26 beam splitter = 0.42 mirror = 0.87 filter = 0.2 window = 0.92 QE = 0.47 (p/4) * ( 1 / f/# )2 = 0.031416 steradians cos(20) = 0.939693 (detector angle to the line of sight) 4K x 4K chip = 16777216 pixels

for a final countrate for the F280N D2 lamp calsystem observing mode of 42 в 1 e 11 в 0.26 в 0.42 в 0.87 в 0.2 в 0.92 в 0.47 в 0.031416 в 0.9397 / 4096 = 60.7 In the initial calsystem SYNPHOT build, using the Excel spreadsheet preliminary lamp spectrum and QE, the resulting countrate was 62.1 e-/s/pix:
calcphot wfc3,uviscalsys,d2,f280n spectrum=uvis.d2lamp.tab form=counts wave=wave.tab Mode = band(wfc3,uviscalsys,d2,f280n) Pivot Equiv Gaussian Wavelength FWHM 2796.218 31.37834 band(wfc3,uviscalsys,d2,f280n) Spectrum: uvis.d2lamp.tab VZERO (COUNTS s^-1 hstarea^-1) 0. 62.12714
2

or 62.1 e-/s/pix, within ~2% of the Excel spreadsheet value. The SYNPHOT path for F280N included the new wfc3 uvis calsystem lightpath table based on the contributions listed in the Excel spreadsheet, not including the lamp, the F280N filter throughput as a function of wavelength (from the pre-existing SYNPHOT system), the Excel-based detector QE table, and a spectrum (uvis.d2lamp.tab) based upon the lamp fluxes in the Excel spreadsheet.

5


Instrument Science Report WFC3 2004-07 All filters Calsystem countrate comparisons were done for all applicable filters by ratioing the countrate predictions from the Excel spreadsheet, using the updated bandwidths and filter throughputs, to the initial SYNPHOT build; the results, as a function of filter central wavelength, are presented in Figure 3. The left plot shows the first ratio results, using the SYNPHOT tables as defined to this point. Ideally, however, the SYNPHOT predictions should match the Excel predictions as closely as possible; to minimize the differences between the two systems, the initial ratios were fit as a function of central wavelength with a second-order polynomial and the slight curvature removed from the countrate ratios by including a correction based upon the fit in the SYNPHOT uvis calsystem table. Ratios of the resulting new SYNPHOT countrates to the Excel countrates are shown in the right plot, along with a fit to the new set of ratios. The inital ratios ranged from 0.89 to 1.5, with a standard deviation of ~0.18; the final ratios ranged from 0.8 to 1.3, with a standard deviation of 0.15. Figure 3: Ratio of SYNPHOT to Excel spreadsheet countrate estimates for the D2 calsystem as a function of central wavelength. Left plot shows ratios using countrates from the first SYNPHOT build, while the right plot shows ratios using the adjusted SYNPHOT calsystem uvis table. p
2 2

1.5

1.5

ratio (SYNPHOT/Excel)

ratio (SYNPHOT/Excel)

1

1

.5

.5

0 2500 3000 wavelength 3500

0 2500 3000 wavelength 3500

Updating CCD QE Curves and Lamp Spectra Once the SYNPHOT and Excel countrates were in reasonable agreement (see Figure 3), the new lamp curves and flight CCD QE curves were swapped into SYNPHOT. The spectra for the new lamps (shown in Figure 1), as measured in the optics lab at GSFC, were first converted from µW/cm2/nm to photons/sec/е, then multiplied by a factor of 5 to account for the GSFC detector intercepting only 20% of the lamp light. For QE response, the detector QE curves measured at the GSFC Detector Characterization Lab (Landsman, 2002) were used. The DCL curves contain a correction for the quantum yield, to account

6


Instrument Science Report WFC3 2004-07 for the tendency for single UV photons shortward of ~3400е to generate more than one electron. In order to estimate calsystem countrates expected at the detector, the quantum yield correction was removed from the QE curves, following the prescription provided by DCL (Hill, 2003); Figure 4 presents the adjusted QE curves used in SYNPHOT. The quantum yield effect is ~1.7 at 2000е, down to 0 at 3396.8 (Hill, 2003). To provide worst-case calsystem countrates, the lowest UV QE chips in each build were used: #178 for build 1 and #50 for the flight backup build 2. Figure 4: Flight detector QE curves; dashed lines are the QE curves from DCL which include a correction for quantum yield while the solid line has the correction removed and was used in the SYNPHOT estimates.
QE CCD178 1 1 QE CCD50

.75

.75

throughput

throughput

.5

.5

.25

.25

0 2000

2500

3000 wavelength

3500

4000

0 2000

2500

3000 wavelength

3500

4000

Countrate Predictions The resulting SYNPHOT countrate predictions for the five new lamps are shown in Figure 5 and tabulated in Appendix A. For completeness, the countrate predictions for 31298 are included as well, but note that due to a variety of factors, including a broken weld, it is no longer a flight candidate. Of all the lamps, both old and new, only the nonflight candidate failed the CEI specification of 16.7 e-/s/pix and then, only in two narrow band filters. The low countrates for this particular lamp may have been due to a conservative extrapolation of the spectrum out to 4000е; the scans of the new, more rugged lamps have data out to that wavelength and no extrapolations of those spectra were necessary. For flight build 1, which includes CCD178, all lamps exceed the CEI requirement by at least a factor of two. For the backup flight build, which includes CCD50, all of the new lamps easily meet the CEI specification and all but 1 lamp meet twice the CEI specification. Lamp 31774, in the F280N filter, falls just slightly below twice the requirement.

7


Instrument Science Report WFC3 2004-07

Figure 5: SYNPHOT countrate predictions for WFC3 UVIS calibration system as a function of filter number. Symbols cross, circle, plus, diamond, and asterisk represent the five new rugged lamps 31774, 31775, 31776, 31886, and 31887; box symbols represent the low flux lamp 31298 from the first batch of lamps (not a flight candidate). Solid line is the CEI specification; dashed line is twice the level of the CEI specification.
104 D2 lamp countrates, CCD178 10
4

D2 lamp countrates, CCD50

103

10

3

e-/sec/pix

e-/sec/pix

102

10

2

101

10

1

0

250

300 filter

350

0

250

300 filter

350

Conclusions
The D2 calibration system component throughput estimates from Ball have been migrated over to SYNPHOT, along with the detector QE curves from the GSFC DCL and the lamp spectra from the GSFC optics lab. Calsystem countrate predictions using the new SYNPHOT tables show that all of the five rugged lamps meet the CEI specification of 16.7 e-/s/ pix in filters blueward and including F395N, in either flight build1 or build2. Furthermore, for build1, all lamps exceed at least twice the CEI specification; for build2, all but one lamp (31774) exceed the CEI requirement by a factor of two in all applicable filters.

Acknowledgements
We thank D.Ebbets for helpful discussions of D2 lamps and the WFC3 calsystem.

8


Instrument Science Report WFC3 2004-07

References
Baggett, S., and Quijada, M., "Lifetime Test of Deuterium Lamp for the WFC3 Calibration Subsystem," ISR 2003-11, Nov 2003. Boucarut, R., Garcia, K., and Kutina, R., "Lifetime Test Plan for the HST Wide Field Camera 3 Calibration Deuterium Lamp," April 2003, GSFC internal document. Hill, R., priv.comm, Nov 2003. Hubble Space Telescope Wide Field 3 Contract End Item Specification, May 2000, HST Library #TM-0277550, section 4.11.4, pg 4-46. Landsman, W., tables of DCL UVIS detector QE curves, priv.comm, June 2002. Lupie, O. and Boucarut, R.,"WFC3 UVIS Filters: Measured Throughput and Comparison to Specifications," WFC3 ISR 2003-02, Feb 2003. Sullivan, J., "UVIS Calibration Subsystem Predicted Flux at UVIS Detector,", Ball SER OPT-60, July 2001.

9


Instrument Science Report WFC3 2004-07

Appendix A.
Table 1. Estimated countrates, in e-/s/pix, for D2 the lowest and highest flux flight candidate lamps 31298 countrates are listed as well though it is no noise calculations, countrates can be corrected for the factor in the last column.
lamp 31298 filter ccd 178 f218w f225w f232n f243n f275w f280n f300x f336w f343n f373n f378n f387n f390m f390w f395n 1968.0 2353.7 73.8 96.8 1868.0 66.1 3376.7 1186.4 520.1 43.8 74.1 13.6 84.2 632.4 15.8 ccd 50 1429.0 1740.5 56.4 76.6 1451.9 50.8 2617.9 890.3 395.1 34.1 58.1 11.0 67.9 489.9 13.2 ccd 178 1235.6 1490.8 48.9 63.9 1222.7 40.3 2258.8 1030.0 511.0 78.2 163.2 46.0 343.8 1701.9 137.9 ccd 50 899.8 1106.0 37.4 50.5 950.0 31.0 1750.2 775.9 388.9 60.9 128.3 37.4 281.6 1399.9 115.9 ccd 178 2132.3 2549.6 80.0 104.7 1999.3 70.0 3658.4 1475.8 702.7 99.9 207.0 57.8 431.4 2168.6 172.0 ccd 50 1548.1 1884.9 61.1 82.7 1553.9 53.9 2835.9 1110.3 534.4 77.7 162.7 47.0 353.3 1781.1 144.6 lamp 31774

calsystem. Listed are the countrates for 31774 and 31887; for comparison, longer flight candidate. For signal to the quantum yield effect by dividing by
lamp 31887 quantum yield factor 1.54 1.50 1.46 1.40 1.27 1.21 1.28 1.03 1.01 1. 1. 1. 1. 1. 1.

10