Open Access Paper
13 August 2014 Modeling instrumental field-dependent aberrations in the NIRC2 instrument on the Keck II telescope
Breann N. Sitarski, Gunther Witzel, Michael P. Fitzgerald, Leo Meyer, Andrea M. Ghez, Randall D. Campbell, Jessica R. Lu, Keith Matthews, Peter Wizinowich, Jim Lyke
Author Affiliations +
Abstract
We present a model of field-dependent aberrations arising in the NIRC2 instrument on the W. M. Keck II telescope. We use high signal-to-noise phase diversity data employing a source in the Nasmyth focal plane to construct a model of the optical path difference as a function of field position and wavelength. With a differential wavefront error of up to 190 nm, this effect is one of the main sources of astrometric and photometric measurement uncertainties. Our tests of temporal stability show sufficient reliability for our measurements over a 20-month period at the field extrema. Additionally, while chromaticity exists, applying a correction for field-dependent aberrations provides overall improvement compared to the existing aberrations present across the field of view.

1.

INTRODUCTION AND GOALS

Anisoplanatism and the field-dependence of the anisoplanatic point spread function (PSF) greatly influences astrometric and photometric accuracy. In crowded fields like the Galactic Center, insufficient knowledge of the PSF variation accounts for the largest error in astrometric measurements (Witzel et al., this volume). The varying morphologies affect not only astrometric analyses, but photometric analyses as well. As in our first publication on the project,1 we divide the definition of the optical transfer function (OTF) into two parts: a static, instrument-based term, and the Adaptive Optics (AO) response to the atmospheric PSF term that takes into account the 00211_psisdg9148_91486T_page_1_1.jpg profile of the atmospheric turbulence.

A static, instrument-based term can still have field dependences and can be further decomposed into on- and off-axis terms.1 The on-axis term is expected to be estimated by other means, such as on-sky phase diversity (Jolissant et al., this volume, and Ragland et al., this volume) or through crowded-field estimation. We present a fiber phase diversity phase map grid with data taken between October 2012 and June 2014 that uses an improved data collection methodology. We study the effects of time variability in our instrumental data as well as the affects of chromaticity in our data and discuss whether we can apply this Fe II-based grid to other data sets, including broadband K′ data.

2.

FIBER PHASE DIVERSITY DATA COLLECTION

We work to measure the differential instrumental aberrations in NIRC2 with phase diversity data taken with a illuminated fiber internal to NIRC2 itself. By measuring the phase of aberrations at a specific position on the detector, we aim to transform phases at specific points to phases at arbitrary field points. Keck Observatory currently uses phase diversity to optimize static aberrations in their image sharpening procedure, which is run prior to every adaptive optics (AO) night at Keck Observatory on the NIRC2 instrument. In order to extract the phase at each point, an image is taken at Δz = -2 mm, -4 mm, and -6 mm, where Δz is difference in defocus term.

Traditionally, image sharpening only needs to focus on low-order modes, and phase diversity collection made use of short exposure times (tint = 0.181 sec, ncoadds = 50, and the CDS readout mode). While this method is efficient in terms of total integration time spent at each field position (72 seconds), it produces a very low signal-to-noise ratio phase diversity map (Figure 1). This is, in large part, due to insufficient integration times for the more out-of-focus images that are used when the phase is extracted.

Figure 1.

A comparison between the old method of collection Phase Diversity data and our new method, which results in a much higher signal-to-noise ratio phase map.

00211_psisdg9148_91486T_page_2_1.jpg

Aiming to obtain higher signal-to-noise ratio phase maps, we developed an improved method of data collection. Instead of keeping the integration time constant for each Δz, we increased the integration times and number of reads and decreased the number of coadds (Table 1). While increasing the overall integration time (213 seconds), this yields a phase diversity map with only 12.5 nm RMS random measurement error in each phase map.

Table 1.

An Improved Method for Phase Diversity data collection

Δz (mm)tint (sec)ncoaddsnreads
00.18150CDS
-21.5208
-47616
-6180264

In order to get the phase at each position, we image sharpen at the center (at pixel [545, 495] measured from zero-based index) then take phase diversity data across the field of view. By extracting phase information at various points, we can measure the field dependence of the aberrations.

2.1

Wavefront Error and Field Variation Parameters

In total, we have taken ~55 hours of fiber phase diversity data across the field of view of NIRC2 in the fashion described in the above section. This yields approximately 150 phase diversity phase maps, some taken at the same position over various times to measure the affects of time variability on the data (see section X). The full grid of phase diversity data is shown in Figure 2, where each phase map is centered at the position it was acquired. All of these data were taken at the Fe II (λ0= 1.65 μm) narrow-band filter, following the image sharpening procedure at Keck Observatory. The effects of chromaticity are described in section X.

Figure 2.

Our full phase diversity data grid, where phase diversity data is taken at various points across the field of view. There are 150 different data sets (some taken at the same position to study time variability affects). All data is taken through the Fe II narrow-band filter on NIRC2 (λ0 = 1.65 μm).

00211_psisdg9148_91486T_page_3_1.jpg

At each field position, we take three phase diversity data sets (including the OPD map, and images at Δz = 0, -2, -4, and -6 mm), allowing us to determine the errors on the Zernike terms. We acquire a center position data set ([x, y] = [545, 495]) at the beginning and end of each data acquisition run which are averaged together. The center position is then subtracted out from all other field points for that specific data acquisition run in order to obtain the differential position with regards to the center position.

The field extrema points are those furthest from the center position-the four corners of the detector. The relative wavefront error between the center reference point and the corner points is 190 nm RMS at Fe II. In fact, the wavefront error difference between adjacent samples can be used to quantify the error associated with the field sampling. The average wavefront error difference is 43 nm RMS at Fe II. The smallest and largest wavefront error differences between adjacent OPD maps is 87 and 15 nm RMS at Fe II, respectively.

Adjacent samples allow us to quantify the error associated with field sampling. Figure 3 shows a map and the connecting points of the adjacent phase maps with their wavefront error; Figure 3 also plots the difference in RMS Wavefront Error as a function of the distance to the adjacent phase map. While the maps with nearby neighbors have fairly low error (below the temporal error), they exhibit quite a lot of scatter. This can be caused by field variation because phase maps near extrema may differ more drastically than those close to the image sharpening position. The average wavefront error difference between adjacent samples is 43 nm RMS, while the smallest/largest WFE difference between adjacent samples is 15 and 87 nm RMS, respectively. These all fall below the average wavefront error difference relative to the center of the field (96.5 nm RMS).

Figure 3.

Grid of adjacent phase maps and the associated wavefront error between the adjacent samples.

00211_psisdg9148_91486T_page_4_1.jpg

3.

WAVELENGTH VARIATION OF PHASE DIVERSITY DATA

We take all of our phase diversity data at the narrow-band Fe II filter (λ0=1.6544 μm; Δλ=0.0256 μm). However, all Galactic Center data is taken through the K′ (λ0=2.124 μm, Δλ=0.351 μm) and L′ (λ0=3.776 μm, Δλ=0.700 μm) filters. Furthermore, all Keck image sharpening data is also taken at Fe II; it is therefore necessary to see if the phase diversity phase map grid outlined in section X is usable at other wavelengths.

We took phase diversity data at several non-thermal narrow- and medium-band filters available on NIRC2: He 1 A (λ0=1.0847 μm; Δλ=0.0182 μm); Pa γ0=1.096 μm; Δλ=0.016 μm); J continuum (λ0=1.2132 μm; Δλ=0.0198 μm); Pa β (λ0=1.2903 μm; Δλ=0.0193 μm); H continuum (λ0=1.5804 μm; Δλ=0.0232 μm); CH4S (λ0=1.5923 μm; Δλ=0.1257 μm); Fe II; CH4L (λ0=1.6809 μm; Δλ=0.1368 μm); He 1 B (λ0=2.0563 μm; Δλ=0.0326 μm); Br γ (λ0=2.1686 μm; Δλ=0.0326 μm); and K continuum (λ0=2.2706 μm; Δλ=0.0296 μm). We acquired the phase diversity using the method described in Section X. However, we optimized integration times to obtain less than 10,000 DN/coadd since nonlinearity becomes an issue with the NIRC2 detector.

Assuming that the phase maps extracted from the phase diversity data are linear combinations of the Zernike polynomials, we extract the defocus term from each phase map, as it is visually apparent that there is some defocus offset (see Figure 4). In plotting this defocus term as a function of wavelength, some sort of dependence of defocus on wavelength becomes apparent (Figure 5). Keck Observatory has implemented a focus compensation for closed loop operation on sky to account for the focus offset. The defocus could arise due to a transmissive slab in a converging beam (e.g., dichroic, window).

Figure 4.

Center phase maps extracted from data taken at various narrow- and medium-band filters available on NIRC2. There is a visible apparent focus offset.

00211_psisdg9148_91486T_page_4_2.jpg

Figure 5.

A plot of the extracted Zernike coefficient for the defocus mode as a function of wavelength. There is a linear trend between the wavelength and magnitude of the defocus term.

00211_psisdg9148_91486T_page_5_1.jpg

4.

COMPATIBILITY OF FE II AND BR γ FIELD OPD MAPS

Since our entire grid of phase diversity data is currently taken at the Fe II filter, we need to determine if we can extend the use of our grid to different wavelengths, primarily since the majority of our Galactic Center science data is taken at the broadband K′ filter (λ0=2.124 μm; Δλ=0.351 μm).27

Comparing the field extrema (the corner s of the detector) allows us to investigate the wavelength dependence of the aberrations. After subtracting out the respective center phase maps for a given wavelength from the extrema phase maps, we compare the Br γ and Fe II optical path differences. A Zernike decomposition of the optical path differences is shown in Figure 6. The Br γ and Fe II are comparable to within ~80 nm RMS. While that is a large value, it overall much less than the 190 nm extrema-to-center difference. Even so, it is still more advantageous to apply this instrumental model with the Fe II data than to not apply the model at all.

Figure 6.

Zernike decomposition of the optical path differences between Brγ and Fe II at three of the field extrema. The lower-order terms clearly dominate the overall wavefront error.

00211_psisdg9148_91486T_page_6_1.jpg

Due to the brightness of the fiber, K′ data was not obtained at the field extrema since we had to subarray down to a 320 × 312 grid around the very center of the NIRC2 detector in order to not saturate the in-focus Δz = 0 mm image. We obtained positions at the extrema of this subarray, and decomposed them in the way described above but at K′ and corresponding field positions at Fe II. The Zernike decomposition is shown in Figure 7. The largest optical path differences is 50.9 nm RMS.

Figure 7.

Zernike decomposition of the optical path differences between K′ and Fe II at positions fairly close to the center of the detector that stayed in the 320 × 312 subarray of the detector. The lower-order terms again dominate the overall wavefront error.

00211_psisdg9148_91486T_page_6_2.jpg

5.

TIME VARIABILITY

The time difference between acquiring many of the field positions, ranging from October 2012 to June 2014, poses a major concern in using a grid of phase diversity data. In order to probe the effects of time variability, we took multiple data sets at the field extrema and the center positions over this time period.

We took the center position twice over a span of five hours to probe the temporal drift on short time periods. Since we take three sets of data (and therefore three phase maps) at each position, we use an average phase map for each time and subtract them to solve for the residual. The temporal drift over 5 hours is therefore 18 nm RMS (see Figure 8). Between October 2012 and June 2014, the temporal drift between the center phase maps is 41.6 nm RMS, obtained using a similar analysis as described above. The measurement repeatability at the center positions is 12.5 nm RMS.

Figure 8.

Variation of the center phase map over 5 hours. The residual seems to consist of high-order noise.

00211_psisdg9148_91486T_page_7_1.jpg

We carried out a similar exercise at the field extrema. Between October 2012 and May 2014, the temporal drift over the field extrema is 53.1 nm RMS, with a measurement repeatability of 15.43 nm RMS. It is noteworthy that the 53.1 nm RMS temporal difference falls significantly below the 190 nm RMS center to corner difference in phase maps (see Figure 9).

Figure 9.

Variation of the one of the extrema points (corner positions) as a function of time. The general morphology of the phase map stays the same, but the residual seems to consist of high-order noise.

00211_psisdg9148_91486T_page_7_2.jpg

However, it appears that much of the error is common mode. The top row of Figure 10 shows phase maps from three different epochs. At each extrema point, we averaged over all the epochs and subtracted that to yield the second row. To determine if there were any common-mode errors in a given epoch, we averaged the results of the second row over a given epoch and subtracted it from row 2. This results in field-dependent wavefront error stability after removing common-mode errors of <20 nm RMS.

Figure 10.

There seems to be common-mode error that exists across the different epochs and different field points. The middle row shows the average of the top row subtracted out; the last row shows the average of all extrema data at the respective epoch subtracted out.

00211_psisdg9148_91486T_page_8_1.jpg

6.

ON-SKY TESTS AND VALIDATION

In Figure 11, we show an on-sky test of the instrumental model that we have developed. The observed on-axis PSF is the guide star for the Natural Guide Star AO system we utilized on the Keck II telescope.8, 9 We transform the observed on-axis PSF to the model shown on the far right of the image by multiplying by the OTF ratio generated from the nearest phase maps of the on-axis and off-axis points. It is clear from the comparison of the middle frame (observed off-axis) and the far-right frame (model) that many complex features are reproduced in our model, as denoted by the arrows.

Figure 11.

Transformation of the on-axis PSF to the model. Complex structure from the observed off-axis PSF is reproduced by our model in this binary star test.

00211_psisdg9148_91486T_page_8_2.jpg

7.

CONCLUSIONS AND FUTURE WORK

We have presented a large phase diversity data collection effort to map the aberrations across the field of view of the NIRC2 detector on the Keck II telescope. This data collection consists of a grid of phase diversity data taken at Fe II (λ0 = 1.65 μm), along with multi-wavelength studies and temporal studies from October 2012 to June 2014.

Our efforts have shown that while time variability and chromaticity have an effect, their are intrinsically less than simply not correcting for the field-dependent aberrations at all arising from the instrumental PSF (see Table 2). In fact, as shown in Figure 5, we may be able to calibrate out the lower-order modes causing a discrepancy between Brγ and Fe II by using a method similar for calibrating out the defocus term.

Table 2.

An Improved Method for Phase Diversity data collection

MeasurementWavefront Error (WFE; nm RMS)
Center to Extrema max WFE Difference190
Temporal drift of center phase map over 5 hours18
Temporal drift of center from October 2012 - June 201441.6
Temporal drift of field extrema from October 2012 - May 201453.1
Measurement repeatability of center phase map12.5
Measurement repeatability of field extrema phase maps15.4
Compatibility of Br γ and Fe II at corners80
Compatibility of K′ and Fe II at K′ corners50.9
Average WFE Difference between adjacent maps43
Smallest/largest WFE difference between adjacent maps15/87
Average WFE difference relative to the center over the field96.5

Future work includes expanding this work to create an interpolated model for each Zernike coefficient across the field of view of the NIRC2 detector. Preliminary results are favorable to this method and would allow for correction not be based on the nearest phase diversity phase map as we are currently implementing. We will combine this work with models of the AO response to the atmospheric PSF for both Laser Guide Star AO and Natural Guide Star AO and apply it to both crowded regions and other science cases (Lu et al., this volume).

8.

ACKNOWLEDGEMENTS

We would like to acknowledge the financial support of the W. M. Keck Foundation for this project. AMG and MM are also supported by NSF AST-0909218. The authors would also like to thank Jay Anderson, Richard Dekaney, Brent Ellerbroek, James Graham, Claire Max, Andrei Tokovinin, and Marcos van Dam for their valuable input to the project. BNS would also like to thank Saundra Albers.

REFERENCES

[1] 

Fitzgerald, M.P., Witzel, G., Britton, M.C., Ghez, A.M., Meyer, L., Sitarski, B.N., Cheng, C., Becklin, E.E., Campbell, R.D., Do, T., Lu, J.R., Matthews, K., Morris, M.R., Neyman, C.R., Tyler, G.A., Wizinowich, P.L., and Yelda, S., , “Modeling anisoplanatism in the Keck II laser guide star AO system,” Proc. SPIE 8447 (July 2012).Fitzgerald, M.P., Witzel, G., Britton, M.C., Ghez, A.M., Meyer, L., Sitarski, B.N., Cheng, C., Becklin, E.E., Campbell, R.D., Do, T., Lu, J.R., Matthews, K., Morris, M.R., Neyman, C.R., Tyler, G.A., Wizinowich, P.L., and Yelda, S., , “Modeling anisoplanatism in the Keck II laser guide star AO system,” Proc. SPIE 8447 (July 2012).

[2] 

Ghez, A.M., Hornstein, S.D., Lu, J.R., Bouchez, A., Le Mignant, D., van Dam, M.A., Wizinowich, P., Matthews, K., Morris, M., Becklin, E.E., Campbell, R.D., Chin, J.C.Y., Hartman, S.K., Johansson, E.M., Lafon, R.E., Stomski, P.J., and Summers, D.M., “The First Laser Guide Star Adaptive Optics Observations of the Galactic Center: Sgr A*’s Infrared Color and the Extended Red Emission in its Vicinity,” 635 1087 –1094 (2005). Google Scholar

[3] 

Ghez, A.M., Salim, S., Weinberg, N.N., Lu, J.R., Do, T., Dunn, J.K., Matthews, K., Morris, M.R., Yelda, S., Becklin, E.E., Kremenek, T., Milosavljevic, M., and Naiman, J., “Measuring Distance and Properties of the Milky Way’s Central Supermassive Black Hole with Stellar Orbits,” 689 1044 –1062 (2008). Google Scholar

[4] 

Lu, J.R., Ghez, A.M., Hornstein, S.D., Morris, M.R., Becklin, E.E., and Matthews, K., “A Disk of Young Stars at the Galactic Center as Determined by Individual Stellar Orbits,” 690 1463 –1487 (2009). Google Scholar

[5] 

Yelda, S., Ghez, A.M., Lu, J.R., Do, T., Meyer, L., and Morris, M.R., , “Adaptive optics observations of the galactic center young stars,” Proc. SPIE 8447 (July 2012).Yelda, S., Ghez, A.M., Lu, J.R., Do, T., Meyer, L., and Morris, M.R., , “Adaptive optics observations of the galactic center young stars,” Proc. SPIE 8447 (July 2012).

[6] 

Yelda, S., Ghez, A.M., Lu, J.R., Do, T., Meyer, L., Morris, M.R., and Matthews, K., “Properties of the Remnant Clockwise Disk of Young Stars in the Galactic Center,” 783 131 (2014). Google Scholar

[7] 

Lu, J.R., Do, T., Ghez, A.M., Morris, M.R., Yelda, S., and Matthews, K., “Stellar Populations in the Central 0.5 pc of the Galaxy. II. The Initial Mass Function,” 764 155 (2013). Google Scholar

[8] 

Wizinowich, P.L., Le Mignant, D., Bouchez, A.H., Campbell, R.D., Chin, J.C.Y., Contos, A.R., van Dam, M.A., Hartman, S.K., Johansson, E.M., Lafon, R.E., Lewis, H., Stomski, P.J., Summers, D.M., Brown, C.G., Danforth, P.M., Max, C.E., and Pennington, D.M., “The W.M. Keck Observatory Laser Guide Star Adaptive Optics System: Overview,” 118 297 –309 (2006). Google Scholar

[9] 

van Dam, M.A., Bouchez, A.H., Le Mignant, D., Johansson, E.M., Wizinowich, P.L., Campbell, R.D., Chin, J.C.Y., Hartman, S.K., Lafon, R.E., Stomski, P.J.Jr., and Summers, D.M., “The W.M. Keck Observatory Laser Guide Star Adaptive Optics System,” Performance Characterization, 118 310 –318 (2006). Google Scholar
© (2014) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.
Breann N. Sitarski, Gunther Witzel, Michael P. Fitzgerald, Leo Meyer, Andrea M. Ghez, Randall D. Campbell, Jessica R. Lu, Keith Matthews, Peter Wizinowich, and Jim Lyke "Modeling instrumental field-dependent aberrations in the NIRC2 instrument on the Keck II telescope", Proc. SPIE 9148, Adaptive Optics Systems IV, 91486T (13 August 2014); https://doi.org/10.1117/12.2057140
Lens.org Logo
CITATIONS
Cited by 5 scholarly publications.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Iron

Wavefronts

Data centers

Adaptive optics

Point spread functions

Instrument modeling

Sensors

Back to Top