Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.stsci.edu/hst/wfc3/documents/ISRs/WFC3-2009-17.pdf Äàòà èçìåíåíèÿ: Thu Nov 19 20:30:45 2009 Äàòà èíäåêñèðîâàíèÿ: Tue Nov 24 12:32:37 2009 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: arp 220

ST-ECF Instrument Science Report WFC3-2009-17

WFC3 SMOV proposal 11552: Calibration of the G141 grism
H. Kuntschner, H. Bushouse, M. KÝmmel, J. R. Walsh October 16, 2009 A

BSTR A C T

Based on the SMOV observations, the performance of the WFC3 near-IR grism G141 has been assessed. The location of the +1st spectral order relative to the source positions in the direct imaging filters F140W and F160W are determined, field dependent trace and wavelength solutions are derived, and the absolute throughput is quantified. Similar information is derived for higher spectral orders and the -1st order. The trace and wavelength solutions were found to be linear functions, varying smoothly over the field of view. Generally, there is good agreement (<= 1 pixel) with the calibration solutions derived during the ground calibrations in TV2 and TV3. For the +1st order aperture corrections are provided as a function of wavelength. The cross-dispersion PSF is evaluated and agrees with the design expectations as well as ground calibration measurements. All calibrations in this ISR are given with respect to source positions derived from F140W images. Source positions obtained from F160W images show only a small systematic offset relative to F140W images of X=0.057±0.003 and Y=-0.045± 0.004. The +1st order of the G141 grism shows a mean dispersion of 46.5 å/pixel with a resolving power of R~120 at 13000 å. The throughput of HST + WFC3 G141 peaks at 14350 å with 47.8% and is above 10% between 10800 and 16900 å.

T h e S p a c e T e le s c o p e E u ro p e a n C o o rd in a tin g F a c ility . A ll R ig h ts R e s e rv e d .


ST-ECF Technical Instrument Report WFC3-2009-17

1. Introduction
The Wide Field Camera 3 (WFC3) is fitted with three grisms for slitless spectroscopy. In the UVIS channel there is one grism, G280, for the near-UV to visible range (200 - 400nm). The NIR channel has two grisms, G102 and G141, for the shorter (800 - 1150nm) and longer NIR wavelengths (1100-1700nm), respectively. The results from ground calibration efforts, including trace, wavelength, flat-field, throughput and aperture correction calibrations are summarized in Kuntschner et al. (2008a,b) and Kuntschner et al. (2009). This report presents the first in-orbit calibrations for the G141 grism derived from SMOV and selected Cycle 17 calibration observations.

2. Observations
The calibrations presented in this report are based on the SMOV program 11552 (PI: Bushouse), which comprises observations of the HST primary flux standard star GD153, and the wavelength calibrator planetary nebula HB12 (= PN G111.802.8) at several positions over the field-of-view (hereafter FoV). Furthermore, we made use of the Cycle 17 calibration program 11937 (PI: Bushouse), which comprises observations of the wavelength calibrator PN Vy2-2 over 9 different field positions. A summary of the pointing locations for all observations is given in Figure 1.
PROP11552 - GD153
1000 800 600 400 200 0 0 1000 800 600 400 200 0 0

PROP11552 - HB12
1000 800 600 400 200 0 0

PROP11937 - Vy2-2

Y [pix]

Y [pix]

200

400 600 X [pix]

800

1000

200

400 600 X [pix]

800

1000

Y [pix]

200

400 600 X [pix]

800

1000

Figure 1: Distribution of target positions as seen on the F140W images for the SMOV proposal 11552 (GD153 and HB12) as well as Cycle 17 calibration proposal 11937 (Vy2-2). For each position one or more G141 grism exposures were acquired. A detailed listing of all exposures used in this report is given in the Appendix. 2


ST-ECF Technical Instrument Report WFC3-2009-17

3. Analysis
In this section we describe the analysis of the calibration data obtained for the G141 grism in proposals 11552 and 11937 yielding calibrations for the trace, the wavelength solution and the absolute throughput. Furthermore, we investigate offsets between the filters F140W and F160W. If not noted otherwise, analysis was carried out on fully processed images (i.e. "flt" files). The files were produced with the calwf3 pipeline, which carries out the basic data reductions including the following tasks: · · · · · · · · · · Flag known bad pixels in the data quality (DQ) array Subtract bias drifts determined from the reference pixels Identify pixels in the initial ("zeroth") read that contain detectable source signal Subtract the zeroth read Calculate a noise model for each pixel and record in the error (ERR) array Subtract dark image Correct for photometric non-linearity Convert the data to signal rates (in DN per second) Perform "up-the-ramp" fitting and cosmic-ray rejection Perform flat-fielding and apply gain calibration (to electrons per second)

The grism images are not flat-fielded with the normal procedure since for grism spectroscopy a given pixel can receive any wavelength for which the combination of camera and grism shows throughput. In our analysis we employ a special flatfield procedure where for each extracted pixel a wavelength dependent flat-field is applied; therefore calwf3 only applies the gain calibration to grism images. The grism flat-field information was determined during ground calibration and is stored in a flat-field cube (see Kuntschner et al 2008b). In Figure 2 we show a typical direct image ­ G141 grism image pair by using observations of the flux standard star GD153.

3


ST-ECF Technical Instrument Report WFC3-2009-17

Figure 2: Example of a typical direct image ­ grism image pair for the G141 grism. The top panel shows a F140W image (iab904meq_flt.fits) of the flux standard GD153 and the bottom panel shows a G141 grism exposure (iab904mgq_flt.fits) for the same region on the sky. The point like source on the left-hand side of the grism image is the zeroth order spectrum; the +1st, +2nd and +3rd orders follow to the right. Faint traces of +1st orders from other sources in the field of view are visible in the grism image. For each image we show the full extent of the detector in the x-direction and about 200 pixels in the y-direction.

3 .1 .

Trace calibration

In order to establish a good in-orbit trace calibration, we utilize the G141 grism exposures of the PN HB12. For the trace calibration, we are not interested in the primary target itself but in the other point sources around the primary target. This field (galactic latitude: -02.8498 degrees) comprises a relatively good compromise between field coverage + object density and thus avoiding too much spectral overlap between different sources in the FoV (see Figure 2). The ST-ECF aXe software package for the reduction of slitless spectroscopy data treats the spectral traces and wavelength solutions defined with respect to the position of the source in the direct image. The centroids of all sources above a given threshold in the F140W images (Xref, Yref) were determined with SExtractor (Bertin &Arnouts, 1996). These positions were assumed not to change between the observations of each direct image grism pair. The spectra of all sources were traced as a function of X=X-Xref in the detector X-direction by measuring the centroids of 7-10 pixel wide bins, using custom-written IDL programs. We found the traces of all orders (excluding the zeroth order) to be well fit by straight lines with standard deviations of 0.1 pixels. The trace definitions are of the form (Y-Yref)=DYDX_0 + DYDX_1*X, where DYDX_0 and DYDX_1 are field dependent and given in the usual format used by the ST-ECF aXe reduction package, e.g., DYDX_1 = a0

4


ST-ECF Technical Instrument Report WFC3-2009-17 +a1*Xref + a2*Yref + ... (see also the aXe manual1 for more details).

Figure 2: Image of the G141 grism exposure (iab908s4q_flt.fits). The primary target HB12 is located in the centre. Most of the other prominent spectra are first orders of different (stellar) objects in the FoV used for the trace calibration. The final, field dependent trace solution for the +1st order of the G141 grism was derived from 162 different traces covering semi-uniformly the FOV. A graphical representation of the measured offsets and slopes as a function of Xref and Yref positions is given in Figure 3. The offset shows a marked trend with Yref position while the slope shows a more complex dependence on Xref and Yref. The offset is well represented with a field dependent fit using only a linear function (rms = 0.08 pixel), while the slope required a fit with quadratic terms (rms = 0.0006). The inorbit calibrations show reasonably good agreement with the ground calibrations obtained in TV2 and TV3 with differences of ~0.3 pixel in the offset and good
1

http://www.stecf.org/software/slitless_software/axe/ 5


ST-ECF Technical Instrument Report WFC3-2009-17 agreement for the slope. The results of our field dependent fits are given in Table 1 and are also shown in Figure 3. For the -1st, +2nd and +3rd orders a similar procedure was carried out, however, due to the reduced field coverage for these orders the field dependent fits were more restricted and the trace prediction accuracy is reduced. Results are summarized in Table 1. For the 0th order we adopt the trace determined in the ground calibrations. The +4th order spectra are visible for sources located at Xref 270, however, this order is not analysed as part of this ISR and will be determined at a later stage.

6


ST-ECF Technical Instrument Report WFC3-2009-17

1000 800

Yref [pix]

600 400 200

200

400 600 Xref [pix]

800

1000

2.0 1.5 1.0 0.5 0.0 200 400 600 Xref [pix] 800 1000

2.0 1.5 1.0 0.5 0.0 200 400 600 Yref [pix] 800 1000

Offset [pix]

0.016 0.014 0.012

Slope

0.008 0.006 0.004 0.002 200 400 600 Xref [pix] 800 1000

Slope

0.010

Offset [pix]

0.016 0.014 0.012 0.010 0.008 0.006 0.004 0.002 200 400 600 Yref [pix] 800 1000

Figure 3: The trace fits for the G141 1st order spectra are shown as a function of Xref and Yref position (diamond symbols). A total of 162 different source positions are used covering well the field-of-view (top left panel). The final field-dependent trace solution fit (see also Table 1) is shown as the grid of black dots.

7


ST-ECF Technical Instrument Report WFC3-2009-17 Table 1: Field dependent trace descriptions for G141. Where available errors are given below the values
Term DYDX_A_0 error DYDX_A_1 error DYDX_B_0 DYDX_B_1 XOFF_B YOFF_B a0 1.96882E+00 8.09111E-02 1.04275E-02 5.94731E-04 6.58628E-01 0.0 -1.92240E+02 0.0 a1(X) +1 9.09159E-05 3.57223E-06 -7.96978E-06 4.34517E-07 1.58398E-04 -2.31440E-03 +2 DYDX_C_0 error DYDX_C_1 error DYDX_D_0 error DYDX_D_1 error DYDX_E_0 error DYDX_E_1 error 2.26624E+00 4.83732E-01 5.39806E-03 1.57094E-03 3.05905E+00 5.02518E-01 5.65708E-03 1.22990E-03 4.19045E+00 1.11680E+00 1.39687E-02 2.44552E-03 -1.78192E-04 1.92052E-05 2.78779E-06 3.53707E-07 0.0 0.0 0.0 0.0 nd st

a2(X) order -1.93260E-03 3.12042E-06 -2.49607E-06 3.57986E-07 0
th

a3(X^2) 1.45963E-09 3.87141E-10 -

a4(X*Y) 1.39757E-08 3.29421E-10 -

a5(Y^2) 4.84940E-10 3.08712E-10 -

order

-1.62512E-03 1.11089E-02 order

-2.71785E-03 1.34415E-05 4.84354E-06 2.47603E-07 +3
rd

-

-

-

order

-5.51206E-03 2.65696E-05 6.87491E-06 3.43741E-07 -1
st

order

-8.53483E-03 1.30599E-04 -1.77629E-05 3.05845E-06

Notes: The 0th order trace solutions are taken from the ground calibrations.

3 .2 .

P o s itio n d e p e n d e n c e o n filte r

The trace and wavelength calibrations reported in this ISR are based on the source positions measured from the standard filter F140W. An alternative filter for the G141 grism would be F160W. If the F160W filter would introduce a shift in source pixel positions with respect to the F140W filter, the trace and wavelength calibrations are directly affected. For the WFC3 G102 grism and its associated filters F098M and F105W, we detected shifts of up to 0.3 pixel (see Kuntschner et al. 2009). We tested for these potential shifts by using the imaging of Proposal 11937 (PN Vy2-2). For each pair of F140W and F160W images taken at the same 8


ST-ECF Technical Instrument Report WFC3-2009-17 POSTARGS, we ran Sextractor and matched the source positions. In total there are 9 image pairs, which yielded 2492 matched object positions after some robust cleaning had been applied. The resulting mean offsets were derived with a biweight mean estimator (IDL biweight_mean.pro) and are given in Table 2. The analysis showed little evidence for a field-position dependence of the shifts as can be seen from Figure 4. Since the shifts are smaller than 1/10 of a pixels it will generally not be required to apply any offsets. Table 2: Source position offsets between F140W and F160W X offset [pixel] Y offset [pixel] F140W ­ F160W 0.057±0.003 -0.045±0.004

9


ST-ECF Technical Instrument Report WFC3-2009-17

F140W - F160W
1.5 1.0

Delta Y [pix]

0.5 0.0 -0.5 -1.0 -1.5 -1.5

-1.0

-0.5 0.0 0.5 Delta X [pix]

1.0

1.5

X offset: 0.057 +/- 0.003
300 250 300 200 150 100 100 50 0 -1.0 0 -0.5 0.0 0.5 Delta X [pix] 1.0 200 400

Y offset: -0.045 +/- 0.004

N

N

-1.0

-0.5 0.0 0.5 Delta Y [pix]

1.0

1.5 1.0

1.5 1.0

Delta X, Y [pix]

0.5 0.0 -0.5 -1.0 -1.5 200 400 600 X [pix] 800 1000

Delta X, Y [pix]

0.5 0.0 -0.5 -1.0 -1.5 200 400 600 Y [pix] 800 1000

Figure 4: Source position shifts between the filters F140W and F160W. A total of 2492 sources were matched in the direct imaging of Proposal 11937. Individual offsets and the mean are shown in the top panel. Histograms versus x and y-axis are shown in the middle panels whereas potential offset dependence with X and Y position is explored in the bottom panels. Solid lines show a robust linear fit, while the dashed lines indicate the mean shift in x (blue) and y (red).

10


ST-ECF Technical Instrument Report WFC3-2009-17

3 .3 .

Wavelength solutions

In order to establish the in-orbit wavelength calibration, the PN HB12 was observed during SMOV (Proposal 11552). Despite the extended nature of the object, the central pixels of the source saturate even in the first read of the direct image, which makes it difficult to determine an accurate source position. The PN Vy2-2 was observed as a wavelength calibrator as part of the Cycle 17 calibration program (Proposal 11937). This PN is somewhat fainter and accurate positions of the source could be derived from the direct imaging. However, also here the first read is slightly saturated and the final positions were derived from the zeroth read using the "*_raw.fits" files. For wavelength calibration one would ideally want point sources with clearly identifiable emission lines over the wavelength range and at the resolution of the grism. While bright, nearby PNs do show a reasonably good distribution of emission lines, they are slightly extended at the spatial resolution of the WFC3 IR camera. Since the PNs Vy2-2 and HB12 have already been used for NICMOS grism calibration (see Pirzkal et al. 2009) and they offer a reasonable balance between brightness, number of emission lines and wavelength coverage, we opted for using those targets in the first wavelength calibrations of the WFC3 IR grisms. Future calibration programs may consider the use of e.g., WR stars or M-dwarfs. Using the trace calibration described in Section 3.1 we extracted the spectra of Vy2-2 at all nine field positions with the help of the aXe software. In this extraction a special configuration file was used which provided spectra with counts versus pixel number along the x-axis. For the G141 grism (+1st order) a total of 7 emission features could be identified and used to establish the wavelength solution (see Table 3). Prior to the fitting, the radial velocity of -71 km/s for Vy2-2 and an air to vacuum conversion was applied to the tabulated wavelength of the emission lines (see Table 3). The ground calibration efforts suggested that a linear wavelength solution is a good representation of the true wavelength solution for this grism (Kuntschner et al. 2008a, b). Therefore we established for each spatial position a linear wavelength solution (see Figure 5) where fits give a typical rms scatter of 0.2-0.3 pixel (see Figure 5). The resulting wavelength solution shows typical errors of 9 å for the zeropoint and 0.07 å/pixel for the dispersion (~46.5 å/pixel). Sometimes the emission lines are affected by e.g., cosmic ray hits, hot/dead pixels etc. and thus good emission line centroids could not be established. Therefore we used a 3-sigma rejection iteration for each linear fit to remove outliers.

11


ST-ECF Technical Instrument Report WFC3-2009-17 Table 3: Rest air wavelength of emission features fitted to Vy2-2 Name He I He I He I H B13 H B12 H B11 H B10 Wavelength [å] 10830.3 11969.1 12784.9
2 2

16109.3 16407.2 16806.5 17362.1

A ground based spectrum of Vy2-2 (Hora et al. 1999) suggests that at the spectral resolution of the G141 grism the effective peak wavelength requires adjustments of -10.4 and +34.6 å, for the He I emission lines at 11969.1 and 12784.9 å. Before determining the dispersion solution these adjustments were applied. 12

2


ST-ECF Technical Instrument Report WFC3-2009-17

Figure 5: Linear fit for the wavelength solution of the G141 grism. The example shows a fit for one pointing (Xref=123.56, Yref=919.38). The top plot shows the fit whereas the bottom panel shows the deviations from the linear fit; the dashed lines indicate ±0.25 pixel deviation from the fit. A similar procedure was used to establish wavelength solutions for the +2nd, +3rd and -1st orders of the G141 grism. Due to the finite detector size and therefore reduced wavelength coverage for the higher orders a reduced number of emission lines was used in these fits. The average dispersion per pixel in the higher orders is higher by the approximate ratio of the orders. It is 23.6, 15.7 and -48.7 å/pix for the +2nd, +3rd, and -1st order, respectively. Following the methodology used for the trace calibrations we established field 13


ST-ECF Technical Instrument Report WFC3-2009-17 dependent (2-dimensional) fits for the wavelength solution (see Figure 6 and Table 5). The fit shows a marked trend of the wavelength zeropoint with Xref position. A trend with Yref position was not significant and thus set to 0.0. The zeropoint measurements show an rms scatter of about 8 å with respect to the fit. For the dispersion we find a marked trend with Yref position and a much weaker, quadratic trend with Xref positions. The dispersion varies approximately from 45.0 to 47.5 å/pixel across the full field of view. The measurements show a rms scatter of 0.06 å/pixel with respect to the fit. Compared to the ground calibration measurements in TV3 we find a significant change as detailed in Table 4 for a source position at Xref=340 and Yref=550. However, over most of the wavelength range of the G141 grism the differences are less than one pixel in dispersion direction. The differences are possibly caused by gravity release. Table 4: Wavelength solution comparison between TV3 and SMOV for a source at Xref=340 and Yref=550. Zeropoint [å] Dispersion [å/pixel] TV3 SMOV 8953.0 8986.1 46.792 46.517

14


ST-ECF Technical Instrument Report WFC3-2009-17

Figure 6: The wavelength solution for the G141 1st order spectra is shown as a function of Xref and Yref position (diamond symbols). The final field-dependent wavelength solution (see also Table 5) is shown as the grid of black dots.

15


ST-ECF Technical Instrument Report WFC3-2009-17

Table 5: Field dependent wavelength solution for G141. Where available errors are given below the values
Term DLDP_A_0 error DLDP_A_1 error DLDP_B_0 DLDP_B_1 DLDP_C_0 error DLDP_C_1 error DLDP_D_0 error DLDP_D_1 error DLDP_E_0 error DLDP_E_1 error a0 8.95431E+03 8.14876E+00 4.51423E+01 6.26774E-02 1.28397E+04 2.32405E+03 4.48233E+03 2.79017E+00 2.28943E+01 9.95420E-03 3.00187E+03 6.60071E+00 1.52552E+01 1.23680E-02 -9.37152E+03 2.56600E+00 -4.68834E+01 1.22200E-02 a1(X) +1 9.35925E-02 1.09748E-02 3.17239E-04 3.98039E-04 +2 7.61896E-02 2.26239E-02 -1.81224E-04 5.51487E-05 1.04205E-01 8.33292E-02 -2.08555E-04 1.41001E-04 0.0 0.0 nd st

a2(Y) order 0.0 2.17055E-03 2.31850E-04 0 order order -2.70153E-03 1.01228E-02 1.40766E-03 2.95185E-05 +3
rd th

a3(X^2) -7.42504E-07 4.45730E-07 -

a4(X*Y) 3.48639E-07 3.20519E-07 3.09213E-07 2.16386E-07 -

order

-1.18134E-03 3.99081E-02 9.55645E-04 6.74659E-05 -1
st

order

1.71214E-01 1.00445E-01 -2.48997E-03 2.29997E-04

In Figure 7 we show the individual emission lines of Vy2-2 after the flux calibration (see Section 3.3) and our field-dependent trace and wavelength solutions have been applied to the data. Overall, the field dependent calibrations provide a satisfactory result. A further test of our in-orbit wavelength calibrations was performed by making use of the flux standard star GD153 (see also Section 3.3). The G141 grism covers the wavelength of the prominent Hydrogen Paschen series lines H P5 (12818.1 nm), which is clearly seen in the white dwarf GD153. Using the direct pixel fitting code ppxf (Cappellari & Emsellem 2004) we fit a model spectrum of GD153 (CALSPEC library; gd153_mod_005.fits) to the observed spectrum. The ppxf IDL routine involves the logarithmic rebinning of the spectrum and we chose a value of 1085 km/s/pixel. We measure a relative velocity between our observations and the model spectrum of 1±70 km/s, which corresponds to 0.0±0.06 velocity pixels. The above numbers indicate an accuracy of the wavelength zeropoint of <0.25 pixel. From the fitting procedure, we can also provide a rough estimate of the spectral 16


ST-ECF Technical Instrument Report WFC3-2009-17 resolution 13000 nm from the agreement at the wavelength of the two Paschen lines and obtain R=116±11 at or a FWHM of 112±11 å. This is lower than what would be inferred cross-dispersion profile (see Figure 13), however, it is broadly in with the 2 * pixel FWHM (2*46.5=93 å).

17


ST-ECF Technical Instrument Report WFC3-2009-17

Figure 7: Selected emission lines in the spectra of Vy2-2. The aXe software was used to extract the spectra utilizing our field-dependent trace and wavelength solutions established in Section 3.1 and 3.2. Shown in color are the 12 individual spectra covering 9 spatial pointings (see also Figure 1). The 4 central spectra were combined and are shown in black. The color-coding is given in the top left panel. For each panel and emission line we give the identification of the lines in the title and the central wavelength and ± 0.5 pixel offsets are given as vertical red lines. The central pixels of the He I line (top left panel) in one spectrum (green, ibbu02tzq_flt.fits) are affected by bad detector pixels and therefore the line appears to have a dip at the center.

18


ST-ECF Technical Instrument Report WFC3-2009-17

3 .4 .

T h r o u g h p u t m e a s u r e m e n ts

Using the trace and wavelength solutions derived in the previous sections, the spectra of the flux standard star GD153 observed in Proposal 11552 were extracted. Note, that we make use of our own flat-field calibrations for the G141 grism, which significantly improve the quality of the spectra. There were four, slightly dithered, G141 exposures taken near the central position of the FoV. A further two exposures cover the top left and bottom right corners of the FoV (see Figure 1 and Table A1). For the flux calibration we combine the +1st order of the four central exposures into a single spectrum with the help of the spectral drizzle option in the aXe software (aXedrizzle). This spectrum is used to establish the flux calibration for the G141 grism. Firstly, the spectrum was converted to units of [e- / å / sec] and divided by a smoothed version of the model spectrum taken from the HST CALSPEC library (gd153_mod_005.fits). Smoothing in a sliding 7 pixel window provides the final sensitivity function. The error is evaluated by taking the standard deviation of the four individual sensitivity curves from the central pointings at each wavelength bin and imposing a minimum error of 1%. Using this new sensitivity curve, Figure 8 shows (from the top) the aXe-extracted spectra in electrons per second per pixel, in flux units and the ratio of the observed flux over the model spectrum versus wavelength. In the top plot of Figure 8, significant variations between spectra can be seen since the spectral dispersion (å/pixel) changes as a function of position within the field-of-view. The overall agreement between the fully calibrated standard star spectra (middle plot) is good for the wavelength range 1050 to 1700 nm. We note, that the two spectra which show up to 4-5% less flux are the ones located at the upper left and lower right positions of the FoV. This is clear evidence for a large-scale illumination variation which still requires correction in the grism flat-field cubes. Observations of the flux standard star GD71 in the Cycle 17 calibration program cover 9 pointings over the FoV and will enable such corrections. Sensitivity curves were also derived for the 0th, +2nd, +3rd and -1st orders. For the +2nd and +3rd orders we combined the 4 central exposures, whereas for the -1st order only one exposure was available. In Figure 9 we summarize the results for all sensitivity curves. Compared to the TV3 results we find up to 10% more throughput in the +1st order. The TV3 results are shown as green lines in Figure 9. The peak throughput of HST + WFC3 G141 is 47.8% at 1435 nm while it is >10% over the wavelength range 1080 to 1690 nm.

19


ST-ECF Technical Instrument Report WFC3-2009-17

WFC3 G141 grism - SMOV
500
Counts [e-/s] Flux [erg s-1 cm-2 A-1] Observed / standard

400 300 200 100 10
-14

10

-15

10-16 1.10 1.05 1.00 0.95 0.90 1000 1200 1400 Wavelength [nm] 1600

Figure 8: The aXe extracted 1-dimensional spectra of the standard star GD153 in units of electrons per second per pixel, fully calibrated flux units and the ratio of the observed flux over the model spectrum versus wavelength (from top to bottom). Observations obtained at different field positions are colored, and the combined (aXedrizzle) central spectrum is shown in black. In the individually extracted spectra there are a number of prominent wavelength steps, which show too low flux. These points are caused by pixels, which are flagged in the original flt files. The aXedrizzle spectrum corrects for these effects and thus does not show any deviant wavelength steps.

20


ST-ECF Technical Instrument Report WFC3-2009-17

10

17

WFC3 G141 grism - SMOV
+1st order +2nd order +3rd order -1st order

Sensitivity [e- s-1 A-1 per erg cm-2 s-1 A-1]

10

16

10

15

10

14

1013 1000

1200

1400 Wavelength [nm]

1600

1800

Figure 9: Sensitivities of the G141 grism for the +1st, +2nd, +3rd and -1st orders. The results were derived from the observations of the flux standard GD153 during SMOV (Proposal 11552). For comparison we show in green the measurements obtained during the TV3 ground calibrations.

21


ST-ECF Technical Instrument Report WFC3-2009-17

WFC3 G141 grism - SMOV
1.000

0.100

Total transmission

+1st order +2nd order +3rd order -1st order 0th order OTA

0.010

0.001 1000 1200 1400 1600 Wavelength [nm] 1800

Figure 10: Total throughput estimates for the +1st, +2nd, +3rd and -1st orders of the G141 grism, derived from the sensitivity curves of this report. The 0th order is taken from the TV3 ground calibrations and only shown for comparison. In Figure 11 we compare the total throughput estimates from the WFC3 IR grisms with "synphot" predictions. The synphot files are based and calibrated to WFC3 imaging from the SMOV observations. The plot shows theoretical values, which would be achieved with WFC3 without any filter or grism inserted. The analysis of the G141 grism is described in detail in Kuntschner et al. (2008a, b). The overall agreement to better than ±5% is good and shows that a consistent first order flux calibration of the IR channel was achieved.

22


ST-ECF Technical Instrument Report WFC3-2009-17

WFC3 IR Throughput "Clear"

0.8 0.6
Synphot G102 G141

0.4

0.2

0.0 1.15 1.10
Synphot / Grism

1.05 1.00 0.95 0.90 0.85 800 1000 1200 1400 Wavelength [nm] 1600 1800

Figure 11: Comparison of total throughput "synphot" predictions with SMOV observations of the WFC3 IR grisms. The upper plot shows theoretical values, which would be achieved with WFC3 without any filter or grism inserted. The lower plot shows the ratio of "synphot" predictions over observations.

3 .5 . A p e r tu r e c o r r e c tio n s a n d c r o s s d is p e r s io n PSF
Using the aXe software (aXedrizzle option) we derived a rectified 2-dimensional image for the central exposures of the flux standard star GD153. The resulting image (called "stamp" image) is wavelength calibrated and image distortions, as recorded in our trace calibrations, are removed. From this image, covering a 220 pixel wide aperture around GD153, various sub-apertures (here used in the sense of 23


ST-ECF Technical Instrument Report WFC3-2009-17 a diameter) are extracted with the help of IDL scripts and compared to the flu determined in the largest aperture. This procedure is carried out in five wavelength bins. The resulting aperture correction values are given in Tables 6 and shown in Figure 12. About 95% of the flux is concentrated in an aperture of 15 pixels diameter or 1.924 arcsec. The aperture corrections derived from in orbit data compare well with TV3 results with a maximum difference of 4%. Table 6: Aperture corrections for WFC3 G141
D ia m e te r [arcsec] 0 .1 2 8 0 .3 8 5 0 .6 4 1 0 .8 9 8 1 .1 5 4 1 .4 1 1 1 .6 6 7 1 .9 2 4 3 .7 1 9 7 .5 6 7 1 2 .9 5 4 2 5 .7 7 9 11300 å 0 0 0 0 0 0 0 0 0 0 0 1 . . . . . . . . . . . . 4 8 8 9 9 9 9 9 9 9 9 0 4 1 9 1 3 5 6 7 9 9 9 0 5 1 2 8 9 3 3 1 1 7 9 0 12300 å 0 0 0 0 0 0 0 0 0 0 0 1 . . . . . . . . . . . . 4 7 8 8 9 9 9 9 9 9 9 0 4 9 7 9 2 3 4 5 8 9 9 0 2 1 5 9 2 7 7 6 2 2 5 0 13300 å 0 0 0 0 0 0 0 0 0 0 0 1 . . . . . . . . . . . . 4 7 8 8 9 9 9 9 9 9 9 0 0 7 7 9 2 3 4 5 8 9 9 0 2 0 1 9 0 7 8 7 6 6 9 0 14300 å 0 0 0 0 0 0 0 0 0 0 1 1 . . . . . . . . . . . . 3 7 8 8 9 9 9 9 9 9 0 0 4 5 6 9 1 3 4 5 8 9 0 0 7 0 6 7 6 4 7 6 8 8 1 0 15300 å 0 0 0 0 0 0 0 0 0 1 1 1 . . . . . . . . . . . . 3 7 8 8 9 9 9 9 9 0 0 0 4 4 6 9 1 3 4 5 9 0 0 0 5 1 1 5 4 2 6 6 1 1 2 0 16300 å 0 0 0 0 0 0 0 0 0 0 0 1 . . . . . . . . . . . . 3 7 8 8 9 9 9 9 9 9 9 0 7 2 5 8 0 2 3 4 7 9 9 0 5 4 0 9 4 2 6 6 7 1 4 0

24


ST-ECF Technical Instrument Report WFC3-2009-17

WFC3/IR G141 SMOV Aper corr center
1.0

Normalized Flux

0.8 0.6 0.4 0.2 0.0
Lambda Lambda Lambda Lambda Lambda Lambda 11300 12300 13300 14300 15300 16300

0.1 1.0 10.0 Extraction aperture [arcsec]
Figure 12: Aperture corrections for WFC3 G141 as a function of aperture (=diameter) and wavelength. For comparison TV3 data at 1330 nm is shown as thick grey line. We also investigated the cross dispersion PSF of the central pointing for the flux standard star GD153. Ground calibration results showed that the G141 grism produces, for a Gaussian approximation, a cross dispersion width (FWHM) of XX to XX pixel (see Kuntschner et al. 2009). We show in Figure 13 the SMOV measurements, which indicate a slightly smaller width in comparison with TV3 results. Sampling effects cause the wave like structure of the measurements and we consider the minima as good representations of the true cross dispersion PSF. The measurements are also in agreement with the expectations from the direct imaging also shown in Figure 13. In conclusion, the G141 grism is in focus and produces spectra of point sources in agreement with design expectations.

25


ST-ECF Technical Instrument Report WFC3-2009-17

2.5

cross disp. PSF FWHM (pix)

2.0

1.5

1.0 1000

1200 1400 1600 Wavelength (nm)

1800

Figure 13: Cross-dispersion profile width as measured in FWHM (pixel) for a Gaussian approximation. Diamonds represent the SMOV data for one of the central pointings for grism G141. The small stars show estimates derived from white light images obtained during the TV3 ground calibrations (see Kuntschner et al. 2009). The triangles connected with a solid line show for comparison measurements from the TV3 direct light images as tabulated by Hartig (2008b).

4. Conclusions
This ISR presented the calibration of the WFC3 G141 grism based on SMOV data and one Cycle 17 wavelength calibration proposal. We have established in-orbit, source position dependent calibrations of the trace and wavelength solution for the +1st order "science" spectra. Both, the local trace and wavelength solution can be approximated with a linear solution. The mean dispersion of the G141 grism is 26


ST-ECF Technical Instrument Report WFC3-2009-17 46.5 å/pixel varying from 45.0 to 47.5 å/pixel across the field of view. The inorbit solutions agree well (<= 1 pixel) with the calibrations established during ground calibrations in TV2 and TV3. We also provide, field dependent calibrations for the +2nd, +3rd and -1st spectral orders, albeit with a less rigorous analysis. We establish flux calibrations for the +1st, +2nd +3rd and -1st orders. Furthermore, we provide aperture corrections for the +1st order as a function of wavelength. The throughput of HST + WGC3 G141 peaks at 14350 å with 47.8% and is above 10% between 10800 and 16900 å.

A p p e n d ix
The datasets acquired during SMOV Proposal 11552 (PI Bushouse) and the Cycle 17 calibration program of Proposal 11937 (PI Bushouse) are listed in Table A1 and A2, respectively. Note, that throughout this report we use the name of "HB12" for the target PN G111.8-02.8. Table A1. Summary of observations for the calibration of the grism G141 in Proposal 11552
FILE i i i i i i i i i i i i i i i i i i i i a a a a a a a a a a a a a a a a a a a a b b b b b b b b b b b b b b b b b b b b 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 04m eq 0 4 m fq 04m gq 04m kq 0 4 m lq 04m m q 04m oq 04m pq 04m qq a4m sq a 4 m tq a4m uq 0 0 0 0 0 0 0 0 8 8 8 8 8 8 8 8 s s s s s s s s 3q 4q 8q 9q dq eq hq iq TARGNAME GD GD GD GD GD GD GD GD GD GD GD GD PN PN PN PN PN PN PN PN 153 153 153 153 153 153 153 153 153 153 153 153 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 . . . . . . . . 8 8 8 8 8 8 8 8 0 0 0 0 0 0 0 0 2 2 2 2 2 2 2 2 . . . . . . . . 8 8 8 8 8 8 8 8 FILT E R F140W F160W G 141 G 141 G 141 G 141 F140W F160W G 141 F140W F160W G 141 F140W G 141 F140W G 141 F140W G 141 F140W G 141 E X P T IM E [seco n d s] 5 .8 6 5 .8 6 1 0 2 .9 3 1 0 2 .9 3 1 0 2 .9 3 2 0 2 .9 3 5 .8 6 5 .8 6 1 0 2 .9 3 5 .8 6 5 .8 6 1 0 2 .9 3 2 2 2 2 2 2 2 2 . . . . . . . . 9 9 9 9 9 9 9 9 3 3 3 3 3 3 3 3 PO STARG 1 [arcsec] -2 0 .0 -2 0 .0 -2 0 .0 -2 0 .5 -2 1 .0 -2 1 .2 2 3 .0 2 3 .0 2 3 .0 -5 0 .0 -5 0 .0 -5 0 .0 2 2 5 5 2 2 2 2 0 0 0 0 0 0 3 3 . . . . . . . . 0 0 0 0 0 0 0 0 PO STARG 2 [arcsec] 0 .0 0 .0 0 .0 0 .5 1 .0 1 .2 -4 5 .0 -4 5 .0 -4 5 .0 4 5 .0 4 5 .0 4 5 .0 0 0 5 5 5 5 5 5 . . . . . . . . 0 0 0 0 0 0 0 0

-

G G G G G G G G

10 10 10 10

4 4 4 4 4 4

27


ST-ECF Technical Instrument Report WFC3-2009-17
i i i i i i i i i i a a a a a a a a a a b b b b b b b b b b 9 9 9 9 9 9 9 9 9 9 0 0 0 0 0 0 0 0 0 0 8 8 8 8 8 8 8 8 8 8 s lq sm q spq sqq t9 q ta q td q te q th q tiq PN PN PN PN PN PN PN PN PN PN G G G G G G G G G G 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 . . . . . . . . . . 8 8 8 8 8 8 8 8 8 8 0 0 0 0 0 0 0 0 0 0 2 2 2 2 2 2 2 2 2 2 . . . . . . . . . . 8 8 8 8 8 8 8 8 8 8 F140W G 141 F140W G 141 F140W G 141 F140W G 141 F140W G 141 2 102 2 102 2 102 2 102 2 102 . . . . . . . . . . 9 9 9 9 9 9 9 9 9 9 3 3 3 3 3 3 3 3 3 3 -5 0 -5 0 23 23 -5 0 -5 0 -2 0 -2 0 23 23 . . . . . . . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 5 5 5 5 5 . . . . . . . . . . 0 0 0 0 0 0 0 0 0 0

-

4 4 4 4 4 4

Table A1. Summary of observations for the calibration of the grism G141 in Proposal 11937
FILE i i i i i i i i i i i i i i i i i i i i i i i i i i i i b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b b u u u u u u u u u u u u u u u u u u u u u u u u u u u u 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 tw q tx q ty q tz q u0q u1q u3q u4q u5q uw q uxq uyq v0q v1q v2q v3q v4q v5q v7q v8q v9q vaq vbq vcq veq v fq vgq vhq TARGNAME V V V V V V V V V V V V V V V V V V V V V V V V V V V V Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 FILT E R F140W F160W G 141 G 141 G 141 G 141 F140W F160W G 141 F140W F160W G 141 F140W F160W G 141 F140W F160W G 141 F140W F160W G 141 F140W F160W G 141 F140W F160W G 141 F140W E X P T IM E [seco n d s] 5 .8 6 5 .8 6 1 0 2 .9 3 1 0 2 .9 3 1 0 2 .9 3 1 0 2 .9 3 5 .8 6 5 .8 6 1 0 2 .9 3 5 .8 6 5 .8 6 1 0 2 .9 3 5 .8 6 5 .8 6 1 0 2 .9 3 5 .8 6 5 .8 6 1 0 2 .9 3 5 .8 6 5 .8 6 1 0 2 .9 3 5 .8 6 5 .8 6 1 0 2 .9 3 5 .8 6 5 .8 6 1 0 2 .9 3 5 .8 6 PO STARG 1 [arcsec] -2 0 .0 -2 0 .0 -2 0 .0 -2 0 .5 -2 1 .0 -2 1 .2 -5 0 .0 -5 0 .0 -5 0 .0 -2 0 .0 -2 0 .0 -2 0 .0 2 3 .0 2 3 .0 2 3 .0 -5 0 .0 -5 0 .0 -5 0 .0 2 3 .0 2 3 .0 2 3 .0 -5 0 .0 -5 0 .0 -5 0 .0 -2 0 .0 -2 0 .0 -2 0 .0 2 3 .0 PO STARG 2 [arcsec] 0 .0 0 .0 0 .0 0 .5 1 .0 1 .2 4 5 .0 4 5 .0 4 5 .0 4 5 .0 4 5 .0 4 5 .0 4 5 .0 4 5 .0 4 5 .0 0 .0 0 .0 0 .0 0 .0 0 .0 0 .0 -4 5 .0 -4 5 .0 -4 5 .0 -4 5 .0 -4 5 .0 -4 5 .0 -4 5 .0

28


ST-ECF Technical Instrument Report WFC3-2009-17
ib b u 0 2 v iq ib b u 0 2 v jq V Y 2 -2 V Y 2 -2 F160W G 141 5 .8 6 1 0 2 .9 3 2 3 .0 2 3 .0 -4 5 .0 -4 5 .0

R e fe r e n c e s
Bertin, E., Arnouts, S. 1996, AAS, 117, 393 Cappellari M., Emsellem E., 2004, PASP, 116, 138 Hora, J. L., Latter, W.B., Deutsch, L. K., 1999, APJS, 124, 195 Kuntschner, H. Bushouse, H., KÝmmel, M., Walsh, J. R., 2009, WFC3 Instrument Science Report, WFC3-2009-XX: The SMOV calibration of the WFC3 G102 grism Kuntschner, H. Bushouse, H., KÝmmel, M., Walsh, J. R., 2009, WFC3 Instrument Science Report, WFC3-2009-04: The cross dispersion profiles of the WFC3 grisms Kuntschner, H. Bushouse, H., Walsh, J. R., KÝmmel, M., 2008a, WFC3 Instrument Science Report, WFC3-2008-16: The TV3 ground calibrations of the WFC3 NIR grisms Kuntschner, H. Bushouse, H., Walsh, J. R., KÝmmel, M., 2008b, WFC3 Instrument Science Report, WFC3-2008-15: The TV2 ground calibrations of the WFC3 NIR grisms Pirzkal, N., Bohlin, R. Thatte, D., 2009, NICMOS Instrument Science Report, NICMOS-2009-006: NICMOS Grism Wavelength Calibration

29