Confidence Level and Sensitivity Limits in High Contrast Imaging
In long adaptive optics corrected exposures, exoplanet detections are currently limited by speckle noise originating from the telescope and instrument optics, and it is expected that such noise will also limit future high-contrast imaging instruments for both ground and space-based telescopes. Previous theoretical analysis have shown that the time intensity variations of a single speckle follows a modified Rician. It is first demonstrated here that for a circular pupil this temporal intensity distribution also represents the speckle spatial intensity distribution at a fix separation from the point spread function center; this fact is demonstrated using numerical simulations for coronagraphic and non-coronagraphic data. The real statistical distribution of the noise needs to be taken into account explicitly when selecting a detection threshold appropriate for some desired confidence level. In this paper, a technique is described to obtain the pixel intensity distribution of an image and its corresponding confidence level as a function of the detection threshold. Using numerical simulations, it is shown that in the presence of speckles noise, a detection threshold up to three times higher is required to obtain a confidence level equivalent to that at for Gaussian noise. The technique is then tested using TRIDENT CFHT and angular differential imaging NIRI Gemini adaptive optics data. It is found that the angular differential imaging technique produces quasi-Gaussian residuals, a remarkable result compared to classical adaptive optic imaging. A power-law is finally derived to predict the confidence level detection threshold when averaging a partially correlated non-Gaussian noise.
Searching for faint point sources around bright objects is a challenging endeavor. The atmosphere (Roddier, 1981; Racine et al, 1999; Macintosh et al, 2005), telescope and instruments optics (Marois et al, 2003, 2005; Hinkley et al, 2007) produce speckles having a range of timescales that limit the direct detection of faint companions. From previous theoretical analysis, it is known that the speckle intensity temporal distribution is a modified Rician (Goodman, 1968; Soummer & Aime, 2004; Fitzgerald & Graham, 2006). If a large number of uncorrelated speckle realizations are coadded, from the central limit theorem, then the final residual speckle noise follows a Gaussian intensity distribution. Since atmospheric turbulence produces random speckles that have a very short correlation time, a Gaussian distributed residual speckle noise is commonly assumed for ground-based adaptive optics (AO) long integrations and a detection threshold of 5 is usually considered.
However, careful residual noise analysis of AO images have demonstrated that long exposures are not limited by random short-lived atmospheric speckles but by quasi-static speckles (Marois et al, 2003; Marois, 2004; Masciadri et al, 2005; Marois et al, 2005, 2006) originating from the telescope and instruments. The speckle noise currently limiting high-contrast ground-based imaging is thus very similar to that limiting space-based observations (Schneider & Silverstone, 2003). The typical lifetime of ground-based quasi-static speckles has been found to be several minutes to hours (Marois et al, 2006; Hinkley et al, 2007); the noise in the combination of several images spanning 1 hr is very similar to that in a single image (see Fig. 1 for an example; acquired with NIRI/Altair at the Gemini telescope). In this case, since the quasi-static speckle noise is well correlated for the entire sequence, the central limit theorem does not apply and the speckle noise in the final combined image will be non-Gaussian. Sensitivity limits calculated assuming Gaussian statistics would have lower confidence levels (CL). Finding a robust technique to estimate proper sensitivity limits is fundamental to analyze adequately the sensitivity of an exoplanet survey as a function of angular separation. The contrast limit reached by a survey plays a central role in Monte Carlo simulations to derive exoplanet frequencies around stars and constrain planet formation scenarios (Metchev, 2006; Carson et al, 2006; Kasper et al, 2007; Lafreniere et al, 2007a). Understanding the residual noise statistical distribution is thus important for future dedicated surveys of next generation AO systems like NICI (Ftaclas et al, 2003), the Gemini Planet Imager (GPI, Macintosh et al, 2006), the VLT SPHERE (Dohlen et al, 2006), and as well as future space observatories.
In this paper, a new technique is presented to estimate sensitivity limits of a noise showing arbitrary statistics using a CL approach. The theory behind speckle statistics is summarized in § 2. Then § 3 presents a technique to derive detection thresholds using the probability density function and associated CLs. The technique is applied to simulated (§ 4) and observational (§ 5) data to confirm the theory and to validate the technique. The effect of averaging a sequence of independent non-Gaussian noise realizations is discussed in § 6.1. Concluding remarks follow in § 7.
2 Speckle Noise Statistics
Following the work of Goodman (1968); Soummer & Aime (2004); Fitzgerald & Graham (2006), the speckle intensity probability density function (PDF) for one location in the image plane and random temporal phase errors can be shown to be a modified Rician (MR) function. At a specific location in the image plane, the MR PDF is a function of the local time-averaged static point spread function intensity and random speckle noise intensities :
where is the point spread function (PSF) intensity () and is the zero-order modified Bessel function of the first kind. At a specific point of the PSF, if , relevant to Airy ring pinned speckles, the associated PDF is a Gaussian-like function showing a bright positive tail, while if , relevant to PSF dark rings or coronagaphic PSFs dominated by second order halo speckles, the noise distribution is exponential. The CL for a given detection threshold is simply obtained by:
where is the mean-subtracted PDF. Fig. 2 illustrates the different possible regimes compared to a Gaussian intensity distribution. For a 5 detection threshold, where is the standard deviation of the noise obtained using the robust_sigma IDL algorithm
In the previous speckle PDF analysis, it was shown that the atmospheric speckle noise PDF is obtained by analyzing the temporal variation at one location of the PSF. For a quasi-static speckle noise, this approach is not adequate since the noise does not vary significantly with time. The quasi-static noise PDF can be derived using a very simple argument. If we consider a PSF produced by a circular aperture and if the PDF is obtained by analyzing pixels inside a narrow annulus centered on the PSF core, azimuthal quasi-static speckle noise variations are produced with the same value of (here, is the unaberrated PSF and it is azimuthally symmetric for a circular aperture). The speckle noise inside a narrow annulus and from a single speckle noise realization thus shows the same PDF as a temporal speckle noise variation from random phase screens at any location inside the annulus.
3 Experimental Derivation of the PDF and CL Curves
A robust technique to derive sensitivity limits can be developed using CLs. The pixel PDF inside a specific region of the image is first obtained and the CL curve is then derived and extrapolated to estimate a local detection threshold. To avoid having too many false positive detections without missing possible faint companions, a CL of (5 if Gaussian) is selected here. The basic steps to derive the PDF, to obtain the CL curves and to estimate the CL detection threshold is summarized in table 1.
The local PDF is obtained by producing an histogram of the pixel intensities inside a specific region of the image after subtraction of the mean pixel intensity over the region and division by the noise RMS of the image region. The CL curve as a function of detection threshold can be easily estimated by integrating the PDF inside the interval (see Eq. 2). Due to the limited number of resolution elements (one for PSFs or a pixel for simulated noise images) in an image, the PDF and CL curves will be known up to a certain detection threshold. In theory, for the CL detection threshold considered here, each area where the PDF needs to be estimated should have several million independent resolution elements. In practice, for images typically containing up to (250,000 resolution elements), the PDF will be known only up to CL for Gaussian noise. A model fit using a analysis or a polynomial fit are required to extrapolate the CL curve and obtain the detection threshold corresponding to a CL. Since the CL curves of various distributions are nearly linear in a semi- vs detection threshold plot (see Fig. 2), we have chosen to use a polynomial fit due to its simplicity of implementation, its execution speed, and its accuracy. Due to non-linear effects for detection thresholds near 0, a linear fit is first performed for detection thresholds above 1.5. If the detection threshold for a CL is below 9, a second order polynomial fit is used instead to better approximate the CL curve for quasi-Gaussian statistics. The CL extrapolation accuracy will be analyzed in the next section.
4 Technique Validation with Simulated Data
In this section, the PDF, the CL curve and the detection threshold of simulated data are obtained.
4.1 Simulated PDFs
Simulated noise images using specific PDFs are used to test the algorithm in recovering the proper detection threshold for known PDFs. To test the effect of the image area size on the CL extrapolation accuracy, images of various sizes are produced following a MR of equal to 0.1, 1 and 10 (see Fig. 2). For each size, 25 independent realizations are computed to derive the extrapolation accuracy. Fig. 3 and Table 2 show the CL extrapolation accuracy for simulated statistical distributions. In general, the algorithm slightly underestimates the detection threshold for exponential statistics by 5%, but usually within the 2 error calculated for each area size. Typically, the bigger the area is, the more accurate is the detection threshold. To achieve a detection threshold accuracy of 10% for a CL detection threshold, each PDF needs to be known up to a CL (10,000 resolution elements per area). For instruments with smaller FOVs (several hundred per several hundred ), an area of (2,500 resolution elements) would deliver a detection threshold accuracy of 15% for all types of distribution. A solution to increase the detection accuracy of small FOVs would be to combine observations of several objects of similar magnitudes and observing conditions to increase the number of independant noise realizations in each area.
4.2 Simulated PSFs
The algorithm is now tested using simulated aberrated PSFs. For PSF observations, determination of the PDF is more complex. The speckle noise amplitude is decreasing with angular separation and the PDF may change with angular separation due to relative importance of random atmospheric speckles, photon, background and read noises. Since the speckle noise amplitude decreases with angular separation, a signal-to-noise ratio image is first obtained by dividing the pixel intensities, at each radius, by the standard deviation of the noise at that radius (estimated using the IDL robust_sigma algorithm). Finally, since CL are extrapolated, a compromise needs to be found between having a good radial sampling of the PDF and having sections of images big enough to adequately determine the PDF.
PSF simulations are performed using Fast Fourier Transforms of complex pixel images with a 512 pixels diameter pupil; the full width half maximum (FWHM) of the PSF is 4 pixels. The pupil has uniform amplitude and includes RMS of phase errors generated using a power-law of index . The PSF images are then trimmed to pixels to avoid FFT aliasing effects. A non-aberrated reference PSF is subtracted to remove the Airy pattern and a signal-to-noise image is calculated.
For simplicity, consider the calculation of the PDF within an annulus centered on the PSF. Since the presence of background or companion point sources inside that annulus could bias the statistics for real data (we are assuming that the background star density is such that only one or a few background objects are detected in the field of view (FOV) around any single target; cases with a high background star density will be discussed in section 6.2), we have chosen to divide the annulus in three azimuthal sections containing 50,000 pixels each ( 10,000 resolution elements, see Fig. 4). The median PDF over the three azimuthal sections is calculated. Given the area of these sections, the PDF will be known down to a CL for Gaussian statistics and the detection threshold will be known to 10% accuracy (see Tab. 2).
To further avoid cases where a point source is located at the border of two sections, the entire procedure is repeated by rotating the sections by 30 and 60 degrees which respect to the PSF center, and the median PDF over the three orientations is finally obtained. This procedure is repeated at different angular separations.
Since the PDF is estimated in large areas that may contain speckles with , or , such technique returns an average PDF weighted by the various speckle noise contributions (pinned/unpinned speckles or Gaussian noises). Simulation with and without a coronagraph (simulated with a Gaussian pupil apodizer having a FWHM equal to a quarter of the pupil diameter) and for , and RMS phase aberration are presented (see Fig. 5 and 6). The algorithm clearly detects the MR distribution expected for a pinned speckle dominated PSF and the unpinned (exponential) speckle dominated coronagraphic PSF.
The non-coronagraphic RMS simulations confirm that speckle noise follows a MR (required detection threshold of for a CL), as expected since pinned speckles are dominant for this case. As the quantity of aberrations increases, the ratio of pinned to non-pinned speckles decreases and the noise becomes exponential. For the Gaussian apodized case, since pinned speckles are strongly attenuated, the halo term dominates and the noise is more exponential. Note that none of these curves are expected to be flat as a function of angular separation since the ratio of pinned to unpinned speckles is varying with angular separation, thus changing the pixel intensity distribution, and some noise is expected from the CL curve extrapolation (see Tab. 2). Another simulation was performed using the RMS case to show that if a constant Gaussian noise (background or read noise) is added to the image, the algorithm correctly detects the change of intensity distribution of the pixels at wide separations (see Fig. 7).
In high-contrast imaging observations, a partially correlated reference star PSF is usually subtracted to remove a fraction of the quasi-static speckle noise. Such reference PSF can be obtained by observing a nearby target, by acquiring the same star at another wavelength (simultaneous spectral differential imaging, Marois et al, 2005) or polarization (Potter et al, 2001), or by building the reference using images acquired with different field angles (angular differential imaging, Marois et al, 2006). Such a PSF subtraction is now simulated to estimate how its affect the PDF. The observed PSF is simulated with a RMS phase aberration , with and without a Gaussian apodizer. The reference PSF is constructed by combining a perfectly correlated phase aberration , where is a constant less than 1, with an uncorrelated part such as:
where the total noise RMS of and are equal and the ratio of the noise RMS of and is equal to 0.1, 1 and 10. Unless the background, read noise, random atmospheric speckles or photon noises are achieved by the reference PSF subtraction, the residual PDF is essentially unchanged for cases with and without a coronagraph (see Fig. 8).
5 Application to Observational Data
The steps required to use the algorithm with observational data are similar to the ones described in Tab. 1 and . Only a few additional reduction steps are necessary. Fig. 9 illustrates the various steps of the technique using Gemini data.
Beside the usual data reduction, deviant pixels, like diffraction from the secondary mirror support, must first be masked. Diffraction from the secondary mirror support usually produces a bright concentrated flux emanating from the PSF core along several azimuthal directions. Since this flux is not produced by quasi-static aberrations and is very localized in the image, if we include these pixels in the PDF, they will produced a bright positive tail in the PDF and the CL detection threshold will be overestimated.
Since the main science goal is to detect point sources, noise filtering is also applied to remove the noise that is not at the spatial scale of point sources. An FWHM median filter is first subtracted from the image to reject large spatial period noises. Then, a FWHM median filter is applied to reject bad/hot pixels and smooth out the noise having spatial period below the resolution limit. The image is finally divided, at each radius, by the standard deviation of the noise (again obtained with the IDL robust_sigma algorithm) to obtain a signal-to-noise ratio image.
The algorithm is first tested using data obtained at the Gemini telescope with the Altair adaptive optics system (Saddlemyer et al, 1998) and the NIRI near-infrared camera (Hodapp et al, 2000). These data are part of the Gemini Deep Planet Survey (Lafreniere et al, 2007a) that uses the angular differential imaging (ADI) technique (Marois et al, 2006; Lafreniere et al, 2007b) to detect faint companions. This technique consists in acquiring a sequence of images with continuous FOV rotation. A reference PSF that does not contain any point sources is first obtained by combining images of the sequence, and the quasi-static speckle noise is then attenuated by subtracting the reference ADI PSF. The data for the star HD97334B (program GN-2005A-Q16), acquired on April 18, 2005 with good seeing conditions (Strehl of at H-band), are presented. These data have been reduced, registered and processed using the pipeline described in Marois et al (2006) with the additional steps mentioned above, i.e. pixel masking and noise filtering & normalizing. Given that NIRI images are pixels and PSFs have 3 pixels per FWHM, we have chosen the same areas as the simulated PSFs mentioned above (see § 4.2) to calculate the PDF. For each region of the image, the pixel intensity histogram is obtained and then integrated to derive the CL curve. The CL is then extrapolated using a polynomial fit and the CL detection threshold is estimated. The derived detection thresholds for a CL are presented (see Fig. 10) for a single PSF, a PSF minus the ADI reference PSF, and the combined ADI-subtracted images.
The algorithm is next tested using observational data obtained at the Canada-France-Hawaii telescope using the PUEO adaptive optics system (Rigaut et al, 1998) and TRIDENT near-infrared camera (Marois et al, 2005). TRIDENT is a triple beam multi-wavelength (1.58, 1.625 and 1.68 with 1% bandpass) imager built following the simultaneous spectral differential imaging technique (Racine et al, 1999; Marois et al, 2000). The technique consists in acquiring several images at different wavelengths and subtracting them to attenuate the speckle noise while retaining most of the flux of nearby companions. These data have been acquired as part of a direct imaging survey of stars confirmed to possess exoplanets from radial velocity analysis. The dataset of the star Ups And, acquired on November 14, 2002, is used here. Seeing conditions were relatively good and the Strehl ratio was of the order of . The data have been reduced by subtracting a dark, dividing by a flat field, correcting for bad/hot pixels, and registering the PSF at the image center. An optimized reference PSF was obtained using a star (Chi And) having a similar spectral type, magnitude and acquired at the same DEC and HA to minimize PSFs evolution from differential atmospheric refraction and flexure effects. Performances with and without the simultaneous reference PSF and the Chi And reference PSF subtractions are analyzed. Due to the limited FOV of TRIDENT, sections of 10,000 pixels (500 resolution elements, given that TRIDENT has 5 pixels per ) are used to derive the PDF. Detection thresholds are thus known to 15% (see Tab. 2).
It is clear that both TRIDENT and Gemini raw PSFs are limited by a non-Gaussian noise. For both TRIDENT/CFHT and NIRI/Gemini images, even after subtraction of a reference PSF, the residuals are still dominated by quasi-static speckles rather than averaged atmospheric speckle, background, read, or photon noises. Only the final combined ADI-subtracted image possessed a clear Gaussian-like noise. Fig. 11 shows a visual example of a 5 detection with and without a Gaussian distributed noise after introducing artificial 5 point sources. For the Gaussian distributed noise, only the artificial point sources are detected with a detection threshold at 5,
6.1 PDF Evolution with Quasi-Static Speckle Averaging
It was shown in § 5 that the ADI technique produces a quasi-Gaussian noise. This is mainly due to the FOV rotation that occurs during the observing sequence; the residual noise is averaged incoherently when combining the images after FOV alignment. From the central limit theorem, it is thus expected that the noise in the final combined image shows quasi-Gaussian statistics.
In this section, simulations are presented to estimate the number of independent speckle noise realizations required to converge to a quasi-Gaussian noise intensity distribution. Random noise images having 10 resolution elements are created following an MR distribution having equal to 0.1, 1 and 10. The PDF and CL curves are calculated for a single realization up to the coaddition of 25 independent realizations (see Fig. 12). Typically, 20 independent realizations are required for the MR distribution to converge, to 20%, to a Gaussian distribution. Fig. 13 shows the detection threshold for a CL as a function of , the number of independent noise realizations
where is the number of acquired images in the sequence, is the integration time per image and the speckle noise decorrelation timescale (the equation is valid if ; if then ). These three curves can be well fit by a simple power-law of the form
where is the detection threshold of a single image for a CL. This equation is valid for all types of statistical distributions studied here. Eq. 6 can be used to predict the detection threshold required for a CL and a statistical distribution with a known instantaneous PDF and speckle noise decorrelation timescale. If we consider the Gemini ADI observation ( for a single ADI-subtracted image, see Fig. 10), for a 70 minute observation sequence with minutes (at 2 or 50, see Marois et al, 2006), it is expected that the detection threshold for a CL of the final combined image () will be at 50, in good agreement with the number derived with real images (, see Fig. 10).
6.2 Targets with a High Background Star Density
In some cases, it is desirable to observe an interesting nearby target situated along the galactic plane or in front of the galactic bulge. In those areas, the high background stellar density implies that numerous background stars will be present in each of the areas defined to derive pixel intensity distributions. The detection threshold obtained will be affected by those stars since the algorithm would consider them as speckles. Several techniques can be used to remove the stars before estimating detection thresholds. A simple solution is to subtract these stars, using a non-saturated image of the primary, and mask any remaining contaminated areas with a not-a-number (NaN) mask. If the observations are obtained using the ADI technique then a star-free residual image can be obtained in the following manner. Prior to combining all the ADI-subtracted images together, instead of rotating them by the angle required to align their field of view, they are rotated by the negative of that angle such that all off-axis sources (the background stars) are eliminated by the median combination. As the amplitude of the rotation between the images is the same as for the “proper” combination, the effect of the median combination on the residual noise is expected to be the same and this star-free residual image can be used to estimate the noise distribution. Of course, the proper ADI combination of images must be used to search for companions. Another approach, still within the ADI framework, is to use the final ADI residual image to subtract the off-axis sources from each non-rotated ADI-subtracted image. Then these source-free images are rotated by the negative of the angle needed to align their field of view, such that their median combination eliminates the off-axis sources subtraction residuals. As for the previous technique, this source-free residual image should have the same residual noise distribution as the proper ADI residual image.
A robust technique was elaborated to estimate sensitivity limits using a CL approach. This technique correctly finds the expected MR intensity distributions of simulated and real PSFs, and properly detects a change of PDF as a function of angular separation. Experiments with simulated and observational data confirm the prediction of the theory that raw PSFs obtained with high-contrast imaging instruments are limited by a non-Gaussian noise. A correction factor (up to 3) needs to be applied to detection limits found assuming Gaussian statistics to obtain the desired CL detection threshold. Properly estimating this effect is important for future high-contrast imaging instruments for both ground- and space-based dedicated missions since a loss of a factor of three in contrast results in less sensitivity to low-mass exoplanets or, if a specific contrast needs to be achieved, integration times need to be at least nine times longer. It was shown that the ADI technique is the only observing strategy currently known that generates, intrinsically, a quasi-Gaussian noise at all separations where sufficient FOV rotation has occurred. A simulation has shown that it takes typically independent speckle noise realizations to produce an average speckle noise that shows quasi-Gaussian statistics. A general power-law is derived to predict the detection threshold required when averaging independent speckle noise realization of known PDFs and decorrelation timescale.
|1||Define an image region|
|2||Subtract the mean intensity|
|3||Divide the pixel intensities by the noise RMS value|
|4||Obtain the pixel intensity histogram (PDF)|
|5||Integrate PDF between d to obtain the corresponding CL curve|
|6||Perform a polynomial fit on the CL curve|
|7||Extrapolate the CL curve to derive the CL detection threshold|
- affiliation: Based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council of Canada, the Institut National des Sciences de l’Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii. Based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the Particle Physics and Astronomy Research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), CNPq (Brazil) and CONICET (Argentina).
- affiliation: Institute of Geophysics and Planetary Physics L-413,
Lawrence Livermore National Laboratory, 7000 East Ave, Livermore, CA 94550
- affiliation: Département de physique and Observatoire du Mont Mégantic, Université de Montréal, C.P. 6128, Succ. Centre-Ville,
Montréal, QC, Canada H3C 3J7
- affiliation: Département de physique and Observatoire du Mont Mégantic, Université de Montréal, C.P. 6128, Succ. Centre-Ville,
Montréal, QC, Canada H3C 3J7
- The robust_sigma algorithm uses the median absolute deviation as a first estimate of the standard deviation and then weight points using Tukey’s Biweight; this algorithm provides another step of robustness to avoid biasing the standard deviation estimate if bad pixels are present.
- Only approximately half of the artificial point sources are detected since the artificial sources, being in intensity, vary in S/N by RMS due to the underlying noise in the image. A detection threshold thus misses/detects % of sources.
- Carson, J. C., Eikenberry, S. S., Smith, J. J., Cordes, J. M. 2006, AJ, 132, 1146
- Dohlen, K. et al. 2006, Proc. SPIE, 6269, 93
- Fitzgerald, M. P., Graham, J. R. 2006, 637, 541
- Ftaclas, C., Martín, E. L., Toomey, D. 2003, Proc. SPIE, 211, 521
- Goodman 1968, Introduction to Fourier Optics (New York: McGraw Hill)
- Hodapp, K., Hora, J., Graves, E., Irwin, E. M., Yamada, H., Douglass, J. W., Young, T. T., Robertson, L. 2000, Proc. SPIE, 4008, 1334
- Hinkley, S., Oppenheimer, B. R., Soummer, R., Sivaramakrishnan, A., Roberts, L. C. Jr., Kuhn, J., Makidon, R. B., Perrin, M. D., Lloyd, J. P., Kratter, K., Brenner, D. 2007, ApJ, 654, 633
- Kasper, M., Apai, D., Janson, M., Brandner, W. 2007, A&A, 472, 321
- Lafrenière et al. 2007a, ApJ, accepted, astroph/0705.4290v2
- Lafrenière, D., Marois, C., Doyon, R., Nadeau, D., Artigau, É 2007b, ApJ, 660,770
- Masciadri, E., Mundt, R., Henning, Th., Alvarez, C., Barrado y Navascus, D. 2005, ApJ, 625, 1004
- Macintosh, B., Poyneer, L., Sivaramakrishnan, A., Marois, C. 2005, Proc. SPIE, 5903, 170
- Macintosh et al. 2006, Proc. SPIE, 6272, 20
- Marois, C. 2004, PhD. thesis, Université de Montréal
- Marois, C., Doyon, R., Racine, R. & Nadeau, D. 2000 PASP, 112, 91
- Marois, C., Doyon, R., Racine, R., Nadeau, D. 2003, Astronomy with High Contrast Imaging: From Planetary Systems to Active Galactic Nuclei (eds. C. Aime and S. Soummer), 8, 233
- Marois, C., Doyon, R., Racine, R., Nadeau, D., Riopel, M., Vallée, P., Lafrenière, D. 2005, PASP, 117,745
- Marois, C., Lafrenière, D., Doyon, R., Macintosh, B., Nadeau, D. 2006, ApJ, 641, 556
- Metchev, S. A. 2006, Ph.D dissertation, California Institute of Technology, United States
- Potter, D., Baudoz, P., Guyon, O., Brandner, W., Close, L., Graves, J. E., Northcott, M. 2001, AAS, 33, 812
- Racine, R., Walker, G. A. H., Nadeau, D., Doyon, R., Marois, C. 1999, PASP, 111, 587
- Rigaut, F., Salmon, D., Arsenault, R., Thomas, J., Lai, O., Rouan, D., Véran, J. P., Gigan, P., Crampton, D., Fletcher, J. M., Stilburn, J., Boyer, C., Jagourel, P. 1998, PASP, 110, 152
- Roddier, F. 1981, Prog. Optics, 19, 281
- Saddlemyer, L. K., Herriot, G., Veran, J.-P., Fletcher, J. M. 1998, Proc. SPIE, 3353, 150
- Schneider, G., Silverstone, M. D. 2003, Proc. SPIE, 4860, 1
- Soummer, R., Aime, C. 2004, Proc. SPIE, 5490, Advancements in Adaptive Optics. Edited by Domenico B. Calia, Brent L. Ellerbroek, and Roberto Ragazzoni, 495