A&A 374, 1161-1168 (2001)
DOI: 10.1051/0004-6361:20010619

The achromatic phase knife coronagraph

L. Abe1 - F. Vakili1 - A. Boccaletti2


1 - Observatoire de la Côte d'Azur, Département Fresnel, CNRS UMR 6528, 06460 Saint Vallier de Thiey, France
2 - Collège de France, 11 place Marcelin Berthelot, 75005 Paris, France

Received 20 November 2000 / Accepted 11 April 2001

Abstract
We present the concept of a nulling coronagraph which theoretically attains the brightness ratio of 108 suitable to directly detect exo-Earths in the vicinity of their parent stars. The concept is derived from the classical Foucault's optics test and Zernike's self-destructive phase contrast principle. We propose an achromatic set-up which should dramatically improve its wide band efficiency. The whole system is simulated against a numerical model which validates a formal analytical approach. The coronagraph is then compared to other concepts and its advantages and shortcomings are outlined, suggesting future laboratory developments and tests on the sky.

Key words: techniques: interferometric - stars: binaries: close - stars: planetary systems


1 Introduction

The amazing growth of exo-planetary detection since the discovery of 51 Peg B (Mayor & Queloz 1995) has triggered considerable efforts for finding new methods to record the direct signal of exo-Earths shrouded by the scattered light of their parent star due to telescope optics (Marcy et al. 2000). Indeed the huge contrast ratio of 109 to 106 (in V and N bands respectively) between a G-star and its orbiting exo-Earth, demands enormous dynamic range imaging possibilites that only nulling techniques, whether interferometric or coronagraphic, could attain. Unlike Lyot's mask, nulling coronagraphs offer both the angular resolution, a few tens of milli-arc-seconds (mas), and the required dynamic range to hunt exo-Earths around a statistically meaningful sample of nearby G-type stars. Pending Bracewell interferometers, e.g. Darwin (Leger et al. 1996) or TPF (Beichman et al. 1999) a large optical telescope, e.g. NGST (Mather 1997), with an embarked coronagraph seems the most likely instrument to collect photons of an exo-planet within 2010-2015 horizon. Following the original concept Gay's Achromatic Interferometric Coronagraph, AIC hereafter (Gay & Rabbia 1996), other designs have been proposed (Roddier & Roddier 1997) and (Rouan et al. 2000) which, unlike AIC, present inherent chromaticity limiting their net nulling efficiency over a wide spectral band.

Following our earlier work on the NGST Exo-planet Finder (Boccaletti et al. 2000), we propose hereafter a coronagraphic design which overcomes the chromatism problem of common coronagraphs, i.e. both the retardation and size of the phase mask. In the next section we outline the general theory of our coronagraph. Section 3 gives a generic optical set-up to overcome the chromatism problem. In Sect. 4 the conceptual design is validated by a number of numerical experiments. Finally we compare the theoretical efficiency of our coronagraph and discuss its limitations and sensitivity to various optical and operational parameters. A mathematical description of the concept is also given in the Appendix.

  \begin{figure}
\par\includegraphics[width=8.8cm,height={!}]{Abe10486f1.eps} \end{figure} Figure 1: Pupil intensity with perfect wavefront and its corresponding Airy pattern ( top left and right). Intensity distribution after the Phase Knife Coronagraph has been applied ( bottom-left): the two thin crescents encircle the pupil area perpendicular to the Knife-Edge direction. "Butterfly shape" of the point spread function of a system where half the amplitude is $\pi $-shifted in the image plane ( bottom-right), and where a Lyot stop has been applied in the conjugate pupil plane of (c).
Open with DEXTER

2 Conceptual design

Our coronagraph is based on Zernike's self-destructive interference concept of phase contrast (Zernike, Nobel Prize 1953) also used by Roddier (Roddier & Roddier 1997) for their Phase Mask coronagraph (PM hereafter) and Rouan (Rouan et al. 2000) for the Four Quadrant Phase Mask (FQ-PM). PM splits the amplitude of the image plane into two parts: a core region with a $0.54\lambda/D$diameter shifted by $\pi $ against the rest of the image. The image is assumed an ideal Airy pattern resulting from the diffraction of a completely flat wavefront across a circular pupil. The destructive interference in the conjugate pupil plane between the two half-amplitudes with opposite signs rejects most of the energy off the pupil area. Applying an apodising mask on this pupil, about 0.8 of its geometrical size, results in a 100 times darker Airy disc in the immediately following image plane. Since the $\pi $-shift and the Airy disc size both depend on $\lambda$, PM's efficiency rapidly degrades when changing the wavelength (Guyon et al. 1999). The FQ-PM is not sensitive to the geometric size of the Airy pattern, but still has a chromatic retardation dependance. A solution for the PM chromaticity problem was originally suggested by Labeyrie (Labeyrie et al. 1999) and can also be applied to FQ-PM (Riaud et al. submitted).

For the present concept, we borrow the principle of Foucault's optics test by replacing the knife-edge by a $\pi $-shifted retarding screen (or mirror), hence the name of Phase-Knife Coronagraph (PKC hereafter). When the phase knife-edge crosses exactly the center of the Airy disc the maximum of a destructive interference occurs in the conjugate pupil image which presents two bright crescents symmetrical to the knife edge direction, as amateur mirror-polishers often observe when performing Foucault's test (Fig. 1, bottom-left). Now if we apply a circular apodizing mask on this pupil, occulting its two lateral bright crescents, we will obtain a pseudo-Airy pattern with a strong cone-like transversal extinction area (Fig. 1, bottom-right). Compared to the original input Airy pattern the intensity distribution undergoes here a noticeable nulling parallel to the phase-knife edge, but is rather ineffective in the perpendicular direction. Note that in this concept the geometrical dependence of the phase-mask is naturally eliminated. However if the image of the exo-planet happened to exactly fall along the phase-knife edge its pseudo-Airy pattern would also vanish, exactly like its parent star. Now if we apply the same operation of the phase-knife edge a second time on the output pseudo-Airy pattern in the direction perpendicular to the first phase-knife we obtain an extremely efficient nulling effect since the two perpendicular phase-knives cancel their mutual limb-brightening effect across the conjugate pupil. In this respect the nulling effect of the PKC is equivalent to the FQ-PM coronagraph (Rouan et al. 2000). However, as it will become clear in the next section, our approach offers a straightforward solution with off-the-shelf optical components to the problem of achromaticity and wide band operation of nulling coronagraphs.

  \begin{figure}
\par\resizebox{8.8cm}{!}{\includegraphics{Abe10486f2r.eps}}\end{figure} Figure 2: Generic 3D optical scheme of the PKC. I1 is the input Airy pattern, DY is a direct vision dispersing prism in the Y direction, I2 corresponds to the first chromatic phase-knife parallel to Y. RY is a second direct-vision prism rotated by 180 degrees with respect to DY which superimposes the dispersed phase-knived airy patterns after DY. The following DX and RX operate exactly the same as DY and RY but orthogonal to them. The final coronagraphic pseudo-Airy pattern is depicted in Fig. 3 bottom-right.
Open with DEXTER


  \begin{figure}
\par\resizebox{18cm}{12cm}{\includegraphics{Abe10486f3r.eps}}\end{figure} Figure 3: Different steps of the phase-knife screen effect: ( top-left) polychromatic Airy pattern (bandwidth: 400-800 nm) (corresponding to step I1 in Fig. 2), ( top-middle) dispersed Airy disc, ( top-right) the polychromatic phase-knife where the optical retardation follows the dispersion law (step I2), ( bottom-left) an intermediate image plane where the Airy discs are de-dispersed (step I3), ( bottom-middle) the polychromatic phase-knife applied in the perpendicular direction (step I4), ( bottom-right) the polychromatic mutually phase-knived pseudo-Airy disc.
Open with DEXTER

3 The achromatic PKC

Wide band operation is crucial for any nulling interferometer or coronagraph. In this respect the AIC coronagraph (Gay & Rabbia 1996) is the most efficient concept with the potential bonus of detecting a companion as close as its parent star's first Airy ring. AIC has two major shortcomings however: it rejects half of the incoming light (due to the beam-splitter) and its coronagraphic image becomes symmetrical (true orientations are lost). On the other hand Roddier & Roddier's PM depends both in retardation and size on the wavelength. It conserves the absolute orientation of the image and is more efficient in terms of collected photons from the planet since it needs a less agressive Lyot pupil diaphragm. Finally the FQ-PM (Rouan et al. 2000) overcomes the geometrical problem of Roddier's phase mask but does not give a straightforward solution to the wavelength dependance of this phase. For the PKC, since light-cancellation is obtained in two orthogonal steps with respect to X-Y directions, each step can be used to disperse the image and apply a wavelength-dependant phase retardation.

This is schematically represented in Fig. 2 where a direct vision prism (DVP hereafter) is used to disperse the input image I1 along the Y axis. Here the phase retardation is produced by putting a glass plate with constant thickness on one half of the dispersed image and a thickness variable plate on the other half of the image. For the latter, the thickness follows the dispersion law (for a prism it would vary as $\lambda^$ for instance). Its gradient remains negligible along the dispersion direction, i.e. a few hundreds of nanometers for a wide band PKC, avoiding any tilt between the monochromatic wavefronts (see Sect. 5.3).

A second DVP flipped by 180 $\hbox{$^\circ$ }$ around the propagation direction with respect to the first DVP re-folds the different phase-knived monochromatic images in I3 (Fig. 2). Then the operation is repeated in the perpendicular direction X and finally the cross-phase-knived coronagraphic image I5 undergoes an apodising Lyot mask (not represented in Fig. 2) in the pupil plane immediately following I5 which masks the limb-brightened edges of this pupil (see Fig. .4 top right and bottom right) to obtain the optimum starlight nulling.

Note that using two succesive DVPs for dispersing and refolding monochromatic phase-knived images is motivated by our concern to concentrate the exo-planet photons on the smallest number of detector pixels to minimize read-out noise. Once detection is positive, one could optionally remove the last DVP in order to record a low dispersion spectrum of the exo-planet for its hypothetical biomarker characterization.

4 Numerical simulations

4.1 Cross dispersion effect

In order to assess the expected performance of PKC we carried out a number of numerical experiments producing the final image (detector) plane. The model makes extensive use of 2D FFTs and direct convolution techniques whenever necessary. We worked on $2048\times2048$ double precision arrays to avoid round-out errors and aliasing effects. Figure 3 represents the different steps of the PKC nulling process. We took 20 "monochromatic" Airy discs scaled to their central wavelength from 400 to 800 nm by steps of 20 nm. These monochromatic Airy discs were co-added to mimic a wide band Airy pattern (Fig. 3 top left and mid). Next, the action of the first chromatic phase-knife was applied to dispersed monochromatic Airy disc amplitudes by assuming a linear shift of the Airy pattern along the dispersion direction. The phase-knived Airy amplitudes were then folded to mimic the action of the de-dispering second DVP. The action of the second chromatic phase-knife stage is mimiced by repeating the first stage in the perpendicular direction. Passing from each plane to its conjugate was performed by FFT and by adjusting the center of the diffraction patterns so as to avoid any numerical phase-artifact. Note that the Lyot stop was only applied before forming the last image by applying a diaphragm 0.9D, D being the diameter of the geometrical pupil. The final coronagraphic image is shown in Fig. 3 (bottom right). We note the cross-shaped blind-zone of PKC where the daugther planet PSF is affected by an intensity loss factor of about 10 (see discussion below). In order to assess the final dynamic range that one could attain with a "perfect" PKC we averaged the radial profile of the intensity in the final plane. Figure 4 depicts this profile where one can notice that even at the distance of the first Airy ring the expected contrast attains 106. This value increases up to 108 at the 5th Airy ring, the magic number of exo-Earth detection limit in the visible.

4.2 Ground-based observation simulation

As pointed out in Sect. 5, phase defects of the wavefront produced by atmospheric seeing and figuring errors of the primary mirror is a major problem for detecting faint companions. Current AO are actually unable to smooth in real time the wavefront down to the desired level of exo-planet detection.

We have carried out a numerical simulation to assess the real performance of PKC under atmospheric seeing. We generated a thousand short exposures considering an 8 m circular pupil including central obscuration and spiders. Compensating the wavefront turbulence with a high order AO, we obtained an average Strehl ratio of about 70%. This is typically the best performance expected for the VLT telescope using NAOS (Rousset et al. 1998). Then, we computed the coronagraphic image for each short exposure as described in Sect. 4. Finally, photon noise and readout noise were added to each frame. The long exposure image is obtained by directly co-adding the short exposures. Such an image is shown in Fig. 5. Companions as faint as $\Delta m=9$ located at the third diffraction ring, around the star, become detectable.

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{Abe10486f4.eps} %
\end{figure} Figure 4: Radiallly averaged profile of a polychromatic Airy pattern (solid line) obtained across a spectral bandwidth of 400 nm in the perfect wavefront case. The dotted line corresponds to the achromatic PKC. The maximum cancellation attains 107 to 109 in a large field of view at a distance of $20\lambda /D$ from the center of the field.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{Abe10486f5r.eps} %
\end{figure} Figure 5: Numerical simulation of PKC on an 8 m ground-based telescope equipped with adaptive optics. Atmospheric seeing and high-order compensation (Strehl ratio ${\simeq } 70\%$) are generated with 1000 independant phase maps to provide distorted Airy patterns. Faint companions have been added to short exposures with magnitude differences of 7, 8 and 9 (marked by arrows clockwise from bottom-right to top-left), with respect to the un-masked star. The image total integration time is 60 s. Photon and read-out noise ( $38\,{\rm e}^{-}$/pixel/frame) are added to the frames. Companions as faint as $\Delta m=9$ are directly imaged with standard integration of short frames.
Open with DEXTER

5 Discussion

For possible applications of the PKC to the real case of space-borne telescopes it is interesting to evaluate its sensitivity to different manufacturing and operating parameters. Hereafter we limit ourselves to a few of them (star- tracking, residual chromatism and the inherent blind zones in PKC's accessible field of view). The effects of light scattering due to the primary mirror bumpiness and the possible angular resolution of the on-axis parent star itself will be assessed in a more detailed study in the future. The following sections describe common problems of PKC and FQ-PM (Rouan et al. 2000), the former being presented as a straightforward solution to the chromatism problem of this kind of phase mask coronagraphs.

5.1 Global efficiency vs. centering

Considering an ideal Airy pattern, the residual halo of the coronagraphic PSF is reduced by 106 at the distance of the first Airy ring for a perfectly on-axis star. Tip-tilt errors represent the most degrading effects for using any coronagraphic instrument. They were simulated numerically for the PKC. The plot in Fig. 6 shows that decentering the star by 1/10th of $\lambda/D$ off-axis already drops the nulling effect to $2\times 10^4$for the second diffraction ring (2.5$\lambda/D$). Indeed space-borne telescopes like the NGST will be affected by much slower tip-tilt drifts, therefore meeting the severe condition given above. For other wave-front phase residuals, PKC is obviously less sensitive to centro-symmetrical aberrations.
  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{Abe10486f6.eps}\end{figure} Figure 6: The evolution of the brightness at the second diffraction ring of the coronagraphic PSF for increasing tilt errors, expressed as a fraction of the Airy radius.
Open with DEXTER

5.2 Detectability field

As explained earlier, the PKC presents four inherently pseudo-blind areas along the two orthogonal directions, corresponding to the relative orientation of the two phase knives. We numerically assessed the nulling efficiency near these critical areas by computing the energy transmission for an object located on the "blind axis" and gradually deviating from it. The effect, normalized on an undisturbed planet brightness is plotted in Fig. 7 where the loss is given in percentage versus shift, expressed in units of Airy radii. It exceeds 0.5 at a distance of 0.4 Airy radius. A 90% energy transmission is achieved over 1.4 unit distance. Note that on these two orthogonal axes the transmission is not zero, and some 10% intensity still pass through.

  \begin{figure}
\par\includegraphics[width=8cm,clip]{Abe10486f7.eps}\end{figure} Figure 7: Blind areas of PKC. Curve depicting the PKC effect for an object located on and near the dark axis: a 90 percent energy transmission is achieved over a 1.4 Airy radius distance perpendicularly to the axis.
Open with DEXTER

5.3 Residual chromatism

In sect. 3 we outlined the generic dispersing-refolding optical set-up to solve the problem of chromatism. After the dispersing step, the prismatic-shape of the phase knife must exactly match the dispersion law since a $\pi $- retardation is only valid for a given wavelength. Indeed the higher the dispersion, the smoother the slope of the prismatic phase-knife given the min/max $\pi $-retardations correctly match bandwidth boundaries, even though local deviations from the ideal $\pi $-retardation may occur. We simulated these aspects by introducing slight local variations following different slopes (linear). The influence of the slope is parametrized by a $\kappa$ factor taken as the ratio of the mean Airy size to the total spectrum extent (e.g. the central monochromatic Airy disc of Fig. 3, top-middle, to the extent of the horizontal spectrum). We considered three situations where the total bandwidth $\Delta\lambda=200$ nm (from 600 to 800 nm), and $\kappa=20, 5$ and 4 represented by the dotted, dash-dotted and short dashed curves of Fig. 8. In this figure the long dashed curve corresponds to a case without local phase variation, and the continuous line is the reference Airy profile without coronagraph. As a result one can see that in order to attain good nulling performances a dispersion over PSF ratio of 10 is sufficient which should easily be achieved with ordinary DVP or diffraction gratings. Note that the nulling performances of Fig. 8 differ from those previously presented (Fig. 4) and is only due to the use of smaller numerical arrays ( $1024\times 1024$).

  \begin{figure}
\par\resizebox{8.7cm}{!}{\includegraphics{Abe10486f8.eps}}\end{figure} Figure 8: Loss of nulling due to chromatism. PSF profiles for different local phase variations due to phase knives chromatic efficiency. From bottom to top, (long dashed) without phase variation ($\kappa =0$), (dotted) $\kappa =20$, (dash dotted) $\kappa =5$, (short dashed) $\kappa =4$ and the reference PSF pattern without coronagraph (continuous).
Open with DEXTER

6 Conclusion

In this paper we described the principles of a nulling coronagraph which operates in an image plane and offers achromatic properties both in terms of image geometry and phase shift. The achromatisation is obtained using common optical components and the actual manufacturing should not be difficult. We showed from numerical experiments that in the ideal case of an aberration-free telescope PKC directly attains the required 109 to 106nulling effect in the visible and IR wavelengths for direct detection of exo-planet photons. One of the major shortcoming of PKC is its partial blind cross-shaped areas, roughly two stripes $\lambda/D$ wide, parallel to the 2 successive phase- knives. If the exo-planet's image falls exactly along one of these stripes its Airy pattern is broken in two petals and attenuated by a factor of ten in intensity (see Fig. 1, bottom-right) which means these areas are semi-blind. This can be overcome provided the telescope can rotate around its optical axis, by say 45 $\hbox{$^\circ$ }$ for instance, and subtracting the two long exposures at each rotation angle with the extra bonus of subtracting the background noise. On the other hand the linear scheme of consecutive direct vision prisms and their relay optics for making PKC achromatic appears at first very sensitive to scattered light through the optical train. A workable PKC will operate off-axis with folded optics and baffling at intermediate image planes where care must be taken for extremely efficient coatings. Residual ghost images and scattered light will inevitably add an incoherent halo on the final coronagraphic image and limit the detectable faint planet within any reasonable integration time. A precise optical design involving ray-tracing calculations and post-processing of the collected data will be discussed in a next paper, including laboratory tests.

Another point which has not been discussed in this paper is the sensitivity of the achromatic PKC to different optical aberrations, specially light scattering due to residual mirror bumpiness at different spatial scales. This makes part of an end-to-end study that we are currently undertaking. It should be noticed however that for PKC, as for the other coronagraphs (AIC, PM and FQ-PM) the most critical issue remains the tip-tilt of the wavefront. The image jitter around the intersection of the two orthogonal phase-knives can dramatically reduce the nulling efficiency of PKC. Another critical issue is the chromatic dependance of the optical retardation along the direction of the phase-knife upon which the dispersed Airy disc of the star is formed. The manufacturing of this retardation, quadratic in the case of DVPs, needs careful design and manufacturing. An alternative solution, at the expense of luminosity, is to use diffraction gratings instead of DVPs. One could also think about a folded PKC using single diffraction gratings. A laboratory bread-board could more quantitatively define the nulling performance of PKC against other coronagraphic concepts. Such experiments are most desirable since PKC, thanks to its achromaticity (versus PM and FQ-PM) and throughput (versus AIC), appears as an attractive solution for the focal instrument of hyper-telescopes using a densified pupil (Pedretti et al. 2000).

Acknowledgements
L. Abe is grateful to R. Krawczyk and Alcatel-Space Industries/Cannes for supporting his Ph.D. fellowship. This paper benefitted from discussions with A. Labeyrie, P. Assus, L. Arnold, Y. Rabbia, E. Thomas and D. Mourard. The authors are indebted to D. Mouillet for his constructive and critical refereeing of the original manuscript.

Appendix A: Mathematical formalism

In order to assess the exact intensity distribution in the different pupil/image planes, in the following we give the mathematical expression of the PKC effect on the amplitude of the diffracted light in these planes. Since the PKC has a privileged direction, namely that of its edge, from a formal point of view, the problem can be studied in one dimension.
  \begin{figure}
\par\mbox{\includegraphics[width=4cm,clip]{Abe10486f9a.eps}\hspac...
...m,clip]{Abe10486f9d.eps} }
\mbox{\hspace*{2cm}{c}\hspace*{4cm}{d}}
\end{figure} Figure A.1: a) Amplitude profile of an Airy pattern, b) the Phase Knife function (PK) applied to the Airy pattern, c) is the Fourier transform of a) and d) the Fourier transform of b). To obtain the resulting amplitude in the conjugate pupil plane c) and d) functions have to be convolved. Graphics are drawn with arbitrary units.

Now, the 2D amplitude distribution of a perfect Airy pattern has the form,

\begin{eqnarray*}a(r) = 2\frac{J_1(r)}{r}
\end{eqnarray*}


where:

\begin{eqnarray*}r = \sqrt{x^2 + y^2}
\end{eqnarray*}


using r, x and y as the coordinates in the image plane.

When we want to $\pi $-shift this amplitude, a(r) is multiplied by the pk(x,y)pk(x,y) = Hx(x,y) - Hx(-x,y)

where Hx(x,y) is the 2D Heaviside pseudo-function along the x axis. Accordingly the resulting amplitude in the conjugate pupil plane can be expressed by Fourier transforming a(r) and pk(x,y) and then convolving them. We note hereafter by small letters the amplitudes in the image plane and in capital letters those on pupil planes.

From the expression of pk(x,y) it can be easily seen that its Fourier transform will reduce to a one-dimensional distribution along the u axis (conjugate of x). Since the Airy function is centro-symmetrical and the pk function is even, the Fourier transform profiles can be written as:

\begin{eqnarray*}A(u) = \Pi_{\rho}(u)
\end{eqnarray*}


and

\begin{eqnarray*}PK(u) = -\frac{u}
\end{eqnarray*}


$\Pi_{\rho}(u)$ is the pupil function profile of width $2\rho$, and PK(u) is a hyperbolic function which presents a singularity at u=0. The convolution of the two functions, noted by $\ast$and using $\tau$ as the variable, leads the amplitude distribution:

 \begin{displaymath}%
F(\tau) = A \ast PK(\tau) = \int \Pi_{\rho}(u)PK(\tau-u){\rm d}u.
\end{displaymath} (A.1)

The calculation of this integral can be wisely split into three intervals. One leads to a null integral since the product of $\Pi_{\rho}$ and PK is odd.
  \begin{figure}
\mbox{\includegraphics[width=4.2cm,clip]{Abe10486f10a.eps}\hspace...
...ip]{Abe10486f10b.eps} }
\mbox{\hspace*{2cm}{a}\hspace*{4cm}{b}}
\par\end{figure} Figure A.2: a) Amplitude profile in the pupil plane as a result of the 1-D convolution of $\tilde(u)$ with A(u) and b) the corresponding intensity profile. It corresponds to the case depicted of Fig. 1. Notice that the curve near the pupil edges is not as deep as the 2-D convolution case (see below Fig. .4c). Plots are given in arbitrary units.

The remaining parts can be classically calculated as the integral of a hyperbolic function, leading to a logarithmic type expression,

\begin{eqnarray*}F(\tau) = {\rm ln}(\rho+\tau) - {\rm ln}(\rho-\tau).
\end{eqnarray*}


We notice that $F(\tau)$ presents two singularities on $\rho$ and $-\rho$, which correspond to the edges of the pupil. This would mean that the intensity at these points is infinite which is physically incorrect. The answer is that the pk(x) distribution has finite boundaries, so we must consider a finite boundaries function, say $\tilde{pk}(x)$ then,

\begin{eqnarray*}\tilde{pk}(x) = H(x) - H_{x+d_(x) - H(-x) + H_{d_(-x)
\end{eqnarray*}


where d0 is an arbitrary size given for the phase knife geometrical extent, which is large compared to the Airy disc. This function and its Fourier transform are plotted in Figs. .1b and .1d respectively, the latter being expressed as,

\begin{eqnarray*}\tilde(u) = i\left(-\frac + 2\cos{(\rho u)}/u\right).
\end{eqnarray*}


Note that $\tilde(u)$ is odd and drops to zero for u=0. The convolution between the two functions $\Pi_{\rho}(u)$ and $\tilde(u)$ gives the amplitude distribution in the pupil plane (Fig. .2).

To assess the amplitude over the whole 2D pupil plane, we need to consider one-dimensional pupil profiles for each v, assuming the dimensions of the phase mask are large enough compared to the airy pattern. The convolution with $\tilde{A}(u)$ (which is also one-dimensional) gives the same result as shown in Fig. .2. This agrees completely with the 2D intensity distribution depicted in Fig. 1, bottom-left.

  \begin{figure}
\par\resizebox{7.5cm}{!}{\includegraphics{Abe10486f11.eps}}\end{figure} Figure A.3: 3-D representation of the equivalent quadrant phase mask. The phase goes alternatively from $-\pi /2$ to $+\pi /2$.


  \begin{figure}
\includegraphics[width=8.8cm,clip]{figA4.eps}\end{figure} Figure A.4: Amplitude of the Fourier transform of the equivalent quadrant phase mask a). This function is odd in a symmetry-4. When convolving this function PK(u) with the pupil function P(u), it contributes to a very low resulting intensity inside the pupil area b) and c). Top-right represents the modulus of the complex amplitude in the pupil plane. d) shows this modulus encircled by the pupil area. These figures were numerically computed using an exact convolution, instead of FFTs.

These expressions can be easily extended to the second $\pi $-shifting step in the orthogonal direction necessary for cross-cancelling the diffracted light in the conjugate pupil. To obtain a strong light cancellation, we have to remove the Lyot stop in the first intermediate pupil plane. Indeed the physical effect of two orthogonal phase knives result in a four quadrant mask, where the phase goes alternatively from 0 and $-\pi$ or where the amplitude is alternatively left unchanged and sign-switched (Fig. .3). The mathematical expression of the Fourier transform of the FQ-PM has been calculated and its 2-D representation is shown in Fig. .4 top left. One can notice that this function is also odd along two orthogonal directions. Since the 2-D mathematical extension of relation .1 has no simple expression and would require further study, we numerically computed the convolution operation using the exact mathematical formulae for the 2-D expressions $\Pi_{\rho}$ and  $\tilde$ functions (see Fig. .4). It turns out to completely agrees with results from FFT based simulations.

References

 


Copyright ESO 2001