1. Introduction
In photonics, one of the most intriguing tasks is to visualize a point-like image of an Earth-type planet in the vicinity of a neighboring star (of Solar type) at a distance on the order of 10 parsecs. Formally, the spatial resolution of a meter class telescope in the optical wavelength range is sufficient by Rayleigh criterion to detect the outstanding peak of an exoplanet point spread function (PSF) at a stellocentric is unnecessary distance of 1...10 λ/D.
However, the luminosity ratio of star/planet is about 10
9 … 10
10. The planet has insufficient contrast as the faint light source relative to the host star light source. We shall consider a two-meter class telescope because the effect of exozodi is surrounding and therefore masking a stellar vicinity [
1,
2]. If the telescope has an ideal optical quality and has a stellar coronagraph mounted after, such a coronagraph instrument can show a planet on an attenuated diffraction background of the host star or, more precisely, on the attenuated background of the stellar PSF wings.
Modern space optics offers diffraction-limited resolution; however, optics is not free from residual aberrations and micro-roughness, which result in a visible
halo effect around the main maximum of the PSF. Certainly that masks the image of the faint planet light source. To move from a theoretical to a practical scope, one has to use precise adaptive optics (AO) in order to eliminate the wavefront (WF) error (WFE). The adaptive optics system aims to measure the WFE and then to compensate for it. One fundamentally analyses the diffraction problem resumes wavefront error in terms of complex amplitude [
3]. One has to compensate not only for the phase component of the WFE, but also for its amplitude component.
The amplitude component of the WFE is caused by some zonal inhomogeneity in transmittance or/and reflection; as well, it is due to Fresnel diffraction on the micro-roughnesses, residual aberrations, and aperture boundaries of every optical element (or optical surface) mounted in the plane not being conjugated with the pupil of the optical system, e.g., a secondary telescope mirror, downstream coronagraph optics, a spider, apertures etc. Moreover, a fringe field effect can be caused in the LC layer that could distort the phase and amplitude wavefront profiles in outgoing light from the LC-SLM itself [
4].
A portrait of the Fresnel diffraction can be primitively illustrated via the Talbot effect for a periodic structure or via sharp boundary oscillating. To compensate for the WFE in the Fresnel zone (not in the Fraunhofer zone), we require an increasing number of controllable pixels or actuators if deformable mirrors (DM) are used.
In the present communication, we show new technique and demonstrate the start-up laboratory experiments aiming at the phase–amplitude correction for the WFE in a stellar interference coronagraph.
Omitting details, we shall address the reader to the published papers on the characteristics and operations of the stellar coronagraph, as constructed by a classical Lyot schematic with and without apodization schemes [
5,
6,
7], and, in particular, by a coronagraph based on destructive interference process [
8,
9].
The working principle of an interfero coronagraph [
10] is to superimpose with antiphase two pupil images being shifted [
8], reversed, or rotated [
11], or differently superposed [
12]. The listed interferometers result in a destructive interference process. Therefore, an interferometer eliminates the starlight being received from the on-axis direction. Planetary (or companion) light does not interfere destructively because of the off-axis tilt (or shift or shear in a pupil). The image component from a planet becomes split into two copies; by superposition, they are spatially separated and do not interfere destructively. Interfero-coronagraphy is advantageous in terms of its broad spectral band achromaticity and, more generally, because of its small inner working angle (IWA). The IWA characteristic can be considered as the spatial resolution of a coronagraph instrument. The achromatic interfero-coronagraph (AIC) is known to have one of the smallest possible IWAs (0.38 λ/D), such as it was initially referred to in [
10] with the fixed angle of a 180-deg. rotational shear. Later, it was redesigned into a modified Sagnac scheme, implementing the common path (or cyclic path or cavity) (CP-AIC) [
13] aiming to relax the mechanic instability.
The severe functional disadvantage of an interfero coronagraph is known as the stellar leakage effect. Due to this effect, the starlight cannot be completely suppressed because the apparent size of the star is not physically infinitesimal [
5]. The observed size of a star is far beyond the spatial resolution of a
single dish optical telescope, but physically, the star size causes a non-fully coherent point-like source; considering spatial coherence, it causes an extended source. Numerous methods to reduce the stellar leakage effect, due to the infinitesimal apparent size of the star, have been proposed from the 180-deg. rotational shift to smaller angles [
11] or to apply an apodization, e.g., a Sonine type [
14].
In the present communication, the proposed method is applicable to the listed possible efforts to reduce the star leakage effect.
2. Method
Precise wavefront correction in phase and amplitude is ultimately required as pre-optics schematics for the coronagraph for its functionality. AO compensates for optical aberrations including optical defects given by the optical surfaces and apertures. Light radiation is collected by a telescope under or above the telluric atmosphere, and it is analyzed after a coronagraph by a field camera, spectroscopy, or different instrument. It is important to note that the measurement of the wavefront has to be organized after the coronagraph. Otherwise, a non-common-path wavefront error (NCP WFE) is caused by a different or incompletely similar optical path to the wavefront sensor (WFS). This generates additional phase and amplitude errors that become magnified [
15] by the coronagraph.
Aiming active adaptive optics to control the WFE in phase and amplitude, we used a phase-only liquid crystal spatial light modulator (LC-SLM) mounted in a specific polarization schematic as shown in
Figure 1 inside the dashed-line box I. In our work, we used the reflective type LC-SLM. In
Figure 1, the LC-SLM is shown in transmittance mode for the purpose of simplification. It is better to mount the LC-SLM by an incident angle less than 15° to minimize depolarization effects, which can reduce the SLM modulation contrast [
16]. In our setup the incident angle was about 4°.
Math algebra of the polarization of Jones vectors and matrices was applied to describe the action of the
phase-only LC-SLM being used as a phase and amplitude WFE corrector. We considered the output of the following polarization scheme: entrance linear polarization with the azimuth orientation at 45° to the
X axis, LC-SLM with the phase modulation axis along
X, and the final linear polarizer with transmission at 45° to
X. See this set, in the dash-line box I in
Figure 1.
Then, the phase–amplitude control of the WF is realized as:
where
=
defines an entrance wave, with the polarization linear state with a 45-deg. azimuth respective
X axis, with the complex amplitude
at entrance, where
denotes the modulus (or amplitude) and
—the phase. At the output of box I in
Figure 1, the complex amplitude is denoted by
.
defines the phase-only LC-SLM (liquid crystal spatial light modulator) with the phase angle modulation along X, which is controlled by an externally applied electric voltage.
= defines the linear polarizer oriented by its transmission axis at 45-deg. to the X axis.
From Equation (1c), it can be seen that the control of α leads to a change not only in the phase but also in the amplitude of the complex field. This useful property is next applied to a phase–amplitude correction of the wavefront.
Presented here, a phase–amplitude correction of the WFE is entirely developed for the optical schemes, which use the principle of a nulling interferometer with rotational shift functioning as a stellar coronagraph [
11,
13]. A shear interferometer has two mutually shifted waves being superposed. This allows two degrees of freedom: two controls (modulations)
and
, at two “coupled” points “A” and “B” of the initial wavefront, which are superimposed. We aim these controls to obtain the mostly accurate “dark field” destructive interference condition.
Let us consider that the shear interferometer superposes two electric fields
and
in points “
A” and “
B”. This is illustrated in
Figure 1; see the dot-line box II. Before the interferometer, at the initially separated WF symmetric points “
A” and “
B”, by the SLM we control the complex amplitudes
and
; see Equation (1d). The destructive interference regime requires both the following conditions (i) and (ii).
(i): A zero-phase difference before the interferometer while we count on an additional phase shift of π implemented inside the nulling interferometer [
11,
13]:
(ii): The equality of the moduli of the complex amplitudes:
Therefore, to satisfy both the phase (2a) and amplitude (2b) conditions, we have to consider the electric fields in two
coupled (by following superposition) points
A and
B in a pupil before 180° RSI. Their waves interfere further in the conjugated pupil plane at the dark port (with implemented anti-phase) of the interferometer at the point denoted as
A +
B (See
Figure 1, dot-line box II).
Now, we recall the polarization scheme in
Figure 1, in dashed-line box I. Then the complex amplitudes of the electric field at points
A and
B with respective
and
modulations via the LC-SLM control are defined according to (1a)–(1d):
We stress here, again, our notations for the complex amplitudes moduli: , and phases , (without tilde) which describe the initial wavefront (before the LC-SLM), while the complex amplitudes moduli , and the phases , (with tilde) denote the corresponding values after the LC-SLM and polarizer (after the polarization scheme ; see Equation (1a).
After passing the interferometer, on its dark port, by coherent wave subtraction, we have the complex field amplitude at point
A +
B by the superposition of waves
,
:
The CCD detects the intensity
of the interference pattern, e.g., in the
A +
B point after ensemble averaging (where 〈
〉 defines the ensemble average over an exposure time):
where the
upper underline denotes the complex conjugate (in matrix form it denotes the Hermitian conjugate).
Substituting (4) into (5), we can express the detected intensity
through the complex amplitude moduli
,
and the phase difference
of the electric fields in the entrance pupil WF and the corresponding SLM modulations (
,
):
2.1. Measurement of the Wavefront
For the sake of brevity, we show here the algorithm of the wavefront measurement to determine amplitude and phase difference , by intensity measurements, which is their function, (6a). We study Equations in (6): if one measures a series (j = (1..3) of three intensities , in the plane conjugated to the pupil (the SLM is set in the pupil) by providing certain SLM modulations as parameters and in the coupled points A and B, he determines the phase difference and the amplitudes and for the next step of WFE correction.
For example, a three-step
j = (1..3) procedure to determine
,
,
detects three intensities by SLM modulations set as
:
Entrance wavefront characteristics
,
,
we determine from solving the system of Equation (7) by their expression as (6b):
Certainly, the wavefront measurement method can be modified and improved in accuracy, e.g., by increasing the number
j of measurements
,
, similar to phase shifting interferometry [
17].
We resume here that by means of Equations (8a)–(8c), it is possible to measure the wavefront error in terms , at symmetric (or coupled by electric fields superposition) points A and B. Therefore, we can restore the spatial distribution of WFE pixel wise.
2.2. Wavefront Correction by Phase-Amplitude Modulation
We intend to show here a method to correct the wavefront error, taking into account three measured wavefront characteristics: , . These measured WFE characteristics at the coupled pupil points A and B are associated with subsequent interference in the A + B point. This correction method is easily extended to cover all the pupil points.
2.2.1. Phase Correction
To present the WFE correction problem, at first, we simplify considering the phase-only WFE: , .
Figure 2 depicts schematically the principle of phase correction in the
coupled points
A and
B. Here, the
and
denote the phases in the
coupled points
A and
B (set diametrically opposite) in the pupil plane where the SLM is installed. After passing through the RSI interferometer, an optically conjugated pupil plane is formed. Here, in the points denoted
A +
B,
L (in the left half-plane) and
A +
B,
R (in the right half-plane), the interference signals intensities of interference pattern
,
become proportional to the cosine of
:
Because the cosine function is even visually through the recorded intensity, there is no difference between the left and right half-planes: , but physically, the relationship (9) turns out to be important to organize the proper wavefront correction.
The measured phase difference
is added then as an additional phase shift to the phase of point
A (on the left) or to the right point
B (but with a negative sign
not shown in the figure). In
Figure 2, next after passing through the interferometer RSI, at the two points
A +
B,
L and
A +
B,
R (shown in green), the corrected phase difference becomes
≡ 0, which means the phase-only aberrations (or wavefront errors) have been successfully corrected.
2.2.2. Phase- and Amplitude Wavefront Error Correction
We analyze here the algorithm to correct a more complex WFE by
(in phase),
(and in amplitude), which in contrast to the phase-only correction (discussed above in
Section 2.2.1), contains an additional amplitude imbalance correction option.
For simplicity, we denote
, see Equation (9a,9b). According to Equations (4) and (5), we search for such modulations
and
that cause a zero-intensity at this point:
. We solve therefore the following equation:
aiming to find two unknowns,
and
, each of which depends on four variables: the moduli of complex amplitudes of fields
,
, and the phases
,
. We reduce four variables to two variables, which are physically relevant in the frame of our model, namely the amplitude ratio:
and the phase difference:
. The
and
are real numbers, and they were measured by Equations (8a)–(8c). Equation (10) has the following solution:
We did not simplify the arctangent arguments in
, Equation (11a), and in
, Equation (11b), and we left them in ratios to use the function
[
18].
3. Simulation
Figure 3 illustrates some numerical verifications to validate the analytical solution found by Equation (11). In the figure, an analytical solution via Equation (11) is shown by the
red cross unnecessary mark and coincides visually with the minimum in intensity in the logarithmic scale
of the interference pattern shown by the background colormap. The latter was evaluated by a trivial enumeration of
in the mesh region
by the values
radians, where the used
and
have been chosen arbitrarily. The red cross mark position in
coordinates mesh assigns the global minimum.
In
Figure 4, we analyze an intensity transmission characteristic
, which shows the attenuation of the signal (from a faint exoplanet) in a single pupil pixel because we corrected the amplitude imbalance of
. In a pixel, the planetary intensity signal transmission
becomes reduced relative to unitary because an amplitude correction attenuates the transmission of both signals from the star and planet. We have considered here the WFE in point
A +
B by
. The
Figure 4 background has the colormap that is the two-dimensional distribution of
being evaluated in the domain of
.
In
Figure 5, in a ±15 λ/D focal image domain we show several simulation graphs to compare the images (1) in non-coronagraphic mode, (2) in coronagraphic mode without any WFE correction, (3) in coronagraphic mode with WFE phase-only correction, and (4) in coronagraphic mode with WFE phase and amplitude correction. These images are resumed by azimuth-averaged stellocentric cross-sections (5).
In simulations in
Figure 5, one can see that, caused by diffraction and scattering on residual aberrations and surface micro-roughness, the stellar background shows the gradual suppression by means of the coronagraph without and with the wavefront control. The wavefront control was performed by correcting the phase-only wavefront errors at first, then correcting both the phase and amplitude WFE. In these simulations, the parameters were phase deviations of about λ/50 rms at λ = 600 nm and amplitude deviations of about 5%.
4. Laboratory Experiment
In the laboratory experiment, we studied the above proposed approach to correct the complex value wavefront error simultaneously in phase and amplitude at the wavelength 632.8 nm. In addition to our previously published efforts [
9] where we performed the phase-only wavefront control, here, we succeeded in both measuring the complex number WFE, and then controlling the WFE in phase and amplitude operating by a spatial light modulator (SLM) being inserted in the optical and polarization scheme, as in
Figure 1. We used a commercially available LC-SLM from HOLOEYE© (Berlin, Germany): model PLUTO-2.1-VIS-016. It has an 8 μm pixel size with a resolution up to 1920 × 1080 active pixels and a response time of ~66 ms [
19]. The LC-SLM was capable of providing a maximum phase shift of about 5.2 π (on the working wavelength 630 nm), but it was calibrated to provide only reduced 2π of maximum phase shift, that was sufficient for the described correction method. The camera used in the setup was a CMOS monochrome camera from Edmund Optics© (Barrington, NJ, USA), model EO-5012M with a pixel size of 2.2 μm and a resolution of up to 2056 × 1920. It was synced with the LC-SLM by a dedicated frame sync signal. Several experimental results are shown in
Figure 6, where we demonstrate the gradually reducing azimuth-averaged cross-sections of the non-coronagraphic PSF, the coronagraphic PSF without WFE correction, and the coronagraphic PSF with WFE correction in phase and amplitude; see panel (a). In
Figure 6 in panels (b) and (c), we show noncorrected fragments of pupil intensities after the coronagraph, where the WFE was noncorrected and corrected. Corresponding histograms are shown in panels (d) and (e). If one analyzes pupil fragments (b) and (c), he finds at first the decreased intensity level; see the colorbar scale on the left. Integrally, this decrease corresponds in numbers to the ratio of the coronagraphic PSF without applied WFE correction (shown in the azimuth-averaged section by the red color graph in panel (a)) to the coronagraphic PSF with WFE correction (shown by the orange color graph). Additionally, the corresponding histogram (in panels (d) and (e)) demonstrate that the intensity deviations decrease after correction. Because of amplitude attenuation, we estimated the virtual planetary signal had its throughput at about 0.92 as averaged over all pixels.
5. Discussion
In the present article, we have continued the studies of our previous communication [
9], where we just roughly sketched some general possibilities about how to correct wavefronts, not only in phase but simultaneously in amplitude. This is realized by the LC-SLM providing general polarization modulation at first, then it is converted in phase and in amplitude wavefront components by the linear polarizer mounted after the LC-SLM. Therefore, in the present communication, we have aimed to find and to test a practical analytical solution for wavefront error active correction both in phase and amplitude. This is shown analytically as well by numerical simulation and in a simple experiment in an optical lab.
Here, our goal was not to reach a high coronagraphic contrast in stellar light suppression, but we studied an appropriate optical architecture to control WFE. At first, we measured both the WFE phase and amplitude components. Then we operated by controlling voltages at an adaptive optics unit (LC-SLM) to correct wavefront error simultaneously in phase and amplitude. The advantages are the following: (i) both the wavefront measurement and the correction have been realized without any mechanical motion in the optical scheme and without its reconfiguration; (ii) such a WFE measurement was performed after the coronagraph, therefore we excluded possible non-common-path aberrations (NCPAs).
The simple approach presented here was found analytically and it was tested numerically by simulations using the Proper [
20]. Finally, we designed an optical experiment in the lab, which has demonstrated the right tendency in deeper signal suppression in coronagraphic mode for an on-axis light source (that imitated starlight) by wavefront correction.
In future, we are optimistic about applying this technique onboard of an orbital telescope for direct exoplanet imaging. LC-SLM technology was
scientifically tested under a thermo-vacuum environment [
21], however not yet under the space requirement. A possible LC-SLM mounting can be in a chamber with an optical window, with stabilized normal temperature and pressure. Radiation tests are not known, but a telescope can be launched in a low orbit beneath the Van Allen radiation belt. Active adaptive optics competes with the static correction, but in space, an orbital telescope is stabilized statically up to certain level. Solar radiation plus heat generation by electronics causes a low speed mechanical instability of the primary mirror geometry and to optics having requirements greater than λ/100 in optical wavelength domain. Therefore, precise adaptive optics with the functions of measurement and the control of the wavefront in phase and simultaneously in amplitude remains advanced.