Spatial linear dark field control: stabilizing deep contrast for exoplanet imaging using bright speckles

Kelsey Miller; Olivier Guyon; Jared Males

doi:10.1117/1.JATIS.3.4.049002

30 October 2017 Spatial linear dark field control: stabilizing deep contrast for exoplanet imaging using bright speckles

Kelsey Miller, Olivier Guyon, Jared Males

Author Affiliations +

Journal of Astronomical Telescopes, Instruments, and Systems, Vol. 3, Issue 4, 049002 (October 2017). https://doi.org/10.1117/1.JATIS.3.4.049002

Abstract

Direct imaging of exoplanets requires establishing and maintaining a high-contrast dark field (DF) within the science image to a high degree of precision ( 10−10). Current approaches aimed at establishing the DF, such as electric field conjugation (EFC), have been demonstrated in the lab and have proven capable of high-contrast DF generation. The same approaches have been considered for the maintenance of the DF as well. However, these methods rely on phase diversity measurements, which require field modulation; this interrupts the DF and consequently competes with the science acquisition. We introduce and demonstrate spatial linear dark field control (LDFC) as an alternative technique by which the high-contrast DF can be maintained without modulation. Once the DF has been established by conventional EFC, spatial LDFC locks the high-contrast state of the DF by operating a closed loop around the linear response of the bright field (BF) to wavefront variations that modify both the BF and the DF. We describe the fundamental operating principles of spatial LDFC and provide numerical simulations of its operation as a DF stabilization technique that is capable of wavefront correction within the DF without interrupting science acquisition.

1. Introduction

In the last two decades, the existence of 3,498 exoplanets has been confirmed,¹ and in the upcoming era of 30-m class ground-based telescopes and new space-based observatories, there is the promise of discovery, even characterization, of many more exoplanets, including potentially Earth-like worlds. With such powerful capabilities on the horizon, it has become imperative to push stellar suppression technology to higher precision in order to directly image these exoplanets. Combined with coronagraphy, current speckle nulling techniques² are capable of creating the dark field (DF) in the science image where light from an exoplanet orbiting its star can be detected. In principle, these methods are capable of generating a DF with $10^{- 10}$ contrast, which would enable the imaging of Earth-like exoplanets,³ but to generate and maintain a DF at such contrast requires high-precision measurements of very small changes in the science image.

The DF is very susceptible to small dynamic aberrations in the beam path. Aberrations produce postcoronagraph stellar light leakage and result in a quasi-static speckle field in the image plane that limits the DF contrast.⁴ Conventional wavefront sensors (WFS) that branch off from the main beam path cannot correct for noncommon path errors, and the resulting quasistatic aberrations can eventually become larger than the residual dynamic wavefront errors after correction by the WFS control loop. For high-contrast imaging, the WFS should ideally be common path with all of the optics seen by the science instrument to allow access to all aberrations created in the science beam. This is achieved by focal plane wavefront sensing (FPWFS), which uses the science detector as the WFS to measure the exact aberrations seen by the main science beam.⁵^,⁶ FPWFS techniques like speckle nulling and electric field conjugation (EFC) that have proven capable of generating a DF with high contrast in the lab have also been under consideration for maintenance of the DF.⁷^–⁹ As a control method, FPWFS presents its own set of challenges given that speckles have a quadratic relationship with aberrations. These techniques also rely on phase diversity measurements of the field at the science detector, which require field modulation and multiple images.¹⁰^,¹¹ This field modulation at the science detector, induced by a deformable mirror (DM), throws stellar light back into the DF and disrupts the science measurement. This interruption, which is required to rebuild the DF every time the contrast degrades, fundamentally limits the integration time that can be spent on any given target. The duration of this interruption to the science acquisition is directly related to the contrast of the DF. For deeper contrast, the required exposure time to sense the speckle field increases; therefore, at the $10^{- 10}$ contrast level, multiple images with long exposure times significantly reduce the amount of time that can be spent on observations. The need for modulation, multiple images, and long exposures, consequently makes the use of current speckle nulling methods and EFC nonideal for continuous maintenance of the DF.

Another technique known as the self-coherent camera (SCC) has been under development as a method for obtaining and maintaining the DF without science acquisition competition.⁶ While SCC does not require modulation, it does still utilize the mixing of some starlight with the DF. Linear dark field control (LDFC) does not require any such mixing of starlight with the DF and offers a potential solution for overcoming the limitations presented by speckle nulling and EFC. To avoid disrupting the science measurement with field modulation to rebuild the DF, LDFC locks the high-contrast state of the field once the DF has been constructed using conventional methods like EFC. Using only one image of the bright field (BF), LDFC freezes the state of the field by sensing and canceling changes in the wavefront that result in speckle formation in the image plane. The ability to maintain the DF with a single image yields a substantial increase in time that can be spent in the observation and analysis of exoplanets and will lead to an overall increase in the number of planets detected and analyzed over the lifetime of an instrument.

2. Concept

LDFC maintains high contrast without modulating the field to update the field estimate as is required when using EFC in closed loop. Conceptually, LDFC is inspired by low-order wavefront sensing (LOWFS), a proven technique designed to actively sense and correct low-order wavefront aberrations using light rejected by a reflective stop placed in a pupil plane¹² or focal plane¹³^,¹⁴ within the coronagraph. Instead of sensing only low-order aberrations using starlight rejected by the coronagraph, LDFC also has access to mid- and high-spatial frequencies. It is a FPWFS that operates a closed loop around starlight outside of the DF.

Spatial LDFC freezes the state of the DF by using state measurements of light spatially outside of the DF (see Fig. 1). This method uses the linear signal from the strongly illuminated BF to measure the change in the image plane intensity and uses that variation to calculate the correction required to return the image to its initial state, thereby stabilizing the DF. Unlike other FPWFS techniques, spatial LDFC does not rely on any induced modulation to derive an estimate of the field to be canceled. Instead, spatial LDFC observes the changes in the image intensity with respect to a reference image taken after the DF has been established by conventional speckle nulling methods. This process requires only the reference image and a single image taken at a later time.

Fig. 1

(a) Spatial LDFC: $3.5 λ / D \times 8.5 λ / D$ DF created in simulation using conventional EFC. (b) Wavefront aberrations produce speckles in the BF and the DF, which degrade the DF. The change in intensity between (b) the aberrated BF and (a) the ideal BF is used to measure and cancel the speckles in the DF.

Without the need for modulation or multiple images, spatial LDFC does not interrupt the science measurement, it decreases the time necessary to return the DF to its initial high-contrast state, and consequently it allows for longer, uninterrupted observing at high contrast. This paper introduces spatial LDFC as a more efficient alternative to conventional speckle nulling methods for stabilizing the DF. The theory behind spatial LDFC is laid out here in Sec. 2, and demonstrations of LDFC’s abilities in simulation are shown in Sec. 3. Further discussion of the limitations and null space of spatial LDFC is laid-out in Sec. 4 with concluding remarks in Sec. 5.

2.1.

Theory

Spatial LDFC relies on the linear response of the BF to wavefront perturbations that affect both the BF and the DF; this linearity allows for a closed loop control algorithm directly relating wavefront perturbations to changes in BF intensity. To derive the source of this linear response, we begin with the relationship between an incident wavefront and the resulting image. The complex amplitude of the incident wavefront in a pupil plane $E_{0}$ is linearly related to the complex amplitude at the image plane $E_{t}$ at a given time $t$ . The same linear relationship is true with respect to $E_{DM}$ , the multiplicative complex amplitude introduced by the DM in a conjugate pupil plane, and the complex amplitude at the image plane $E_{t}$ .

When the changes in optical path length (OPL) induced by the DM are very small such that $OPL ≪ 1$ , the resulting field $E_{t}$ at a given time $t$ in the image plane can be written as the sum of the initial pupil plane field $E_{0}$ and the small changes in complex amplitude induced in a conjugate pupil plane by the DM¹⁰

Eq. (1)

E_{t} \approx E_{0} + E_{DM} .

The resulting intensity in the image plane at time $t$ is then given by

Eq. (2)

I_{t} = {| E_{t} |}^{2} .

The total image plane intensity can be written as a sum of three terms: the intensity contribution from the initial pupil field: ${| E_{0} |}^{2}$ , the resulting intensity due to phase perturbations induced by the DM: ${| E_{DM} |}^{2}$ , and the inner product of the initial pupil field and the DM contribution to the complex amplitude

Eq. (3)

I_{t} \approx {| E_{0} |}^{2} + {| E_{DM} |}^{2} + 2 ⟨ E_{0}, E_{DM} ⟩ .

In the DF, the contribution of the initial field to the total intensity is very small, and the total intensity is dominated by the contribution of the DM such that

{| E_{DM} |}^{2} ≫ {| E_{0} |}^{2}

, thereby leading to a quadratic dependence of the DF on the DM input. However, in the BF, the contribution of the initial field to the total intensity dominates the contribution of the DM

Eq. (4)

{| E_{0} |}^{2} ≫ {| E_{DM} |}^{2} .

The intensity of the BF at the image plane at time $t$ can, therefore, be written as a linear function of the complex amplitude contribution of the DM

Eq. (5)

I_{t} \approx 2 ⟨ E_{0}, E_{DM} ⟩ + {| E_{0} |}^{2} .

In Eq. (5), the term

{| E_{0} |}^{2}

is the reference image

I_{ref}

taken after the DF has been established. The BF signal used by spatial LDFC to drive the DF back to its initial state is simply the difference between this reference and a single image

I_{t}

taken at time

t

Eq. (6)

Δ I_{t} = I_{t} - I_{ref} \approx 2 ⟨ E_{0}, E_{DM} ⟩ .

This linear response of the BF to field perturbations controlled by the DM is shown in Fig. 2 alongside the quadratic response of the DF to the same DM perturbation. In this figure, the BF and DF response to the DM field contribution

E_{DM}

is shown in a simulated PSF with a DF established by conventional EFC. A model of a MEMS DM was used to create the DF and then perturb the input wavefront by inducing a positive and negative delay in the optical path with a single actuator. This was done for a range of actuator amplitudes from

- 0.075

to

+ 0.075 μ m

.¹⁵ The resulting intensity response of pixels located in the DF [shown in Fig. 2(a)] is governed by Eq. (3) with the expected quadratic dependence on the field perturbation. The intensity response of the BF [shown in Fig. 2(b)] reveals the predicted linear dependence on the DM-induced field perturbation given by Eq. (6).

Fig. 2

The response of the BF and DF to the same range of DM “poke” amplitudes from $- 0.075$ to $+ 0.075 μ m$ on a single DM actuator. A demonstration of the (a) quadratic response of the DF and (b) the linear response of the BF to the same range of DM poke amplitudes. Each curve in the DF plot is the response of a single pixel in the DF between $7 λ / D$ and $8 λ / D$ , and each curve in the BF plot is the response of a single pixel in the BF between $10.5 λ / D$ and $11.5 λ / D$ .

In closed loop, $Δ I_{t}$ is small, and the linear approximation holds. However, even when initially closing the loop where $Δ I_{t}$ is larger, strict linearity is not required, only a monotonic trend. In instances both of strict linearity or monotonicity, the BF response allows for the construction of a linear servo driven solely by changes in the BF intensity. Unlike EFC and speckle nulling, which use modulation to provide an absolute field measurement, LDFC relies on relative measurements of BF intensity variation in the science image. These intensity variations are used to track and cancel changes in the wavefront that modify both the BF and DF, thereby stabilizing the DF without any disruptions to the science measurement.

2.2.

Calibration

Given the linear relationship between BF intensity and wavefront, using LDFC to stabilize the DF contrast is faster and more robust than using EFC. EFC requires multiple images to estimate the field, whereas each iteration of LDFC requires only one image at the science detector to determine how the field has changed with respect to the initial EFC-derived state. Since this image does not require field probing, which breaks the science measurement, the LDFC servo operates with a 100% duty cycle. Furthermore, LDFC does not rely on complex field estimates that require a model-based complex phase response matrix that is difficult to measure and verify; instead, LDFC relies only on a $DM \to image$ calibration that links a set of DM shapes, or basis functions, to changes in intensity in the science image.¹⁶

For this simulation, the DM influence functions were chosen as the basis functions. The calibration between the image and the basis set is obtained by building a response matrix $M$ whose columns relate the application of each individual influence function to the responding intensity variation at the science detector. Though modal control¹⁷ does offer performance benefits, especially when it maps with the expected temporal evolution of the wavefront error, such modal control tuning has not been explored at this time. For this work, application of the influence function basis set involved the actuation or “poking” of one of the $k$ actuators on the DM that lie within the illuminated system pupil. To begin building $M$ , the ideal reference image $I_{ref}$ with dimensions [ $n_{pix} \times n_{pix}$ ] is recorded after the DF has been established using EFC. To fill each of the $k$ columns in $M$ , a single actuator is poked, the resulting perturbed image $I_{k}$ is measured, and the unperturbed reference image $I_{ref}$ is subtracted off to yield the change in intensity. This is done for all $k$ actuators. The resulting response matrix $M$ has the dimensions [ $n_{pix}^{2} \times k$ ]

Eq. (7)

M [:, k] = I_{k} - I_{ref} .

The matrix

M

records the intensity change of both the BF and DF pixels with respect to each actuator poke. However, spatial LDFC uses only the BF pixels, which respond linearly to wavefront perturbations. The selection of these BF pixels relies on multiple parameters including background flux, flux per speckle, detector efficiency, and SNR. Based on these requirements, a threshold is applied to the initial EFC image

I_{ref}

, which selects only the

n

pixels with intensities greater than or equal to the threshold. The result is an image

I_{ref, n}

that records the initial EFC-state of only the BF pixels. An example of this BF reference image and the corresponding pixel map can be seen in Fig. 3 for a contrast threshold of

10^{- 4.5}

.

Fig. 3

Images of the applied (a) BF pixel mask and (b) the $\log_{10}$ masked reference image ( $I_{ref, n}$ ). The binary mask passes only the pixels at or above the contrast threshold (shown in white). In this image, and for the following demonstrations, the contrast threshold was $10^{- 4.5}$ , and the outer diameter of the masked control area was set to be the control radius of the active area on the DM.

To build the BF response matrix $M$ , the full response matrix is filtered to include only the $n$ BF pixels with intensities above the threshold such that $M = M_{n}$ . This filtered response matrix $M$ is used throughout the operation of spatial LDFC.

2.3.

Closed-Loop Implementation

To implement LDFC in closed-loop, an image $I_{t}$ is taken at time $t$ and the same $n$ BF pixels that pass the threshold are recorded in the BF image $I_{t, n}$ . The BF reference image $I_{ref, n}$ is then subtracted from the new BF image to track the changes that occurred in the BF with respect to the initial EFC BF reference (see Fig. 4)

Eq. (8)

Δ I_{t, n} = I_{t, n} - I_{ref, n} .

This BF intensity change

Δ I_{t, n}

is fit to the pseudoinverse of the BF response matrix

M

, also known as the control matrix, to calculate the DM shape that returns the field to its initial EFC reference state. The DM shape is represented by a vector of individual actuator amplitudes

u_{t}

Eq. (9)

u_{t} = - {(M^{T} M)}^{- 1} M^{T} Δ I_{t, n} .

This pseudoinverse of

M

is implemented by using singular value decomposition (SVD) and applying a threshold to filter out the modes that are not properly sensed by LDFC. For this simulation, the threshold value was chosen based on simulation performance, resulting in the inclusion of 286 out of an initial 398 modes in the pseudo-inversion process. A plot of the singular values of

M

is shown below in Fig. 5.

Fig. 4

The BF pixels that are used in the (a) reference and (b) aberrated images to measure the (c) intensity change $Δ I_{t, n}$ that drives the spatial LDFC control loop to stabilize the DF. In all three images, the $3.5 λ / D \times 8.5 λ / D$ DF can be seen to the right of the PSF core.

Fig. 5

Singular values of the spatial LDFC response matrix $M$ showing the applied SVD threshold (black) as well as the modes that were used in the inversion (blue) and the modes that were discarded (red). Out of 398 total modes, 286 were used for the inversion of $M$ in the following simulations.

In the following simulations, the response matrix $M$ and subsequent control matrix were measured once and applied in closed loop with an initial gain of 0.6. Once the DF contrast converged to $10^{- 7.9}$ , the gain was lowered to 0.1 to maintain the correction. The ensuing process of taking an image, calculating the intensity change of the BF from its reference state, and updating the DM was iterated on to actively freeze the science image field in its initial EFC state.

3. Simulation Results

To demonstrate spatial LDFC’s ability to maintain the high-contrast DF, a 6.5-m telescope system was constructed in simulation, which included a single DM and Lyot coronagraph that removed approximately two orders of magnitudes of stellar light from the final image. The system entrance pupil was a circular 6.5-m diameter, centrally obscured mask with a 30% central obscuration, and 2% spiders [see Fig. 6(a)]. The Lyot coronagraph consisted of a Lyot stop undersized by 1% and a focal plane mask with a diameter of $2.44 λ / D$ . For the system’s DM, a model of a Boston Micromachines 1K DM was defined using 1024 actuators sharing a common Gaussian influence function and 15% interactuator coupling. The diameter of the illuminated pupil projected onto the DM was 6.5 mm, covering $\sim 21$ actuators and lending an outer working angle (OWA), or control radius of $10.5 λ / D$ . Sampling at the science detector was $0.24 λ / D$ per pixel. The source was a magnitude 5 star with sensing done at $λ = 550 nm$ (V band) with 10% bandwidth. The total flux at the entrance pupil was $1.82 \times 10^{9} photons / s$ , and this rate was used to embed photon noise in all of the $I_{t}$ images in Eq. (8). All of these test parameters are listed in Table 1.

Fig. 6

(a) Telescope pupil and standard implementation of EFC using a DM with 398 illuminated actuators to create (b) a $3.5 λ / D \times 8.5 λ / D$ DF centered at $6.75 λ / D$ from the center of the stellar PSF with $10^{- 7.94}$ average contrast.

Table 1

Simulated system parameters used in the following spatial LDFC demonstrations.

Stellar magnitude	5
Total flux	$1.82 \times 10^{9} photons / s$
Noise included	Photon noise
Exposure time	5 s
Source wavelength	550 nm, V band
Source bandwidth	10%
Telescope diameter	6.5 m
Sampling at detector	$0.24 λ / D$ per pixel
# DM actuators used	398 (21 in diameter)
# BF pixels used	4535
BF contrast threshold	$10^{- 4.5}$
Inner working angle (IWA)	$2.44 λ / D$
Outer working angle (OWA)	$10.5 λ / D$

To build the DF, a standard implementation of EFC¹¹ was used to suppress the stellar light to an average contrast floor of $10^{- 7.94}$ within a $3.5 λ / D \times 8.5 λ / D$ region centered at $6.75 λ / D$ from the PSF core [shown in Fig. 6(b)]. This DF was the ideal reference state for LDFC to maintain, and the intensity image $I_{ref}$ was saved as the reference image to be used in the LDFC servo to return the DF to its EFC-derived state.

With the DF established, the spatial LDFC algorithm was implemented as described in Sec. 2.3 to maintain the DF in the presence of two separate injected phase aberrations. In the first case, a single speckle pair was induced in the image plane by applying a sine wave phase perturbation in the pupil. For the second case, a random Kolmogorov phase screen was introduced in the pupil creating multiple speckles in the image plane. In both cases, the same optical system, source, and threshold values were kept constant as were all other simulation parameters. The following sections present spatial LDFC’s response to these two cases.

3.1.

Sine Wave Phase Perturbation

After the DF was constructed, a spatial sine wave phase perturbation with 6 cycles/aperture was introduced into the pupil plane, forming a speckle at $\pm 6 λ / D$ : one speckle within the DF and one speckle within the BF. The sine perturbation was given a 1-nm peak-to-valley (P-V) amplitude in phase, creating a speckle pair with a maximum magnitude of $10^{- 5.0}$ and an average aberrated DF contrast of $10^{- 6.90}$ . The LDFC control loop was run with a gain of 0.6 until the average DF contrast reached $10^{- 7.9}$ , at which point the gain was reduced to 0.1 to maintain the correction. The LDFC control loop was allowed to run for 50 iterations for this demonstration with convergence occurring after six iterations. The results are shown in Table 2 and in Figs. 7 Fig. 8 Fig. 9–10. It should be noted that, in Figs. 9 and 10(b), the LDFC-corrected DF contrast occasionally drops below the initial EFC contrast level. This effect is due to noise fluctuations.

Table 2

Performance with a sine wave phase: initial EFC DF average contrast, magnitude of the injected speckle, average contrast of the aberrated DF, average contrast of the DF after LDFC, total change in contrast for one full LDFC loop, and the number of iterations to converge to the EFC contrast floor.

EFC DF contrast	$10^{- 7.94}$
Speckle magnitude	$10^{- 5.0}$
Avg. DF contrast with speckle	$10^{- 6.90}$
LDFC DF contrast	$10^{- 7.94}$
$Δ Contrast$	$10^{- 1.04}$
# Iterations to converge	6

Fig. 7

The injected 6 cycles/aperture sine wave phase perturbation with (a) a P-V amplitude of 1 nm (RMS: 0.4 nm), (b) the DM response derived by LDFC, and (c) the final residual wavefront error after six iterations. In this case, the residual WFE is dominated by a mode with a frequency of $\sim 14 cycles / aperture$ which falls beyond the spatial frequency limit ( $10.5 cycles / aperture$ ) the DM can correct. This residual WFE is due to the Gaussian shape of the DM’s influence functions which cannot perfectly fit the injected sine wave perturbation, thereby leaving a residual sinusoidal pattern. Scale is given in nm.

Fig. 8

(a) The aberrated PSF with a single $10^{- 5.0}$ magnitude speckle in the DF with $10^{- 6.90}$ average contrast and a matching speckle in the BF, (b) the final LDFC-corrected DF with $10^{- 7.94}$ average DF contrast, and (c) the reference EFC-derived DF also with $10^{- 7.94}$ average DF contrast. Scale is $\log_{10}$ contrast.

Fig. 9

Evolution of the DF over the six iterations [seen in Fig. 10(b)] to converge from a degraded DF average contrast of $10^{- 6.90}$ with a $10^{- 5.0}$ magnitude speckle to the LDFC-corrected DF with $10^{- 7.94}$ average contrast. The ideal DF is shown in the final frame for reference. Scale is $\log_{10}$ contrast.

Fig. 10

Performance of the spatial LDFC servo with a sinusoidal phase perturbation. $Gain = 0.6$ until the DF contrast reached $10^{- 7.9}$ . The gain was lowered to 0.1 for the remaining iterations. (a) Average contrast across the full DF for the pre-EFC PSF, DF post-EFC (blue), DF post-EFC with injected speckle (red), and the corrected DF post-LDFC (black). (b) Average DF contrast (black) over 50 iterations showing convergence to the initial EFC contrast (blue) after six iterations.

3.2.

Kolmogorov Phase Perturbation

In the first case, the injected aberration created a single speckle in the DF and a corresponding speckle in the BF. To demonstrate spatial LDFC’s ability to suppress multiple speckles, a Kolmogorov phase aberration was generated in the pupil plane instead of a sinusoidal phase perturbation (see Fig. 11). The phase perturbation was given a P-V amplitude of 20.5 nm, creating an aberrated DF with an average contrast of $10^{- 6.51}$ . The LDFC control loop was again run with a gain of 0.6 until the average DF contrast reached $10^{- 7.9}$ , at which point the gain was reduced to 0.1 to maintain the correction. The LDFC control loop was allowed to run for 50 iterations for this demonstration with convergence occurring after six iterations. The results are shown in Table 3 and shown in Figs. 11 Fig. 12 Fig. 13–14. As in the previous single speckle demonstration, the LDFC-corrected DF contrast occasionally drops below the initial EFC contrast level in Fig. 14(b). This effect is due to noise fluctuations.

Table 3

Performance with Kolmogorov phase: initial EFC DF average contrast, magnitude of the injected speckle, average contrast of the aberrated DF, average contrast of the DF after LDFC, total change in contrast for one full LDFC loop, and the number of iterations to converge to the EFC contrast floor.

EFC DF contrast	$10^{- 7.94}$
Avg. DF contrast with aberration	$10^{- 6.51}$
LDFC DF contrast	$10^{- 7.94}$
$Δ Contrast$	$10^{- 1.43}$
# Iterations to converge	6

Fig. 11

The injected Kolmogorov phase perturbation with (a) a P-V amplitude of 20.5 nm (RMS: 5 nm), (b) the DM response derived by LDFC, and (c) the final residual wavefront error (RMS: 0.39 nm) after six iterations. Scale is given in nm.

Fig. 12

(a) The aberrated PSF with multiple speckles in the DF and average DF contrast of $10^{- 6.51}$ , (b) the final LDFC-corrected DF with $10^{- 7.94}$ average DF contrast, and (c) the reference EFC-derived DF with $10^{- 7.94}$ average DF contrast. Scale is $\log_{10}$ contrast.

Fig. 13

Evolution of the DF over the six iterations [seen in Fig. 14(b)] to converge from a degraded DF average contrast of $10^{- 6.51}$ to the LDFC-corrected DF with $10^{- 7.94}$ average contrast. The ideal DF is shown in the final frame for reference. Scale is $\log_{10}$ contrast.

Fig. 14

Performance of the spatial LDFC servo with a Kolmogorov phase perturbation. $Gain = 0.6$ until the DF contrast reached $10^{- 7.9}$ . The gain was lowered to 0.1 for the remaining iterations. (a) Average contrast across the full DF for the pre-EFC PSF, DF post-EFC (blue), DF post-EFC with injected speckle (red), and the corrected DF post-LDFC (black). (b) Average DF contrast (black) over 50 LDFC iterations showing convergence to the initial EFC contrast (blue) after six iterations.

4. Discussion

We have demonstrated here that spatial LDFC is capable of locking the DF contrast at its ideal EFC state using only the BF response to a perturbation in the optical path. However, there are limitations to spatial LDFC and a potential null space, which need to be explored. These issues and some potential solutions are addressed below.

4.1.

Limitations of Spatial Linear Dark Field Control

One significant limiting factor for spatial LDFC is DF symmetry. This technique requires access to a BF that is located spatially opposite the DF. Due to this requirement, spatial LDFC is expected to work only with a nonsymmetric DF. However, in the case of a symmetric DF, spectral LDFC offers a possible solution (see Sec. 4.2). Since spectral LDFC relies on speckles that are located spatially within the DF but outside of the control bandwidth, it is not affected by the lack of a BF spatially opposite the DF. In the case of a much larger DF than the one presented here, spatial LDFC is predicted to still be capable of stabilizing the DF, but it cannot use BF speckles at spatial frequencies higher than those present in the DF to do so.

In future work, spatial LDFC’s null space must also be further explored. For this technique, the null space will consist of wavefront errors that affect the DF without changing the BF. One potential example of this null space is the formation of a speckle on a single side of the focal plane due to the combination of phase and amplitude sine wave aberrations. In such a case, if the speckle falls inside the DF, the BF will not see any modulation and will, therefore, be unable to sense and correct the aberration. A second potential null space example would consist of an incident phase aberration sine wave with a phase that creates a BF speckle with a phase that is 90 deg from the local BF phase. This case would not create a linear signal and would, therefore, not be corrected by LDFC. To date, the simulations presented here have not revealed a null space, but it can exist under system-specific conditions. It should also be noted that this paper has specifically explored a system in which aberrations were introduced and corrected by the same DM; in real systems, there will be aberrations that occur outside the DM-conjugate pupil plane and subsequently do not correspond exactly to DM authority. Such cases will require further investigation as spatial LDFC continues to develop.

4.2.

Addressing the Null Space with Spectral Linear Dark Field Control

A potential solution for overcoming spatial LDFC’s null space is to operate a separate version of LDFC simultaneously. This second version, known as spectral LDFC, freezes the state of the DF within the control bandwidth by using state measurements of light outside of the control bandwidth. This method exploits the fixed wavelength relationships that exist between speckles at different wavelengths that were generated by the same aberration. To first order, this fixed relationship scales the speckle separation linearly with wavelength and scales the complex amplitude inversely with wavelength. The complex amplitude speckle field may also interfere with static chromatic coronagraph residuals due to the coronagraph’s finite design bandwidth. These relationships between out-of-band and in-band light allow for the state of the DF within the control band to be monitored and maintained by measurements made of speckles located outside of the spectral control band.¹⁸ Since spectral and spatial LDFC rely on a BF signal from separate dimensions, the null spaces of the two forms of LDFC are not expected to overlap. For this reason, concurrent operation of spectral and spatial LDFC can provide a powerful tool for compensating for the separate null spaces of both techniques.

As an example of the BF signal that can be used by spectral LDFC, Fig. 15 shows the DF created using a phase induced amplitude apodization (PIAA)¹⁹ coronagraph at JPL’s high-contrast imaging testbed (HCIT). Speckles within the DF are shown at multiple wavelengths, both in-band (the science image) and out-of-band (the signal used by spectral LDFC).

Fig. 15

Spectral LDFC example: the $7 λ / D \times 8 λ / D$ DF from JPL’s HCIT created using a PIAA coronagraph. The DF is shown in four individual spectral channels: two channels within the control bandwidth and two out-of-band channels with speckles used to maintain the DF state within the control bandwidth. Also shown is the average DF over the full 10% bandwidth centered at $λ = 800 nm$ .²⁰

While spectral LDFC was not in operation when this image was taken, this is a clear demonstration of a case in which the in-band DF contrast could be maintained by sensing the speckles that are outside the control bandwidth and applying the appropriate wavelength-scaled correction to cancel the in-band speckles. Further development and analysis of this form of LDFC can be found in an upcoming paper by Guyon et al.¹⁸

5. Conclusion

In summary, spatial LDFC acts as an extension of EFC by operating as a servo that can maintain high contrast in the DF during science exposures. Using changes in the BF to provide updates on the state of the field within the DF, spatial LDFC is able to lock the state of the DF after it is established by EFC without relying on field modulation which interrupts the science acquisition and fundamentally limits the exposure time. The substantial increase in uninterrupted observation time spatial LDFC offers makes it a more efficient method than EFC for maintaining deep contrast and will lead to an overall increase in the number of planets detected and analyzed over the lifetime of an instrument. Here, we have introduced the mathematical principles behind spatial LDFC and provided demonstrations of its capabilities through numerical simulation; future work will include further analysis of this technique’s abilities as well as laboratory demonstrations.

Acknowledgments

This work was funded by the NASA Exoplanet Exploration Program Segmented Coronagraph Design and Analysis study effort. The authors would like to acknowledge Dr. Sandrine Thomas for her assistance in implementing EFC for this work. For their continued support, the authors also acknowledge fellow members of the OG and JRM research group: Justin Knight, Jennifer Lumbres, Alexander Rodack, and Lauren Schatz.

References

1.

California Institute of Technology, “,” Confirmed Planets, NASA Exoplanet Science Institute: NASA exoplanet archive“2017). Google Scholar

2.

M. Bottom et al., “Speckle nulling wavefront control for Palomar and Keck,” Proc. SPIE, 9909 990955 (2016). http://dx.doi.org/10.1117/12.2233025 PSISDG 0277-786X Google Scholar

3.

J. T. Trauger and W. A. Traub, “A laboratory demonstration of the capability to image an earth-like extrasolar planet,” Nature, 446 771 –773 (2007). http://dx.doi.org/10.1038/nature05729 Google Scholar

4.

W. Traub, B. Oppenheimer, “Direct imaging of exoplanets,” Exoplanets, University of Arizona Press, Tucson (2010). Google Scholar

5.

O. Guyon, “Limits of adaptive optics for high-contrast imaging,” Astrophys. J., 629 (1), 592 –614 (2005). http://dx.doi.org/10.1086/apj.2005.629.issue-1 Google Scholar

6.

J. R. Delorme et al., “Focal plane wavefront sensor achromatization: the multireference self-coherent camera,” Astron. Astrophys., 588 A136 (2016). http://dx.doi.org/10.1051/0004-6361/201527657 Google Scholar

7.

E. Cady et al., “Electric field conjugation with the project 1640 coronagraph,” Proc. SPIE, 8864 88640K (2013). http://dx.doi.org/10.1117/12.2024635 Google Scholar

8.

J.-B. Ruffio, “Non common path aberrations correction,” (2014). Google Scholar

9.

J. Krist, B. Nemati and B. Mennesson, “Numerical modeling of the proposed WFIRST-AFTA coronagraphs and their predicted performances,” J. Astron. Telesc. Instrum. Syst., 2 011003 (2016). http://dx.doi.org/10.1117/1.JATIS.2.1.011003 Google Scholar

10.

A. Give’on et al., “Broadband wavefront correction algorithm for high-contrast imaging systems,” Proc. SPIE, 6691 66910A (2007). http://dx.doi.org/10.1117/12.733122 Google Scholar

11.

T. D. Groff et al., “Methods and limitations of focal plane sensing, estimation, and control in high-contrast imaging,” J. Astron. Telesc. Instrum. Syst., 2 (1), 011009 (2015). http://dx.doi.org/10.1117/1.JATIS.2.1.011009 Google Scholar

12.

G. Singh et al., “On-sky demonstration of low-order wavefront sensing and control with focal plane phase mask coronagraphs,” Publ. Astron. Soc. Pac., 127 (955), 857 –869 (2015). http://dx.doi.org/10.1086/682726 PASPAU 0004-6280 Google Scholar

13.

O. Guyon, T. Matsuo and R. Angel, “Coronagraphic low-order wave-front sensor: Principle and application to a phase-induced amplitude coronagraph,” Astrophys. J., 693 (1), 75 –84 (2009). http://dx.doi.org/10.1088/0004-637X/693/1/75 Google Scholar

14.

J. K. Wallace et al., “The Gemini Planet Imager calibration wavefront sensor instrument,” Proc. SPIE, 7736 77365D (2010). http://dx.doi.org/10.1117/12.858269 PSISDG 0277-786X Google Scholar

15.

K. Miller and O. Guyon, “Linear dark field control: simulation for implementation and testing on the UA wavefront control testbed,” Proc. SPIE, 9909 99094G (2016). http://dx.doi.org/10.1117/12.2232120 PSISDG 0277-786X Google Scholar

16.

O. Guyon et al., “Wavefront sensing, control, and PSF calibration,” in NASA–Jet Propulsion Laboratory (JPL) Low-Order Wavefront Sensing and Control (LOWFS/C) Conf., (2015). Google Scholar

17.

L. A. Poyneer and J.-P. Véran, “Optimal modal Fourier-transform wavefront control,” J. Opt. Soc. Am. A, 22 1515 –1526 (2005). http://dx.doi.org/10.1364/JOSAA.22.001515 JOAOD6 0740-3232 Google Scholar

18.

O. Guyon et al., “Spectral LDFC: stabilizing deep contrast for exoplanet imaging using out-of-band light,” (2017). Google Scholar

19.

O. Guyon et al., “The phase-induced amplitude coronagraph (PIAA),” Proc. Int. Astron. Union, 1 (C200), 385 –392 (2005). http://dx.doi.org/10.1017/S174392130600963X 1743-9213 Google Scholar

20.

O. Guyon et al., “Phase-induced amplitude apodization (PIAA) technology development, milestone 3: 10% bandpass contrast demonstration,” (2014). Google Scholar

Biography

Kelsey Miller is a PhD student at the University of Arizona, the College of Optical Sciences. She has a bachelor’s degree in physics with a focus in astrophysics from the New Mexico Institute of Mining and Technology and a master’s degree in optical sciences from the University of Arizona. Her doctoral research is in high-contrast exoplanet imaging with a focus on focal plane wavefront sensing and control systems for extreme adaptive optics systems.

Biographies for the other authors are not available.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Kelsey Miller, Olivier Guyon, and Jared Males "Spatial linear dark field control: stabilizing deep contrast for exoplanet imaging using bright speckles," Journal of Astronomical Telescopes, Instruments, and Systems 3(4), 049002 (30 October 2017). https://doi.org/10.1117/1.JATIS.3.4.049002

Received: 4 April 2017; Accepted: 3 October 2017; Published: 30 October 2017

Access the abstract

JOURNAL ARTICLE
11 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 42 scholarly publications.

Explore citations on Lens.org

KEYWORDS

Speckle

Exoplanets

Wavefronts

Modulation

Actuators

Coronagraphy

Phase shift keying

1.

Introduction

2.

Concept