UNSWIRF: A Tunable Imaging Spectrometer for the Near-Infrared

Stuart D. Ryder , Yin-Sheng Sun , Michael C. B. Ashley , Michael G. Burton , Lori E. Allen , John W. V. Storey, PASA, 15 (2), 228
The html and gzipped postscript versions of this paper are in preprint form.
To access the final published version, download the pdf file
.

Next Section: Results
Title/Abstract Page: UNSWIRF: A Tunable Imaging
Previous Section: Overview of the Instrument
Contents Page: Volume 15, Number 2

 

Observing with UNSWIRF

Wavelength Calibration

The first step in commissioning UNSWIRF was to calibrate the relationship between etalon spacing d and peak transmitted wavelength tex2html_wrap_inline921. This is complicated by a number of factors:

  • the actual order n being passed by the blocking filter is not necessarily known in advance;
  • the relationship between etalon spacing d, and the analog-to-digital units Z (-2048 < Z < +2047) employed by the control computer and CS-100 must first be calibrated;
  • because of the wide tuning range demanded, and the multiple layer coatings involved, the effective depth where reflection occurs within the etalon coatings will change with wavelength, and so too will the apparent physical spacing of the plates. This manifests itself as a change in tex2html_wrap_inline1093 with respect to an assumed value of the physical plate spacing (or equivalently, by a non-integral order number n). Let us define tex2html_wrap_inline1097 such that when the angle of incidence tex2html_wrap_inline1099, and the interplate medium is air at 1 atmosphere pressure, equation 1 becomes
     equation121
    and thus
     equation125
    Thus, measurements of the FSR over the full range in wavelength and for multiple orders allows us to derive the order number n and the wavelength variation of apparent plate separation tex2html_wrap_inline1097.
We now outline a novel calibration procedure that deals with these issues, and replicates the actual observing configuration as closely as possible, without the need for external calibration in the laboratory.

First, the etalon and all blocking filters are withdrawn from the light path, and IRIS configured with the narrow (tex2html_wrap_inline1105) slit and its H+K band échelle grism. Images of the emission-line spectra produced by four separate lamps (Ar, Kr, Hg, and Xe) are then taken in order to wavelength-calibrate the four complete échelle orders (covering the range tex2html_wrap_inline1109m after straightening) produced by this grism. Next, the Fabry-Perot is inserted into the beam, the plate spacing set to Z=0, and a continuum source (such as a quartz-iodine lamp) used to illuminate the system. Since there are no blocking filters in place, all orders passed by the etalon are imaged, resulting in the ``picket-fence'' appearance of Figure 1. Since wavelength as a function of position on the array is already known, the position of each peak, and their separations (i.e., tex2html_wrap_inline1093) can be determined. For improved accuracy, these measurements are repeated with etalon settings tex2html_wrap_inline1115, which causes the etalon orders to shift position on the array, as illustrated in Figure 1.

  figure142
Figure: The ``picket fence'' of orders produced using the H+K échelle and narrow slit of IRIS, a QI continuum lamp, and no blocking filters with UNSWIRF is shown for etalon Z=200 (top left), 400 (top right), 600 (bottom left), and 800 ( bottom right). Note also how the spacing between adjacent orders of the etalon (tex2html_wrap_inline1093) increases with wavelength (eqtn 2). The wavelength ranges covered by each complete échelle order are (top to bottom, and from left to right) 1.44-1.70 tex2html_wrap_inline897m, 1.62-1.91 tex2html_wrap_inline897m, 1.86-2.18 tex2html_wrap_inline897m, and 2.17-2.54 tex2html_wrap_inline897m.

As Figure 2 shows, the tex2html_wrap_inline1093 of UNSWIRF is indeed more complex than equation 2 would suggest for a simple dielectric coating. From equations (1) and (2), it follows that for normal incidence in air


 equation162

The upper dashed line in Figure 2 corresponds to tex2html_wrap_inline1093 for a constant physical etalon spacing d of 52.0 tex2html_wrap_inline897m, while the lower dashed line is for tex2html_wrap_inline1139m. For tex2html_wrap_inline1141m, the effective spacing between the plates (tex2html_wrap_inline1143) is 52.0 tex2html_wrap_inline897m, but this grows rapidly at around 2 tex2html_wrap_inline897m due to the nature of the coatings, and is more like tex2html_wrap_inline1149m out to the long wavelength cutoff of UNSWIRF.

  figure171
Figure 2: Variation of Free Spectral Range tex2html_wrap_inline1093 with wavelength tex2html_wrap_inline921. The upper dashed line indicates a constant physical etalon spacing of 52.0 tex2html_wrap_inline897m, while the lower dashed line indicates a constant physical etalon spacing of 61.0 tex2html_wrap_inline897m.

The same data, analysed using equation 5, indicate that the change in spacing resulting from one step in Z is


equation180

Having determined tex2html_wrap_inline1097 over the full wavelength range of UNSWIRF, the accuracy of this calibration has been tested by comparing the predictions of equation 5 for the wavelengths of the arc-lines measured earlier with their known wavelengths. We find an r.m.s. accuracy of 0.04 nm, or tex2html_wrap_inline1163% of the instrumental resolution.

 

Parallelism, Resolving Power, and Finesse

Having calibrated the relationship between d, Z, and tex2html_wrap_inline921, it now becomes possible to set the etalon spacing for the peak transmission of any required wavelength. The next step is to make the plates as parallel as possible, and maintain this parallelism over the full wavelength range, given the residual non-flatness of the plates and/or their coatings. This can most easily be done by scanning in wavelength across an unresolved emission line (e.g., from a discharge lamp or from the OH airglow) and building up a ``cube'' of images, with the third axis representing etalon spacing Z. For each spatial pixel (x,y) in the cube, the spectrum of intensity I(Z) is analysed with the UNSWFIT routine (Section 3.4) to determine the position tex2html_wrap_inline1177 and intensity of the emission peak. Owing to the dependence on tex2html_wrap_inline919 of the condition for peak transmission (equation 1), tex2html_wrap_inline1177 should be a maximum near the centre of the array, and decrease towards the edges, as shown in Figure 3. From analysing such images, we find that the etalon surfaces are flat to tex2html_wrap_inline1183 at a wavelength of 1.65 tex2html_wrap_inline897m, well within the specification.

The way we have chosen to monitor the parallelism is to fit a tangent plane to this surface, and then use the measured x- and y-slopes to correct the parallelism settings passed to the etalon from the CS-100. The parallelism has been found to be weakly, but repeatably, dependent on Z (and thus on tex2html_wrap_inline921, due to the fact that the effective reflection at different wavelengths comes from different depths within the coatings, and the coating thicknesses vary slightly), and this is now accounted for by the observing control software. After optimising the parallelism, we have calculated the resolving power and effective finesse tex2html_wrap_inline931 (equation 4) of UNSWIRF using a series of discharge lamp lines over the available wavelength range, as tabulated in Table 3. The high reflectivity R of the plates will make tex2html_wrap_inline929 the dominant contributor to tex2html_wrap_inline931. In addition, the throughput of the etalon at each wavelength (except for 2.334 tex2html_wrap_inline897m) has been measured by comparing the peak intensities obtained with the etalon in and then out of the beam.

  figure197
Figure 3: UNSWIRF scan of the Ar 1.6520 tex2html_wrap_inline897m line, incrementing the etalon spacing by 2Z units each time. Etalon spacing increases from left to right, and from top to bottom.

   

Line

n tex2html_wrap_inline1211 tex2html_wrap_inline927 tex2html_wrap_inline931 Throughput
Wavelength (tex2html_wrap_inline897m)tex2html_wrap_inline965 (nm) (%)
Ar 1.6437 64 0.378 4348 68 53
Ar 1.6520 65 0.354 4666 72 62
Kr 2.1165 51 0.386 5483 107 56
Kr 2.1903 47 0.437 5011 107 70
Ar 2.2077 50 0.470 4697 94 43
Kr 2.2486 48 0.409 5497 115 43
Ar 2.3133 47 0.560 4130 88 15
Kr 2.3340 46 0.515 4532 99 tex2html_wrap_inline1221
Table 3: Etalon parameters as a function of wavelength

tex2html_wrap_inline965Line wavelengths in air.

The shift in peak transmitted wavelength for the Kr 2.1165 tex2html_wrap_inline897m line, going from the centre to the edge of the etalon, is quite small compared with the instrumental resolution. As can be seen from Figure 4, the shift is tex2html_wrap_inline1227 over the inner 90 pixel diameter, and still tex2html_wrap_inline1229 over the entire usable field of view. In fact, owing to possible non-uniform illumination of the etalon by the discharge lamp, Figure 4 may slightly overestimate this shift. Thus, UNSWIRF is virtually monochromatic, and can in principle be used as a pure tunable line imaging filter, provided velocity gradients and dispersions are small (<50 km tex2html_wrap_inline893), and the line centre wavelength is known in advance. Otherwise, more extensive scanning in wavelength will be necessary (but this of course furnishes, as a spinoff, the velocity field). Although the ability of UNSWIRF to resolve lines is limited by its instrumental profile width (tex2html_wrap_inline1235 km tex2html_wrap_inline893, depending on parallelism), we have found that the profile fitting allows us to measure shifts in the position of the line peak equivalent to velocity changes of <10 km tex2html_wrap_inline893, depending on signal-to-noise of the data, profile shape, and plate parallelism.

  figure226
Figure 4: Contour plot of the shift in transmitted wavelength with position on the IRIS array for the Kr 2.1165 tex2html_wrap_inline897m line, as a fraction of the instrumental profile width tex2html_wrap_inline1211. The heavy black line around the edge marks the unvignetted field of view of the etalon. Beginning from the top right edge, the contours mark a shift of -80%, -40%, -20%, -10%, and -5% of tex2html_wrap_inline1211 relative to the mean peak wavelength at the array centre.

Observing strategies

Observing with UNSWIRF is much like normal narrow-band imaging in the near-infrared. As recommended by the IRIS manual (Allen 1993), readout method 4 is employed, which breaks up each exposure into a series of tex2html_wrap_inline1259 Non-Destructive Readouts (NDRs), enabling on-the-fly bias correction and linearisation, and yielding the lowest possible read-noise (typically tex2html_wrap_inline1261 etex2html_wrap_inline1263). To guarantee background-limited performance, each exposure at a given etalon Z setting is normally 120 s in duration at K (180 s at H), broken up into 12 NDRs. In order to reduce overheads, a complete scan in Z on the object is normally done before moving the telescope to an offset sky position tex2html_wrap_inline1273 away, and repeating the sequence. Except when the scan crosses a strong, and rapidly varying, OH airglow line (usually more of a problem in H-band than K-band), sky subtraction is found to be quite adequate, even when sky frames are taken 15-20 minutes after the matching object frame.

An existing procedure for commanding telescope spatial offsets, written in the AAO DRAMA environment (Bailey et al. 1995), has been enhanced with the ability to request etalon spacing and parallelism changes from the CS-100. To guard against possible drifts in the UNSWIRF etalon parallelism (usually in response to changes in the ambient temperature and/or humidity), the parallelism is normally checked immediately prior to each night of observing, by scanning a calibration lamp line close to the region of interest.

Because of the monochromatic nature of the infrared radiation reaching the IRIS array from UNSWIRF, it is essential that matching sky exposures and dome flatfields be obtained for all of the etalon Z settings used on an object, as otherwise severe fringing can result. Similarly, bright spectroscopic standard stars need to be observed once, and preferably twice at these same Z settings. Although accuracy of the photometry is usually limited by the sky and the array to tex2html_wrap_inline1285% at best, it is necessary to determine the intensity scaling of the continuum images relative to the line peak in order to ensure proper continuum subtraction.

 

Data Analysis

A schematic of the basic data reduction procedure for UNSWIRF is shown in Figure 5, and begins with subtraction of a matching sky frame, followed by division by a normalised matching dome flat. Monochromatic imaging of simple sources then requires just a scaling and subtraction of the off-line frame. For more complex sources, a ``cube'' is constructed from a sequence of such monochromatic images at a constant Z interval, aligned to a common spatial frame defined by field stars. The moments of this cube (integrated line intensity I(x,y), Z position of the line peak tex2html_wrap_inline1177, and the line width tex2html_wrap_inline1295) are extracted by fitting a Lorentzian to the spectrum at each spatial pixel. In general, the line profile is usually too noisy to allow three free parameters (the base level having already been set to tex2html_wrap_inline1297 by the off-line subtraction). Since in most cases the emission-line profile will be unresolved by UNSWIRF, the line width tex2html_wrap_inline1295 can be assumed to be the same as the instrumental profile width, as mapped by the calibration lamp line scans, leaving only two free parameters in the fitting. Finally, all pixels in the vignetted corner regions of the moment maps, as well as any pixels which fall below a specified intensity threshold, or for which the fitted tex2html_wrap_inline1177 lies outside the actual Z range scanned (assuming that the observations did adequately span the line of interest) are blanked out.

  figure245
Figure 5: Schematic diagram of data reduction steps, for monochromatic imaging, using on- and off-line images of the Htex2html_wrap_inline1069 emission around OMC-1.

In order to streamline the processing of UNSWIRF data, a suite of programs has been written using the IRAFgif environment. A listing of these programs and their functions is given in Table 4. With the exception of UNSWFIT, these programs are scripts written in the IRAF Command Language (CL) which execute a series of existing IRAF routines. The UNSWFIT task is a purpose-written SPP (Subset Pre-Processor) program that uses a Lorentzian-fitting algorithm supplied by F. Valdes.

   

Task name Purpose
UNSWBLANK Execute UNSWFIT, then blank incongruous pixels in output maps.
UNSWCAL Convert intensity from (etex2html_wrap_inline1263) stex2html_wrap_inline941 to ergs cmtex2html_wrap_inline1019 stex2html_wrap_inline941 pixeltex2html_wrap_inline941.
UNSWCUBE Sky-subtract, flatfield, rotate, clean, align, continuum-subtract
and stack a series of consecutive images into a datacube.
UNSWDISP Display a ``movie'' of the datacube planes.
UNSWFIT Fit Lorentzian profiles to each datacube pixel, output maps of intensity,
wavelength shift, and profile width.
UNSWFLATS Produce and label flatfields.
UNSWLIN Convert data from ADUs to electrons.
UNSWMASK Mask a map using the same blanking as another map.
UNSWMERGE Sort and stack a series of processed images into a datacube, averaging
repeat data where available.
UNSWPHOT Carry out aperture photometry on a sequence of standard star images.
UNSWPROC Sky-subtract, flatfield, rotate, and clean a series of consecutive images.
UNSWSLOPE Execute UNSWFIT, fit tangent plane to wavelength shift map, and compute
parallelism corrections.
UNSWSPEC Plot a spectrum of intensity vs Z, averaged over a range in x and y.
UNSWVEL Correct wavelength shift map for instrumental shift, convert from Z to km tex2html_wrap_inline893.
Table 4: Summary of UNSWIRF data reduction tasks


Next Section: Results
Title/Abstract Page: UNSWIRF: A Tunable Imaging
Previous Section: Overview of the Instrument
Contents Page: Volume 15, Number 2

Welcome... About Electronic PASA... Instructions to Authors
ASA Home Page... CSIRO Publishing PASA
Browse Articles HOME Search Articles
© Copyright Astronomical Society of Australia 1997
ASKAP
Public