Research Note
Speckle temporal stability in XAO coronagraphic images
II. Refine model for quasistatic speckle temporal evolution for VLT/SPHERE
^{1} Laboratoire Lagrange, UMR 7293, Université de Nice SophiaAntipolis, CNRS, Observatoire de la Côte dÁzur, Bd. de l’Observatoire, 06304 Nice, France
email: patrice.martinez@oca.eu
^{2} European Southern Observatory, KarlSchwarzschildStraße 2, 85748 Garching, Germany
^{3} UJF – Grenoble 1/CNRS – INSU, Institut de Planétologie et d’Astrophysique de Grenoble UMR 5274, 38041 Grenoble, France
^{4} Laboratoire d’Astrophysique de Marseille, UMR 7326, CNRS, Université de Provence, 38 rue Frédéric JoliotCurie, 13388 Marseille Cedex 13, France
^{5} Office National d’Études et de Recherches Aérospatiales (ONERA), Optics Department, BP 72, 92322 Chatillon Cedex, France
Received: 29 November 2012
Accepted: 4 March 2013
Context. Observing sequences have shown that the major noise source limitation in highcontrast imaging is the presence of quasistatic speckles. The timescale on which quasistatic speckles evolve is determined by various factors, mechanical or thermal deformations, among others.
Aims. Understanding these timevariable instrumental speckles and, especially, their interaction with other aberrations, referred to as the pinning effect, is paramount for the search for faint stellar companions. The temporal evolution of quasistatic speckles is, for instance, required for quantifying the gain expected when using angular differential imaging (ADI) and to determining the interval on which speckle nulling techniques must be carried out.
Methods. Following an early analysis of a time series of adaptively corrected, coronagraphic images obtained in a laboratory condition with the highorder test bench (HOT) at ESO Headquarters, we confirm our results with new measurements carried out with the SPHERE instrument during its final test phase in Europe. The analysis of the residual speckle pattern in both direct and differential coronagraphic images enables the characterization of the temporal stability of quasistatic speckles. Data were obtained in a thermally actively controlled environment reproducing realistic conditions encountered at the telescope.
Results. The temporal evolution of the quasistatic wavefront error exhibits a linear power law, which can be used to model quasistatic speckle evolution in the context of forthcoming highcontrast imaging instruments, with implications for instrumentation (design, observing strategies, data reduction). Such a model can be used for instance to derive the timescale on which noncommon path aberrations must be sensed and corrected. We found in our data that quasistatic wavefront error increases with ~0.7 Å per minute.
Key words: instrumentation: high angular resolution / techniques: high angular resolution / telescopes
© ESO, 2013
1. Introduction
Observing sequences have shown that the major noise source limitation in highcontrast imaging is the presence of instrumental speckles and, more precisely, quasistatic speckles (Marois et al. 2003; Boccaletti et al. 2003, 2004; Hinkley et al. 2007). Speckle noise comes from wavefront errors caused by various independent sources, and evolves on different timescales depending on their nature. The first class of speckle to overcome comes from the large, dynamical wavefront error that the atmosphere generates, but realtime adaptive optics systems measure and correct it down to fundamental limitations. The fastvarying speckle noise floor left uncorrected by such systems would average out over time, because it consists in a random and uncorrelated noise for which the intensity variance converges to a zero contribution for an infinitely long exposure. For a non photon noiselimited observing run, speckle noise associated to wavefront aberrations introduced in the optical train must be tackled. In this context, instrumental speckles can be divided into two different types: (i) longtimescale wavefront errors present in the optical train (e.g., optical quality, misalignment errors) that generate static speckles that constitute a deterministic contribution to the noise variance; and (ii) slowly varying instrumental wavefront aberrations (amplitude and phase errors) from various causes, among others mechanical or thermal deformations. The latest have evolved on a shorter timescale than longlived aberrations, and are the socalled quasistatic speckles.
Instrumental speckles average to form a fixed pattern, which can be calibrated to a certain extent. A deterministic contribution to the noise variance, such as static speckles can be calibrated easily, while using a reference image, timevariable quasistatic noise can be subtracted, but its temporal evolution ultimately limits this possibility. In particular it is understood that some timescales have more of an impact than others. This is especially true because quasistatic speckles interact with other aberrations, referred to as the pinning effect, or speckle cross terms. The timescale of quasistatic speckles evolution is essential for understanding and predicting the performance of the new generation of instruments such as SPHERE (Beuzit et al. 2008), GPI (Macintosh et al. 2008), HiCIAO (Hodapp et al. 2008), and Project 1640 (Hinkley et al. 2011). The temporal evolution of these quasistatic speckles is in particular needed to quantify the gain expected with angular differential imaging (ADI, Marois et al. 2006), as well as to determine the timescale on which speckle nulling techniques should be carried out. For instance, a typical hourlong ADI observing sequence provides a partial selfcalibration of the residuals after a rotation of ~ 1λ/D at a given angular separation, which generally requires fewer than a few minutes (e.g., 5 to 7 mn at 1″ on a 8m class telescope for stars near the meridian in Hband), though it depends on wavelength, telescope latitude, and object declination. Residual speckles with decorrelation times faster than the time needed to obtain the ADI reference image cannot then be removed, while quasistatic speckles associated with larger timescales can usually be subtracted.
In this context, several authors have investigated the decorrelation timescale of quasistatic residuals in the particular context of ADI but at moderate 20−40% Strehl levels (Marois et al. 2006; Lafrenière et al. 2007; Hinkley et al. 2007), while in a former paper (Martinez et al. 2012, hereafter Paper I), we explored the realm of very high Strehl ratios (extreme adaptive optic systems, XAO). In Paper I, we analyzed a time series of adaptively corrected, coronagraphic images recorded in the laboratory with the highorder test bench (HOT), a versatile highcontrast imaging, adaptive optics bench developed at ESO. We showed that quasistatic aberrations exhibit a linear power law with time and are interacting through the pinning effect with static speckles. We examined and discussed this effect using the statistical model of the noise variance in the highcontrast imaging proposed by Soummer et al. (2007). In particular, we found that quasistatic speckles, fastevolving on the level of a few angstroms to nanometers on a timescale of a few seconds, explained the evolution of our sensitivity through amplification of the systematics. It is believed that this effect is a consequence of thermal and mechanical instabilities of the optical bench. The HOT bench is indeed localized in a classical laboratory setting, and was not initially designed to guarantee stability or mechanical stability at the level that would be required/expected for an actual highcontrast imaging instrument. Indeed Paper I emphasizes the importance of such stability for the next generation of highcontrast instruments, but the estimates found in this former analysis (amplitude, and slope of the temporal evolution) could not fairly be considered as representative of a realistic situation in order to affect operational aspects or designs of current and future real instruments.
In this paper, we confirm the results presented in Paper I with more representative measurements. We analyze a time series of adaptively corrected, coronagraphic images with the SPHERE instrument in a thermally and actively controlled environment that reproduces the realistic conditions encountered at the telescope. In this context, we propose a refined model of quasistatic speckle evolution that can be used for forthcoming highcontrast imaging instrument classes which SPHERE represents.
The paper reads as follows. In Sect. 2, we briefly recall the formalism of the statistical model of the noise variance in the highcontrast imaging context proposed by Soummer et al. (2007) and used in Paper I to discuss our former results. In Sect. 3, the experimental conditions are described, and in Sect. 4 we analyze and discuss the results. Finally, in Sect. 5, we draw conclusions.
2. Speckle noise and dynamical range
In Paper I, we examined and discussed our data using the statistical model of the noise variance in highcontrast imaging proposed by Soummer et al. (2007). Following the same formalism as for Paper I, we briefly recall here the main equations.
Soummer et al. (2007) proposes the following analytical expression for the variance of the intensity, including speckle and photon noise in the presence of static, quasistatic, and quickly varying aberrations, in the context of a propagation through a coronagraph: (1)where I denotes the intensity, is the variance of the photon noise, and N is the ratio of fastspeckle and slowspeckle lifetimes. The intensity produced by the deterministic part of the wavefront, including static aberrations, is denoted by I_{c}, while the I_{s} terms corresponds to the halo produced by random intensity variations, i.e. atmospheric (I_{s1}) and quasistatic contributions (I_{s2}). In this generalized expression of the variance, several contributions can be identified by order of appearance: (1/) the atmospheric halo; (2/) the quasistatic halo; (3/) the atmospheric pinning term, the speckle pinning of the static aberrations by the fastevolving atmospheric speckles; (4/) the speckle pinning of the static by quasistatic speckles; and finally (5/) the speckle pinning of the atmospheric speckles by quasistatic speckles. Equation (1) provides useful insight into the impact of quasistatic speckles and their interactions through the pinning phenomenon with other aberrations present in a real instrument.
As for Paper I, we focus our analysis on contribution (4/) to this noise budget in XAO coronagraphic images. In particular, we are interested in the speckle pinning of the static by the quasistatic speckles when no atmospheric contribution is present (dynamical speckles), i.e., when the contribution of I_{s1} to Eq. (1) can be neglected. Furthermore, our study concerns a situation where the photon noise is not limiting, so that the contribution from the noise variance can be neglected. In such conditions, Eq. (1) can be simplified such that (2)and the present study focuses on the effect of the cross term I_{c}I_{s2}. Since we can fairly assume that I_{s2} ≪ I_{c}, I_{s2} can be neglected except in the crossterm, and the noise variance in the raw coronagraphic image finally becomes (3)Ultimately, very deep dynamic range imager such as SPHERE, or GPI, aim to calibrate static speckles (I_{c}) such that I_{c} ≈ I_{s2}, which would largely reduce speckle pinning. Speckle nulling techniques to correct for remnant quasistatic aberrations would be a must for accessing deeper contrast level, nonetheless, the temporal characteristic of I_{s2} remains a key parameter in these circumstances (Eq. (2)).
From Eq. (3) a breakdown of this pinning effect can be carried out at the level of differential images. Raw coronagraphic images are dominated by static speckle noise. This means that the interaction between the quasistatic terms of Eq. (3), being timedependent, and static terms, assumed to be timeindependent, can be studied through differential imaging from a time series of raw coronagraphic images, which simply refers to the difference in intensity between an image recorded at time t_{0} + Δt and the reference image registered at t_{0}. In this situation, a similar expression of the noise variance for the differential images (σ_{DI}) can be derived as the difference of Eq. (3) evaluated at t_{0} + Δt, to that of the reference, at t_{0}, and reads as (4)where ΔI_{s2} represents the quasistatic evolution between the two successive images. Therefore, the quasistatic contribution can be expressed as (5)A general expression of the speckle intensity (Racine et al. 1999) is (6)where S can be related to the wavefront error φ using Maréchal’s approximation (Born & Wolf 1993) (7)and the contribution from static speckles to the wavefront error (φ_{s2} ≪ φ_{c}, and φ_{s1} neglected) can be expressed as (8)Similarly, the contribution from quasistatic speckles to the wavefront error in the differential images can be expressed as (9)Using Eqs. (8) and (9), the analysis of both direct and differential images allows to characterize the temporal properties of static and quasistatic aberrations. We note that both Eqs. (8) and (9) are approximated expressions that rely on simple rules of thumbs and assumptions (e.g., general expression of speckle intensity, lowphase aberration regime) to provide wavefront estimation per Fourier component.
3. Experimental conditions
3.1. Optical setup
SPHERE, which stands for SpectroPolarimetric Highcontrast Exoplanet REsearch, is a secondgeneration instrument for the Very Large Telescope, aiming at direct detection and spectral characterization of extrasolar planets. The instrument is now nearing completion in its final integration stage in Europe, before shipping to Chile. Now in its final test phase, it offers a unique opportunity to carry out research on static and quasistatic speckles. SPHERE is a unique instrument, with first light in 2013, including a powerful extreme adaptive optics system (SAXO), an infrared differential imaging camera (IRDIS), an integral field spectrograph (IFS), and a visible differential polarimeter (ZIMPOL). The time series of adaptively corrected, coronagraphic images that will be discussed through this paper have been obtained using SPHERE with the IRDIS instrument, so that for the sake of clarity, only the systems that have been used and are relevant will be described further.
The SPHERE adaptive optics for exoplanet observation (SAXO) uses a 41 × 41 actuator deformable mirror (DM) of 180 mm diameter with interactuator stroke >±1 μm, a maximum stroke >±3.5 μm, and twoaxis tiptilt mirror (TTM) with ±0.5 mas resolution. The wavefront sensor is a 40 × 40 subapertures ShackHartmann sensor equipped with a spatial filter for aliasing minimization. During the test, no dynamical turbulence was present in the system, except the low internal turbulence. Optical elements are installed inside an hermetical enclosure. As a result, SAXO is used to correct for internal turbulence, static aberrations, and guarantee image and pupil stability. As image and pupil stability are essential in highcontrast imaging, differential image movements due to thermomechanical effects are measured in realtime using an auxiliary NIR tiptilt sensor located close to the coronagraph focus plan, and corrected with a differential tiptilt mirror in the wavefront sensor arm. Similarly, pupil runout is accounted for and corrected by a pupil tiptilt mirror at the entrance of the instrument. The nearIR Strehl ratio is >95%. In particular, non common path aberrations (from the DM to the detector) are measured offline using a phase diversity algorithm and compensated for by reference slopes adjustment. In these conditions, static wavefront errors left uncorrected in the system are estimated at the level of ~6 nm rms.
SPHERE offers various coronagraphic possibilities (classical Lyot coronagraphs, apodized Lyot coronagraphs – ALC –, and achromatic four quadrants phase masks). During the experiment, a 5.2λ/D ALC was used.
The infrared dual beam imaging and spectroscopy (IRDIS) subsystem includes a spectral range from 950 to 2320 nm and an image scale of 12.25 mas per pixel (Nyquist sampling at 950 nm). The field of view is greater than 11″ square, with a 2k × 2k HawaiiIIRG detector. The main mode of IRDIS is the dualband imaging (DBI), providing images in two neighboring spectral channels with minimized differential aberrations. Ten different filter couples are defined corresponding to different spectral features in modeled exoplanet spectra. During the experiment the narrow Hband couple filters have been used (centered around 1593 nm, and 1667 nm, R = 30). The DBI mode is not of the interest of the present study since we are interested in the temporal evolution of aberrations; nevertheless, it allows simultaneous characterization of quasistatic speckles at two different wavelengths for comparison/confirmation purposes.
Finally, the SPHERE instrument is installed on an active tripod damping system and fully covered by an hermetic enclosure.
Fig. 1 Coronagraphic images recorded at t_{0}, t_{0} + 10 mn, and t_{0} + 100 mn. The Strehl ratio is ~95%. The arbitrary color scale and dynamic range (identical for the three images) were chosen to enhance the contrast for the sake of clarity. The field covered in the images is ~2″ square. 

Open with DEXTER 
Fig. 2 Differential coronagraphic images. Left: difference of the t_{0} + 10 mn image to the reference t_{0}. Right: difference of the t_{0} + 100 mn image to the reference t_{0}. The increase in the strength of residuals can be qualitatively observed here (the intensity range is similar in both images). 

Open with DEXTER 
3.2. Thermal functioning
The SPHERE instrument must conform to several environmental specifications. In particular, it will operate under a temperature range from 5 to 18 °C, while the highest gradient of temperature at the VLT is − 0.9 [entity ! #x20 ! ] °C/h and − 1.4 [entity ! #x20 ! ] °C/h respectively for 80% and 95% of the nights. To validate that SPHERE is compliant with such conditions, cold tests have been set up for functional and performance evaluation at different ambient temperatures within the operational range and with realistic transient conditions reflecting situations encountered at the VLT. The whole SPHERE instrument was installed in a cold tent (150 m^{3}) and cooled down with an efficient air conditioning system. The cooling system allows reducing the temperature inside the tent by ~10 °C compared to the outside (in the integration hall). Numerous temperature probes have been installed at several critical locations in the instrument to accurately monitor the evolution inside/outside the SPHERE enclosure.
3.3. Observing run and data reduction
The observing run consists in recording a time series of AO closeloop IRDIS coronagraphic images. Each coronagraphic image of the time series corresponds to a series of 3 s short exposure images, averaged out over ~3 mn. During the experiment, the temperature outside the enclosure (inside the cold tent) was monitored and programed to decrease from 15 °C with a rate of 2 °C per hour. The total duration of the experiment was 100 mn. The gradient inside the enclosure was estimated at the level of 0.2 °C/h. The datareduction process corrects for bad pixels and background, and normalizes the images with respect to exposure time and flux. Before subtraction, a fine, subpixel correction of the residual tiptilt component is performed on the raw images. The analysis was done on the two Hband filters data set available and converge to similar results, so that only one data set is presented here. Depending on the nature of the image analyzed, we applied different metrics. The contrast refers to the ratio of intensity in the raw coronagraphic image, averaged azimuthally at a given angular separation, to the peak intensity of the direct flux. When studying a differential image, implying the subtraction of a reference image, the average contrast is no longer suitable. The detectability is used then, which stands for the azimuthal standard deviation measured in a ring of width λ/D. It quantifies the ability to distinguish a companion at a given angular distance.
4. Analysis and interpretation
Figure 1 presents three coronagraphic images extracted from the time series: the coronagraphic image recorded at t_{0} (the reference), t_{0}+10 mn, and t_{0}+100 mn. In Fig. 2, the corresponding differential images are shown: the difference of the t_{0} + 10 mn image to the reference t_{0} (left), and the difference of the t_{0} + 100 mn image to the reference t_{0} (right). The corresponding profiles for both raw and differential images are presented in Fig. 3.
Qualitatively, the raw coronagraphic images (Fig. 1) demonstrate starlight attenuation, and exhibit static speckles with lower intensity in the AOcorrection domain (from the second or third wing of the coronagraphic image to the rise in the AO cutoff frequency). A radial trend in speckle intensity is observable in the image: speckles closer to the center of the image are brighter. The central part of the coronagraphic image is dominated by diffraction residuals. The AO cutoff frequency can be readily seen in the image owing to the slope of intensity in the speckle field at 0.8″ (20λ/D) from the center, as expected. Outside the inner domain defined by the AO cutoff frequency, the AO system cannot measure or correct the corresponding spatial frequencies. Various bright spots are observable at ~40λ/D (twice the AO cutoff frequency) and correspond to the inter actuator pitch spatial frequency.
Fig. 3 Contrast profiles of a time series of coronagraphic images (top) and detectability (1σ) of the differential images (subtraction of the time series of coronagraphic images to the reference one, bottom). 

Open with DEXTER 
As observed in Fig. 1 and presented in Fig. 3 (top curves), raw coronagraphic images are dominated by the static contribution, for which contrast profiles are stable over time at any angular separation. Raw coronagraphic images are dominated by static speckle noise, since no evolution in the speckle field is visually detectable over the three images presented that cover the temporal duration of the experiment. This is consistent with the fact that the interaction between the quasistatic terms of Eq. (3), being timedependent, and static terms, assumed to be timeindependent, can be studied through differential imaging from a time series of raw coronagraphic images, and not directly at the level of raw images.
The static wavefront error amplitude has been evaluated on the raw coronagraphic images using Eq. (8), and converges to the value of ~5 nm rms, though it might be slightly underestimated (see Sect. 3.1). This is basically a factor of 10 lower than the estimation obtained on the HOT bench, and it contributes to the increase in dynamic range of the instrument to that of the performance obtained with the HOT bench in both classical and differential imaging (virtually by an order of magnitude). In contrast to raw images, an evolution in the speckle field is seen in the differential images. The profiles presented in Fig. 3 (bottom curves) clearly indicate that the detectability level degrades with time. These profiles demonstrate that raw coronagraphic images evolve temporally, since they are less and less correlated with the reference over time, even though such an evolution cannot be readily seen in raw images. Furthermore, this degradation of the detectability is effective at all angular separations. This result is compliant with observations formally reported in Paper I. From the time series of differential images, using Eq. (9), we derive the quasistatic wavefront error contribution per Fourier component, of the pinning effect at several angular separations.
Fig. 4 Time variability of wavefront error due to quasistatic speckles, evaluated at various angular separations (observational data). 

Open with DEXTER 
Figure 4 shows the temporal evolution of φ_{s2} at 3, and 10 λ/D, i.e. the quasistatic wavefront error (rms) as a function of time at several angular separations. It clearly indicates that φ_{s2} depends on time and increases with time, justifying the constant degradation of detectability observed in Fig. 3 (bottom curves). This is true for all angular separations, and the shorter the separation the higher the amplitude. Quantitatively, the level of quasistatic wavefront error that limits the sensitivity in the differential images (~10^{7} to ~10^{8}, 1σ) is found to be in the regime of ~0.1 Å to 1 Å. The power law of the temporal evolution of the quasistatic wavefront error is derived and exhibits a similar tendency regardless of the angular position in the field. It can be fitted by a linear function of time. The parameters for the linear fits are presented in the legend of Fig. 4, and can be used to extrapolate a model for quasistatic speckle evolution in the context of highcontrast imaging instruments. From the linear fits derived at various angular separations, and considering that the difference in amplitude is not significant, we attempt to generalize the expression of the power law for any angular position in the field, which reads as the following approximation: (10)where t is the time in minutes, and Δφ_{s2} is expressed in nm rms.
In Paper I, the parameters found for the fit were 0.250 for the value at the origin, and a slope of 0.012. This means that SPHERE is definitely more stable in term of thermomechanical variations than the HOT bench, though this was expected (a slope reduced by at least an order of magnitude, since data were not recorded a transient regime as performed here, but in stabilized temperature environment). While scaling the fit parameters to model quasistatic temporal evolution observed with the HOT bench to real instrument was highly non trivial as it significantly depends on the instrument environment (temperature or pressure changes, mechanical flexures, etc.), we believe that the estimates found with SPHERE and presented in Eq. (10) are representative of the new highcontrast imaging instrument generation that SPHERE represents, and can be used as such. In our data we found that quasistatic wavefront error increases with ~0.7 Å per minute. Quasistatic speckles, fastevolving on the level of a few angstroms on a timescale of a few seconds, explained the evolution of our sensitivity through amplification of the systematics. It is believed that this effect is a consequence of thermal and mechanical instabilities (pressure changes, mechanical flexures, etc.) of the optical bench. This is qualitatively supported by the fact that outside the AO control radius, a characteristic butterfly shape drawn by the speckle pattern in the differential image at t_{0} + 100 mn is observable (Fig. 2, right). This indicates that a beam shift is at work, though it is difficult to obtained quantitative information on it.
Fig. 5 Quasistatic aberrations power spectral density on various timescales (using Eq. (5)). 

Open with DEXTER 
I functions described in Sect. 2 essentially represent power spectral densities (PSD). While I_{c} symbolizes the PSD of the static wavefront, ΔI_{s2} stands for the PSD of the differential aberrations, which is calculated by Eq. (5). Plotting Eq. (5) gives access to the quasistatic aberration PSD as a function of time (Fig. 5). From the DSPs presented in Fig. 5, we can observe that (1/) at very low frequencies (from 0.1 to 3 cycles per pupil) the PSD is temporally roughly stable and exhibits a f^{0} power law, where f is the spatial frequency, while it must be pointed out that in this frequency domain, Eq. (5) might not be entirely valid; (2/) the PSD exhibits a f^{5} power law at low frequencies (from 3 to 8 cycles per pupil), (3/) the PSD exhibits a f^{0} power law at intermediate frequencies (from 8 to 20 cycles per pupil), then is essentially dominated by a white noise; (4/) the power law is again f^{5} at high frequencies (outside the AO control domain, from 20 to 30 cycles per pupil, noisedominated at higher spatial frequencies). From Fig. 5 it is difficult to extract any other unambiguous information, or initiate preliminary explanations, such as the f^{5} power law contradicting the generally adopted and standard PSD slope of f^{2} for surface optics (static aberrations).
5. Conclusion
In this paper, we have confirmed the results formerly reported in Paper I with new measurements, and derived a practical and generalized expression to model quasistatic speckles temporal evolution for any angular position in the field (Eq. (10)), in the context of the forthcoming highcontrast planet imagers. Quasistatic aberrations observed in a time series of extreme adaptive opticscorrected coronagraphic images exhibit a linear power law with time and are interacting through the pinning effect with static aberrations. We examined and discussed this effect using the statistical model of the noise variance in highcontrast imaging proposed by Soummer et al. (2007), where the effect of
pinning quasistatic to static speckles as described by the expression for the variance (Eq. (9)) is found to reflect the situation in our data set fairly well.
We found that quasistatic wavefront error increases with a rate of ~0.7 Å per minute. In this context, Eq. (10) is meant to feed complex system analysis for timescale specifications and affect assessments in the context of calibration/operational strategies for highcontrast imaging class instruments. In addition, Eq. (5) provide useful insight into the PSD of quasistatic aberrations in the system and, at first glance to identify what moves in the system. The proposed model of quasistatic speckles temporal evolution can be used to derive timescales for calibration/operational aspects, such as ADI, or non common path aberrations correction.
The case considered in this present paper represents a static system subject to a thermal gradient. The foreseen impact of rotating components in the instrument, such as, the atmospheric dispersion compensators (ADCs, not seen by the AO system), or the derotator (seen by AO system) will be treated in a separated study.
Acknowledgments
SPHERE is an instrument designed and built by a consortium consisting of IPAG, MPIA, LAM, LESIA, Laboratoire Lagrange, INAF, Observatoire de Genève, ETH, NOVA, ONERA, and ASTRON in collaboration with ESO.
References
 Beuzit, J.L., Feldt, M., Dohlen, K., et al. 2008, in SPIE Conf. Ser., 7014 [Google Scholar]
 Boccaletti, A., Chauvin, G., Lagrange, A.M., & Marchis, F. 2003, A&A, 410, 283 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Boccaletti, A., Riaud, P., Baudoz, P., et al. 2004, PASP, 116, 1061 [NASA ADS] [CrossRef] [Google Scholar]
 Born, M., & Wolf, E. 1993, Principles of Optics (Cambridge University Press) [Google Scholar]
 Hinkley, S., Oppenheimer, B. R., Soummer, R., et al. 2007, ApJ, 654, 633 [NASA ADS] [CrossRef] [Google Scholar]
 Hinkley, S., Oppenheimer, B. R., Zimmerman, N., et al. 2011, PASP, 123, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Hodapp, K. W., Suzuki, R., Tamura, M., et al. 2008, in SPIE Conf. Ser., 7014 [Google Scholar]
 Lafrenière, D., Marois, C., Doyon, R., Nadeau, D., & Artigau, É. 2007, ApJ, 660, 770 [NASA ADS] [CrossRef] [Google Scholar]
 Macintosh, B. A., Graham, J. R., Palmer, D. W., et al. 2008, in SPIE Conf. Ser., 7015 [Google Scholar]
 Marois, C., Nadeau, D., Doyon, R., Racine, R., & Walker, G. A. H. 2003, in Brown Dwarfs, ed. E. Martín, IAU Symp., 211, 275 [NASA ADS] [Google Scholar]
 Marois, C., Lafrenière, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556 [NASA ADS] [CrossRef] [Google Scholar]
 Martinez, P., Loose, C., Aller Carpentier, E., & Kasper, M. 2012, A&A, 541, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Racine, R., Walker, G. A. H., Nadeau, D., Doyon, R., & Marois, C. 1999, PASP, 111, 587 [NASA ADS] [CrossRef] [Google Scholar]
 Soummer, R., Ferrari, A., Aime, C., & Jolissaint, L. 2007, ApJ, 669, 642 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Coronagraphic images recorded at t_{0}, t_{0} + 10 mn, and t_{0} + 100 mn. The Strehl ratio is ~95%. The arbitrary color scale and dynamic range (identical for the three images) were chosen to enhance the contrast for the sake of clarity. The field covered in the images is ~2″ square. 

Open with DEXTER  
In the text 
Fig. 2 Differential coronagraphic images. Left: difference of the t_{0} + 10 mn image to the reference t_{0}. Right: difference of the t_{0} + 100 mn image to the reference t_{0}. The increase in the strength of residuals can be qualitatively observed here (the intensity range is similar in both images). 

Open with DEXTER  
In the text 
Fig. 3 Contrast profiles of a time series of coronagraphic images (top) and detectability (1σ) of the differential images (subtraction of the time series of coronagraphic images to the reference one, bottom). 

Open with DEXTER  
In the text 
Fig. 4 Time variability of wavefront error due to quasistatic speckles, evaluated at various angular separations (observational data). 

Open with DEXTER  
In the text 
Fig. 5 Quasistatic aberrations power spectral density on various timescales (using Eq. (5)). 

Open with DEXTER  
In the text 