Selfcoherent camera as a focal plane wavefront sensor: simulations
R. Galicher^{1,2}  P. Baudoz^{1,2}  G. Rousset^{1,2}  J. Totems^{2,3}  M. Mas^{1,2}
1  LESIA, Observatoire de Paris, CNRS, University
Pierre et Marie Curie Paris 6 and University Denis Diderot Paris 7, 5
place Jules Janssen, 92195 Meudon, France
2  Groupement d'Intérêt Scientifique Partenariat Haute
Résolution Angulaire Sol Espace (PHASE) between ONERA,
Observatoire de Paris, CNRS and University Denis Diderot Paris 7, France
3 
Onera / DOTA  Chemin de la Hunière  91761 Palaiseau Cedex, France
Received 16 July 2009 / Accepted 22 September 2009
Abstract
Context. Direct detection of exoplanets requires high
dynamic range imaging. Coronagraphs could be the solution, but their
performance in space is limited by wavefront errors (manufacturing
errors on optics, temperature variations, etc.), which create
quasistatic stellar speckles in the final image.
Aims. Several solutions have been suggested for tackling this
speckle noise. Differential imaging techniques substract a reference
image to the coronagraphic residue in a postprocessing imaging. Other
techniques attempt to actively correct wavefront errors using a
deformable mirror. In that case, wavefront aberrations have to be
measured in the science image to extremely high accuracy.
Methods. We propose the selfcoherent camera sequentially used
as a focalplane wavefront sensor for active correction and
differential imaging. For both uses, stellar speckles are spatially
encoded in the science image so that differential aberrations are
strongly minimized. The encoding is based on the principle of light
incoherence between the hosting star and its environment.
Results. In this paper, we first discuss one intrinsic
limitation of deformable mirrors. Then, several parameters of the
selfcoherent camera are studied in detail. We also propose an easy and
robust design to associate the selfcoherent camera with a coronagraph
that uses a Lyot stop. Finally, we discuss the case of the association
with a fourquadrant phase mask and numerically demonstrate that such a
device enables detection of Earthlike planets under realistic
conditions.
Conclusions. The parametric study of the technique lets us
believe it can be implemented quite easily in future instruments
dedicated to direct imaging of exoplanets.
Key words: instrumentation: high angular resolution  instrumentation: interferometers  instrumentation: adaptive optics  techniques: image processing  techniques: high angular resolution
1 Introduction
Exoplanets are typically 10^{7} to 10^{10} fainter than their host and are often located within a fraction of an arcsecond from their star. Numerous coronagraphs have been proposed to reduce the overwhelming light of a star to obtain a direct imaging of extrasolar planets (Guyon et al. 2005; Rouan et al. 2000; Mawet et al. 2005). Several of them provide observations (Boccaletti et al. 2004; Schneider et al. 1998a). But performance is limited by wavefront errors in the upstream beam for all these coronagraphs and the final focal plane image shows stellar speckles. The effect of most of these aberrations can be corrected by adaptive optics (AO) or eXtreme AO (XAO, Vérinaud et al. 2008) but the uncorrected part generates quasistatic speckles, which limit the image contrast (Macintosh et al. 2005; Cavarroc et al. 2006). To reduce this speckle noise, differential imaging techniques attempt to subtract a reference image of the stellar speckles from the science image (star plus companion).
Several ways are used to measure that reference: spectral characteristics (Marois et al. 2000,2004; Racine et al. 1999), polarization states (Baba & Murakami 2003; Stam et al. 2004), differential rotation in image (Marois et al. 2006; Schneider et al 1998b), or incoherence between stellar and companion lights (Guyon 2004). The selfcoherent camera with which we work uses the last property (Baudoz et al. 2006; Galicher & Baudoz 2007). But before using one of these a posteriori techniques, we may actively correct quasistatic wavefront errors so that a first speckle reduction is achieved and differential imaging techniques have less to do. Because of the low level of aberrations that must be achieved (a few nanometers), the wavefront sensor of that loop has to measure for phase and amplitude errors in the final science image to avoid differential errors introduced by classical wavefront sensor (ShackHartmann for example, Shack & Platt 1971). Codona & Angel (2004) suggest using the incoherence between companion and stellar lights and use a modified MachZender interferometer to encoded the stellar speckles.
The instrument we propose, the self coherent camera (SCC), is based on the same property and uses Fizeau interferences. We insist on how SCC can be used both as a wavefront sensor for active correction (called step A in this paper) and as a differential imaging technique (step B). In Galicher et al. (2008), we numerically demonstrated that, applying step B after step A, a selfcoherent camera associated with a deformable mirror and a perfect coronagraph detects earths (contrast of ) from space in a few hours under realistic assumptions (zodiacal light, photon noise, readout noise, phase errors of 20 nm rms, 20% bandwidth and effective bandwidth, 8 m telescope with a 25% throughput and nm). Here, we report results from a parametric study of the SCC, and propose an easy and robust design to associate it with coronagraphs that use a Lyot stop. We explain the performance in the case of a fourquadrant phase mask coronagraph (FQPM, Rouan et al. 2000). Section 2 recalls the principle of the technique and presents the estimators of the pupil complex amplitude (phase and amplitude errors, step A) and companion images (step B). Section 3 provides the assumptions and criterions used for the parametric studies. Section 4 is a general study (no SCC) of one intrinsic limitation for deformable mirrors and presents the best contrast they can provide. The signaltonoise ratios on both SCC estimators (wavefront and companion) are developed in Sect. 5. Section 6 estimates the required stability of the reference beam. We report the effects of amplitude aberrations on the SCC performance in Sect. 7. The last two sections are the most important in the paper. The first one is the study of the chromatism on the SCC performance when a perfect coronagraph is used (Sect. 8). In the second one (Sect. 9), we propose a device to associate the SCC and any coronagraph using a Lyot stop. We discuss the case of a FQPM coronagraph: earths are detected using such a coronagraph with the SCC just by adding a small hole to the Lyot stop.
Figure 1: Selfcoherent camera principle schematics. 

Open with DEXTER 
2 Selfcoherent camera principle
The goal of the selfcoherent camera is the measurement (step A) of phase and amplitude aberrations in the pupil upstream the coronagraph and the speckle field estimation (step B) without introducing any noncommon path errors. To do so, we use spatial interferences in the science image to encode the stellar speckles that are directly linked to the wavefront aberrations. This section briefly recalls how to use the selfcoherent camera to measure wavefront errors and also to reduce the speckle noise in the image.
More details can be found in Galicher et al. (2008); Baudoz et al. (2006). As in these previous papers, we consider hereafter only space observations. Figure 1 presents the instrument schematics. A deformable mirror, located in a plane conjugated to the entrance pupil, reflects the beam coming in from the telescope. We then split the beam. In the image channel, the beam goes through a coronagraph. In the reference channel, we suppress all the companion light and extract a beam containing only light from the host star (cf. Sect. 9). Finally, we recombine the two beams in a Fizeau scheme to obtain spatial fringes in the science image on the detector. Phase and amplitude aberrations give a coronagraphic residue in the last focal plane. And the reference channel induces spatial interferences on these stellar speckles (cf. Fig. 2), whereas it does not have any impact on a possible companion image since companion light is not coherent with star light. The stellar speckles are thus encoded (modulated), whereas the companion image is not.
Figure 2: Science image of the selfcoherent camera. Stellar speckles of the coronagraphic residue are fringed, hence spatially encoded. No companion is present. 

Open with DEXTER 
 Step A, wavefront estimation and correction: estimate phase and amplitude errors from the focal plane image and correct for them using a deformable mirror (correction loop without any differential errors, Sect. 2.2).
 Step B, companion detection: record the science image when the best correction is achieved and postprocess that image to overcome the DM limitation (Sect. 2.3).
2.1 Notations
Before giving the expression of the wavefront estimator, we summarize the notations of this paper in Table 1.
Table 1: Table of the notations.
Their significance is given in the text when necessary.2.2 Step A: correction loop
2.2.1 Estimation of the aberrated wavefront
In Galicher et al. (2008), we proposed an estimator of phase errors in the pupil upstream of the coronagraph. Here, we give the SCC estimator of both phase and amplitude aberrations, in other words, the pupil complex amplitude.
That amplitude is extracted from the science image I (Fig. 2). Using notations from Table 1, we can write the interferential image I as
where the wavelength belongs to . The first three terms are the coronagraphic residue (stellar speckles to be reduced), the reference image, and the companion image (to be detected). They constitute the unmodulated part of the interferential image. The last term traduces the spatial speckle modulation: interferences with a interfringe. The required pupil separation (Fig. 1) not to overlap the peaks in the spatial frequency plane (inverse Fourier transform of the image I) is at least:
We use a margin (factor 1.1).
We follow Galicher et al. (2008) to estimate wavefront errors, and we extract the modulated part
of
because it contains a linear combination of
(what we look for) and
.
To do it, we select one of the two lateral peaks of the inverse Fourier transform of I and apply a Fourier transform to it:
We then assume a small effective bandwidth (cf. Sect. 8) and estimate the pupil complex amplitude downstream the coronagraph:
where is the inverse Fourier transform, u the cooordinate in the correlation plane of I, and F is applied to partially compensate for chromatic effects (see Sect. 8.2). Assuming small phase () and amplitude (a) aberrations and a perfect coronagraph, we find the complex amplitude downstream of the coronagraph (see Eqs. (6) to (8) in Galicher et al. 2008, for details):
(5) 
where is the star flux upstream of the coronagraph. From the two last equations, we deduce
We can thus estimate wavefront errors upstream of the coronagraph from I_{}, extracted from the science image. The other terms can be estimated as follows. The chromatic factor is computed, and its expression is given in Sect. 8.2. The constant is uniform (Table 2). We can precisely estimate it from flux calibration of the space telescope. Finally, the reference amplitude in the focal plane has to be recorded before the beginning of the loop and has to be stable. We point out that such a stabilization can be achieved (Sect. 6).
Once wavefront aberrations are estimated, we correct for them with a deformable mirror as explained in the next section.
2.2.2 Correction of the aberrated wavefront
In the whole paper, we assume that we use only one deformable mirror with actuators to correct for the wavefront errors. If only phase aberrations exist, it is possible to clean all the centered area in the science image. If both phase and amplitude aberrations are present, we can only clean half of that region, and we adopt the method proposed by Bordé & Traub (2006) to make the phase error screen hermitian. As the estimation is not perfect (small aberration linearization, noises, reference division) and the speckles are not static but quasistatic, we need to iterate a few times.
To achieve highcontrast imaging, we need to measure wavefront aberrations with high accuracy and to drive the deformable mirror faster than the aberration evolution time. As the required time to work out the SCC estimator is relatively short (three fast Fourier transform), the correction speed is determined by the integration time needed to reach a reasonable signaltonoise ratio for the wavefront measure. There is then a compromise between the achievable image contrast and the pointed stellar flux. That compromise is roughly the same regardless of the speckle calibration technique. Recent industrial studies for a spacebased 1.5 m telescope estimate that stabilities can be expected as low as 1 pm per hour (Guyon 2009).
2.3 Step B: companion estimation
Once the best correction is achieved (last iteration of
step A), we apply the postproccessing algorithm described in Baudoz et al. (2006) and Galicher & Baudoz (2007) to suppress the largest part of the speckle residue. We estimate the companion image
using
where I_{} is defined by Eq. (3) and is the unmodulated part of the science image I (Eq. (1)), which can be written as
To extract it, we select the central peak of the inverse Fourier transform of the science image (Eq. (1)) and apply a Fourier transform to that selection. Finally, is the reference image. It has to be estimated as explained in the introduction of Sect. 6.
3 Assumptions and criteria
This section introduces the assumptions of our numerical studies and the criteria for optimizing all the parameters of the instrument.
3.1 Assumptions
In the whole paper, we assume spatial observations (no dynamic aberrations) and achromatic coronagraphs without defects (perfect coronagraph or FQPM). The power spectral density (PSD) of static phase errors in the instrument upstream of the coronagraph varies as f^{3}, where f is the spatial frequency, which corresponds to typical VLT mirror aberrations (Bordé & Traub 2006). Amplitude aberration PSD is flat or evolves as f^{3}.
Table 2: Simulation assumptions for each section.
The reference complex amplitude in the detector plane is the Fourier transform of: the pupil upstream of the coronagraph (same wavefront errors) densified to obtain a pupil diameter instead of D, if a perfect coronagraph is used. The flux is set to verify the condition of Sect. 5.1 in the corrected area.
 a pupil of diameter extracted from the Lyot stop plane (phase and amplitude errors depend on the coronagraph) if a FQPM coronagraph is used (see Sect. 9).
To simulate the science image I (cf. Eq. (1)) in polychromatic light, we sum 5 monochromatic images corresponding to wavelengths uniformly distributed in the considered bandpass . Shannon sampling is imposed for fringes , and follows Eq. (2). When photon noise is simulated, we consider a G2 star at 10 pc observed by a space telescope with of throughput and a bandwidth of at nm. We indicate the exposure time and the telescope diameter case by case. Table 2 regroups the other assumptions, which change from one section to another. Only one parameter generally takes realistic values so that its impact can be determined independently to the other parameters, set to ideal values. At the end of Sect. 9, all parameters are set to realistic values to determine the SCCFQPM performance.
3.2 Criteria
To optimize the different parameters of the selfcoherent camera, we define two criteria.
3.2.1 Averaged contrast
The first criterion, called C_{1}, gives the averaged contrast achieved in the corrected area of the SCC science image. This area varies as a function of three parameters : deformable mirror size, polychromatism and amplitude aberrations. Considering an deformable mirror:
 in monochromatic light and without amplitude aberrations, is the centered region because we reduce the corrected area by shrinking it by a factor 1.05 (Sect. 3.1) and do not account for borders (factor 0.9) of that area where the detection is prevented by the diffracted light of the uncorrected speckles.
 in polychromatic light and without amplitude aberrations, following results of Sect. 8.1, we reduce to its intersection with the line of width and of same direction as fringes.
where I_{0} is the maximum intensity of the star image without coronagraph or SCC, and I is given by Eq. (1). The smaller C_{1}, the better the correction and the fainter companions can be detected. Since C_{1} is an average over the corrected area, it does not give any information on preferential positions in the image for the detection.
As explained in Sect. 2.2.2, several iterations are required to reach very highcontrast imaging. That is why, when C_{1} is considered, we study its convergence speed as a function of the number of correction steps.
3.2.2 5 detection
The second quantity we use to compare the different configurations is the 5
detection
,
which we define as
(10) 
where is the azimuthal standard deviation of the considered image at radial separation . As for C_{1}, the smaller , the better the correction and the fainter companions can be detected. The metric is a function of the angular separation in the considered image. It is then possible to determine preferential angular separations for the detection. However, does not give any information on favored directions. If is used, we plot it against the angular separation.
4 Intrinsic deformable mirror limitation
To derive SCC performance, we needed to know the intrinsic limitation of deformable mirrors that we use. In that section, we set out that intrinsic limitation (no SCC) if no apodization is used. We considered a single telescope associated with a perfect coronagraph and a deformable mirror. All the assumptions are recalled in Sect. 3.1. We simulated the coronagraphic image, without SCC, with a full DM correction (perfect estimation of the phase errors). We ploted against the angular separation in Fig. 3 for three different DM sizes and the same phase screen. The residue after the perfect coronagraph (no DM) is overplotted.
Figure 3: 5detection versus angular separation in the image provided by a perfect coronagraph associated to a single telescope in monochromatic light without DM (solid black line) or with DM: , and from top to bottom. The initial phase error rms amplitude is 20 nm with a power spectral density in f^{3}. Residual phase error levels are given for each DM in nm. 

Open with DEXTER 
As expected, the cutoff frequency equals , and higher frequencies are not corrected: no change for with or without correction further than . For , the detection limit depends on the diffracted energy of uncorrected speckles located outside the corrected area ( ). The larger , the larger and the lower the level of uncorrected speckles. That is why, for a fixed , the correction is better with a DM than with a .
To improve the detection limit we can increase the number of actuators (degrees of freedom of the correction). However, nowadays, manufacturing issues limit that number to for MicroElectroMechanical Systems (MEMS). For the EELT planet finder, EPICS, a device is under study but has not been developed yet. Another way to improve the detection limit is to minimize wavefront errors but it would be very optimistic to assume less than 20 nm rms for telescopes larger than a few meters. Another solution to rule out the detection limit is apodize the pupil so that the uncorrected speckles would spread less light into the corrected area. Give'on et al. (2006) assume a Kaiser function to apodize the pupil and achieve a 10^{10} contrast (1 detection) at in the focal plane. But such an apodization profile is are not easy to produce. With an amplitude mask (Kasdin et al. 2005; Vanderbei et al 2003), its throughput would be about , which significantly increases the exposure time for detecting faint companions. Another apodization technique, called phaseinduced amplitude apodization (Guyon et al. 2005), could be used because its throughput is about . But manufacturing problems limit its performance, even if solutions are under study (Pluzhnik et al. 2006).
Finally, if no apodization is used, the best DM correction shown in Fig. 3
 cannot be overcome by any device used to estimate wavefront errors because it is an intrinsic limitation depending on the diffraction of uncorrected speckles;
 comes from a speckle residue so that the sole way to improve it is a postcorrection of the image by differential imaging techniques. One example is the selfcoherent camera (step B, Sect. 2.3).
5 Signaltonoise ratio
We study the impact of the reference flux on the signaltonoise ratio of wavefront error estimation (step A, Sect. 5.1) and companion detection (step B, Sect. 5.2).
5.1 Step A: wavefront error measurement
The purpose of step A is to measure wavefront errors directly from the noisy science image I. The interesting quantity is then and, here, we work out the signaltonoise ratio of that measure considering the photon noise.
To estimate
,
we use Eq. (3) and obtain in monochromatic light
where the n index refers to the noisy quantities and where I_{} is the Fourier transform of one of the two lateral correlation peaks of the science image I. Assuming is almost uniform in the image (see Sect. 6), Eq. (11) implies
where we have replaced by the intensity of the reference channel on the detector . Considering photon noise, variance of can be approximated by (see Eq. (42) of Appendix A):
where is the number of pixels per resolution element in the science image. From Eqs. (12) and (13), if the measure is done with a signaltonoise ratio higher than , we can write
(14) 
Supplanting I_{} by (Eq. (11)) and by , we find:
(15) 
We want to estimate for wavefront errors (directly linked to ) and we know the reference intensity . In the case of an association with a FQPM (Sect. 9.3.1), before beginning the correction loop. Under that condition, . Using the Shannon fringe sampling ( ) and the minimum separation between pupils of image and reference channels (Eq. (2)), the signaltonoise ratio at the beginning of the loop roughly equals . After a few correction steps, we hope to achieve . In that case, . Finally, wavefront estimation (step A) can be done with a good accuracy as soon as is greater than a few photons per pixel. The sole problem would be a zero value of the reference flux because the speckles would not be encoded and would be zero (see Sect. 6).
5.2 Step B: companion detection
Step B corresponds to the companion estimation from the science image I of the last iteration of the correction (end of step A). It would be quite exhausting to study the exact propagation of noises through all steps of the algorithm that we use to compute the estimator of Eq. (7). We may present such a study in a future paper but it is not the purpose of the present one. Here, we express the dependence of the signaltonoise ratio of the estimator on fluxes in the image channel ( , speckle noise) and in the reference channel ( ).
We first assume uncorrelated noises for the unmodulated (
)
and modulated (
)
parts of the science image I (not rigorously exact). Considering the recorded reference intensity
is almost uniform on the detector (see Sect. 6), we can write the variance of the estimator
(Eq. (7)):
We restrict the covariance of to its variance (see Eq. (40) of Appendix A) and we suppose the distribution of the noise on still follows a Poissonian law. We can then write the fourth momentum of :
Finally, using Eqs. (42) and (43) of Appendix A (variances of and ) and Eqs. (16) and (17), we obtain an approximation of the noise on :
where I roughly equals to . We insist that this expression is not exact, but we did numerically check that it roughly approximates the real noise. We can then derive the ideal flux regime ( and ) to optimize the speckle substraction when applying estimator of Eq. (7). To do it, we assume the companion intensity is low compared to the coronagraphic residue ( ). If this is not the case, the companion is detected and the SCC not useful.
Weak fluxes
If
photon per pixel, the variance of the noise on
reduces to
(19) 
and the signaltonoise ratio of the companion detection is
(20) 
which is not interesting.
Dominating coronagraphic residue
If
,
and
,
Eq. (18) becomes
(21) 
If , noise on the companion estimator is greater than the initial speckle noise and the SCC is useless. If , speckle noise is lightly reduced and the SCC acts positively.
Dominating reference image
If
and
,
we simplify Eq. (18) as
(22) 
The referenceimage photon noise dominates the estimator noise. As it is greater than the speckle residue that we want to reduce ( ), this case is not interesting again.
Strong fluxes for both channels
If
and
photon per pixel, the variance of the noise on the companion estimator
(Eq. (18)) reduces to
(23) 
which means that it is proportional to the photon noise of the coronagraphic residue (of which the variance is ): the speckle noise is suppressed. The signaltonoise ratio of the companion detection is then
(24) 
We have assumed to verify the Shannon sampling of fringes ( ) and the minimum separation between image and reference pupils (Eq. (2)). The factor 3 reveals that the noise is smoothed and averaged when we select the correlation peaks and work out and by fast Fourier transform.
Finally, the best way to estimate the companion image from Eq. (7) with a good signaltonoise ratio is to set strong fluxes in both image and reference channels ( and larger than 1 photon per pixel on the detector).
For both wavefront estimation (step A) and companion estimation (step B), the range of accepted reference fluxes is broad and we propose to work with around a few tens of photons per pixel in the science image.
6 Reference beam
Like all differential imaging techniques, SCC needs the recording of a reference whose the complex amplitude is used to estimate wavefront errors (Eq. (6)). A first problem (zero division) appears at the positions where is zero because the condition of Sect. 5.1 is not verified and the corresponding speckles are not encoded well. The second problem is that it is impossible to simultaneously measure and the science image. We then propose to record the reference complex amplitude before beginning the loop and we call it . But in that case we need that recording to be stable in time. In Sect. 6.1, we present a way to answer both problems at the same time. Then, we study the impact of errors on : spatial drift 6.1 and optical path difference variation (Sect. 6.2) between the recording of and the beginning of the loop.
It is important to distinguish the reference amplitude used to simulate the SCC science fringed image (see Sect. 3.1) and the reference amplitude used in the estimator of Eq. (6). The second one is an estimation of the a priori unknown first one. In this paper, we use the diffracted complex amplitude by a diameter pupil free from any aberrations for . Flux and position of are set using a recorded image of the reference channel before the beginning of the correction loop.
6.1 Reference pupil diameter and spatial drifts of the image
The solution proposed to stabilize the reference channel and avoid the zero divisions is to use a small diameter for the reference pupil. In that way, , roughly equal to the central part of the Airy complex amplitude, has very low sensitivity to wavefront variations in the reference channel, and its first dark ring is pushed away from the center of the image. At this point, we would like to set the smallest that we could. However, in a real setup, the light entering the reference channel is extracted from the coronagraph rejected light, and its flux has a finite value. A minimum value for is then required to verify the condition of Sect. 5.1 and a tradeoff has to be derived for each setup. In the case of an association with a FQPM coronagraph (Rouan et al. 2000), we find that a reasonable value of is between 10 to 30 in function of wavefront error level (Sect. 9.3), and we often use in the paper.
We now quantify the impact of errors on the recorded reference amplitude . We first consider tiptilt variations in the reference channel that induce spatial drifts of the reference image on the detector between and . With and under the assumptions mentioned in Sect. 3.1, we find that for a (respectively 6) drift of the reference image, the correction is effective and roughly converges to the DM limitation (twice the DM limitation). The specification is then not critical, and in Sect. 9, we propose a compact, robust, and very simple device where the reference is stable enough in time.
6.2 Optical path difference
Another potential error on the estimated reference amplitude is the variation of the optical path difference (OPD) between the image and the reference channels. If the OPD varies, the center of the fringe pattern is shifted on the detector and, if we do not account for it, the wavefront error estimation is degraded. To quantify that limitation, we consider that the recorded reference corresponds to a zero OPD and we simulate correction loops for several nonzero OPD under assumptions of Sect. 3.1. Figure 4 gives the convergence of the averaged contrast C_{1} in the corrected area.
Figure 4: Averaged contrast C_{1} in the corrected area against correction steps for different optical path differences (OPD) between reference and image channels. We assume a deformable mirror. 

Open with DEXTER 
7 Amplitude aberrations
In that section, we look at the impact of amplitude errors in the pupil plane upstream of the coronagraph. We simulated the correction loop for several amplitude aberration levels ( to rms). With a sole DM, amplitude aberrations induce a reduction of the corrected area by a factor 2 and speckles of the uncorrected halfarea diffract their light into the corrected halfarea as seen in Fig. 5  amplitude aberrations and other assumptions detailed in Table 2.
Figure 5: SCC Science image after 4 steps of correction. Phase errors are set to 20 nm rms and amplitude aberrations to rms before correction. 

Open with DEXTER 
Figure 6: Detection at against the angular separation in the image at the 4th iteration of the correcion loop. We recall the coronagraphic residue without any correction (top black solid line). 

Open with DEXTER 
8 Chromatism impact
The point spread function (PSF) size and the interfringe are both proportional to wavelength (cf. Eq. (1)). In white light (small ), the science image is the superposition of all the monochromatic images over the considered bandwidth and fringes are blurred as seen in Fig. 7.
Figure 7: Science images in monochromatic ( left) and polychromatic ( , right) lights. These images correspond to the iteration 0 (no wavefront correction). 

Open with DEXTER 
8.1 Field limitation
We call
and
the largest and the smallest wavelengths of the spectral bandwidth .
The mean interfringe is
(see Eq. (2)). We assume the white fringe  null OPD  is in the center of the image. We derive the distance
where the fringe systems for
and
are shifted by half an interfringe:
(25) 
with . We deduce and for :
In the perpendicular fringe direction, fringes become blurred as from . In the fringe direction, fringes are not blurred since OPD between the two channels is the same for all wavelengths of the bandwidth. In that way, speckles of the science image are correctly encoded in the fringe direction whereas they are not in the perpendicular fringe direction as from . This effect is visible in the Fig. 8 where we show the 10th iteration of the correction for a deformable mirror and two bandwidths. For these images, we do not apply the chromatic factor (Eq. (28)) in the wavefront estimator (Eq. (6)). Even if fringes are not blurred in the fringe direction, the speckles are dispersed (PSF size is wavelength dependent), and wavefront errors are not estimated perfectly. Hereafter, we study the combined impacts of both effects (fringe and speckle dispersions).
Figure 8: Science image for the 10th iteration of the correction for two bandwidths: ( left) and 100 ( right). We assume a deformable mirror. The fringe direction is from bottom left to top right. 

Open with DEXTER 
(27) 
where is the separation between the reference and the image pupils (cf. Table 1). If we consider (Eq. (2)), we find , 75, and 150, respectively, for a , , and deformable mirror.
8.2 Chromatic factor
The wavefront error estimator given by Eq. (6) requires knowledge of the chromatic factor F that we optimize in the current section. In Galicher et al. (2008), we used
derived from the model of light propagation through SCC:
However, (close to cardinal sin function) has zero values in the corrected area as soon as ( ). Speckles located at the zeros of are not encoded well and can be neither estimated nor corrected. This limitation is similar to the zero values of the reference amplitude (Sect. 6.1). Using the estimator of Eq. (6), pupil spatial frequencies corresponding to separations higher than (blurred speckles) cannot be recovered. But we can restrain their impact on the corrected area, and we test five chromatic factors:
F_{1}  =  
F_{2}  =  (29)  
F_{3}  =  
F_{4}  =  1  
F_{5}  = 
We attempt to correct for unblurred speckles (Sect. 8.1) and minimize the impact of the uncorrected speckles (blurred ones). To compare the different chromatic factors, we consider assumptions given in Table 2 and look at the evolution of the averaged contrast in the corrected area C_{1} during the correction for . We find that setting the factor F to provides the best results (no divergence). That factor reduces the impact of the uncorrected speckles because it multiplies their intensity by a low value in the estimator of Eq. (6). It acts as a regularization of the uncorrected spatial frequencies. We use that chromatic factor in all the paper from now. A more sophisticated estimator (regularization of a minimization without linearization of wavefront errors and without assumptions on ) could certainly be developed but has not been studied yet.
8.3 Correction level
The convergence speed of the correction loop slightly decreases when the light becomes more and more chromatic: 3 steps in monochromatic light and 10 steps for for a . At the same time, the averaged correction gets worst as seen in Fig. 9 where the criterion C_{1} at the 10th loop iteration is plotted against the spectral bandwidth ( to ) under the assumptions of Table 2.
Figure 9: Averaged contrast C_{1} in the corrected area at the 10th iteration of the correction versus the spectral resolution for three deformable mirrors. 

Open with DEXTER 
From these results, it is clear that chromatism strongly limits SCC performance and that the selfcoherent camera cannot be used with a classical bandwidth . Working with a narrow bandwidth is not the solution because of photon noise. We foresee two possibilities. We can develop software solutions and modify the wavefront estimator of Eq. (6) to account for polychromatic dispersions of speckles and fringes (regularization of a minimization without linearization of aberrations and without assumptions on ). Hardware solutions are also conceivable. We could associate the selfcoherent camera with an integral field spectrometer (IFS) at modest resolution ( ). We could estimate wavefront errors for each wavelength channel. We could also develop a new algorithm to process the data of all the channels at the same time to optimize the estimation. This solution will be studied in future work. Another hardware solution is a Wynne compensator that we describe in Sect. 8.4.
8.4 Wynne compensator
The image widths (about and ) and the fringe period ( ) are proportional to wavelength, which is the reason for the SCC chromatism limitation (Sects. 8.1 to 8.3). In the context of speckle interferometry (Labeyrie 1970), Wynne (1979) proposed a device to correct for such a spectral dependence over a wide spectral range: the Wynne compensator (Fig. 10).
Figure 10: Scheme of a Wynne Compensator composed by two triplet lenses made with two kinds of glasses (called 1 and 2). Dispersion is linear with the wavelength and the outgoing beam is collimated. 

Open with DEXTER 
Such a Wynne compensator can be associated to the SCC by adding it just before the recombining optic. Since magnification is proportional to the wavelength, D, , and are proportional to and we obtain nonblurred fringes all over the detector (corresponding to the monochromatic case). A simulation is shown in Fig. 11 for a spectral resolution of at 650 nm. Since only magnification is important, we do not consider any coronagraph and the incoming wavefront is assumed to be aberrationfree.
Figure 11: Science image without ( left) or with (right) Wynne compensator. The initial spectral resolution is 6.5 at 650 nm. For the corrected image, the effective bandwidth is about ( right). We use log scale. 

Open with DEXTER 
Finally, even if chromatism seems to be a hard point of the selfcoherent camera technique, several solutions are conceivable: more sophisticated estimators, hardware solutions (association with a Wynne compensator or an integral field spectrometer at modest resolution).
9 SCC and real coronagraphs
The previous sections present a parametric study of the selfcoherent camera without mentioning any concrete setup and assuming a perfect coronagraph as in Galicher et al. (2008). In Baudoz et al. (2006) and Galicher & Baudoz (2007), we proposed a selfcoherent camera device built as an interferometer: a beamsplitter to create the reference channel, a pinhole to filter the reference beam, a delay line to ensure a null optical path difference, and a lens to recombine image and reference beams. The disadvantage of that device is the delay line that has to be controlled with very high accuracy in real time (Sect. 6.2). In Sects. 9.1 and 9.2, we propose a new robust design for associating the selfcoherent camera with a coronagraph that has a Lyot stop plane. We describe in detail the case of a FQPM (Rouan et al. 2000) in Sect. 9.3.
9.1 How to use the Lyot stop plane
Using a coronagraph that needs a Lyot stop is very interesting because it can be easily associated with the SCC. Such a coronagraph rejects only the stellar light outside the pupil  no companion light  so that the light stopped by the Lyot stop comes only from the hosting star and can be used to create the reference channel. We propose adding a small noncentered pupil to the classical Lyot stop as shown in the schematics of Fig. 12 for a FQPM coronagraph.
Figure 12: Schematics of the selfcoherent camera associated to a coronagraph that uses a Lyot stop plane. The pupil intensity distribution in the Lyot stop plane is given for a FQPM coronagraph. 

Open with DEXTER 
9.2 Wavefront estimation
The complex amplitude estimated by Eq. (4) is the pupil amplitude in the Lyot stop plane downstream the coronagraph. The second step is determine the complex amplitude upstream of the coronagraph where the deformable mirror stands. If the coronagraph is perfect and wavefront aberrations are small, we can use the estimator of Eq. (6). If the coronagraph is not perfect, we have to model the light propagation through it.
Calling M the mask function in the focal plane and L the classical Lyot stop (sole image channel), we have
where * denotes the convolution. We deduce within the L area:
where (1/M)_{0} is the inverse of the mask function M where M is not zero and equals 0 elsewhere. This expresses that we cannot estimate the spatial frequencies for which the mask M has stopped the energy (ie. M=0). We finally estimate the complex amplitude upstream of the coronagraph from Eqs. (4) and (31). We notice that phase masks are not bounded by this limitation since they do not block light.
9.3 SCC and fourquadrant phase mask
The FQPM coronagraph uses a Lyot stop and can be associated with the selfcoherent camera in a device, called SCCFQPM. The FQPM has been described in detail in previous papers (Riaud et al. 2001,2003; Rouan et al. 2000; Boccaletti et al. 2004). We recall the coronagraphic mask M induces a phase shift on two quadrants of a diagonal and no phase shift in the two other quadrants. In Fig. 12, we explain where the image and reference pupils are picked in the Lyot stop of an SCCFQPM. In Sect. 9.3.1, we first give the quantity of energy that goes through the reference pupil to determine the optimized ratio of pupil diameters (cf. Sect. 6.1). We then study the optical path difference between reference and image channels (cf. Sect. 6.2) in Sect. 9.3.1. Finally, we examine the performance of the SCCFQPM under realistic conditions in Sect. 9.3.2. In the whole section, we assume a perfect achromatic FQPM: infinitely thin transitions and perfect phase shifts.9.3.1 Implementation
Reference flux
As explained in Sect. 6.1, the ratio between image and reference pupil diameters
has to be optimized for each SCC design. In Table 3, we give the averaged intensities of the image and reference PSFs in the corrected area of a
DM for different
and phaseerror levels. All fluxes are expressed in photons per pixel.
To establish the table, we consider the device of Fig. 12 working under assumptions of Table 2. We put the reference pupil behind one of the quadrants of the FQPM (see Sect. 9.3.1), and we set the stellar flux using assumptions of Sect. 3.1, a 20 s exposure time, and 4 m diameter telescope.
Table 3: Averaged fluxes of image ( ) and reference ( ) channels in the corrected area of a DM for different (first raw) and phase error levels (second raw).
Coronagraphic residue in the image channel depends roughly quadratically on the phase error level as expected because the coronagraphic pupil is approximately . In contrast, the rejected star energy remains slightly constant in the reference pupil ( ). We notice also that roughly evolves as the inverse square of . We propose setting to 20 so that the condition of a few photons per detector pixel coming in from the reference channel (see Sect. 5.1) can be respected with exposure times that are not too long (speckle lifetime). It is important to notice that this choice is not strict. To strictly set , we have to know the lifetime of the speckles to be corrected. This would impose a maximum exposure time to record the interferential images I (aberrations have to be static to be estimated from I). From that maximum time, we could choose the apropriate value to let the reference flux verify the condition of Sect. 5.1. To make that choice, we may notice again that higher values of are interesting to ensure stability of the reference image and to limit the presence of low values of in the corrected area (cf. Sect. 6.1).
Optical path difference
We have shown in Sect. 6.2 that we need to accurately control the OPD between reference and image channels (accuracy and stability of
). In the case of an SCCFQPM, if the reference pupil is picked behind one of the
phase shift quadrants of the focal mask, the optical path difference in the Lyot stop plane in monochromatic light is .
Fringes are shifted by half an interfringe, and we have to multiply by 1 the estimator of Eq. (31).
If the reference pupil is picked behind another quadrant, the OPD
is null and we can use directly the estimator of Eq. (31). The directions to avoid are the transition direction of the FQPM for which the OPD is not just 0 or
and above all is not uniform over the reference pupil so that it is not
easy to account for it in the estimator. In a real experiment, we
propose to pick the reference pupil behind one of the four quadrants
avoiding the transition directions and to calibrate the sign of the
estimator. We are confident that the OPD is stable in time in the
compact device proposed in Fig. (12) and laboratory demonstrations are in progress to confirm it.
9.3.2 Performances
Impact of chromatism
If we directly apply the wavefront estimator of Eq. (31), the correction loop becomes slower than in the perfect coronagraph case (Sect. 8.3)
because the FQPM model is not perfect. We then employ a more
complex model accounting for the exact FQPM impact on the
first 999 Zernike polynomials (Noll 1976). Using that model, we plot the averaged contrast C_{1} in the corrected area of a SCCFQPM versus the correction step for several spectral bandwidths (Fig. 13). Assumptions are given in Sect. 3.1 and Table 2.
Figure 13: Averaged contrast C_{1} in the corrected area of a SCCFQPM versus the number of correction iterations for several bandwidths . The dotanddash line and the dashed line represent respectively the best DM correction and the monochromatic case with a perfect coronagraph. 

Open with DEXTER 
SCCFQPM detections
In the last section, we check for the detection efficiency of SCCFQPM under realistic assumptions detailed in Sect. 3.1 and Table 2
for spatial observations. The variance of the whole phase error
is 20 nm rms, and we set the initial astigmatism defects
in the FQPM transition direction
to 1 nm rms (levels before correction). We account
for amplitude aberrations of rms. We simulated twelve
earths, located symmetrically with respect to the image center at 3, 6, 9, 12, 15, and
.
For the 4 m telescope that we consider (throughput of ),
these separations correspond to 0.12 (1.24),
0.25 (2.48), 0.37 (3.71), 0.50 (4.95), 0.62 (6.19),
and 0.74 (7.43) arcsec (AU). One
superearth is located at
(8.25 AU) and a 10^{9} Jupiterlike planet sits at
(15 AU).
All spectra are flat and the hosting star is a G2 star
at 10 pc. We closed the loop to correct for wavefront
errors (step A, Sect. 2.2)
during 7 h25 min, which corresponds to 6
iterations. After step A, we apply the companion estimator and
obtain the image presented in Fig. 14. We consider photon noise but no readout noise.
Figure 14: Detections of earths (white circles), superearth (blue circle) and Jupiter (red circle) with a SCCFQPM downstream a 4 m space telescope and a deformable mirror. The image is the result of a wavefront correction loop (step A, 7 h 25 min) and a companion estimation (step B). 

Open with DEXTER 
Table 4: Comparison between the contrasts and angular separations measured in the SCC image and the simulated values.
Finally, planets as faint as earths are detectable by SCCFQPM in a few hours from space under realistic assumptions.
10 Conclusion
In Sect. 4, we provided the intrinsic limitation for deformable mirrors controlled by the algorithm of Bordé & Traub (2006), under realistic yet optimistic assumptions (no dead actuators, continuous face sheet). It is important to keep in mind that this limitation does not depend on the technique used to estimate for wavefront errors since we assumed a perfect estimation. One way to improve the DM best contrast could be an apodization of the pupil so that the uncorrected speckles would diffract their light in a more restricted area. However, all the techniques proposed to apodize a pupil (Pluzhnik et al. 2006; Guyon et al. 2005; Kasdin et al. 2005; Vanderbei et al 2003) come with throughput problems or manufacturing limitations.In Sects. 5 to 8, we gave the results of the parametric study of a selfcoherent camera (SCC) used as a focal plane wavefront sensor and associated with a perfect coronagraph and a deformable mirror. Several points do not seem to be critical for the technique: reference flux (Sect. 5), error on the exact position of the reference image, and diameter of the reference pupil (Sect. 6.1). On the contrary, two points are more critical:
 optical path difference between the reference and the image channels. We have to know and control this OPD with an accuracy of about . If we associate the SCC with a coronagraph using a Lyot stop as described in Sect. 9, the hardwar optical path difference is always zero and we only have to control the end of the setup. Only common optics are used in this setup. We could put the device in a closed box to avoid differential air variations. We plan to check for the level of the variations in the optical path difference in the device of Fig. 12 in a laboratory experiment.
 chromatism. As shown in Sect. 8, chromatism is the most critical point in the technique. The main consequence is the reduction of the corrected area, in other words, the field of view of the image as shown in Fig. 8. The uncorrected speckles spread light and limit the contrast of the detection. We are studying software and hardware solutions to minimize that effect. For example, for the former we could develop a more sophisticated wavefront estimator using a minimization with regularization terms to account for the spectral dispersion of the speckles ( ) and of the fringes (). Hardware solutions are certainly more appropriate. We presented the Wynne compensator in Sect. 8.4. According to numerical simulations, it would enable working with a classical bandpass (15%, ) in the visible light with a DM. We will test such a Wynne compensator in a laboratory experiment. A second hardware solution to overcome the chromatism limitation could be the association of the SCC with an integral field spectrometer at modest spectral resolution ( ). This solution is very attractive because it would directly provide a companion spectra, but no work has been done on it yet.
Section 9 presented a very simple and robust design that associates the selfcoherent camera with any coronagraph which uses a Lyot stop. The most interesting point is that the optical path difference between the two channels is constant per construction. In Sect. 9.3, we studied in detail the case of the association of the SCC with a FQPM coronagraph (SCCFQPM) and we showed in Sect. 9.3.2 that the performance is very attractive and comparable to the case of a perfect coronagraph that is presented in Galicher et al. (2008). Detections of earths, superearth, and 10^{9} Jupiter under realistic assumptions are numerically demonstrated for an SCCFQPM in polychromatic light ( ) using a Wynne compensator (reducing the effective bandwidth to ) in 7 h 25 min with a 4 m space telescope.
The next steps are laboratory demonstrations of both SCC capabilities: focalplane wavefront estimation (step A, Sect. 2.2) and companion estimation by differential imaging (step B, Sect. 2.3). We will also attempt to overcome the poor estimate of the astigmatism in the FQPM transition direction. New algorithms have already been developed for using a DM interaction matrix including the impact of the whole instrument: DM, coronagraph, and SCC. It will be tested in a laboratory experiment very soon. We also study the SCC association with other coronagraphs like an annular groove phase mask (Mawet et al. 2005).
AcknowledgementsWe thank Rémi Soummer for private communications about his paper ``Fast computation of Lyotstyle coronagraph propagation'' (Soummer et al. 2007b), which was very useful for simulating the polychromatic images and the different values of the parameter.
Appendix A
In that appendix, we present how photon noise and readout noise
propagate through the numerical algorithms providing the wavefront
estimation (Eq. (6)) and the companion estimation (Sect. 5.2). We call
the noisy intensity of the recorded sience image:
where is the noiseless intensity at position of the interferential image (Eq. (1)) and the noise. For photon noise (Poissonian distribution), the mean of is zero and its variance at the position is . For a uniform readout noise (Gaussian distribution), the mean of is null and its variance is over the whole image. For both photon noise and readout noise, we assume the noise is not correlated between pixels located at different positions and (the spatial covariance of is zero):
where is the Kronecker symbol and E the statistical mean. The selfcoherent camera technique uses , the inverse Fourier transform of the science image . Let be the noise on at the position u in the Fourier plane, such that:
Using Eqs. (32) and (34), the monochromatic case at and the Fourier transform properties, we find that the mean of is zero and its spatial covariance is
From Eq. (33), we determine the covariance of in the case of photon noise
and in the case of read out noise:
The spatial covariance of the inverse Fourier transform of photon noise (Poissonian distribution) is not reduced to its variance whereas it is for readout noise (Gaussian distribution).
During the SCC image processing, to estimate
,
we isolate one of the two lateral peaks of
(Baudoz et al. 2006; Galicher & Baudoz 2007) using a circular binary mask of diameter
(see Sect. 6.1 for that approximation) and we apply a Fourier transform. The noiseless part of that Fourier transform is called
(Eq. (3)). The noisy part is
(38) 
with the noise at the position. The mean of is zero because . From Eq. (36), we find that in the case of a photon noise:
where represents the integral over a pupil of diameter D. Calling where J_{1} is the Bessel function of the first kind corresponding to a diameter circular aperture, we can write
We deduce from Eq. (40) the variance of in the case of photon noise:
Practically speaking, we record the interferential image I on a finite number of pixels and numerically process data. This meansthat the Fourier transform of Eq. (39) is a fast Fourier transform and the width of is pixels. Finally, the convolution in the expression of roughly corresponds to averaging the noise over pixels and
A similar result is found for the noise on the unmodulated part of the image (Eq. (8)):
The sole difference is the size of the selecting binary mask in the correlation plane.
References
 Aime, C. 2005, PASP, 117, 1012 [NASA ADS] [CrossRef] [Google Scholar]
 Aime, C., Soummer, R., & Ferrari, A. 2002, A&A, 389, 334 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baba, N., & Murakami, N. 2003, PASP, 115, 1363 [NASA ADS] [CrossRef] [Google Scholar]
 Baudoz, P., Boccaletti, A., Baudrand, J., et al. 2006, Proc. IAU Colloquium 200, 553 [Google Scholar]
 Boccaletti, A., Riaud, P., Baudoz, P., et al. 2004, PASP, 116, 1061 [NASA ADS] [CrossRef] [Google Scholar]
 Bordé, P., & Traub, W. A. 2006, ApJ, 638, 488 [NASA ADS] [CrossRef] [Google Scholar]
 Cavarroc, C., Boccaletti, A., Baudoz, P., Fusco, T., & Rouan, D. 2006, A&A, 447, 397 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Codona, J., & Angel, R. 2004, ApJ, 604, L117 [NASA ADS] [CrossRef] [Google Scholar]
 Galicher, R., & Baudoz, P. 2007, C. R. Physique 8, 333 [NASA ADS] [CrossRef] [Google Scholar]
 Galicher, R., Baudoz, P., & Rousset, G. 2008, A&A [Google Scholar]
 Give'On, A., Kasdin, N. J., Vanderbei, R. J., et al. 2006, JOSA A, 23, 1063 [NASA ADS] [CrossRef] [Google Scholar]
 Goodman 2000, Statistical Optics (Wiley Classics Library) [Google Scholar]
 Guyon, O. 2004, ApJ, 615, 562 [NASA ADS] [CrossRef] [Google Scholar]
 Guyon, O. 2009, Priv. Comm. [Google Scholar]
 Guyon, O., Pluzhnik, E. A. Galicher, R., et al. 2005, ApJ, 622, 744 [NASA ADS] [CrossRef] [Google Scholar]
 Huang, L., Rao, C., & Jiang, W. 2008, OSAJ, 16, 108 [Google Scholar]
 Kasdin, N. J., Vanderbei, R. J., Littman, M. G., et al. 2005, ApOpt, 44, 1117 [NASA ADS] [Google Scholar]
 Kuchner, M. J., & Traub, W. A. 2002, ApJ, 570, 900 [NASA ADS] [CrossRef] [Google Scholar]
 Labeyrie, A. 1970, A&A, 6, 85 [NASA ADS] [Google Scholar]
 Macintosh, B., Poyneer, L., Sivaramakrishnan, A., et al. 2005, SPIE Conf., 5903, 170 [NASA ADS] [Google Scholar]
 Marois, C., Doyon, R., Racine, R., et al. 2000, PASP, 112, 91 [NASA ADS] [CrossRef] [Google Scholar]
 Marois, C., Racine, R., Doyon, R., Lafrenière, D., & Nadeau, D. 2004, ApJ, 615, L61 [NASA ADS] [CrossRef] [Google Scholar]
 Marois, C., Lafrenière, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556 [NASA ADS] [CrossRef] [Google Scholar]
 Mawet, D., Riaud, P., Absil, O., et al. 2005, ApJ, 633, 1191 [NASA ADS] [CrossRef] [Google Scholar]
 Noll, R. J. 1976, JOSA, 66, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Pluzhnik, E. A., Guyon, O., Ridgway, S. T., et al. 2006, ApJ, 644, 1246 [NASA ADS] [CrossRef] [Google Scholar]
 Rabbia, Y., Gay, J., & Rivet, J.P. 2008, CRAS, 8, 385 [Google Scholar]
 Racine, R., Walker, G. A. H., Nadeau, D., Doyon, R., & Marois, C. 1999, PASP, 111, 587 [NASA ADS] [CrossRef] [Google Scholar]
 Riaud, P., Boccaletti, A., Baudrand, J., et al. 2001, PASP, 113, 1145 [NASA ADS] [CrossRef] [Google Scholar]
 Riaud, P., Boccaletti, A., Baudrand, J., et al. 2003, PASP, 115, 712 [NASA ADS] [CrossRef] [Google Scholar]
 Roddier, F., & Roddier, C. 1997, PASP, 109, 815 [NASA ADS] [CrossRef] [Google Scholar]
 Rouan, D., Riaud, P., Boccaletti, A., Clénet, Y., & Labeyrie, A. 2000, PASP, 112, 1479 [NASA ADS] [CrossRef] [Google Scholar]
 Schneider, G., Becklin, E. E., Smith, B. A., et al. 1998a, BASS, 30, 1363 [Google Scholar]
 Schneider, G., Thompson, R. I., Smith, B. A., et al. 1998b, SPIE Conf., 3356, 222 [NASA ADS] [CrossRef] [Google Scholar]
 Shack, R. B., & Platt, B. C. 1971, JOSA, 61, 656 [Google Scholar]
 Soummer, R., Ferrari, A., Aime, C. et al. 2007a, ApJ, 669, 642 [NASA ADS] [CrossRef] [Google Scholar]
 Soummer, R., Pueyo, L., Sivaramakrishnan, A., et al. 2007b, Optic Express, 15, 15935 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Stam, D. M., Hovenier, J. W., & Waters, L. B. F. M. 2004, A&A, 428, 663 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Trauger, J. T., & Traub, W. A. 2007, Nature, 446, 771 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Vérinaud, C., Korkiakoski, V., Yaitskova, N., et al. 2008, SPIE Conf., 7014 [Google Scholar]
 Vanderbei, R. J., Spergel, D. N., & Kasdin., N. J. 2003, ApJ, 599, 686 [NASA ADS] [CrossRef] [Google Scholar]
 Wynne, C. 1979, Opt. Comm., 28, 21 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: Table of the notations.
Table 2: Simulation assumptions for each section.
Table 3: Averaged fluxes of image ( ) and reference ( ) channels in the corrected area of a DM for different (first raw) and phase error levels (second raw).
Table 4: Comparison between the contrasts and angular separations measured in the SCC image and the simulated values.
All Figures
Figure 1: Selfcoherent camera principle schematics. 

Open with DEXTER  
In the text 
Figure 2: Science image of the selfcoherent camera. Stellar speckles of the coronagraphic residue are fringed, hence spatially encoded. No companion is present. 

Open with DEXTER  
In the text 
Figure 3: 5detection versus angular separation in the image provided by a perfect coronagraph associated to a single telescope in monochromatic light without DM (solid black line) or with DM: , and from top to bottom. The initial phase error rms amplitude is 20 nm with a power spectral density in f^{3}. Residual phase error levels are given for each DM in nm. 

Open with DEXTER  
In the text 
Figure 4: Averaged contrast C_{1} in the corrected area against correction steps for different optical path differences (OPD) between reference and image channels. We assume a deformable mirror. 

Open with DEXTER  
In the text 
Figure 5: SCC Science image after 4 steps of correction. Phase errors are set to 20 nm rms and amplitude aberrations to rms before correction. 

Open with DEXTER  
In the text 
Figure 6: Detection at against the angular separation in the image at the 4th iteration of the correcion loop. We recall the coronagraphic residue without any correction (top black solid line). 

Open with DEXTER  
In the text 
Figure 7: Science images in monochromatic ( left) and polychromatic ( , right) lights. These images correspond to the iteration 0 (no wavefront correction). 

Open with DEXTER  
In the text 
Figure 8: Science image for the 10th iteration of the correction for two bandwidths: ( left) and 100 ( right). We assume a deformable mirror. The fringe direction is from bottom left to top right. 

Open with DEXTER  
In the text 
Figure 9: Averaged contrast C_{1} in the corrected area at the 10th iteration of the correction versus the spectral resolution for three deformable mirrors. 

Open with DEXTER  
In the text 
Figure 10: Scheme of a Wynne Compensator composed by two triplet lenses made with two kinds of glasses (called 1 and 2). Dispersion is linear with the wavelength and the outgoing beam is collimated. 

Open with DEXTER  
In the text 
Figure 11: Science image without ( left) or with (right) Wynne compensator. The initial spectral resolution is 6.5 at 650 nm. For the corrected image, the effective bandwidth is about ( right). We use log scale. 

Open with DEXTER  
In the text 
Figure 12: Schematics of the selfcoherent camera associated to a coronagraph that uses a Lyot stop plane. The pupil intensity distribution in the Lyot stop plane is given for a FQPM coronagraph. 

Open with DEXTER  
In the text 
Figure 13: Averaged contrast C_{1} in the corrected area of a SCCFQPM versus the number of correction iterations for several bandwidths . The dotanddash line and the dashed line represent respectively the best DM correction and the monochromatic case with a perfect coronagraph. 

Open with DEXTER  
In the text 
Figure 14: Detections of earths (white circles), superearth (blue circle) and Jupiter (red circle) with a SCCFQPM downstream a 4 m space telescope and a deformable mirror. The image is the result of a wavefront correction loop (step A, 7 h 25 min) and a companion estimation (step B). 

Open with DEXTER  
In the text 
Copyright ESO 2010