Evaluation of HD 80606 Observing and Image Analysis Strategies
- Defocusing is the Winner -
This web page describes an evalutaion of several observing and image analysis strategies for the exoplanet HD 80606b. The usual procedures of maintaining sharp focus and using a photometry aperture radius ~4 times FWHM was compared to an alternative observing strategy, defocusing, which for my telescope system produced superior results. This is due in part to the brightness of HD 80606 and also to the fact that its equally bright orbital companion, HD 80607, is located only 20 "arc away.

Introduction

HD 80606 has an exoplanet that orbits with a 111-day period in a very eccentric orbit (0.93). The orbit is inclined close to 90 degrees because when it is behind the star it undergoes a "secondary" transit, as observed with the Spitzer Space Telescope by Laughlin et al (2009). It is not known if the planet transits in front of the star, for a primary transit, because that part of its orbit is much farther from the star. An opportunity for observing a possible primary transit exists on Feb 13/14 UT of 2009. This web page has been created for offering suggestions on the best strategy for conducting these observations and for processing the images to produce a light curve. 

Hardware

My telescope is a 11-inch aperture Celestron (CPC1100) mounted on an equatorial wedge. I use a focal reducer and SBIG AO-7 image stabilizer in front of a SBIG ST-8XE CCD camera. The AO-7 employs a tip/tilt mirror that can be operated at several Hz, and when mirror motion exceeds a user-specified limit a correcting nudge signal is sent to the telescope drive motors. The polar axis has been aligned to the north celestial pole with an accuracy of ~2 'arc. This assures that when the AO-7 keeps an autoguided star fixed to a pixel location on the autoguider chip the star field does not rotate about that star and simultaneously move the star field across the main CCD chip. This arrangement assures that the star field is fixed with respect to the CCD pixel field to within a few pixels during the course of a many hour observing session. The image scale is 1.21 "arc/pixel, and the FOV is 31 x 21 'arc. I use a wireless focuser (Craycroft type) for adjusting the position of the optical back-end (the mirror is not moved). MaxIm DL 5.03 is used to control the telescope, wireless focuser, AO-7 and CCD camera.

Each observer's hardware is different, and atmospheric seeing is different each night, so the following is presented with the purpose of illustrating concepts that can be used to guide observers in choosing optimum observing and image analysis strategies.

Review of Usual Strategy for Exoplanet Transit Light Curves

Normally, transiting exoplanet light curves are produced using sharply-focused images and a photometry aperture radius ~ 4 x FWHM. Slightly smaller and larger photometry apertures can sometimes produce small improvements (e.g., ApertureRadius/FWHM within range 3 to 5). However, when a star is close to the target star it may not be possbile to employ this standard aperture setting. Exposure times are chosen that yield maximum counts for stars of interest that are slightly less than the linearity limit (52,000 counts for my CCD) at times of best atmospheric seeing. Here's an image of the HD 80606 star field.


Figure 1. Star field surrounding HD 80606 (indicated) and it's orbital companion star, HD 80607 (to its left). FOV = 31 x 21 'arc, north up, east left. FWHM ~ 4.2 "arc.


Figure 2. Zoom of previous image showing a photometry aperture centered on HD 80606 and HD 80607 20.5 "arc to the east. FOV = 9.3 x 8.2 'arc. The aperture radius is 9.7 "arc (8 pixels), which is ~2 x FWHM for most images of this observing session.

Fig. 2 gives the impression that it is possible to employ a small photometry aperture for stars as close as 20.5 "arc (17 pixels). However, the brightness/contrast setting for this image gives a false impression of each star's PSF (point-spread function). Here's a different brightness/contrast choice.


Figure 3.  Zoom showing a photometry aperture radius of 10.9 "arc (9 pixels), which places a boundary of the signal aperture half-way between the two stars.

Fig. 3 shows a signal aperture choice that should be least affected by atmospheric seeing changes that cause the PSF of the nearby star, HD 80607, to "spill over" into the signal aperture while the same seeing change causes the target star's PSF to "spill out" of the signal aperture by approximately the same amount. For this observing session, and for exposure times of 30 seconds (the longest exposure time that still prevents saturation), PSFs have FWHM values that are typically 4.5 "arc. This aperture corresponds to the ratio ApertureRadius /  FWHM  = 2.0, which is not within the customary 3 to 5 range where best results are usually obtained. This suggests that using this aperture it may not be possible to obtain a high precision light curve.

Three Strategy Candidates

In anticipation of problems with the usual observing strategies I note three approaches to observing HD 80606:

    1) Choose a photometry aperture radius of 1/2 the separation of HD 80606 and HD 80607 (i.e., 1/2 of 20.5 "arc ~ 10 "arc), and hope this is adequate (e.g., Fig.3),
    2) Choose a photometry aperture that includes both stars, and accept the added noise from the larger number of pixels within the signal aperture (e.g., Fig. 8, below), and
    3) Defocus all images the same amount and use a large photometry aperture to include both stars (Fig. 10, below).

Note that a defocused image can accomodate a much longer exposure time than a sharply focused image, thus improving the duty cycle (defined as "time exposing an image divided by interval between images"). For example, one of my image sets had an exposure time of 10 seconds (sharply focused) and the defocused images had an exposure itme of 60 seconds. Since my overhead is 12 seconds (8 seconds for image download and 4 seconds for the AO-7 to adjust tracking), the corresponding duty cycles are 45% and 83%. The number of photons collected per unit time are in the same ratio, 45 and 83, so the defocused image set yields almost twice as many photons per unit observing time as the 10-second image set. Duty cycle advantages for defocusing are greater for bright stars since brightness calls for short exposure times for sharply focused images (to avoid saturation).

Strategy 1: Sharp focus, small aperture

This strategy is as close to what normally is done as possible; only the photometry aperture is smaller than normal because of the nearby star HD 80607. The exposure time was 30 seconds, which kept all stars in the FOV just below saturation. The following image shows my choice for reference stars.


Figure 4. Sharply focused and 30-second exposure star field with 10 reference stars indicated.

This strategy was employed for an hour, and the next figure shows processed light curve.


Figure 5. Light curve of HD 80606 for the sharply focused, 30-second exposure image set.

The single image RMS is equivalent to a normalized 2-minute exposure having RMS = 3.5 mmag per image. Taking the 71% duty cycle into account these data correspond to achieving a precision of 4.15 mmag per 2 minutes of observing time. This is greater than I normally achieve (edven for 11th magnitude stars), so this appears to support my expectation that the precision suffers from having to use ApertureRadius/FWHM = 2 instead of the normally optimum 4.

It occured to me that I could improve precision if shaper images could be achieved. One way to improve sharpness is to use shorter exposures. So I took another hour of 10-second exposures. 


Figure 6. Shorter exposure image (10 seconds instead of 30 seconds) showing the same photometry aperture setting as before (9x20x12 pixels). 


Figure 7.
Light curve for above set of images.

The precision per image, normalized to a 2-minute exposure, is 2.06 mmag, which is the kind of improvement that what I was hoping for. However, 10-second imaging has a duty cycle of 45%, and adjusting for this yields a precision of 3.07 mmag per 2 minutes of observing time. This is still significantly better than the precision of the 4.15 mmag per 2 minutes of observing time calculated from the 30-second exposure image set.

Strategy 2: Sharp focus, large aperture

The 10-second images were processed using an aperture that included the flux fromboth HD 80606 and HD 80607, as shown in the next image.


Figure 8. Using a large aperture for the same images as in the previous analysis (10-second exposures). Aperture dimensions are 28x8x16 pixels.


Figure 9. Light curve from 10-second, sharply focused images but using a large photometry aperture that includes both stars. 

The large aperture led to a worsened RMS per image. After correcting for duty cycle the precision is 4.61 mmag per 2-minutes of observing time.

Strategy 3: Defocused, large aperture

The following defocused image is from a 1-hour observing session immediately following the previous image set.


Figure 10. Defocused image showing the same large photometry aperture used with the sharply focused images. Exposure time is 60 seconds.


Figure 11. Light curve from defocused images and large photometry aperture that included both stars.

The 2-minute equivalent RMS per image is 2.02 mmag, and correcting for duty cycle yields 2.22 mmag per 2 minutes of observing time. This is better than any of the other observing and analysis strategies.

Why? 

I used to think that the only reason for defocusing was to improve duty cycle. But compare the "sharply focused image set with large aperture" with "defocused image set with large aperture" (same apertures), and the RMS per 2-minute equivalent image, unadjusted for duty cycle, is better for the defocused image than for the sharp image. The two precisions are 3.09 mmag per image and 2.02 mmag per image. Since duty cycle has nothing to do with these two precision calculations there's something going on that has to do with the spreading out of where photo-electrons are produced.

Could it be that when most of the star's photo-electrons originate from only a few pixels that the readout noise per pixel, or the thermal noise per pixel, has a greater percentage effect on the total flux than when many pixels are involved? When most of the photo-electrons come from just a few pixels could small errors in the flat field contribute a greater error in the total flux than when these flat field errors are averaged-out over more pixels? Could someone help me understnad this?

Conclusion

It's not necessary to understand why a defocused observing session produces better results than a sharply-focused observing session, for a specific target, when all we want to know is how to achieve the best quality light curve. I concluse that for my hardware, and this target, and the atmospheric seeing that existed for this test, it is better to observe defocused!

Since other observers will have different image scales and different atmospheric seeing, if the best quality light curve is to be achieved on February 13/14 it will be necessary for each observer to perform a simple version of what I described on this web page.

____________________________________________________________________

WebMaster: B. GaryNothing on this web page is copyrighted. This site opened:  2009.01.30 Last Update:  2009.01.30