Issue 
A&A
Volume 557, September 2013



Article Number  A9  
Number of page(s)  13  
Section  Astronomical instrumentation  
DOI  https://doi.org/10.1051/00046361/201321706  
Published online  13 August 2013 
Estimation and correction of wavefront aberrations using the selfcoherent camera: laboratory results
^{1} LESIA, Observatoire de Paris, CNRS, UPMC Paris 6 and Denis Diderot Paris 7, Meudon, France
email: johan.mazoyer@obspm.fr
^{2} Laboratoire d’Astrophysique de Marseille, CNRS, AixMarseille Univ., Marseille, France
Received: 15 April 2013
Accepted: 20 June 2013
Context. Direct imaging of exoplanets requires very high contrast levels, which are obtained using coronagraphs. But residual quasistatic aberrations create speckles in the focal plane downstream of the coronagraph which mask the planet. This problem appears in groundbased instruments as well as in spacebased telescopes.
Aims. An active correction of these wavefront errors using a deformable mirror upstream of the coronagraph is mandatory, but conventional adaptive optics are limited by differential path aberrations. Dedicated techniques have to be implemented to measure phase and amplitude errors directly in the science focal plane.
Methods. First, we propose a method for estimating phase and amplitude aberrations upstream of a coronagraph from the speckle complex field in the downstream focal plane. Then, we present the selfcoherent camera, which uses the coherence of light to spatially encode the focal plane speckles and retrieve the associated complex field. This enabled us to estimate and compensate in a closed loop for the aberrations upstream of the coronagraph. We conducted numerical simulations as well as laboratory tests using a fourquadrant phase mask and a 32 × 32 actuator deformable mirror.
Results. We demonstrated in the laboratory our capability to achieve a stable closed loop and compensated for phase and amplitude quasistatic aberrations. We determined the bestsuited parameter values to implement our technique. Contrasts better than 10^{6} between 2 and 12 λ/D and even 3 × 10^{7} (rms) between 7 and 11 λ/D were reached in the focal plane. It seems that the contrast level is mainly limited by amplitude defects created by the surface of the deformable mirror and by the dynamic of the detector.
Conclusions. These results are promising for a future application to a dedicated space mission for exoplanet characterization. A number of possible improvements have been identified.
Key words: instrumentation: high angular resolution / instrumentation: adaptive optics / techniques: high angular resolution
© ESO, 2013
1. Introduction
Direct imaging is crucial to increase our knowledge of extrasolar planetary systems. On the one hand, it can detect longorbit planets that are inaccessible for other methods (transits, radial velocities). On the other hand, it allows the full spectroscopic characterization of the surface and atmosphere of exoplanets. In a few favorable cases, direct imaging has already enabled the detection of exoplanets (Kalas et al. 2008; Lagrange et al. 2009) and even of planetary systems (Marois et al. 2008, 2010). However, the main difficulties of this method are the high contrast and small separation between the star and its planet. Indeed, a contrast level of 10^{10} has to be reached within a separation of ~0.1′′ or lower to allow the detection of rocky planets.
To reduce the star light in the focal plane of a telescope, several coronagraphs have been developed, such as the fourquadrant phase mask (FQPM) coronograph (Rouan et al. 2000), the vortex coronograph (Mawet et al. 2005) and the phaseinduced amplitude apodization coronograph (Guyon et al. 2005). However, the performance of these instruments is drastically limited by phase and amplitude errors. Indeed, these wavefront aberrations induce stellar speckles in the image, which are leaks of the star light in the focal plane downstream of the coronagraph. When classical adaptive optics (AO) systems correct for most of the dynamic wavefront errors that are caused by to atmosphere, they use a dedicated optical channel for the wavefront sensing. Thus, they cannot detect quasistatic noncommon path aberrations (NCPA) created in the differential optical paths by the instrument optics themselves. These NCPA have to be compensated for using dedicated techniques, for groundbased telescopes as well as for spacebased instruments.
Two strategies have been implemented to overcome the quasistatic speckle limitation. First, one can use differential imaging techniques to calibrate the speckle noise in the focal plane. Theses methods can use either the spectral signature and polarization state of the planet or differential rotation in the image (Marois et al. 2004, 2006). Second, even before applying these postprocessing techniques, an active suppression of speckles (Malbet et al. 1995) has to be implemented to reach very high contrasts. It uses a deformable mirror (DM) controlled by a specific wavefront sensor that is immune against NCPA. The techniques developed for this purpose include dedicated instrumental designs (Guyon et al. 2009; Wallace et al. 2010), or creating of known phases on the DM (Bordé & Traub 2006; Give’on et al. 2007) to estimate the complex speckle field.
Groundbased instruments that combine these two strategies are currently being developed, such as SPHERE (Beuzit et al. 2008) and GPI (Macintosh et al. 2008), to detect young Jupiterlike planets with an expected contrast performance of 10^{6} at 0.5′′. Better contrasts might be achieved to reach the rocky planet level with instruments using dedicated active correction techniques embedded in space telescopes (Trauger & Traub 2007).
In this context, we study a technique of wavefront sensing in the focal plane that allows an active correction in a closed loop. This paper has two main objectives. First we give an overview of how the amplitude and phase errors upstream of a coronagraph can be retrieved from the complex amplitude of the speckle field (Sect. 2) and how they can be compensated for using a DM (Sect. 3). In Sect. 4, we introduce the selfcoherent camera (SCC; Baudoz et al. 2006; Galicher et al. 2008). This instrument uses the coherence of the stellar light to generate Fizeau fringes in the focal plane and spatially encode the speckles. Using both the aberration estimator and the SCC, we are able to correct phase and amplitude aberrations. The second objective of the paper is a laboratory demonstration of the active correction and an experimental parametric study of the SCC (Sect. 5).
Fig. 1
Principle of a coronagraph (top). Aberrations in the entrance pupil plane induce speckles in the focal plane (bottom). 
2. Wavefront estimator in the focal plane of a coronagraph
In this section, we aim to prove that one can retrieve the wavefront upstream of the coronagraph using the measured complex amplitude of the electric field in the focal plane downstream of the coronagraph. We assume in the whole section that we can measure this complex amplitude without error using an undetermined method. We describe one type of this method (the SCC) in Sect. 4. In Sect. 2.1, we express the complex electric field that is associated to the speckles as a function of the wavefront errors in the pupil upstream of a phase mask coronagraph. From this expression, we propose an estimator of the wavefront errors from the speckle electric field (Sect. 2.2) and analyze its accuracy for an FQPM (Sect. 2.3).
2.1. Expression of the complex amplitude of speckles in the focal plane as a function of the initial wavefront
We consider here a model of a phase mask coronagraph using Fourier optics. Figure 1 (top) presents the principle of a coronagraph. We assume that the star is a spatially unresolved monochromatic source centered on the optical axis. The stellar light moves through the entrance pupil P. Behind this pupil, the beam is focused on the mask M in the focal plane, which diffracts the light. Hence, the non aberrated part of the stellar light is rejected outside of the imaged pupil in the next pupil plane and is stopped by the Lyot stop diaphragm L. The aberrated part of the beam goes through the Lyot stop, producing speckles on the detector in the final focal plane (Fig. 1, bottom).
We note whith α and φ the amplitude and phase aberrations in the entrance pupil plane and define the complex wavefront aberrations Φ as (1)The complex amplitude of the star in the entrance pupil plane can be written as (2)where ψ_{0} is the mean amplitude of the field over the pupil P, ξ the coordinate in the pupil plane and λ the wavelength. We assume that the aberrations are small and defined in the pupil P (PΦ = Φ), thus (3)The complex amplitude of the electric field behind the coronagraphic mask M in the first focal plane is (4)where ℱ is the Fourier transform. Using Eq. (3), we can write the electric field before the Lyot stop (5)where ∗ is the convolution product. We call Φ_{M} the aberrated part of the field after the coronagraph: (6)After the Lyot stop L, the electric field ψ_{S} is (7)We assume a coronagraph for which the non aberrated part of the electric field is null inside the imaged pupil. This property of the perfect coronagraph (Cavarroc et al. 2006) has also been demonstrated analytically for several phase coronagraphs such as FQPM coronagraphs (Abe et al. 2003) and vortex coronagraphs (Mawet et al. 2005). The remaining part Φ_{L} of the normalized electric field after the Lyot stop reads (8)In the final focal plane, the complex amplitude A_{S} is (9)This complex amplitude is directly related to the wavefront aberrations in the entrance pupil. If one can measure A_{S}, we can invert Eq. (9) and retrieve the complex wavefront errors Φ in the entrance pupil.
2.2. Wavefront estimator
We still assume in this section that an undefined method provides access to A_{S}. Using this complex amplitude A_{S} as the measurement, we therefore propose the following estimator Φ_{est} for the wavefront: (10)This estimator can be used for any phase mask coronograph (for which M is nonzero over the full focal plane). To justify the pertinence of this estimator, we can rewrite it using the variables of our model. Using Eq. (9), in a noisefree measurement case, this estimator reads (11)Theoretically, if no Lyot stop is applied (L = 1), Eq. (11) becomes Φ_{est} = PΦ = Φ. We propose this estimator based on the assumption that most of the information about the aberrations is not diffracted outside of the imaged pupil by the coronagraphic mask. Therefore, using this assumption, we intuit that for L ≃ P, we still have Φ_{est} ≃ Φ. This assumption is verified by the simulation in Sect. 2.3, and by the experiment described in Sect. 5.
Fig. 2
Simulations of an aberrated phase in the entrance pupil plane (φ), and the real part of the field in the next pupil plane before (φ_{M}) or after (φ_{L}) the Lyot pupil. We also show the estimate (φ_{est}) and the difference between φ_{est} and φ. L = P in this case. 
For a symmetrical phase mask such as the FQPM, either ℱ^{1} [M] and are real. Thus, in the estimator we can separate the real (φ_{est}) and imaginary part (α_{est}) of the estimator in Eq. (11): (12)This relation ensures that within the limits of our model, this estimator independently provides estimates of the phase and amplitude aberrations.
2.3. Performance of the estimation
In this section we test the accuracy of the estimation φ_{est} for a phase aberration φ and no amplitude aberrations (α = 0). In the following numerical simulations, we assumed an FQPM coronagraph. It induces a phase shift of π in two quadrants with respect to the two others quadrants. We simulated FQPM coronagraphs in this paper using the method described in Mas et al. (2012). This coronagraph is completely insensitive to some aberrations, for instance to one of the astigmatism aberrations (Galicher 2009; Galicher et al. 2010). Because these aberrations introduce no aberration inside the Lyot pupil, we are unable to estimate them. We assumed an initial phase with aberrations of 30 nm root mean square (rms) over the pupil at λ = 635 nm, with a power spectral density (PSD) in f^{2}, where f is the spatial frequency.
In these simulations, we studied two cases. First, we used a Lyot pupil of the same diameter as the entrance pupil (L = P). Then, we studied the case of a reduced Lyot (D_{L} < D_{P}, where D_{L} and D_{P} are the diameters of the Lyot and entrance pupil, respectively).
2.3.1. Case L = P
Figure 2 shows the effect of phaseonly aberrations φ in different planes of the coronagraph. Starting from the left, we represent the initial phase φ, the real part of the amplitude due to the aberrations φ_{M} (φ_{L}) before the Lyot stop (after the Lyot stop), derived from Eq. (6) (Eq. (8)) for phaseonly aberrations. The last two images are the estimator φ_{est} and the difference between the estimate and the entrance phase aberrations (φ − φ_{est}).
The estimate φ_{est} is very close to the initial phase φ. For initial phase aberrations of 30 nm rms, the difference φ − φ_{est} presents a level of 10 nm rms in the entire pupil. The vertical and horizontal structures in this difference are due to the cutoff by the Lyot stop of the light diffracted by the FQPM (the light removed between φ_{M} and φ_{L}), which leads to an imperfect estimate of the defects on the pupil edges. Aberrations to which the FQPM coronagraph is not sensitive (such as astigmatism) are also present in this difference.
Assuming a perfect DM, we can directly subtract φ_{est} from φ in the entrance pupil. Then, we can estimate the residual error once again, and iterate the process. The aberrations in the Lyot pupil φ_{L} converge toward zero (0.2 nm in ten iterations). This is important because these aberrations are directly linked to the speckle intensity in the focal plane downstream of the coronagraph. However, the difference φ − φ_{est} does not converge toward zero in the entrance pupil. The fact that φ_{L} converges toward zero proves that the residual phase is only composed of aberrations unseen by the FQPM.
Fig. 3
Simulations of an aberrated phase with a localized default in the entrance pupil plane (φ), and in the next pupil plane (φ_{M}, φ_{L}). We represent the entrance pupil size by a dark ring around φ_{L}. We show the estimate (φ_{est}) and the difference between φ_{est} and φ in the last two images. D_{L} = 95%D_{P} in this case. 
2.3.2. Case D_{L}<D_{P}
In a more realistic case, we aim to remove all the light diffracted by the coronagraphic mask, even for unavoidable misalignments of the Lyot stop. For this reason, the Lyot stop is often chosen to be slightly smaller than the imaged pupil. We consider here a Lyot stop pupil L_{95%} whose diameter is D_{L} = 95%D_{P}. In a first part, we show below that phase defects at the edge of the entrance pupil can be partially retrieved, then we study the convergence of the estimator in this case.
As in Fig. 2, Fig. 3 corresponds to the simulation of the consecutive steps of the model (φ, φ_{M} before Lyot stops, φ_{L} after Lyot stop, then estimated phase φ_{est} and difference with initial phase). We added a small localized phase default, indicated by the black arrow, inside the entrance pupil P, but outside of the Lyot stop L_{95%} (Fig. 3, left). Around φ_{L}, the complex amplitude after the Lyot stop, we drew a circle corresponding to the entrance pupil, slightly larger than L_{95%}. For an FQPM, the additional defect is diffracted in the Lyot stop plane (φ_{M}). After applying the Lyot stop of 95% (φ_{L}), most of the default disappears, but we can still see its signature. As the estimator φ_{est} deconvolves by the phase mask, it partially retrieves the default, as seen in the estimate (indicated by the black arrow). In the error (φ − φ_{est}), we notice a remarkable cross issued from this defect, which is due to the information lost during the filtering by the Lyot stop.
The wavefront estimation is limited when compared to the D_{L} = D_{P} case. Because of the light filtered by the Lyot stop, some information about the wavefront aberrations close to the border of the entrance pupil is inevitably lost. Due to these unseen aberrations, φ − φ_{est} does not converge toward zero. However, the residual aberrations in the Lyot pupil φ_{L} still converge toward zero, practically as quick as in the L = P case.
By this first rough analysis, we see that we can estimate phase aberrations upstream of a coronagraph using the complex amplitude A_{S} of the speckle field and Eq. (10). The same conclusion can be drawn for amplitude aberrations and, because the estimation is linear, for a complex entrance wavefront. In the next section, we demonstrate that one can compensate for the wavefront errors in the entrance pupil.
3. Entrance pupil wavefront correction
In this section, we use the estimator Φ_{est} to numerically simulate the correction of phase and amplitude aberrations in a closed loop. In the loop, we can remove constant factors in the estimator, which can be adjusted with a gain g(13)We still assume that we have a perfect sensor that measures the complex amplitude A_{S} in the focal plane downstream of the coronagraph. We used a deformable mirror (DM) of N × N actuators upstream of the coronagraph, in the entrance pupil plane. We started with phaseonly correction. We explain how to correct for the effects of phase and amplitude aberrations with only one DM in Sect. 3.3.
We define the correction iterative loop by the expression of the residual phase φ^{j + 1} at iteration j + 1: (14)where is the shape of the DM at iteration j, is the incremental command of the DM actuator i at iteration j + 1, and f_{i} is the DM influence function, i.e., the WF deformation when only poking the actuator i. Note that φ^{j} is the phase to be estimated at iteration j + 1. The objective is now to determine the command vector from the phase estimator .
3.1. Wavefront aberration minimization
To derivate the command vector , we minimize the distance between the measurements and the measurements that accounts for the parameters to be estimated. Ideally, we would like to find the { k_{i} } minimizing the distance between the residual phase and the DM shape: (15)As presented in the previous section, a possible estimator of φ_{j} is given by Eq. (13), allowing us to compute from (directly related to φ^{j}). So we minimize (16)where D is the interaction matrix because we are using a linear model. This matrix is calibrated offline directly using the wavefront sensor (Boyer et al. 1990). As in the conventional leastsquares approach, we derived the pseudo inverse of D, denoted , by the singular value decomposition (SVD) method. Therefore, the command vector solution of Eq. (16) is given by (17)Equation (17) is applicable for different estimators (only changes). To create the interaction matrix, we poked one by one the actuators while the others remain flat, as shown in Fig. 4. Each estimated phase vector obtained hence gives the column of the interaction matrix corresponding to the moved actuator. The influence function (which we simulated as a Gaussian function) is at the left, the estimator given by Eq. (13) at the right. At the center, we also plot another estimator of the wavefront that does not include the deconvolution by the coronagraph mask. It is defined by (18)The chosen estimator applied to the influence function must be as spatially localized as possible: we have to filter the noise and it is far more efficient if the relevant information is gathered around one point. For this reason, it is preferable to use Φ_{est} (Fig. 4, right) instead of Φ_{est,2} (Fig. 4, center). We notice in Fig. 4 that even after deconvolution by the FQPM, the estimate (right) shows a negative cross centered on the poked actuator, whereas it is not present in the initial phase (left). This artifact is generated to the transitions of the quadrants, which diffract the light outside of the Lyot stop. This cross can be a problem for two reasons. First, it is difficult to properly retrieve it in a noisy image. Then, because it enhances the crosstalk between the actuator estimates, it may lead to unstable corrections (see Sect. 4.4).
Fig. 4
Simulations of the influence function f_{i} in the pupil plane (left), and the effect using the two different estimators: we deconvolve by the mask (Φ_{est}, right) or not (Φ_{est,2}, center). 
For these reasons, we chose to create a synthetic interaction matrix, i.e., use a slightly different model for the estimation.
3.2. Synthetic interaction matrix
One interest of interaction matrices is to calibrate the misregistration between the DM and the wavefront sensor when considering a complex optical system. It also allows us to calibrate the shape and magnitude of each actuator response. Since the cross in Fig. 4 (right) is 40 times less intense than the poke actuator in the center, it may be only partially retrieved in noisy images, which may lead to an unstable loop. To avoid this problem, we decided to build a synthetic interaction matrix based on the measured position and shape of each actuator. Because we considered an iterative measurement and correction loop, we finally discarded the magnitude calibration of each actuator, which lead to a slight increase of the required iteration number for convergence. We decided to only adjust the misregistration and estimation shape of the actuator set in the output pupil on the measured interaction matrix, as seen by the sensor.
Out of the N × N actuators of the square array in the DM, we chose to limit the response adjustment on the 12 × 12 actuators centered on the entrance pupil. These actuators were alternately pushed and pulled using known electrical voltages and we recorded A_{S} for both positions. Assuming the complex amplitude A_{S} is a linear function of the wavefront errors in the entrance pupil plane and that other aberrations of the optical path remain unchanged between two consecutive movements (the same pushed and pulled actuator), the difference between these two movements leads to the estimate of the influence function of a single poked actuator using Eq. (13). For each estimate, we adjusted a Gaussian function defined by its width and position in the output pupil. From these 144 Gaussian fits, we can build the actuator grid as observed in the plane, where aberrations are estimated and determine the interactuator distance in each direction and the orientation of this grid. We also determined the median width of the adjusted Gaussian functions and computed a synthetic Gaussian function, which was translated onto the adjusted actuator grid to create a new set of N × N synthetic estimates .
From these synthetic estimates, we built the synthetic interaction matrix D^{synth}. Some of the actuators are outside of the entrance pupil, and their impact inside the pupil is negligible. We excluded these actuators from D^{synth}. For any wavefront estimate, the distance to be minimized is now (19)The solution is given by the pseudoinverse of the interaction matrix D^{synth} using the SVD method.
3.3. Phase and amplitude correction using one DM
As explained in Bordé & Traub (2006), a complex wavefront Φ = φ + iα can be corrected for on half of the focal plane with only one DM. The idea is to apply a real phase on the DM to correct for the phase and amplitude on half of the focal plane. Because the Fourier transform of a Hermitian function is real, we define as (20)where is the complex conjugated of A_{S} and introduce it into the estimator of Eq. (13). The resulting estimated wavefront is real, which allows its correction with only one DM.
Now we have a solution to correct for the wavefront aberrations upstream of the coronagraph when the complex amplitude in focal plane is known. We introduce in Sect. 4 a technique to measure A_{S}: the selfcoherent camera.
4. Selfcoherent camera: a complex amplitude sensor in focal plane
The selfcoherent camera (SCC) is an instrument that allows complex electric field estimations in the focal plane.
4.1. SCC principle
Figure 5 (top) is a schematic representation of the SCC combined with a focal phase mask coronagraph and a DM. We added a small pupil R, called reference pupil, in the Lyot stop plane of a classical coronagraph. R selects part of the stellar light that is diffracted by the focal coronagraphic mask. The two beams are recombined in the focal plane, forming Fizeau fringes, which spatially modulate the speckles. In the following, we call the SCC image the image of the encoded speckles (Fig. 5, bottom). In this section, we briefly demonstrate that this spatial modulation allows us to retrieve the complex amplitude A_{S}. A more complete description of the instrument can be found in Baudoz et al. (2006), Galicher et al. (2008, 2010).
Fig. 5
Principle of the SCC combined with a coronagraph and a DM (top). A small hole is added in the Lyot stop plane to create a reference channel. In the final focal plane (bottom), the SCC image is formed by speckles encoded with Fizeau fringes. 
The electric field ψ in the pupil plane (Eq. (7)) after the modified Lyot stop is (21)where ξ_{0} is the separation between the two pupils in the Lyot stop, and δ is the Kronecker delta. ψ can also be written as (22)where ψ_{S} is the complex amplitude in the Lyot stop, defined in Eq. (8), and ψ_{R} is the complex amplitude in the reference pupil. We denote with A_{R} its Fourier transform, the complex amplitude in the focal plane, of the light issued from the reference pupil. In monochromatic light, the intensity I = ℱ [ψ] ^{2} recorded on the detector in the final focal plane can then be written as (23)where A^{∗} is the conjugate of A and x the coordinate in the focal plane. The two first terms are the intensities issued from Lyot and reference pupils, and provide access only to the square modulus of the complex amplitudes. The two correlation terms that create the fringes directly depend on A_{S} and A_{R}.
When an offaxis source (planet) is in the field of view, its light is not diffracted by the coronagraphic mask. Thus, it does not go through the reference pupil. Because the lights of the offaxis and inaxis sources are not coherent, the offaxis light amplitude in the focal plane does not appear in the correlation terms (i.e., its image is not fringed).
4.2. Complex amplitude of the speckle field
In this section, we demonstrate that we can use the SCC image to estimate the complex amplitude of the speckle field. We first apply a numerical inverse Fourier transform to the recorded SCC image (Eq. (23)), (24)where I_{S} =  A_{S}  ^{2} and I_{R} =  A_{R}  ^{2} are the intensities of the speckles and reference pupil, and u is the coordinate in the Fourier plane.
Fig. 6
Correlation peaks in the Fourier transform of the focal plane. The inverse Fourier transform of I_{S} + I_{R} is circled in blue. The inverse Fourier transform of is circled in red. 
ℱ^{1} [I] is composed of three peaks centered at u = [− ξ_{0}/λ,0, + ξ_{0}/λ] (Fig. 6). We denote with D_{L} the diameter of the Lyot pupil and with D_{L}/γ the diameter of the reference pupil (γ > 1). The central peak is the sum of the autocorrelation of the Lyot and reference pupils ℱ^{1} [I_{S} + I_{R}]. Its radius is D_{L} because we assume γ > 1. The lateral peaks of the correlation (ℱ^{1} [I_{−}] and ℱ^{1} [I_{+}] hereafter) have a radius (D_{L} + D_{L}/γ)/2. Thus the three peaks do not overlap only if (Galicher et al. 2010) (25)which puts a condition on the smallest pupil separation. The lateral peaks are conjugated and contain information only on the complex amplitude of the stellar speckles that are spatially modulated on the detector. When we shift one of these lateral peaks to the center of the correlation plane (u = 0), its expression can be derived from Eq. (24): (26)Assuming γ ≫ 1, we can consider that the complex amplitude in the reference pupil is uniform and that is the complex amplitude of an Airy pattern. Therefore, knowing A_{R}, we can to retrieve the complex amplitude A_{S} in the focal plane using the SCC (where is not zero): (27)
4.3. SCC wavefront estimation
Equation (13) shows how to estimate the wavefront upstream of a coronagraph is estimated using the complex amplitude of the speckle field in the focal plane. Combining Eqs. (10) and (27), we have an estimator of the wavefront aberrations Φ as a function of I_{−}: (28)This estimator is only limited in frequency by the size of the reference pupil. Indeed, where the reference flux is null, the speckles are not fringed and their estimate cannot be achieved. Small reference pupils produce large point spread functions (i.e., with a first dark ring at large separation) and allow estimating A_{S} in a large area of the focal plane. The influence of the reference pupil size is detailed in Sect. 5.5.
Figure 7 summarizes the steps followed to estimate the phase and amplitude aberrations with the SCC. From the fringed focal plane, we used a Fourier transform to retrieve I_{−}, from which we deduced the complex amplitude of the speckle field A_{S}. Using the estimator, we measured the phase and amplitude aberrations.
Fig. 7
Steps followed to estimate the phase and amplitude from SCC images. 
4.4. Correction loop
We can use this wavefront estimator to control a DM and correct for the speckle field in the focal plane as explained in Sect. 3. The DM has a finite number of degrees of freedom and thus can only correct for the focal plane in a limited zone. If the reference pupil is small enough (γ ≫ 1), the point spread function (PSF)  A_{R}  ^{2} is uniform over the correction zone (). We discuss this assumption in Sect. 5.5. Under this assumption, Eq. (28) becomes (29)As described in Sect. 3, we removed the constant terms in the estimation and put them into the gain g. From Φ_{est}, we created a synthetic matrix, as explained in Sect. 3.2. Similarly, the other estimator Φ_{est,2} introduced in Eq. (18), becomes (30)Using the interaction matrices deduced from these estimators (Φ_{est}, Φ_{est,2}) and the synthetic one, we studied the correction loop. We simulated a DM with 27 actuators across the entrance pupil. To build these matrices, we only selected the actuators with a high influence in the pupil (633 actuators were selected for this number of actuators in the pupil). Lyot stop and entrance pupil have the same radius, and we chose γ = 40 for the reference pupil size.
In Fig. 8, we plot the singular values (SV), normalized to their highest values, derived from the inversion of the matrices D obtained using the estimators Φ_{est} and Φ_{est,2} and of the interaction matrix built from . As already underlined, the cross in Φ_{est} or Φ_{est,2} (Fig. 4, center and right) correlates the estimates of different actuators and therefore leads to lower SV (up to five times lower for the lowest SV). When inverted in , low SV lead to higher values (in absolute values) and amplify the noise in Eq. (17). Applied to noisy data, such matrices may lead to an unstable correction. Even in a noisefree case, simulations of the correction with the three methods and the same number of actuators used (633) showed that only the synthetic matrix leads to a stable correction.
Fig. 8
Singular values, normalized to their respective maximum, issued from the inversion of the interaction matrices D, obtained using the two estimators Φ_{est} (red, dotted) and Φ_{est,2} (black, solid) and the synthetic matrix (blue, dashed) for γ = 40. 
4.5. Optical path difference
Between the Lyot stop and the detector, the beam is split into two paths (image and reference), which encounter different areas in the optics. Thus, differential aberrations exist Galicher et al. (2010). However, because the reference pupil is small (γ ≫ 1), the main aberration is an optical path difference (OPD) between the two channels. In this section we study how this OPD impacts the SCC performance.
4.5.1. Influence of an OPD on the correction
Given an OPD d_{op}, we can define phase difference φ_{op} = 2πd_{op}/λ. This phase difference modifies the I_{−} originally defined in Eq. (26): (31)The phase and amplitude estimate φ_{est,op} and α_{est,op} can be expressed as a function of the estimates made without an OPD (d_{op} = 0): (32)and thus (33)Hence, even phaseonly aberrations (such as the movements of the DM) have an influence on the estimated amplitude (i.e., the imaginary part of the estimator Φ_{est}) for a nonzero OPD. In this section, we make two assumptions. First, that the DM is perfect and we can correct for any desired phase in the pupil plane. Second, that the only error in the estimator is due to the OPD: if d_{op} = 0, the estimator retrieves the exact phase and amplitude (Φ_{est} = φ + iα).
We started the loop with a phase φ^{0} and an amplitude α. After j iterations the phase in the pupil plane φ^{j} is the difference between the previous phase φ^{j − 1}, and the estimate of this previous phase . Under the previous assumptions, we have α_{est} = α and φ_{est} = φ^{j − 1}, and (34)Because the OPD biases the estimation, the correction introduces an error at each iteration. This sequence converges if  1 − cos(φ_{op})  < 1. This assumption (− π/2 < φ_{op} < π/2) is always satisfied in real cases. Its limit φ_{end} satisfies the equation (35)Therefore, for a nonzero OPD and a phaseonly correction, the SCC correction converges, but the errors on the final phase depend on the uncorrected amplitude aberrations α.
To estimate the OPD effect on the level of the focal plane intensity, we considered the complex amplitude in the focal plane as a linear function of the phase and amplitude aberrations in the entrance pupil plane. We can thus evaluate the energy in the focal plane as a linear function of  φ  ^{2} +  α  ^{2}. Without an OPD, a perfect phaseonly correction would leave a level of speckles only dependent on the entrance amplitude aberrations  α  ^{2}. With an OPD, this level is slightly higher:  α  ^{2}(1 + tan(φ_{op})^{2}). For a realistic phase difference of 0.1 radians, the difference in intensity in the speckle field between the case with and without an OPD would be 1%. The impact on the correction is only weak.
The problem occurs when we try to correct phase and amplitude at the same time with one DM. We explained how to do this in Sect. 3.3. For φ_{op} ≠ 0, numerical simulations as well as tests on an optical bench show that the correction is unstable: at each iteration, we raised the phase aberrations by trying to correct for the amplitude aberrations and vice versa. Thus, we need an estimate of the OPD to stabilize the correction.
4.5.2. Estimation and correction of the OPD
In the construction of the synthetic matrix, (Sect. 3.2), we studied the difference of two SCC images produced by wavefronts that only differ by a movement of an actuator. Because the DM is in the pupil plane, the estimator applied to this difference is real for of an OPD equal to zero. However, for a nonzero OPD d_{op} and using Eq. (31) with α = 0, we deduce . For each of the 12 × 12 actuators used to build the synthetic matrix, the arctangent of the ratio of the imaginary part on the real part of leads to an estimate of the OPD. Due to the noise in the image, small differences in the OPD estimate can appear from one actuator to another. Calculating the median of the estimated OPDs, we obtain the measured phase difference . We modified I_{−} accounting for this OPD and our estimator (Eq. (29)) becomes (36)We use this new estimator from now on.
The OPD variations during the correction are a problem that has to be carefully considered for a telescope application. In the current installation (bench under a hood, room temperature stabilized) these variations are much slower than the time of a correction loop. Moreover, one can change the value of φ_{op} directly during the correction to compensate for slight changes. However, in an operational instrument, this problem will be taken into account by design to comply with the stability requirements (Macintosh et al. 2008).
5. Correction in the focal plane using the selfcoherent camera: laboratory performance
5.1. Laboratory test bench
We tested the SCC on a laboratory bench at the Observatoire de Paris. A thorough description of this optical bench is given in Mas et al. (2010). We briefly present the main components used in the experiments of the current paper:

1.
a quasimonochromatic laser diode emittingat 635 nm;

2.
a tiptilt mirror built at LESIA, used to center the beam on the coronagraphic mask (Mas et al. 2012). The tiptilt mirror can also be used in the closedloop as an offload for the DM;

3.
a Boston Micromachines DM of 32 × 32 actuators on a square array. Each actuator has a size of 300 μm. We currently use an entrance pupil of 8.1 mm and thus 27 actuators across the pupil;

4.
an FQPM optimized for 635 nm;

5.
a Lyot stop with a diameter of 8 mm for an entrance pupil of 8.1 mm (98.7% filtering) and in the same plane, reference pupils of variable diameters: 0.3 mm (γ = 26.6), 0.35 mm (γ = 22.8), 0.4 mm (γ = 20) and 0.5 mm (γ = 16) and 0.8 mm (γ = 10);

6.
a CCD camera of 400 × 400 pixels with a readout noise of 16 electrons/pixel and a full well capacity of 13 000 electrons/pixel.
We used the Labview software to control the bench and the DM and applied the closedloop correction at 20 Hz.
5.2. Dark holes
Fig. 9
Dark holes recorded on the laboratory bench for correction with two different sizes of square mask S_{q}: K_{Sq} = 20.8λ/D_{L} (left) and K_{Sq} = 24.5λ/D_{L} (center). The dark hole recorded on the laboratory bench for a correction in phase and amplitude with a square mask of size K_{Sq} = 24.5λ/D_{L} (right). These images use a different intensity scale but the same space scale. 
Owing to the limited number of actuators on the DM, only spatial frequencies lower than the DM cutoff can be corrected for. For a given diameter D_{L} of the Lyot pupil, the highest frequency attainable for a N × N actuators DM (N actuators across the pupil diameter) is Nλ/(2D_{L}) in one of the principal directions of the mirror and in the diagonal. The largest correction zone, called dark hole () in Malbet et al. (1995) is the zone in the image plane. During the numerical process of the SCC image (Fig. 7), we can decide to reduce the correction to a smaller zone than the one allowed by the number of actuators of the DM. This can be implemented in the SCC correction by multiplying I − by a square mask S_{q}. Modifying Eq. (36), the estimation becomes (37)where S_{q} equals 1 on a square area of K_{Sq}λ/D_{L} × K_{Sq}λ/D_{L} in the center of the image and 0 everywhere else.
Using an SCC with a reference pupil of 0.5 mm (γ = 16), we applied Eq. (37) to estimate the upstream wavefront. We used a square zone to restrain the correction zone to 24.5λ/D_{L} to optimize the correction of the DM. We built a synthetic interaction matrix as described in Sect. 3.2. The pseudo inverse of D^{synth} was used to control the DM in a closed loop using Eq. (14). The correction loop was closed at 20 Hz for the laboratory conditions and ran for a number of iterations large enough (j > 10) for the DM to converge to a stable shape. We recorded focal plane images during the control loop. The typical result obtained on the optical bench for this reference and square zone sizes and for phaseonly correction is shown in Fig. 9 (center). We also show an image of a obtained with a correction with a square zone of size K_{Sq} = 20.8λ/D_{L} (Fig. 9, left). A specific study of the size of the correction zone is made in Sect. 5.4. In Fig. 9, dark zones represent low intensities. The eight bright peaks at the edges are caused by high spatial frequencies due to the printthrough of the actuators on the DM surface. These peaks are uncorrectable by nature, but probably do not strongly alter the correction because they are situated at more than 20λ/D_{L} from the center.
As explained in Sect. 3.1, the correction of phase and amplitude with only one DM is possible by replacing A_{S} by in Eq. (13). With Eq. (20), we similarly define the hermitian function from I_{−}. Using Eq. (27) and the assumption that is an Airy pattern, a phase and amplitude correction is therefore possible by replacing I_{−} by in Eq. (37). This correction allows one to go deeper in contrast but limits the largest possible dark hole to half of the focal plane: . On this half plane, we can also choose to reduce the correction to a smaller zone. A resulting dark hole is presented in Fig. 9 (right) for K_{Sq} = 20.8λ/D_{L}.
5.3. SCC performance
In this section, we present contrast results obtained on the laboratory bench for phaseonly correction and for amplitude and phase correction. We used a reference pupil of 0.5 mm (γ = 16) to estimate the upstream wavefront and a square zone of size K_{Sq} = 24.5λ/D_{L} to optimize the correction of the DM.
5.3.1. Phaseonly correction
Fig. 10
Radial profiles of the azimuthal standard deviation (in rms) of the intensities in the focal plane typically obtained with this method for phaseonly correction, for simulation (blue dashed line) and laboratory bench result (red solid line), for γ = 16 and a square zone of size K_{Sq} = 24.5λ/D_{L}. We also plot in this graph the simulation of the focal plane obtained using the amplitude aberrations recorded and no phase aberrations (black dashdotted line). 
The speckles near the FQPM transitions are brighter than those in other parts of the . Moreover, the contrast in these region is not relevant, because the image of a planet located on a transition would be distorted and strongly attenuated. Therefore, for phaseonly correction, we chose to measure the radial profile of the SCC image only on the points (x,y) which verify (38)We calculated the profiles by normalizing the intensities by the highest value of the PSF measured through the Lyot pupil and without coronagraphic mask. In practice, we moved the source away from the coronagraph transitions to measure this PSF. In the following figures, the distances to the center are measured in λ/D_{L}. Figure 10 shows the radial profile of the azimuthal standard deviation of the intensities obtained in phaseonly correction in the focal plane zone described in Eq. (38). The detection level reaches a contrast level of 10^{6} between 6 and 12 λ/D_{L} and 3 × 10^{7} at 11 λ/D_{L}. As shown in Fig. 9 (center), speckles are still present in the dark area. Since we only corrected for the phase, we can suspect amplitude effects.
Fig. 11
Pupil illumination recorded on the laboratory bench. 
To estimate the amplitude aberration level, we recorded the pupil illumination on the optical bench without coronagraph, shown in Fig. 11. The amplitude defect level is estimated to be about 10% RMS in intensity. The period of the actuator pitch clearly appears in this pupil image. Due to vignetting effects by the focal coronagraphic mask, these highfrequency structures of the DM surface create illumination variations across the pupil. The first effect of these highfrequency aberrations are bright speckles outside the corrected zone (mostly on the eight bright peaks). The second effect is more critical for our purpose. Because the level of highfrequency amplitude errors varies across the pupil, it creates lowfrequency amplitude aberrations, which induce bright speckles in the center of the correction zone.
To compare the level of the recorded speckles with the one expected using amplitude and phase errors, we simulated the expected focal plane image. We used the amplitude aberrations deduced from the intensity measurement on the laboratory bench (Fig. 11). From these amplitude aberrations, we first simulated the the focal plane without phase errors (just amplitude errors). The profile of this focal plane is plotted in Fig. 10 with a black dotdashed line. We then simulated a phaseonly correction, assuming initial phase aberrations of 16 nm RMS over the pupil, and a power spectral density (PSD) in f^{2} where f is the spatial frequency. These simulation results (blue dashed line) are compared to the experimental measurement (red line) in Fig. 10. The level and shape of the two curves are very similar. They show the same structure around 27λ/D_{L}, due to the eight bright peaks created by amplitude aberrations. These curves inside the match the simulation of the focal plane without amplitude aberrations. It seems that in phaseonly correction, we corrected all phase aberrations and that we are only limited by amplitude errors.
5.3.2. Phase and amplitude correction
The simulation without amplitude errors (only phase aberrations) shows that a contrast level of 10^{10} can be reached, as previously shown in Galicher et al. (2010). Since the amplitude errors set the limits of our phaseonly corrections, we aim to correct both phase and amplitude at the same time. However, with only one DM, the corrected zone is smaller by half, as shown in Fig. 9 (right). Therefore, the radial profile measurement zone becomes (39)The results for this correction are plotted in Fig. 12 as a dashed blue line for the simulations and as a red line for the laboratory bench results. When correcting for the phase and amplitude aberrations, we obtain contrasts better than 10^{6} between 2λ/D_{L} and 12λ/D_{L}, and better than 3 × 10^{7} between 7λ/D_{L} and 11λ/D_{L}. This is an improvement compared to the phaseonly correction. The simulated profiles match the laboratory results from 0 to 8 λ/D_{L} and outside of the .
Between 8 and 12 λ/D_{L}, the experimental correction shows a plateau at 3 × 10^{7}, while the simulation correction goes deeper. This plateau is a distinctive feature of a limitation caused by the low dynamic range of the detector (our CCD camera has a full well capacity of 13 000 electrons/pixels for a readout noise of 16 electrons/pixels). This is confirmed by the last images of the loop which show speckle levels below the readout noise between 8 and 12 λ/D_{L}: the speckles beyond the readout noise are not visible and thus beyond correction. However, this problem can be solved by using a detector with a better dynamic range.
The number of incoming photons from the observed source is a critical problem of any specklecorrection technique: the speckles can only be corrected for to a certain level of contrast if the source is bright enough for them to be detected above photon and detector noise at these levels. Although we can correct in a closed loop at 20 Hz in the laboratory, the correction rate in a real telescope observation will be limited by the shortest exposure time necessary. This shortest exposure time depends on several parameters such as stellar magnitude, observational wavelengths, telescope diameter, or dynamic range of the camera.
The contrast level in the numerical simulation is limited to 10^{7}. This is due to the highamplitude defects (10% in intensity) introduced by the DM in the pupil. Indeed, the bright speckles of the uncorrected halfarea diffract their light into the corrected halfarea. This limit, independent of the estimation method (Give’on et al. 2006; Galicher et al. 2010), may be lowered by the introduction of a second DM on the optical bench (Pueyo et al. 2010).
Fig. 12
Radial profiles of the azimuthal standard deviation (in rms) of the intensities in the focal plane typically obtained with this method for phase and amplitude correction, for simulation (blue dashed line) and laboratory bench result (red solid line), for γ = 16 and a square zone of size K_{Sq} = 24.5λ/D_{L}. 
In the next sections (Sects. 5.4 and 5.5), we study the influence of different parameters on the SCC performance.
5.4. Size of the corrected zone
In this section, we compare the performance for different sizes of the square zone S_{q}. Using the modified estimator introduced in Eq. (37), and for different square zone sizes K_{Sq}, we experimentally closed the loop and recorded images after convergence. In these tests, we used N = 27 actuators across the pupil diameter and γ = 16, with phaseonly correction. As explained in Sect. 5.2, for this number of actuators, we have (as D_{L}/D_{P} = 8/8.1). We tested the case (K_{Sq} = ∞) and three others: K_{Sq} = 26.4λ/D_{L}, which is only slightly smaller than size of the largest and two smaller square zones (K_{Sq} = 20.8λ/D_{L} and K_{Sq} = 24.5λ/D_{L}). The images obtained in the last two cases can be seen in Fig. 9: K_{Sq} = 20.8λ/D_{L} (left) and K_{Sq} = 24.5λ/D_{L} (center).
Figure 13 presents the radial profiles of the focal planes obtained on the laboratory bench, normalized by the highest value of the PSF obtained without coronagraphic mask.
Fig. 13
Experimental radial profile comparison of dark holes obtained on the test bench without square zone (directly using the estimator described in Eq. (36)) (red, solid) and with square zones of different side lengths: K_{Sq} = 26.4λ/D_{L} (blue, dotted), K_{Sq} = 24.5λ/D_{L} (green, dashed) and K_{Sq} = 20.8λ/D_{L} (black, dotdashed). These phaseonly corrections were achieved with a γ = 16 reference pupil. The intensities are normalized by the highest value of the PSF obtained without coronagraphic mask. 
The red, solid curve shows the result for K_{Sq} = ∞, without square zone. The blue dotted line represents the result of a square mask of size K_{Sq} = 26.4λ/D_{L}, which is only slightly smaller than the actual cutoff frequency of the DM. In this case, we prevented the correction of speckles outside of the and obtained a great improvement inside the (0 to 13.5 λ/D_{L}) and a small depreciation outside (13.5 to 15.5 λ/D_{L}). Using a smaller correction zone (K_{Sq} = 24.5λ/D_{L} green dashed line) still improves the correction but to the detriment of the size of the (the contrast starts to rise around 12λ/D_{L}). Finally, we see that a smaller square zone (K_{Sq} = 20.8λ/D_{L}, black, dotdashed) produces a smaller but not shallower .
Going from K_{Sq} = ∞ to K_{Sq} = 24.5λ/D_{L}, the contrast in the progressively deepens. This is because correcting fewer of the highest frequencies with a constant number of actuators, we free degrees of freedom. However, for K_{Sq} < 25.5λ/D_{L}, the contrast level does not improve because we reach the level of the speckles created by the amplitude aberrations. Additional shrinking would only reduce the size of the . Thus, the reduction of the corrected zone in the wavefront estimation greatly improves the correction performance (up to a factor 10) with only a small reduction of the size. This effect was described in Bordé & Traub (2006) using 1D simulations.
It is important to note that this improvement does not come from the phenomenon of aliasing in the estimation (Poyneer & Macintosh 2004). Indeed, only the correction is enhanced by this process, because the estimation remained unchanged. The wavefront estimation with the SCC is only limited in frequency by the size of the reference PSF: we can estimate speckles as long as the reference flux is not null, i.e., as long as the speckles are fringed. In most cases (see next section), the first dark ring of the reference PSF is larger than the correction zone and the frequencies inside the PSF’s first dark ring are well estimated.
5.5. Size of the reference pupil
Fig. 14
Experimental radial profiles of the PSFs for reference pupils from γ = 10 to γ = 22.8 recorded on the optical bench. The distance to the center is in λ/D_{L}. These reference PSFs are normalized by the highest value of the Lyot PSF obtained without coronagraphic mask. The vertical line correspond to the frequency cutoff for N = 27 actuators in the entrance pupil (). 
In this section, we study the effect of the size of the reference pupil on the performance of the SCC. In the previous sections, we used two assumptions on the size of the reference pupil. First, in Sect. 4.2, we assumed a reference pupil small enough to consider that the influence of the aberrations inside such a reference pupil is negligible. Simulations showed that even for small γ, the level of aberrations in the reference pupil is very low and uncorrelated to the level of aberrations in the entrance pupil. Second, in Sect. 4, we assumed a reference pupil small enough to consider constant over the correction zone in the focal plane. As previously mentioned, the highest frequency attainable by the DM is . Using the first assumption, is a perfect PSF whose first dark ring is located at 1.22λγ/D_{L}. Thus, is roughly constant over the if (40)For N = 27 actuators in the entrance pupil, Eq. (40) reads γ > 15.6. In Fig. 14, we plot the radial profiles of  A_{R}  ^{2} recorded on the optical bench for γ from 10 to 22.8. We observe a wide range of intensity levels for different reference pupils (from 10^{6} for γ = 10 to 3 × 10^{8} for γ = 22.8). A reference pupil with γ = 10 (blue, solid) does not satisfy Eq. (40), and the first ring of its PSF is inside the correction zone (vertical orange dashed line). We test this case independently in Sect. 5.5.2. The other reference pupils are studied in Sect. 5.5.1.
5.5.1. Impact of small reference pupils
The size of the reference pupil can influence the correction in two different ways: it changes the signaltonoise ratio (S/N) on the fringes and modifies the flatness of the reference PSF over the correction zone. We develop these effects in this order in this section.
The S/N on the fringes is critical, because I_{−} can only be retrieved with wellcontrasted fringes. The S/N is directly related to the reference pupil size. Using Eq. (23), we deduce that the peaktopeak amplitude of the fringes in the focal plane is 2  A_{S}   A_{R} . Thus, if  A_{S}  and  A_{R}  are expressed in photons, and assuming only photon and readout noise, the S/N can be written as (41)where σ_{cam} is the standard deviation of the detector noise in photons. A higher S/N allows a better estimate of the speckle complex amplitude and thus, a better correction of the aberrations. One can notice that this S/N can be simplified depending on the relative values of its different terms. We quickly study the following cases:

if  A_{S}  ≈  A_{R}  ≪ σ_{cam}, S/N → 0. In this case, the correction is limited by the dynamic range.

if σ_{cam} ≪  A_{S}  and  A_{R}  ≪  A_{S} , S/N ~ 2  A_{R} . Initial case, at the beginning of the correction, when the Lyot pupil is a lot brighter than the reference pupil. The S/N is only a function of  A_{R} .

if σ_{cam} ≪  A_{R}  and  A_{S}  ≪  A_{R} , S/N ~ 2  A_{S} . The S/N is decreasing with deepening correction. The reference brightness is not important.
Equation (41) shows that this S/N is an increasing function of  A_{R} , but for deep corrections ( A_{S}  ≪  A_{R} ), the impact of the size of the reference is probably very weak.
The second effect is due to the assumption of a constant reference PSF over the correction zone. Variations of in the correction zone distort the wavefront estimation. This effect advocates for small reference pupils (large γ): a reference pupil of γ = 16 generates an that varies from 1 to 0.03 inside a correction zone of 27 × 27λ/D_{L}. For this reference pupil, the fringe intensity is weaker at the edges of the . Therefore, the estimate is less accurate at these locations.
Using simulation tools, where we can change the camera and photon noise easily, we were able to isolate these two different effects and analyzed their influence on the performance of the instrument separately. A more detailed study has previously been presented in Mazoyer et al. (2012).
Here, we experimentally tested the influence of the reference size. We used 27 actuators across the pupil diameter and K_{Sq} = 24.5λ/D_{L} with phaseonly correction. Figure 15 shows the radial profiles of the SCC image in RMS obtained on the laboratory bench for different reference pupils (γ = 16 and γ = 22.8), normalized by the highest value of the PSF obtained without coronagraphic mask. These results show that a large reference pupil (γ = 16) is preferable, even at the edge of the , where the reference PSF for γ = 16 is fainter than the reference PSF for γ = 22.8. Comparing the contrast levels obtained in this figure with those in Fig. 14 for all reference pupils, we deduce that we are still in the case  A_{R}  ≪  A_{S} . Deeper corrections would normally depend less on the size of the reference pupil.
Fig. 15
Radial profiles obtained on the laboratory bench for two different reference pupils (γ = 16 and γ = 22.8) with N = 27 actuators and K_{Sq} = 24.5λ/D_{L} for the size of the corrected zone. These contrasts are normalized by the highest value of the PSF obtained without coronagraphic mask. 
When we use the SCC as a planet finder there is another impact to consider: detection is possible only if the planet intensity is higher than the photon noise of the reference pupil. This effect advocates for small reference pupils. A tradeoff study of the reference size is needed depending on expected planet intensity and the actual contrast that can be achieved. A more complete study of the noise in the SCC estimation is given in Galicher et al. (2010).
5.5.2. Effect of large reference pupils
Fig. 16
PSF of the 0.8 mm reference pupil (γ = 10) (right). From this PSF we constructed the sign mask (center). The white zones are uniform and equal to 1 and the black zones are equal to − 1. Multiplying I_{−} by this mask, the correction can be achieved (right) for this reference pupil. 
In this section, we experimentally prove that we can still achieve a correction inside the using a reference pupil that does not satisfy Eq. (40) by modifying the phase estimator. This correction has previously been simulated in Galicher et al. (2010). is still considered as the complex amplitude of a perfect PSF, but we cannot consider it uniform anymore over the . First, the speckles in the first dark ring of this PSF are not fringed, because the reference PSF intensity is null at this location. The wavefront errors that produce these speckles are not estimated and are thus not corrected for. Second, the sign of and changes between the first and the second dark ring (i.e., between 1.22 and 2.23λ/D_{L}). These speckles are fringed and we can estimate the wavefront errors that produce them when we consider the sign change. Hence, when Eq. (40) is not satisfied, instead of constant, we assume constant and change the sign of over the correction zone. We now estimate (42)where , is the sign of the real part of . This function is represented in Fig. 16 (center).
In practice, to achieve the correction with this reference pupil, we multiplied I_{−} by the mask in Fig. 16 (center), where the white zones (the black zones) are constant and equal to 1 (− 1). To build this mask, we recorded the reference PSF (Fig. 16, left). From this PSF, we were able to find the dark rings of the complex amplitude. We were able to build the sign of the real part of the complex amplitude.
The tests on the optical bench were conducted using the 0.8 mm reference pupil (γ = 10) and the process described in Sect. 5.3. We used no square zone. The resulting is presented in Fig. 16 (right). We distinctly see the first reference ring at 1.22λ/D_{R}. As expected, the speckles on this ring are not corrected for, because they are not fringed. Nevertheless, apart from this ring, the whole is corrected. Although correction with a large reference pupil is possible, the level of speckle suppression is much lower (better contrast) than with smaller reference pupils (higher γ), because the speckles of the uncorrected dark ring diffract their light into the corrected zone (Galicher et al. 2010; Give’on et al. 2006).
We showed in Sect. 5.5.1 that the SCC used with a reference pupil that obeys Eq. (40) shows a better performance. However, some cases (many aberrations due to an unknown initial position of the DM, for example) may require the use of large reference pupils that produce highly contrasted fringes even with very aberrated wavefronts. The correction can then be initiated by correcting for low spatial frequencies (usually dominating the wavefront errors). Finally, the large reference is replaced with a smaller reference (which satisfies Eq. (40)) to correct higher frequencies and reach better contrast levels.
6. Conclusion
In Sect. 2.1, we used Fourier optics to model the propagation of light through a coronagraph. We then proposed a method for estimating phase and amplitude aberrations in the entrance pupil from the complex electric field measured in the focal plane after a fourquadrant phase mask coronagraph. We used this model to correct phase and amplitude aberrations in a closed loop using a DM in the pupil plane, even for a Lyot pupil smaller than the entrance pupil.
We implemented this technique, associated with a selfcoherent camera as a focal plane wavefront sensor. We corrected for phase and amplitude aberrations in a closed loop which led to speckle suppression in the central area of the focal plane (called dark hole).
We tested these methods on a laboratory bench where we were able to close the loop and obtain a stable correction at 20 Hz. When correcting for phase aberrations only, we obtained contrast levels (rms) better than 10^{6} between 6 and 12 λ/D_{L} and 3 × 10^{7} at 11 λ/D_{L}. We proved that we corrected for most phase aberrations in the dark hole and that the contrast is limited by high amplitude aberrations (10% rms in intensity) induced by the DM. When correcting for the phase and amplitude aberrations using one DM, we obtained contrast level better than 10^{6} between 2 λ/D_{L} and 12 λ/D_{L}, and better than 3 × 10^{7} between 7 λ/D_{L} and 11 λ/D_{L}. The simulation performance was limited by the diffraction of the speckles of the uncorrected area in the focal plane created by the amplitude defects. In addition, in laboratory tests, the contrast is currently limited by the camera dynamics in the aberration estimation.
We experimentally proved that a small shrinking of the size of the correction zone can improve the contrast the contrast up to a factor 10. We analyzed the influence of the reference pupil radius on the performance of the SCC and proved that the reference of γ = 16 (the larger reference pupil possible with a nonzero reference flux inside the correction zone) provides the best correction in our case.
To enhance the performance of the selfcoherent camera even more, we plan several improvements. First, one can directly minimize A_{S}, the speckle complex field measured by the SCC and not the phase estimated in the pupil plane. This approach has started to show good results (Baudoz et al. 2012) for the simultaneous correction of amplitude and phase. The correction for the amplitude errors can probably also be improved by the use of two DMs. Moreover, solutions are considered to use the SCC with wider spectral bandwidths. First tests in polychromatic light have already been conducted and show promising results (Baudoz et al. 2012). A preliminary study of these effects has been published (Galicher et al. 2010). A forthcoming paper will present a new version of the SCC that will probably overcome the current chromatic limitation.
Acknowledgments
J. Mazoyer is grateful to the Centre National d’Études Spatiales (CNES, Toulouse, France) and Astrium (Toulouse, France) for supporting his PhD fellowship. SCC development is supported by CNES (Toulouse, France).
References
 Abe, L., Domiciano de Souza, Jr., A., Vakili, F., & Gay, J. 2003, A&A, 400, 385 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baudoz, P., Boccaletti, A., Baudrand, J., & Rouan, D. 2006, in IAU Colloq. 2006, Direct Imaging of Exoplanets: Science & Techniques, eds. C. Aime, & F. Vakili, 553 [Google Scholar]
 Baudoz, P., Mazoyer, J., Mas, M., Galicher, R., & Rousset, G. 2012, in SPIE Conf. Ser., 8446 [Google Scholar]
 Beuzit, J.L., Feldt, M., Dohlen, K., et al. 2008, in SPIE Conf. Ser., 7014 [Google Scholar]
 Bordé, P. J., & Traub, W. A. 2006, ApJ, 638, 488 [NASA ADS] [CrossRef] [Google Scholar]
 Boyer, C., Michau, V., & Rousset, G. 1990, in SPIE Conf. Ser. 1237, ed. J. B. Breckinridge, 406 [Google Scholar]
 Cavarroc, C., Boccaletti, A., Baudoz, P., et al. 2006, in SPIE Conf. Ser., 6271 [Google Scholar]
 Galicher, R. 2009, Ph.D. Thesis, Université Denis Diderot Paris 7, France [Google Scholar]
 Galicher, R., Baudoz, P., & Rousset, G. 2008, A&A, 488, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Galicher, R., Baudoz, P., Rousset, G., Totems, J., & Mas, M. 2010, A&A, 509, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Give’on, A., Kasdin, N. J., Vanderbei, R. J., & Avitzour, Y. 2006, J. Opt. Soc. Am. A, 23, 1063 [NASA ADS] [CrossRef] [Google Scholar]
 Give’on, A., Belikov, R., Shaklan, S., & Kasdin, J. 2007, Opt. Express, 15, 12338 [NASA ADS] [CrossRef] [Google Scholar]
 Guyon, O., Pluzhnik, E. A., Galicher, R., et al. 2005, ApJ, 622, 744 [NASA ADS] [CrossRef] [Google Scholar]
 Guyon, O., Matsuo, T., & Angel, R. 2009, ApJ, 693, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Kalas, P., Graham, J. R., Chiang, E., et al. 2008, Science, 322, 1345 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Lagrange, A., Kasper, M., Boccaletti, A., et al. 2009, A&A, 506, 927 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Macintosh, B. A., Graham, J. R., Palmer, D. W., et al. 2008, in SPIE Conf. Ser., 7015 [Google Scholar]
 Malbet, F., Yu, J. W., & Shao, M. 1995, PASP, 107, 386 [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] [EDP Sciences] [Google Scholar]
 Marois, C., Macintosh, B., Barman, T., et al. 2008, Science, 322, 1348 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Marois, C., Zuckerman, B., Konopacky, Q. M., Macintosh, B., & Barman, T. 2010, Nature, 468, 1080 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Mas, M., Baudoz, P., Rousset, G., Galicher, R., & Baudrand, J. 2010, in SPIE Conf. Ser., 7735 [Google Scholar]
 Mas, M., Baudoz, P., Rousset, G., & Galicher, R. 2012, A&A, 539, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mawet, D., Riaud, P., Absil, O., & Surdej, J. 2005, ApJ, 633, 1191 [NASA ADS] [CrossRef] [Google Scholar]
 Mazoyer, J., Baudoz, P., Mas, M., Rousset, G., & Galicher, R. 2012, in SPIE Conf. Ser., 8442 [Google Scholar]
 Poyneer, L. A., & Macintosh, B. 2004, J. Opt. Soc. Am. A, 21, 810 [NASA ADS] [CrossRef] [Google Scholar]
 Pueyo, L., Shaklan, S. B., Give’On, A., et al. 2010, in Adaptative Optics for Extremely Large Telescopes [Google Scholar]
 Rouan, D., Riaud, P., Boccaletti, A., Clénet, Y., & Labeyrie, A. 2000, PASP, 112, 1479 [NASA ADS] [CrossRef] [Google Scholar]
 Trauger, J. T., & Traub, W. A. 2007, Nature, 446, 771 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Wallace, J. K., Burruss, R. S., Bartos, R. D., et al. 2010, in SPIE Conf. Ser., 7736 [Google Scholar]
All Figures
Fig. 1
Principle of a coronagraph (top). Aberrations in the entrance pupil plane induce speckles in the focal plane (bottom). 

In the text 
Fig. 2
Simulations of an aberrated phase in the entrance pupil plane (φ), and the real part of the field in the next pupil plane before (φ_{M}) or after (φ_{L}) the Lyot pupil. We also show the estimate (φ_{est}) and the difference between φ_{est} and φ. L = P in this case. 

In the text 
Fig. 3
Simulations of an aberrated phase with a localized default in the entrance pupil plane (φ), and in the next pupil plane (φ_{M}, φ_{L}). We represent the entrance pupil size by a dark ring around φ_{L}. We show the estimate (φ_{est}) and the difference between φ_{est} and φ in the last two images. D_{L} = 95%D_{P} in this case. 

In the text 
Fig. 4
Simulations of the influence function f_{i} in the pupil plane (left), and the effect using the two different estimators: we deconvolve by the mask (Φ_{est}, right) or not (Φ_{est,2}, center). 

In the text 
Fig. 5
Principle of the SCC combined with a coronagraph and a DM (top). A small hole is added in the Lyot stop plane to create a reference channel. In the final focal plane (bottom), the SCC image is formed by speckles encoded with Fizeau fringes. 

In the text 
Fig. 6
Correlation peaks in the Fourier transform of the focal plane. The inverse Fourier transform of I_{S} + I_{R} is circled in blue. The inverse Fourier transform of is circled in red. 

In the text 
Fig. 7
Steps followed to estimate the phase and amplitude from SCC images. 

In the text 
Fig. 8
Singular values, normalized to their respective maximum, issued from the inversion of the interaction matrices D, obtained using the two estimators Φ_{est} (red, dotted) and Φ_{est,2} (black, solid) and the synthetic matrix (blue, dashed) for γ = 40. 

In the text 
Fig. 9
Dark holes recorded on the laboratory bench for correction with two different sizes of square mask S_{q}: K_{Sq} = 20.8λ/D_{L} (left) and K_{Sq} = 24.5λ/D_{L} (center). The dark hole recorded on the laboratory bench for a correction in phase and amplitude with a square mask of size K_{Sq} = 24.5λ/D_{L} (right). These images use a different intensity scale but the same space scale. 

In the text 
Fig. 10
Radial profiles of the azimuthal standard deviation (in rms) of the intensities in the focal plane typically obtained with this method for phaseonly correction, for simulation (blue dashed line) and laboratory bench result (red solid line), for γ = 16 and a square zone of size K_{Sq} = 24.5λ/D_{L}. We also plot in this graph the simulation of the focal plane obtained using the amplitude aberrations recorded and no phase aberrations (black dashdotted line). 

In the text 
Fig. 11
Pupil illumination recorded on the laboratory bench. 

In the text 
Fig. 12
Radial profiles of the azimuthal standard deviation (in rms) of the intensities in the focal plane typically obtained with this method for phase and amplitude correction, for simulation (blue dashed line) and laboratory bench result (red solid line), for γ = 16 and a square zone of size K_{Sq} = 24.5λ/D_{L}. 

In the text 
Fig. 13
Experimental radial profile comparison of dark holes obtained on the test bench without square zone (directly using the estimator described in Eq. (36)) (red, solid) and with square zones of different side lengths: K_{Sq} = 26.4λ/D_{L} (blue, dotted), K_{Sq} = 24.5λ/D_{L} (green, dashed) and K_{Sq} = 20.8λ/D_{L} (black, dotdashed). These phaseonly corrections were achieved with a γ = 16 reference pupil. The intensities are normalized by the highest value of the PSF obtained without coronagraphic mask. 

In the text 
Fig. 14
Experimental radial profiles of the PSFs for reference pupils from γ = 10 to γ = 22.8 recorded on the optical bench. The distance to the center is in λ/D_{L}. These reference PSFs are normalized by the highest value of the Lyot PSF obtained without coronagraphic mask. The vertical line correspond to the frequency cutoff for N = 27 actuators in the entrance pupil (). 

In the text 
Fig. 15
Radial profiles obtained on the laboratory bench for two different reference pupils (γ = 16 and γ = 22.8) with N = 27 actuators and K_{Sq} = 24.5λ/D_{L} for the size of the corrected zone. These contrasts are normalized by the highest value of the PSF obtained without coronagraphic mask. 

In the text 
Fig. 16
PSF of the 0.8 mm reference pupil (γ = 10) (right). From this PSF we constructed the sign mask (center). The white zones are uniform and equal to 1 and the black zones are equal to − 1. Multiplying I_{−} by this mask, the correction can be achieved (right) for this reference pupil. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.