Issue 
A&A
Volume 663, July 2022



Article Number  A49  
Number of page(s)  15  
Section  Astronomical instrumentation  
DOI  https://doi.org/10.1051/00046361/202243145  
Published online  13 July 2022 
Loworder wavefront control using a Zernike sensor through Lyot coronagraphs for exoplanet imaging
Blind stabilization of an image dark hole
^{1}
Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS,
Laboratoire Lagrange,
06108
Nice, France
email: raphael.pourcelot@oca.eu
^{2}
Space Telescope Science Institute,
3700 San Martin Drive,
Baltimore, MD
21218, USA
^{3}
Aix Marseille Université, CNRS, CNES, LAM (Laboratoire d’Astrophysique de Marseille)
UMR 7326,
13388
Marseille, France
^{4}
ONERA, The French Aerospace Lab, BP72,
29 avenue de la Division Leclerc,
92322
Châtillon Cedex, France
^{5}
Thales Alenia Space,
5 Allée des Gabians – B.P. 99 –
06156
Cannes la Bocca Cedex – France
^{6}
Hexagon Federal,
Chantilly,
VA 20151, USA
^{7}
Department of Mechanical and Aerospace Engineering, Princeton University,
Princeton,
NJ 08540, USA
^{8}
The Institute of Optics, University of Rochester,
Rochester,
NY 14627, USA
^{9}
NASA Goddard Space Flight Center,
Greenbelt,
MD 20771, USA
Received:
17
January
2022
Accepted:
13
April
2022
Context. The combination of large segmented space telescopes, coronagraphy, and wavefront control methods is a promising solution for producing a dark hole (DH) region in the coronagraphic image of an observed star in order to study planetary companions. The thermal and mechanical evolution of such a highcontrast instrumental setup leads to wavefront drifts that degrade the DH contrast during the observing time, thus limiting the ability to retrieve planetary signals.
Aims. Lyotstyle coronagraphs are starlightsuppression systems that remove the central part of the image for an unresolved observed star, that is, the point spread function, with an opaque focal plane mask (FPM). When implemented with a flat mirror containing an etched pinhole, the mask rejects part of the starlight through the pinhole which can be used to retrieve information about loworder aberrations.
Methods. We propose an active control scheme using a Zernike wavefront sensor (ZWFS) to analyze the light rejected by the FPM, control loworder aberrations, and stabilize the DH contrast. We first present the concept formalism and then describe how we characterized the sensor behavior in simulations and in the laboratory. We performed experimental tests to validate a wavefront control loop using a ZWFS on the HiCAT testbed.
Results. By controlling the first 11 Zernike modes, we show a decrease in the standard deviation of the wavefront error by a factor of up to 9 between open and closedloop operations using the ZWFS. In the presence of wavefront perturbations, we show the ability of this control loop to stabilize a DH contrast around 7 × 10^{−8} with a standard deviation of 7 × 10^{−9}.
Conclusions. Active control with a ZWFS proves to be a promising solution in Lyot coronagraphs with an FPMfiltered beam for controlling and stabilizing loworder wavefront aberrations and DH contrast for exoplanet imaging with future space missions.
Key words: instrumentation: high angular resolution / methods: data analysis / telescopes
© R. Pourcelot et al. 2022
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the SubscribetoOpen model. Subscribe to A&A to support open access publication.
1 Introduction
With about 5000 known planets orbiting nearby stars^{1}, the search to find an Earth analog around a Sunlike star has received increasing attention over the past few years. Still, to assess the resemblance of a given exoplanet candidate with Earth, characteristics such as atmospheric composition and physical parameters need to be retrieved. Direct imaging and spectroscopy of exoplanets represent a promising approach to collect photons from a planetary companion around its parent star and perform spectral characterization of the planet. However, this proves extremely challenging as exoplanets, and in particular Earth twins, are up to 10^{10} times fainter than their host star at angular separations shorter than 100 mas in the visible wavelength range.
Since the mid 2010s, circumstellar disks and planetary companions have been routinely observed around nearby stars with several stateoftheart exoplanet imagers on eightmeterclass groundbased telescopes (Macintosh et al. 2014; Jovanovic et al. 2015; Beuzit et al. 2019), with large surveys being carried out (e.g., Nielsen et al. 2019; Vigan et al. 2021). Such instruments have so far enabled the observation of young gas giant planets with a planettostar flux ratio, or contrast, of down to 10^{−6} at angular separations of as small as 200 mas from their host star in the near infrared (NIR; e.g., Macintosh et al. 2015; Chauvin et al. 2017; Keppler et al. 2018; Currie et al. 2020). The next generation of large groundbased facilities includes the extremely large telescopes (ELTs, Simard et al. 2016; Fanson et al. 2018; Ramsay et al. 2020) which are designed to help us to further understand the formation and evolution of planetary systems. Equipped with instruments to image small gaseous or large telluric planets and spectrally analyze their atmosphere, they are expected to achieve contrasts of down to 10^{−8} at distances shorter than 50 mas in the NIR (Fitzgerald et al. 2019; Guyon et al. 2019; Kasper et al. 2021). In addition, future large observatories in space with highcontrast capabilities are envisioned options to observe small rocky planets with 10^{−10} contrast at less than 50 mas in reflected light in the visible wavelength range.
Following this path, in its 2020 Astronomy and Astrophysics Decadal Survey, the United States National Academies of Science, Engineering, and Medicine (NASEM) recommends the development of a space telescope of at least 6 m in diameter, large enough to provide the sensitivity and resolution required for the detection and characterization of Earthlike exoplanets (National Academies of Sciences, Engineering, & Medicine 2021). Such a large telescope is expected to be deployable and have a segmented primary mirror to fit current and expected launch vehicles. This mission will build on the recent studies from two NASA mission concepts, the Habitable Exoplanet Imaging Mission (HabEx, Gaudi et al. 2020) and the Large Ultraviolet Optical Infrared surveyor (LUVOIR, The LUVOIR Team 2019). While HabEx has an approach based on a 4mdiameter monolithicmirror telescope, LUVOIR relies on a concept with a segmented primary mirror with a diameter of up to 15m.
To overcome the large flux ratio between an observed star and its planet, several teams are currently investigating different schemes of highcontrast instruments with coronagraphic capabilities. Preliminary coronagraph designs, simulations, and experiments show 10^{−9}–10^{−10} raw contrast capabilities, with the assumption of very low wavefront errors or sufficient photon flux for wavefront correction loops (e.g., Ruane et al. 2018; N’Diaye et al. 2016; Por 2020; Sirbu et al. 2021; Seo et al. 2019). Carried out to refine the stability specifications for LUVOIR, the ULTRA study presents a total wavefront error budget of 158 pm root mean square (RMS) on all timescales and for aberrations with all spatial frequency contents (Coyle et al. 2019). To meet these constraints, active optics with deformable mirrors (DMs) provide a promising solution for two main reasons: (i) They are able to produce a highcontrast region, or dark hole (DH), in the coronagraphic image of an observed star that can be used to enhance the planet signal, and (ii) they can maintain this DH during observations by correcting for the aberrations induced by the thermal and mechanical evolution of the observatory and measured by wavefront sensors (WFSs) in real time.
The firststep active optics system leads to the generation of a DH with an optimal wavefront shape to reach the targeted contrast. However, this wavefront shape experiences small perturbations due to mechanical flexures and thermal variations during the observing time. These perturbations lead to aberrations that degrade the achieved DH contrast. Several secondstep wavefrontsensing schemes and control solutions have been proposed in recent years to stabilize the DH contrast during the observation, specific to each coronagraph design (e.g., Jovanovic et al. 2018). Among these wavefrontsensing solutions, the Zernike wavefront sensor (ZWFS, e.g., Zernike 1934; Bloemhof & Wallace 2003; Dohlen 2004; Wallace et al. 2011; N’Diaye et al. 2013b) has proven to be one of the most efficient sensors in terms of photon usage (Guyon 2005; Chambouleyron et al. 2021), and is capable of measuring aberrations down to the picometer level (Ruane et al. 2020; Steeves et al. 2020). In the context of the Nancy Grace Roman Space Telescope (RST, Spergel et al. 2015; Krist et al. 2015), this solution was adopted in the presence of a monolithic primary mirror to control lineofsight jitterpointing errors, thermal focus drifts, and other loworder modes. On RST, the ZWFS uses the light rejected in reflection by the focal plane mask (FPM) of the Hybrid Lyot Coronagraph or the Shaped Pupil Coronagraph (Shi et al. 2016, 2018).
In the case of segmented aperture telescopes, such as LUVOIR, recent developments have considered Lyottype coronagraphs (LCs) as promising solutions for starlight suppression. LCs use a FPM that can be implemented using a flat mirror with a pinhole to reject starlight. Such a configuration was studied in the Lyot project (Sivaramakrishnan et al. 2007) and was successfully implemented on the exoplanet imagers Palomar P1640 (Oppenheimer et al. 2012) and the Gemini Planet Imager for exoplanet observations (Macintosh et al. 2014). Taking advantage of the available photons from the FPMfiltered beam, we explore the use of a secondstep active control loop using a ZWFS to control loworder aberrations. The aim is to maintain the wavefront shape at a point that yields a minimum DH contrast.
In this paper, we present the concept of a ZWFS in a LC FPMfiltered beam, the formalism to describe the optical propagation, and the analytical phase reconstruction of the loworder aberrations with our ZWFSbased approach. Using the Highcontrast imager for Complex Aperture Telescopes testbed (HiCAT, N’Diaye et al. 2013a, 2014, 2015; Leboulleux et al. 2016; Soummer et al. 2018; Moriarty et al. 2018), we implemented and tested our approach for experimental validations in the presence of a segmented aperture. We studied the ZWFS response for the first 11 Zernike modes and with different DM offsets, which correspond to the commands that are sent to the DMs to either maximize the optical wavefront quality or produce an image DH on HiCAT. Finally, we show the stability results for the wavefront error and DH contrast that are achieved with our concept in the presence of drifts introduced on DMs or with various configurations of internal air turbulence on the bench. Specific measurements of segment piston, tip, and tilt errors are not considered in this study and are left for a future analysis.
2 Zernike Wavefront Sensor for a FPMfiltered Beam
In the context of LCs, we explore strategies to use the light rejected by their FPM with a ZWFS to control loworder aberrations. Figure 1 shows a schematic diagram of this configuration. The complex amplitude of the electric field ψ_{0} in the pupil plane 0 at the entrance pupil of the system is spatially filtered by the FPM in the downstream focal plane: we study the light that goes through the mask to the reimaged pupil plane A. From this plane onward, the propagation is identical to the standard ZWFS described by N’Diaye et al. (2013b). The filtering step induced by the FPM modifies the electric field, and this impact is detailed in the following section. We then present the corresponding analytical phase reconstruction and the closedloop strategy to control loworder wavefront errors.
2.1 Formalism of the ZWFS in a Filtered Beam
We present the formalism that relates the electric fields in the entrance pupil and on the ZWFS camera. The formalism is detailed in two parts: first we describe the interference induced by the Zernike phase mask, and then we provide details of the effect of the FPM filtering.
Derived from the phase contrast technique proposed by Zernike (1934), the ZWFS uses a phaseshifting mask that is located in a focal plane in which the point spread function (PSF) of an unresolved source is formed. Its typical size is one resolution element ~λ/D, where λ and D denote the source wavelength and the aperture diameter, respectively. This mask, usually manufactured as a glass plate with a small phase dimple, produces interference between the central core of the PSF that goes through the dimple and the rest of the light that passes outside the dimple. The light propagation through the ZWFS can be described using the formalism presented in N’Diaye et al. (2013b): we refer to ψ_{A} which is the complex amplitude of the electric field in the relay pupil plane A between the FPM and the Zernike phase mask. This complex amplitude can be decomposed into a pupil amplitude P and phase φ such that (1)
In the downstream focal plane B in which the Zernike mask is inserted, the complex amplitude ψ_{B} can be expressed as (2)
where denotes the Fourier transform of f, and t the transmission of the Zernike mask. If M_{Z} is the tophat function describing the mask dimple with a value of 0 outside the mask and 1 inside, and θ is the phase shift introduced by the mask, t can be expressed as (3)
For the standard ZWFS, θ is usually equal to π/2 but the formalism remains valid for any value that is not a multiple of 2π (absence of a phase mask shift). Finally, the complex amplitude ψ_{C} on the ZWFS camera in the reimaged pupil plane C can be written as (4) (5) (6)
where ⊗ is the convolution product and , which is the wave diffracted by the Zernike mask. The intensity I_{C} = ψ_{C}^{2} on the detector can then be expressed as a function of φ: (7)
This expression does not make any assumption on P: the geometry and the amplitude of the pupil can be arbitrary and therefore our formalism is also valid for a pupil containing a grayscale apodization. Filtering the beam using a FPM does not change the relation between the phase in the intermediate pupil A and the intensity on the detector.
We define φ_{0} as the phase of the wave in the entrance pupil plane 0 and ψ_{0} as its electric field: (8)
The next step is to find a relation between the phase error φ_{0} in the entrance pupil plane 0 and φ, the phase in plane A. In the general case, this link can be established through the electric fields ψ_{0} and ψ_{A}: (9)
where M_{FPM} is the tophat function of diameter d_{FPM} that describes the FPM pinhole. In the general case, the relation between φ_{0} and φ cannot be established in a straightforward way and more complex strategies have to be used, such as a datadriven calibration with a Jacobian matrix, as detailed in Sect. 2.3.
Another aspect of this description deals with the properties of b in the presence of the FPM, which can be implemented in reflection (e.g., a mirror with a pinhole) or in transmission. By replacing ψ_{A} in the definition of b by Eq. (9), we obtain: (10)
In this equation, is the Fourier transform of M_{Z} × M_{FPM}, with × the elementwise multiplication. For most applications, the Zernike phase dimple diameter d_{Z} is smaller than the FPM pinhole diameter in angular units: d_{Z} ≈ 1 λ/D and d_{FPM} is larger than 4 λ/D. In this case, M_{Z} × M_{FPM} is equal to M_{Z}. As a result, can be simplified as . The wave diffracted by the mask can therefore be given by (11)
In the regime of small aberrations (φ_{0} ≪ 1 rad), the central core of the PSF can be considered unaberrated, meaning that b is constant, as described in N’Diaye et al. (2013b). With the expression of ψ_{0} in Eq. (8), b can be written as (12)
Under these simplifying assumptions, the term b hence only depends on the entrance pupil aperture and the Zernike phase mask, but not on the FPM. This aspect proves useful for the reconstruction method detailed in the following section.
Fig. 1 Schematic diagram of the ZWFS downstream of a binary amplitude FPM. The light in the entrance pupil comes from the telescope and deformable mirrors in plane 0. It is focused on the FPM, which is reimaged on the Zernike mask in plane B. The intermediate pupil plane A is not required. 
2.2 Analytical Reconstruction
We now present the inversion of Eq. (7) to retrieve φ from the measurement I_{C}.
In the smallaberrations regime (φ ≪ 1 rad), we can use a secondorder Taylor expansion of φ as detailed in N’Diaye et al. (2013b) to derive the following expression: (13)
This quadratic equation is presented as a good approximation in the [−π/4, π/4] range for φ.
Assuming ρ and b are nonnull, this seconddegree polynomial equation can be solved by computing the discriminant Δ with (14)
Finally, φ can be retrieved using (15)
In the general case (without any assumption on φ), it is possible to derive the exact expression for the phase φ from Eq. (7), yielding (16)
Using either the secondorder approximation or the exact inversion, the phase reconstruction always requires an estimation of P and b. The latter can be computed using a numerical propagation of the entrance pupil plane amplitude P_{0}. To eliminate systematic errors using an analytical input, b can also be computed using the pupil amplitude in plane A P, as b is not modified by the FPM filtering. A measurement of p^{2} can be recorded on the ZWFS camera by simply moving the Zernike mask out of the beam. In most optical systems, the assumption is valid as the electric field in the entrance pupil is positive or null for all kinds of pupils and apodizations. However, in the presence of FPM filtering, the pupil amplitude P is not necessarily positive or null in the pupil plane, because the FPM filtering induces πphase shifts in some areas of the pupil measurement. To still make use of the empirical estimation of P, an analytical pupil sign mask input has to be provided which allows us to determine the sign of P all over the pupil. Such a mask can be obtained by numerically propagating P_{0} through the optical system.
Figure 2 shows an example of the image I_{C} recorded on the HiCAT bench, normalized by P^{2}, and the reconstructed wavefront φ. Normalizing the image removes some artifacts coming from pupil illumination discrepancies. The full phase reconstruction requires the use of Eqs. (15) or (16), including the estimation of b.
Throughout this paper, we use the opensource code pyZELDA (Vigan & N’Diaye 2018) with the updated formalism of N’Diaye et al. (2013b) presented here to take into account the FPM filtering. Our code performs phase reconstruction using Eq. (15), as we are working with small phase errors.
Fig. 2 Example of intensity map as measured by the ZWFS (left), normalized by the P^{2} measurement (middle), and the reconstructed phase map using the secondorder analytical reconstructor (right). 
2.3 Closedloop Sensing and Control
We assume that a DM is located in an upstream pupil plane conjugated with plane 0 to inject Zernike modes in our system, and to control wavefront errors.
The phase reconstruction using our analytical method based on phase conjugation can be used as a wavefront error estimator in closedloop for wavefront stabilization (Pourcelot et al. 2021). Here we work with the approach based on the use of a Jacobian matrix J, often referred to as an interaction matrix. This proves suitable for enhanced controlloop stability, in particular for long experiments in which improper alignment between the Zernike camera and the DM may cause divergence in the phase reconstruction process. The Jacobian matrix approach is a powerful way to relate the measured filtered aberrations to the standardized unfiltered Zernike modes. In addition, its singular value decomposition (SVD) gives the singular modes which can easily be filtered to select the modes to control. The phase reconstruction with our analytical method is used as an additional diagnostic tool, in particular to estimate the total wavefront error (WFE) in a reallife system.
The datadriven calibration for this Jacobian matrix is performed by poking the DM actuators with m modes from a suitable basis. A DM probe introducing an optical path difference (OPD) value p is used and the corresponding image I_{C} in its vectorized version is recorded within a preestablished region of interest (ROI) with n pixels on the ZWFS camera. For each mode i, the intensity responses to a positive and a negative probe are recorded, yielding and measurements for I_{C}. The matrix J is therefore constructed by concatenating the vectors for the m modes. As this matrix is rectangular, we perform a SVD to produce its pseudoinverse C, the control matrix. Using the SVD, the Jacobian matrix is expressed as (17)
where S is a diagonal matrix whose values are the singular values of J. The matrices U and V are square unitary matrices of respective dimensions m × m and n × n. The matrix C is then obtained by computing (18)
where S^{+} denotes the pseudoinverse of S, with S^{+} = (S^{T} S)^{−1} S.
In this expression, the diagonal matrix S^{+} contains the inverse coefficients of the nonnull S coefficients. The smallest singular values of S correspond to poorly sensed modes or noisy artifacts that require large strokes on the DM to be controlled. They can be set to zero without excessive loss of information. The matrix S^{+} can therefore be replaced by its filtered version for which k represents, in the list of singular values, the rank of the last value that is left nonnulled. Correspondingly, C is replaced by C_{k} by substituting S^{+} with . This filtering helps to stabilize the control loop by avoiding the introduction of correction errors.
For a given image I_{C} within the ROI, we obtain the list of coefficients with (19)
These coefficients are associated with the k modes and allow us to produce a phase map to apply to the DM shape for modal control. In a reallife system, a control loop gain of smaller than 1 is used to ensure its convergence towards a null wavefront error with the control algorithm.
Fig. 3 Semiunfolded layout of the HiCAT testbed. The apodizer labeled in the top right corner of the diagram is currently replaced by a flat mirror. 
3 ZWFS Characterization
We characterize the behavior of the ZWFS in the presence of a FPMfiltered beam using the setup that was recently installed on HiCAT.
3.1 Experimental Setup on HiCAT
HiCAT is a highcontrast testbed for developing and maturing exoplanet imaging technologies for future large segmentedaperture telescopes. Its goal is to demonstrate a complete systemic approach, combining a segmented aperture, coronagraphy, and wavefront sensing and control to achieve a 10^{−8} contrast in air and in broadband light in the visible wavelength range. This testbed is installed in a cleanroom environment and inside an enclosure that is stabilized in humidity and temperature by a valvecontrolled air flow. Tuning it allows us to change the air turbulence strength. Several light sources are available on HiCAT and for our experiments we work in monochromatic light at the wavelength λ_{0} = 640 nm.
The current design of the testbed is presented in the partially unfolded diagram in Fig. 3. Among different components, HiCAT includes an IrisAO PTT111L deformable mirror with 37 segments to mimic a segmented aperture. By default, the IrisAO uses a calibrated flat map to minimize the WFE due to segmentation. The current starlight suppression system is a classical Lyot coronagraph (CLC) composed of a reflective FPM made up of a flat mirror with a central pinhole, and a Lyot Stop (LS). The imaging camera, a ZWO ASI178, is used to display the coronagraphic image of the source with a sampling of 10.5 pix per λ_{0}/D_{pup}, with D_{pup} denoting the entrance pupil diameter. A dedicated pupil plane is also available to insert an apodizer and implement the Apodized Pupil Lyot Coronagraph (APLC, e.g., Soummer 2005; N’Diaye et al. 2016), but in our experiment this component is replaced with a flat mirror. The FPM pinhole diameter is d_{FPM} = 455 µm, leading to an angular size of 8.52 λ_{0}/D_{pup}, with an entrance pupil diameter D_{pup} of 19.55 mm. The beam transmitted through the FPM pinhole is caught by offaxis parabolas (OAPs) to feed several sensors: a target acquisition (TA) camera in a focal plane for FPM centering, a phase retrieval (PR) camera, and a ZWFS for loworder aberration measurements.
The ZWFS mask consists of a fused silica substrate with an ionetched dimple in the center which is aligned with the source image in the associated focal plane. The Zernike dimple has a depth of 280 nm and a diameter of 54.3 µm, leading to an angular diameter d_{Z} = 1.02 λ_{0}/D_{pup}. The mask was manufactured by SILIOS Technologies using photolithography with a reactive ionetching process, as described in N’Diaye et al. (2010). The ZWFS camera is a ZWO ASI178 with 3096 × 2080 pixels and a pitch of 2.4 µm. The speed limit of the camera is 60 Hz in fullframe operations. As the source provides enough photons, the typical exposure time is 0.5 ms in our experiments.
To shape the wavefront, HiCAT uses two continuous Boston kiloDMs: DM1 in a pupil plane and DM2 out of a pupil plane to allow both phase and amplitude control. A prior calibration of HiCAT was performed to experimentally produce DM flat commands (or ‘flat maps’) that minimize the number of static aberrations on the testbed (Brady et al. 2018). HiCAT also uses both of these DMs to further reject starlight and enhance the contrast in a given region of the coronagraphic image using DH optimization algorithms. Optimal DM commands for DH generation are obtained with the stroke minimization (SM) algorithm (Pueyo et al. 2009; Mazoyer et al. 2018a,b). An example of typical DM commands producing a DH with this algorithm is shown in Fig. 4, with deflections ranging from approximately −50 nm to +50 nm. These commands produce a DH ranging from 7.6 to 12.8 λ_{0}/D_{pup} in the focal plane, as operated on HiCAT as of December 2021. When these commands were generated, they produced a DH with a 2.5 × 10^{−8} average contrast. The aim of our current work is to reduce the inner working angle (IWA) while keeping the same contrast level, but this is out of the scope of the present paper. Throughout this paper and as handled on HiCAT, we present DM commands as the OPD that we aim to introduce. The DM commands in volts are computed using the conversion factors from nm to volts that were measured with a Fizeau interferometer. These gains were computed using the voltage ratio between the flat map and a pure defocus map. These defocus map voltages have been calibrated with the Fizeau interferometer to produce a pure defocus. The expected and introduced OPD will mainly differ because of actuator coupling and nonlinearity for large strokes.
For our experiments, we operate our wavefront control loop using the ZWFS with two different DM settings as a starting point: (i) with the DM flat maps providing the optimal optical quality of the testbed and (ii) with the DM maps as provided by SM, producing a coronagraphic image DH. Both settings will allow us to monitor the wavefront and DH stabilization with our active control loop. When using DM SM maps, we apply a setandforget operation mode: we first run the SM algorithm to obtain the DM commands and then run the ZWFS experiment. The delay between the SM and ZWFS experiments is no longer than a few minutes so as to prevent excessive evolution of the testbed conditions in this time. The SM and ZWFS operations have so far been performed separately and successively. Further developments will enable us to combine SM control and ZWFS closedloop control simultaneously.
In the rest of the paper, J is computed using n at around 24 000 pixels, with the exact number depending on the selected ROI. The DM probes that introduce Zernike modes in OPD are applied on top of the flat maps in case (i) or on top of the SM commands in case (ii), described immediately above. From empirical tests to optimize the control, we use p = 20 nm RMS, m = 12 and we only filter out the smallest singular value to keep k = 11 modes. We use the Noll convention (Noll 1976) for Zernike modes, poking the first 12 modes from piston to vertical secondary astigmatism. The DM response to piston, which is not a pure piston, is poorly sensed by the ZWFS. Its sensitivity to a given mode can be computed as the variance of the coefficients of the corresponding column in J. In our experiments, the ZWFS sensitivity to piston is 18 times smaller than the averaged sensitivity to the other 11 modes. Similarly, it is five times smaller than the sensitivity to trefoil, which is the secondleastsensed mode. For better loop stability, this requires filtering at least one mode in the inversion of J. Piston is only kept in the probes within the control code for simplicity in operations. In our experiment, the requested pokes are pure Zernike modes. Their difference with the introduced modes is calibrated by J. Further studies will be required to assess the maximum number of Zernike modes to be controlled or to address alternative modes that may be more suitable for this control. Finally, for the control loop, the implemented controller is an integrator, with an integral gain of 1 and a proportional gain of 0.15.
Fig. 4 Example of DM commands on both continuous deformable mirrors produced by the SM algorithm that can generate a DH on HiCAT. The DM commands in units of OPD is an approximation of the introduced deflection. The black circle represents the active area of the continuous DMs. For scaling comparison, IrisAO segments are plotted in transparency. The segment footprints are not realistic here as the DM2 is out of the pupil plane. 
3.2 Modal Response Curves
Figure 5 presents the ZWFS response curves for different commands with Zernike modes when operated with a Jacobian matrix to investigate the validity range of this datadriven calibration. These response curves are obtained by multiplying the intensity measurement on the ZWFS camera with the control matrix for different aberration amplitudes of the introduced probes. We reiterate that the calibration probes used to build the matrix J are performed using Zernike modes with p = 20 nm RMS.
In Fig. 5, we compare the response curves to loworder Zernike modes in simulation and experiment for two base DM settings: the flat maps and the SM maps. The Zernike modes are generated in terms of requested OPD. Because of the linear approximation in the conversion from nm to volts and the unknown exact DM response, the exact introduced OPD is unknown. The ZWFS is calibrated with a Jacobian matrix built around the DM flat shapes. Except for tip and tilt aberrations, which we discuss later on, the simulation and experiment results show good agreement for aberration amplitudes smaller than 15nm RMS. The slope of the curves is very close to the unity response. For larger aberrations, the experimental curves show a mode underestimation of typically 10 nm RMS for an input of 40 nm RMS. This effect most likely comes from the quadratic deflectiontovoltage relation of the DMs (Morgan et al. 2019), which is not taken into account when applying DM commands on HiCAT. For openloop operations like the ones shown here, we reach the limit of this DM control that assumes a perfect DM response. This approximation remains nevertheless sufficient for closedloop control as detailed in Sect. 2.3.
In the scope of future operations with changing DM offsets, an important point is to analyze the ZWFS behavior when we add an offset to the DM commands that were used during the Jacobian matrix calibration. For most modes, and as predicted by simulations, the sensitivity of the sensor is decreased when we perform the response measurement with the SM commands presented in Fig. 4. Only the defocus and the coma show a slight increase. As for the tip and tilt modes, they present a larger loss in sensitivity even at small aberrations. This originates from the combination of two factors. First, the ZWFS robustness to these modes is known to be reduced compared to higher order modes, as they represent a ZWFS mask misalignment with respect to the PSF. Second, the testbed drifts due to the delay between the SM command computation and the measurement of the response are likely to induce these tips and tilts. As a result, applying the SM commands that are computed around a slightly different testbed state will lead to a more visible effect on tip and tilt than on other modes.
For most science cases, the ZWFS will be operating in a closed loop with aberration values smaller than 20nm RMS. In this regime, the sensor shows a behavior that should not lead to the degradation of wavefront corrections, as its response remains linear.
Fig. 5 Response curves of the ZWFS to loworder Zernike modes, in simulation (dotted lines) and from HiCAT experiments (solid lines with cross markers). The ZWFS has been calibrated around the DM flat maps that optimize wavefront quality. Response curves are measured both around these flat maps (blue curves) and around SM commands (red curves). 
3.3 Cross Talk
In addition to sensitivity, an important aspect of the sensor is crosstalk between sensed modes: specifically, it unclear whether or not some modes trigger a response on other modes. Figure 6 shows the mode coupling matrix of the ZWFS for the first ten measured Zernike modes. While we have k = 11, only the first ten modes are shown for the sake of simplicity. These measurements come from the same data set as Fig. 5 with p = 20nm RMS. For both experiments, the global matrix shape is diagonal, which is consistent with the fact that the probes are identical to those used for the construction of J. Outside the diagonal and for both experimental cases, defocus is measured as a combination of tip, tilt, and both comas. We attribute this effect to a slight misalignment in the testbed that breaks the symmetry of the optical system, probably in the relay optics between the FPM and the ZWFS. For instance, a slight offaxis defocus could produce this effect. As the nonlinear response of the DMs to voltage is not taken into account, modal crosstalk can also be expected when adding an offset to the commands. Further tests are required to precisely identify the origin of this effect.
The change of DM shape offsets maintains the similarity of the coupling matrix between the two experiments, which remains mostly diagonal. The main differences are created by the responses of the tip and tilt. These two modes show specifically a stronger coupling with both coma modes. The defocus also shows a coupling with both astigmatism modes that was not present with the offset of the flat maps. Overall, these limited changes due to DM offset modifications will unlikely impact the stability of the control loop, as we show in the following section.
4 Stability Results
In this section, we investigate the wavefront and DH contrast stability with the ZWFS loop for different conditions of wavefront perturbations.
4.1 WFE Calibration and Stability
HiCAT is operating in air and naturally suffers from air turbulence, in particular in the optical path before the FPM. The air turbulence contribution is mainly composed of loworder aberrations, which allows us to test the ZWFS loop on an unknown perturbation and correct for them with DM1. With the air valve, we can adjust the strength of the dry air flow and open or close the testbed enclosure in order to generate different levels of air turbulence.
We study the wavefront statistics between an openloop configuration with only WFE monitoring and a closedloop configuration that includes WFE correction. We explore four different configurations, referring to the amount of turbulence. The first one, “closed”, corresponds to the case with the air valve completely closed, with no air flow, which causes the bench temperature and humidity to drift, but reduces the turbulence to a minimum. The second one, “low”, is related to a minimum air flow through the valve to stabilize the bench in humidity and temperature, which is the standard way of operation on HiCAT. The third one, “medium”, leads to the case with a small opening between the HiCAT enclosure covers and increased dry air flow. The last one, “high”, corresponds to the case with the removal of one lateral panel of the enclosure and dry air flow at a maximum. In this latter configuration, the PSF typically jitters around values of ±0.4 λ_{0}/D_{pup} at frequencies typically below 10 Hz. We note that these medium and high turbulence modes are never used in normal operations for HiCAT, and are only used here to increase the perturbation and study the performance of our closedloop ZWFS in suboptimal conditions.
In closed loop, the control is performed with a prebuilt Jacobian matrix, with m = 12 Zernike mode probes, with k = 11 modes being kept in the inversion. To avoid the effect of fitting error estimation on the whole WFE in the pupil as much as possible, the overall standard deviation of the wavefront is computed using the analytical phase reconstruction described in Sect. 2.2. Figure 7 shows the wavefront error estimates for the four different configurations of air turbulence on the testbed in open and closed loop. The closed loop runs at a typical frequency of 20–30 Hz, limited by the current computer calculation speed. The differences in perturbation between all the situations are clearly visible in both open loop before 600 s, and in closed loop after 600 s. For a closed air valve, a drift is present in open loop, degrading the tip mode; see top plot. This effect is due to a thermal drift of the testbed: the dry air is used to keep the temperature and humidity steady. With restored air flow, this effect is not visible in the other configurations. The wavefront fluctuations become larger as the injected air turbulence increases.
Table 1 gives the temporal standard deviation σ for tip, tilt, and defocus during the four experiments. From our experience on the testbed, tip is the most perturbed mode by air turbulence, but the behavior of other modes is very similar. The air turbulence degrades σ_{tip} by a factor of 3.5, from 1.49 to 5.29 nm RMS, between the low and high turbulence configurations. In all four configurations, closing the loop improves the residual tip, with the gain on σ_{tip} ranging from 1.7 with the high turbulence to 9 without any air flow. In the most stable configuration, with no air flux, σ_{tip} drops from 2.17 to 0.24 nm RMS (and typically from ±0.2 λ_{0}/D_{pup} to ±0.05 λ_{0}/D_{pup}), reaching subnanometric stability with our control loop. Improvements are also noticed with the other modes, confirming the ability of our control loop to correct for tip, tilt, and defocus and enhance the stability of the first Zernike modes in all configurations. A similar behavior is also observed with higher order modes as illustrated in the following.
We also performed a temporal spectral analysis of our results; see Fig. 8. The plot shows the openloop (blue) and closedloop (red) power spectrum densities (PSDs) as a function of the temporal frequency ƒ for the first three Zernike modes in the four different air turbulence regimes. All openloop experiments show a characteristic decrease in the PSD, which is consistent with the ƒ^{−5/3} Kolmogorov statistics before reaching the noise floor of the higher temporal frequencies. The turbulence regime alters the frequency cutoff of this behavior from around 0.2 Hz for the closed valve case to several Hertz for the high turbulence regime. The PSDs in the openloop operation do deviate from this ƒ^{−5/3} behavior in the low temporal frequency regime, most likely because of the bench enclosure that limits the turbulence developments at low spatial and temporal frequencies.
Still, all of the four closedloop cases show an improvement of the temporal PSD in the lowfrequency regime by two orders of magnitude from open to closed loop. For the closed and low air valve cases, the closedloop PSDs are almost flat over all frequencies and reach the noise floor, meaning that the correction loop is close to the optimal correction in these conditions. The correction could be further improved by enhancing the signaltonoise ratio in the measurements. This could be done by binning the images on the ZWFS camera, currently of 224 × 224 pixels in size after 4 × 4 binning, or by using temporal averages of several frames, at the cost of a lower loop speed, which is currently around 20–30 Hz. For medium and highturbulence regimes, the closedloop PSDs show an increase from the lowest frequencies to 1 Hz at which the blue and red curves meet, showing that the loop struggles to correct for the aberrations with frequencies higher than 0.1 Hz. A larger gain, a faster loop speed, or an improved control method such as a proportionalintegralderivative (PID) controller or predictive control could be useful for improving the stability in these conditions. Nevertheless, this high level of air turbulence is never met in standard HiCAT operation. Overall, our control loop using ZWFS improves the temporal PSDs of the wavefront errors in all conditions, showing that it is a suitable solution for controlling and stabilizing loworder modes for highcontrast instruments with LCs.
Fig. 6 Mode coupling matrix of the ZWFS to Zernike modes using a Jacobian matrix calibrated with Zernike modes around flat DMs. Introduced and measured modes are given on the horizontal and vertical axes. Left: with flat DMs. Right: after introduction of SM offset on the DMs. Modes are introduced with 20 nm RMS in each case. 
Temporal WFE standard deviation in open and closed loop in nanometers, and the factor of improvement between both configurations for the tip, tilt, and defocus modes and for each air turbulence configuration.
Fig. 7 Temporal evolution of the wavefront error with different levels of air turbulence. Global WFE measured by secondorder reconstructor is plotted in magenta. Tip, tilt, and defocus coefficients estimated by the inverse Jacobian matrix are plotted in blue, orange, and green, respectively. The plots are ranked from top to bottom from no external air flow to maximum air flow. The experiment consisted of 15 minutes of open loop, followed by 15 minutes of closed loop, under natural drifting conditions. Nominal operations for HiCAT correspond to “closed” and “low”, and the “medium/high” configurations are only used to study the performance of the ZWFS in a suboptimal regime. 
4.2 Contrast Stability
The ultimate goal of correction for loworder aberrations is to help maintain the optimal contrast in the DH. We explore the contrast evolution in a setandforget mode: we use the SM algorithm to compute optimal DM commands that allow for a DH at a few 10^{−8} contrast, apply these commands on the DMs, and let the DH contrast drift, with or without the ZWFS loop running. With the current large IWA (7.6 λ_{0}/D_{pup}), the DH contrast is dominated by mid to highorder aberrations. The impact of the air turbulence and the ZWFS correction is therefore challenging to observe. To artificially increase the contrast drift over the experiment duration of 30 min, we introduce Brownian noise on DM1. This random walk is performed every other second by drawing random Zernike coefficients with a uniform law in the range [−0.5, 0.5] nm RMS for the first ten Zernike modes. The average WFE of each individual perturbation is 0.79 nm RMS and the perturbations are introduced on DM1 in a cumulative way. We monitor the wavefront errors and the DH contrast first in open loop and then in close loop for 900 s each. The upper plot in Fig. 9 shows examples of wavefront maps measured with our ZWFS during the experiment. A clear degradation of the wavefront error can be observed during the openloop experiment with a residual WFE of 3.7, 16.5, and 19.2nm RMS at 0, 400, and 800 s respectively. After closing the loop at 900 s, the wavefront errors show a pattern very similar to the initial one at Os, with a WFE of 3.7 nm RMS at 1780 s, showing a good recovery of the initial wavefront state.
The evolution of the corresponding first ten Zernike coefficients estimated with the command matrix C are presented in the upper plot of Fig. 10 for the calibration around SM commands. As expected, the coefficients show improved stability after closing the loop. The temporal standard deviations are given for the first ten modes in Table 2. Once the loop is closed, the noisiest mode is tilt with σ_{tilt} = 760 pm, the mean value of σ being 380 pm between all the modes. Overall, the wavefront error from the ten Zernike coefficients drops from 10.4 to 1.4 nm RMS, showing an improvement by a factor of 7.6. This improvement is confirmed by the analysis of the temporal PSDs for the first ten Zernike coefficients; see Fig. A.1. A gain of a few orders of magnitude is observed in the lowtemporal frequency regime from open to closedloop operation.
The bottom plot in Fig. 9 displays focalplane DH images corresponding to the wavefront measured by our ZWFS. The contrast degradation is visible at 400 and 800 s, but the ZWFS closed loop restores a DH qualitatively similar to the initial one. The bottom plot in Fig. 10 shows the contrast evolution for two different calibrations: one where J was computed around the SM commands (solid line) and one where it was computed around flat commands (dashed line). However, the SM commands applied during both experiments were the SM commands at all times.
Regarding the contrast, we reiterate that our goal is to study the maintenance of a precomputed DH with a given contrast and shape: the ZWFS control loop is not designed to improve the ultimate contrast or DH digging capabilities. As of December 2021, the capabilities of HiCAT allow the generation of a DH with an average contrast of ~2.5 × 10^{−8} while SM is running. Once we interrupt the SM loop, the environmental testbed conditions and lack of any further action on the testbed mean that the contrast will tend to remain at twice this value, around ~5 × 10^{−8}, over a few tens of minutes. As such, we are first and foremost limited by the contrast in the starting DH. After introducing some known perturbations that lead to a contrast degradation, we aim to recover the initial ~5 × 10^{−8} contrast thanks to our ZWFSbased control loop.
The spatially averaged contrast values in the DH are given in Table 3. While in open loop, the contrast decreases by a factor of 10 from 5 × 10^{−8} to 5 × 10^{−7} because of the introduced wavefront drifts during the experiment; closing the loop brings the contrast back to approximately its initial value of around 5–6 × 10^{−8}. In this experiment, the average closedloop contrast for the ZWFS calibration around SM commands and flat maps are 7.6 × 10^{−8} and 6.2 × 10^{−8}, respectively. We achieve a contrast standard deviation σ_{contгast} that goes from 2.2 × 10^{−7} in open loop to around 7 × 10^{−9} in closed loop for both calibrations, which corresponds to an improvement by a factor of 30. We identify at least two reasons to explain the small difference between the initial (5 × 10^{−8}) and recovered (6.2 and 7.5 × 10^{−8}) contrasts. First, as perturbations are still introduced during the closedloop operation, the lag in the control loop will cause some residual WFE. While our loop is running at 20–30 Hz, a faster loop will help to reduce the lag and this residual effect. Second, while we correct for loworder modes, there are also mid and highorder aberration drifts that remain uncontrolled here, degrading the contrast DH with speckles. Although the spatially filtered ZWFS approach will only correct for loworder aberrations, it can be used in conjunction with other highorder techniques to make the most efficient use of the rejected stellar photons.
In Fig. 10 (bottom), the difference between the two curves is due to realization effects: the SM commands calculated at a given time do not provide optimal DH contrast at a different point in time. The achieved contrast also greatly depends on the bench evolution between the moment the SM commands are computed and the acquisition of the reference measurement for the ZWFS. This is clearly visible in Fig. 9: the initial DH, at 0 s, shows some speckles that are usually removed by SM. Still, the contrast degradation at 800 s due to the wavefront drift is canceled after the closedloop control.
Overall, we show that the stabilization of loworder aberrations on HiCAT is possible down to a subnanometric level and we demonstrate the efficiency of our control loop with ZWFS to prevent loworder WFE from degrading the contrast during observations. This correction has been performed in a blind way without interactions between the different control loops: our control loop based on the ZWFS uses no information from the imaging camera and it therefore stabilizes the image DH without seeing it. Further work involving synchronization with SM and running parallel loops on HiCAT will be investigated to stabilize the testbed at its best contrast performance.
Fig. 8 Temporal PSDs for the tip, tilt, and defocus in different air turbulence conditions, in openloop (blue) and closedloop (red) operations. 
Fig. 9 Examples of wavefront differences from reference as reconstructed by the ZWFS (top) and of DH images (bottom) during the contrast stabilization experiment with the calibration around flat maps. The iterations at 0, 400, and 800 s are from the openloop experiment. The loop is closed after 900 s. In this experiment, we are adding “artificial drifts” in the form of loworder aberration random walk added to the continuous DM. At the end of the openloop phase (900 s), the added aberration is about 20 nm RMS. At the end of the closedloop phase (1800 s), the added aberration is about 25 nm RMS. This experiment also corresponds to the nominal turbulence regime (“low”). 
Fig. 10 Wavefront evolution and contrast in the DH as a function of time in openloop (first 900 s) and closedloop (last 900 s) configurations under artificial drift conditions. An additional random perturbation of ten Zernike modes combined with 0.79 nm RMS on average is applied every 2 s. This amounts to a total aberration of about 25 nm RMS at the end of the 30min experiment. Top: Zernike coefficients estimated with the inverse Jacobian matrix calibrated around SM commands. These have been artificially shifted along the vertical axis to enhance readability: all of them actually oscillate around zero in closed loop. Bottom: mean contrast in the DH in the case of a Jacobian matrix calibration performed around SM commands (continuous line) and around flat DM maps (dashed line). 
Openloop and closedloop temporal wavefront error standard deviations in nanometers for the first ten Zernike modes controlled during the DH stabilization experiment, and the gain between both operations.
Statistics of the contrast spatial average in the DH, ranging from 7.6 λ_{0}/D_{pup} to 12.8 λ_{0}/D_{pup}, during the stabilization experiments for flat and DH calibrations.
5 Conclusion and Perspectives
Here, we validate an approach to control loworder wavefront errors in Lyottype coronagraphs using a ZWFS with the light rejected by the coronagraphic FPM. In this work, the existing formalism of the ZWFS is extended for wavefront reconstruction with a FPMfiltered beam.
To validate our approach experimentally, we used the HiCAT testbed in such a configuration. We characterized the ZWFS response to Zernike modes when calibrated with a modal Jacobian matrix. The ZWFS shows a good response to loworder Zernike modes and limited crosstalk between modes. Changing the offsets on the DMs to generate a DH with a 2.5 × 10^{−8} contrast and a 7.6 λ_{0}/D_{pup} IWA in the focal plane barely impacts the ZWFS response. With deeper contrasts and smaller IWAs, the DH DM commands may have larger strokes. Additional tests are still necessary to extend the validity of our findings.
We validated the modal control of loworder aberrations on HiCAT by closing the loop between the ZWFS and DM1. In different air turbulence regimes, we show a clear wavefront correction and stabilization down to a temporal standard deviation of a few hundred picometers for the tip, tilt, and defocus modes. In the presence of artificial random drifts introduced on DM1, our control loop is able to blindly stabilize the DH contrast at a level of 7 × 10^{−8} with a standard deviation of 7 × 10^{−9}.
Further investigations are required to demonstrate full operational capability. In this paper, only 11 Zernike modes were controlled. In the case of HiCAT, with a FPM mask cutoff radius of 4.2λ_{0}/D_{pup}, more Zernike modes could be measured. Regarding wavefront stability, controlling more modes is possible in general. However, these higher order modes include a higher spatial frequency content that is filtered by the FPM and therefore not measurable by this ZWFS. Unconstrained by the control loop, these higher order components add aberrations that create light leakage in the DH, which is no longer stable. An improved modal control would have to be used with modes that can better describe the system, as envisaged by Laginja et al. (2022), or with KarhunenLoeve modes. Furthermore, we currently use a monochromatic light source with enough photons for a high signaltonoise ratio. A more realistic approach would be to validate the ZWFS in broadband light and in a photonstarved regime.
The experiments were conducted on HiCAT with a segmented aperture that was flattened throughout with a calibration map. A crucial step will be to study the response of the ZWFS to segment piston, tip, and tilt aberrations, and correction in combination with the loworder modes.
Finally, we worked in a setandforget setting for the DM SM commands. The next step will consist of operating ZWFS while running SM and possibly other experiments in parallel. Our wavefront control loop currently runs at 30 Hz, limited by the camera readout speed and the computer processing time. The speed could be improved with a faster camera and software improvements.
This first demonstration of a control loop with ZWFS in the rejected light of a CLC and in the presence of a segmented aperture proves very promising to prepare for space coronagraphs to observe exoEarth candidates in line with the recent recommendations of the 2020 Decadal Survey for a future segmented space observatory with highcontrast imaging and spectroscopy.
Acknowledgements
R.P. acknowledges PhD scholarship funding from Région ProvenceAlpesCôte d’Azur and Thales Alenia Space. The authors are especially thankful to the extended HiCAT team (over 50 people) who have worked over the past several years to develop this testbed. This work was supported in part by the National Aeronautics and Space Administration under Grant 80NSSC19K0120 issued through the Strategic Astrophysics Technology/Technology Demonstration for Exoplanet Missions Program (SATTDEM; PI: R. Soummer). E.H.P. was supported by the NASA Hubble Fellowship grant HSTHF251467.001A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS526555.
Appendix A Temporal PSD for the loworder Modes During the DH Maintenance Experiment
Figure A.1 shows the temporal PSDs corresponding to the experiment presented in Sect. 4.2. The impact of the random walk is visible in the openloop operation (blue curves) when comparing the PSDs of the tip, tilt, and defocus in this plot and in Fig. 8 with the “low” air configuration. Observed at the lowest temporal frequencies, the PSD peaks reach 10^{−1}nm^{2}/Hz with the air turbulence only and about 1 nm^{2}/Hz with the additional perturbations introduced on DM1. Overall, with the added perturbations, all the modes in Figure A.1 present larger PSD values than the three modes shown in Fig. 8, at all temporal frequencies, as expected.
Fig. A.1 Temporal PSDs for ten Zernike modes from tip to spherical aberration, in openloop (blue) and closedloop (red) operations. The data are related to the experiment presented in Sect. 4.2. 
In closedloop operations, the final performance in Figure A.1 and 8 is similar with and without perturbations. Most modes present PSDs close to the noise floor around 10^{−5} nm^{2}/Hz. The cutoff frequency of the control loop appears to be the same, at around 1 Hz. As a result, Fig. A.1 shows a gain of up to almost five orders of magnitude in the PSDs at low temporal frequencies from open to closedloop operation. Tip, coma 90, trefoil 30, and spherical aberrations present a slight bump around 10^{−1} Hz, which is possibly due to the fact that perturbations are introduced every other second. This analysis in the presence of perturbations validates the ability of our control loop with ZWFS to recover the aberration level that is observed in the absence of introduced perturbations.
References
 Beuzit, J.L., Vigan, A., Mouillet, D., et al. 2019, A&A, 631, A155 [Google Scholar]
 Bloemhof, E.E., & Wallace, J.K. 2003, SPIE Conf. Ser., 5169, 309 [NASA ADS] [Google Scholar]
 Brady, G.R., Moriarty, C., Petrone, P., et al. 2018, SPIE Conf. Ser., 10698, 106986I [NASA ADS] [Google Scholar]
 Chambouleyron, V., Fauvarque, O., Sauvage, J.F., et al. 2021, A&A, 650, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chauvin, G., Desidera, S., Lagrange, A.M., et al. 2017, A&A, 605, L9 [Google Scholar]
 Coyle, L.E., Knight, J.S., Pueyo, L., et al. 2019, SPIE Conf. Ser., 11115, 111150R [NASA ADS] [Google Scholar]
 Currie, T., Guyon, O., Lozi, J., et al. 2020, SPIE Conf. Ser., 11448, 114487H [NASA ADS] [Google Scholar]
 Dohlen, K. 2004, EAS Pub. Ser., 12, 33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fanson, J., McCarthy, P.J., Bernstein, R., et al. 2018, SPIE Conf. Ser., 10700, 1070012 [NASA ADS] [Google Scholar]
 Fitzgerald, M., Bailey, V., Baranec, C., et al. 2019, BAAS, 51, 251 [NASA ADS] [Google Scholar]
 Gaudi, B.S., Seager, S., Mennesson, B., et al. 2020, ArXiv eprints [arXiv:2881.86683] [Google Scholar]
 Guyon, O. 2005, ApJ, 629, 592 [NASA ADS] [CrossRef] [Google Scholar]
 Guyon, O., Bottom, M., Chun, M., et al. 2019, BAAS, 51, 203 [Google Scholar]
 Jovanovic, N., Martinache, F., Guyon, O., et al. 2015, PASP, 127, 890 [Google Scholar]
 Jovanovic, N., Absil, O., Baudoz, P., et al. 2018, SPIE Conf. Ser., 10703, 107031U [NASA ADS] [Google Scholar]
 Kasper, M., Cerpa Urra, N., Pathak, P., et al. 2021, The Messenger, 182, 38 [Google Scholar]
 Keppler, M., Benisty, M., Müller, A., et al. 2018, A&A, 617, A44 [Google Scholar]
 Krist, J., Nemati, B., Zhou, H., & Sidick, E. 2015, SPIE Conf. Ser., 9605, 960505 [NASA ADS] [Google Scholar]
 Laginja, I., Sauvage, J.F., Mugnier, L.M., et al. 2022, A&A, 658, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Leboulleux, L., N'Diaye, M., Riggs, A.J.E., et al. 2016, SPIE Conf. Ser., 9904, 99043C [NASA ADS] [Google Scholar]
 Macintosh, B., Graham, J.R., Ingraham, P., et al. 2014, Proc. Natl. Acad. Sci., 111, 12661 [NASA ADS] [CrossRef] [Google Scholar]
 Macintosh, B., Graham, J.R., Barman, T., et al. 2015, Science, 350, 64 [NASA ADS] [CrossRef] [Google Scholar]
 Mazoyer, J., Pueyo, L., N'Diaye, M., et al. 2018a, AJ, 155, 7 [Google Scholar]
 Mazoyer, J., Pueyo, L., N'Diaye, M., et al. 2018b, AJ, 155, 8 [Google Scholar]
 Morgan, R., Douglas, E., Allan, G.W., et al. 2019, Micromachinesm, 10, 366 [NASA ADS] [CrossRef] [Google Scholar]
 Moriarty, C., Brooks, K., Soummer, R., et al. 2018, SPIE Conf. Ser., 10698, 1069853 [NASA ADS] [Google Scholar]
 National Academies of Sciences, Engineering, & Medicine 2021, Pathways to Discovery in Astronomy and Astrophysics for the 2020s (Washington, DC: The National Academies Press) [Google Scholar]
 N'Diaye, M., Dohlen, K., Cuevas, S., et al. 2010, A&A, 509, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 N'Diaye, M., Choquet, E., Pueyo, L., et al. 2013a, SPIE Conf. Ser., 8864, 88641K [NASA ADS] [Google Scholar]
 N'Diaye, M., Dohlen, K., Fusco, T., & Paul, B. 2013b, A&A, 555, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 N'Diaye, M., Choquet, E., Pueyo, L., et al. 2014, AAS Meeting Abs., 223, 149.13 [Google Scholar]
 N'Diaye, M., Mazoyer, J., Choquet, É., et al. 2015, SPIE Conf. Ser., 9605, 96050I [NASA ADS] [Google Scholar]
 N'Diaye, M., Vigan, A., Dohlen, K., et al. 2016, A&A, 592, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nielsen, E.L., De Rosa, R.J., Macintosh, B., et al. 2019, AJ, 158, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Noll, R.J. 1976, J. Opt. Soc. Am., 207 [NASA ADS] [CrossRef] [Google Scholar]
 Oppenheimer, B.R., Beichman, C., Brenner, D., et al. 2012, SPIE Conf. Ser., 8447, 844720 [NASA ADS] [Google Scholar]
 Por, E.H. 2020, ApJ, 888, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Pourcelot, R., Vigan, A., Dohlen, K., et al. 2021, A&A, 649, A170 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pueyo, L., Kay, J., Kasdin, N.J., et al. 2009, Appl. Opt., 48, 6296 [NASA ADS] [CrossRef] [Google Scholar]
 Ramsay, S., Amico, P., Bezawada, N., et al. 2020, SPIE Conf. Ser., 11447, 114471U [NASA ADS] [Google Scholar]
 Ruane, G., Riggs, A., Coker, C.T., et al. 2018, SPIE Conf. Ser., 10698, 106984U [NASA ADS] [Google Scholar]
 Ruane, G., Wallace, J.K., Steeves, J., et al. 2020, J. Astron. Teles. Instrum. Syst., 6, 045005 [NASA ADS] [Google Scholar]
 Seo, B.J., Patterson, K., Balasubramanian, K., et al. 2019, SPIE Conf. Ser., 11117, 111171V [NASA ADS] [Google Scholar]
 Shi, F., Balasubramanian, K., Hein, R., et al. 2016, J. Astrono. Teles. Instrum. Syst., 2, 011021 [NASA ADS] [CrossRef] [Google Scholar]
 Shi, F., Seo, B.J., Cady, E., et al. 2018, SPIE Conf. Ser., 10698, 106982O [NASA ADS] [Google Scholar]
 Simard, L., Ellerbroek, B., Bhatia, R., Radovan, M., & Chisholm, E. 2016, SPIE Conf. Ser., 9908, 99081V [NASA ADS] [Google Scholar]
 Sirbu, D., Marx, D., Belikov, R., et al. 2021, SPIE Conf. Ser., 11823, 118230R [NASA ADS] [Google Scholar]
 Sivaramakrishnan, A., Oppenheimer, B.R., Soummer, R., et al. 2007, in In the Spirit of Bernard Lyot: The Direct Detection of Planets and Circumstellar Disks in the 21st Century, ed. P. Kalas (Berkeley, CA, USA: University of California), 21 [Google Scholar]
 Soummer, R. 2005, ApJ, 618, L161 [Google Scholar]
 Soummer, R., Brady, G.R., Brooks, K., et al. 2018, SPIE Conf. Ser., 10698, 106981O [NASA ADS] [Google Scholar]
 Spergel, D., Gehrels, N., Baltay, C., et al. 2015, ArXiv eprints [arXiv:1583.83757] [Google Scholar]
 Steeves, J., Wallace, J.K., Kettenbeil, C., & Jewell, J. 2020, Optica, 7, 1267 [NASA ADS] [CrossRef] [Google Scholar]
 The LUVOIR Team. 2019, ArXiv eprints [arXiv:1912.86219] [Google Scholar]
 Vigan, A., & N'Diaye, M. 2018, Astrophysics Source Code Library [record ascl:1806.003] [Google Scholar]
 Vigan, A., Fontanive, C., Meyer, M., et al. 2021, A&A, 651, A72 [EDP Sciences] [Google Scholar]
 Wallace, J.K., Rao, S., JensenClem, R.M., & Serabyn, G. 2011, SPIE Conf. Ser., 8126, 81260F [NASA ADS] [Google Scholar]
 Zernike, F. 1934, MNRAS, 94, 377 [Google Scholar]
All Tables
Temporal WFE standard deviation in open and closed loop in nanometers, and the factor of improvement between both configurations for the tip, tilt, and defocus modes and for each air turbulence configuration.
Openloop and closedloop temporal wavefront error standard deviations in nanometers for the first ten Zernike modes controlled during the DH stabilization experiment, and the gain between both operations.
Statistics of the contrast spatial average in the DH, ranging from 7.6 λ_{0}/D_{pup} to 12.8 λ_{0}/D_{pup}, during the stabilization experiments for flat and DH calibrations.
All Figures
Fig. 1 Schematic diagram of the ZWFS downstream of a binary amplitude FPM. The light in the entrance pupil comes from the telescope and deformable mirrors in plane 0. It is focused on the FPM, which is reimaged on the Zernike mask in plane B. The intermediate pupil plane A is not required. 

In the text 
Fig. 2 Example of intensity map as measured by the ZWFS (left), normalized by the P^{2} measurement (middle), and the reconstructed phase map using the secondorder analytical reconstructor (right). 

In the text 
Fig. 3 Semiunfolded layout of the HiCAT testbed. The apodizer labeled in the top right corner of the diagram is currently replaced by a flat mirror. 

In the text 
Fig. 4 Example of DM commands on both continuous deformable mirrors produced by the SM algorithm that can generate a DH on HiCAT. The DM commands in units of OPD is an approximation of the introduced deflection. The black circle represents the active area of the continuous DMs. For scaling comparison, IrisAO segments are plotted in transparency. The segment footprints are not realistic here as the DM2 is out of the pupil plane. 

In the text 
Fig. 5 Response curves of the ZWFS to loworder Zernike modes, in simulation (dotted lines) and from HiCAT experiments (solid lines with cross markers). The ZWFS has been calibrated around the DM flat maps that optimize wavefront quality. Response curves are measured both around these flat maps (blue curves) and around SM commands (red curves). 

In the text 
Fig. 6 Mode coupling matrix of the ZWFS to Zernike modes using a Jacobian matrix calibrated with Zernike modes around flat DMs. Introduced and measured modes are given on the horizontal and vertical axes. Left: with flat DMs. Right: after introduction of SM offset on the DMs. Modes are introduced with 20 nm RMS in each case. 

In the text 
Fig. 7 Temporal evolution of the wavefront error with different levels of air turbulence. Global WFE measured by secondorder reconstructor is plotted in magenta. Tip, tilt, and defocus coefficients estimated by the inverse Jacobian matrix are plotted in blue, orange, and green, respectively. The plots are ranked from top to bottom from no external air flow to maximum air flow. The experiment consisted of 15 minutes of open loop, followed by 15 minutes of closed loop, under natural drifting conditions. Nominal operations for HiCAT correspond to “closed” and “low”, and the “medium/high” configurations are only used to study the performance of the ZWFS in a suboptimal regime. 

In the text 
Fig. 8 Temporal PSDs for the tip, tilt, and defocus in different air turbulence conditions, in openloop (blue) and closedloop (red) operations. 

In the text 
Fig. 9 Examples of wavefront differences from reference as reconstructed by the ZWFS (top) and of DH images (bottom) during the contrast stabilization experiment with the calibration around flat maps. The iterations at 0, 400, and 800 s are from the openloop experiment. The loop is closed after 900 s. In this experiment, we are adding “artificial drifts” in the form of loworder aberration random walk added to the continuous DM. At the end of the openloop phase (900 s), the added aberration is about 20 nm RMS. At the end of the closedloop phase (1800 s), the added aberration is about 25 nm RMS. This experiment also corresponds to the nominal turbulence regime (“low”). 

In the text 
Fig. 10 Wavefront evolution and contrast in the DH as a function of time in openloop (first 900 s) and closedloop (last 900 s) configurations under artificial drift conditions. An additional random perturbation of ten Zernike modes combined with 0.79 nm RMS on average is applied every 2 s. This amounts to a total aberration of about 25 nm RMS at the end of the 30min experiment. Top: Zernike coefficients estimated with the inverse Jacobian matrix calibrated around SM commands. These have been artificially shifted along the vertical axis to enhance readability: all of them actually oscillate around zero in closed loop. Bottom: mean contrast in the DH in the case of a Jacobian matrix calibration performed around SM commands (continuous line) and around flat DM maps (dashed line). 

In the text 
Fig. A.1 Temporal PSDs for ten Zernike modes from tip to spherical aberration, in openloop (blue) and closedloop (red) operations. The data are related to the experiment presented in Sect. 4.2. 

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.