Issue 
A&A
Volume 566, June 2014



Article Number  A31  
Number of page(s)  13  
Section  Astronomical instrumentation  
DOI  https://doi.org/10.1051/00046361/201323258  
Published online  04 June 2014 
Apodized phase mask coronagraphs for arbitrary apertures
II. Comprehensive review of solutions for the vortex coronagraph
^{1}
Univ. Grenoble Alpes, IPAG,
38000
Grenoble,
France
email:
alexis.carlotti@obs.ujfgrenoble.fr
^{2}
CNRS, IPAG, 38000
Grenoble,
France
^{3}
Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD
21218,
USA
^{4}
European Southern Observatory, Alonso de Cordova 3107, Vitacura, Santiago, Chile
Received:
16
December
2013
Accepted:
7
April
2014
Context. With a clear circular aperture, the vortex coronagraph perfectly cancels an onaxis point source and offers a 0.9 or 1.75λ/D inner working angle for topological charge 2 or 4, respectively. Current and nearfuture large telescopes are onaxis, however, and the diffraction effects of the central obscuration, and the secondary supports are strong enough to prevent the detection of companions 10^{3}–10^{5} as bright as, or fainter than, their host star.
Aims. Recent advances show that a ring apodizer can restore the performance of this coronagraph by compensating for the diffraction effects of a circular central obscuration in a 1D modeling of the pupil. Our aim is to extend this work and design optimal apodizers for arbitrary apertures in 2D in order to tackle the diffraction effects of the spiders and other noncircular artifacts in the pupil.
Methods. We fold this analytical result into a numerical optimization scheme that yields hybrid coronagraph designs that combine the advantages of the vortex coronagraph (small in IWA) and of shaped pupils coronagraphs (robustness to central obscuration and pupil asymmetric structures). The transmission of the apodizer is maximized, while constraints are set on the extremum values of the electric field that is computed in chosen regions of the Lyot plane through closed form expressions derived even for topological charges. Optimal apodizers are computed for topological charges 2 and 4 vortex coronagraphs and for telescope apertures with 10–30% central obscurations and 0%, 0.5%, and 1% thick spiders.
Results. We put the results of our numerical optimizations in perspective with the analytical solutions and show that our apodizations converge to the ring apodizations. We then characterize the impacts of the obscuration ratio and the thickness of the spiders on the throughput and the IWA. For the apodized charge2 vortex coronagraph the throughputs are slightly below those of the ring apodized vortex coronagraph, and the inner working angle is mostly unaffected by the apodization. The throughputs of the apodizers for the charge4 vortex coronagraph are higher than those of the ring apodized vortex coronagraph. This effect increases with the obscuration ratio, though the inner working angle does, too, and it ranges between 2 and 3λ/D.
Conclusions. The results presented in this paper show that high contrast at small inner working angles can be obtained with a vortex coronagraph for onaxis telescopes, in spite of the presence of a secondary mirror and its secondary support structures.
Key words: instrumentation: high angular resolution / techniques: high angular resolution / methods: analytical / methods: numerical
© ESO, 2014
1. Introduction
The spectral characterization of Earthlike planets around M, F, G, K stars at a few tens of parsecs from our Sun requires a 10^{7} to 10^{10} contrast at a few tens of milliarcseconds (mas) from the host star in 20% bandwidths. Selfluminous planets, which are only tens to hundred million years young, require a less demanding contrast to be imaged than older planets. Since the observation of the former must be preferably done in the infrared part of the spectrum, the chromatic scaling of the pointspread function (PSF) of the instrument partially compensates for this advantage, however.
Onaxis 30–40 m extremely large telescopes equipped with nextgeneration coronagraphs may provide the contrast, the resolution, and the large number of photons that are mandatory for groundbased observations. Few coronagraphs, however, can efficiently cope with the diffraction effects of the central obscuration and the secondary supports.
Adapting a coronagraph to an arbitrary aperture has been an intense research topic in the past few years. Soummer et al. (2009, 2011) explain how the spheroidal prolate apodization of an apodized pupil Lyot coronagraph (APLC) can be adapted iteratively to a given aperture. Pueyo & Norman (2013) study a similar problem, but explicitly constrains the contrast in chosen regions of the image plane, while in the previous case the high contrast was adjusted by varying the radius of the Lyot mask. Nonetheless, APLCs suffer from a rather large IWA.
Shaped pupils, which were initially optimized for highcontrast imaging in one dimension (Spergel & Kasdin 2001; Vanderbei et al. 2003; Kasdin et al. 2007), can be numerically optimized in two dimensions (2D) for any telescope aperture (Carlotti et al. 2011; Vanderbei 2012). Their versatility, robustness, and achromaticity make them good candidates for compact coronagraphic instruments: unlike APLCs, shaped pupils do not rely on a Lyot mask or a Lyot stop to create a high contrast, although a field stop is probably mandatory, given the limited dynamic range of detectors. Like APLCs, 2Dshaped pupil coronagraph have relatively large IWA, usually 3–5 λ/D for apertures with obscurations of 10–30% and for a 10^{7}−10^{10} contrast. The smallest IWA are usually obtained at the expense of the size of the discovery space.
All these coronagraphs require an extreme adaptive optics system to correct for the phase and the amplitude aberrations of the wavefront. This system can be composed of a single deformable mirror (DM), but corrections occur on only one side of the image plane (Malbet et al. 1995). As proposed in Pueyo & Kasdin (2007) and demonstrated in Pueyo et al. (2009), a system of two DMs in nonconjugate planes makes it possible to create symmetric dark holes in broadband. Such a system was recently test in JPL’s highcontrast imaging testbed where it was used to obtain a 4 × 10^{9} in monochromatic light (Riggs et al. 2013).
As suggested by Kasdin et al. (2011), it appears that the wavefront control system can be actively used to create high contrast with arbitrary apertures, thus relaxing the specifications of the coronagraph. Pueyo & Norman (2013) show how to use two DMs to spatially redistribute the energy density in the pupil plane so as to artificially decrease the size of the spiders and, for segmented telescopes, the gaps between the segments. Carlotti et al. (2013) detail methods to optimize stroke commands to be sent to a DM to create a 10^{6}−10^{7} contrast with a centrally obscured aperture.
Pupil mapping is the core principle of the phase induced amplitude apodization (PIAA, Guyon 2003) technique. Combined with a complex amplitude focal plane mask (PIAACMC, Guyon et al. 2010), it offers a very promising solution to the problem of small inner working angles as it results in much smaller IWA (down to 0.64λ/D), while offering about 50% throughput. The currently investigated manufacturing process of the focal plane mask involved in this instrumental concept remains challenging, however (Newman et al. 2013).
With a clear circular aperture of diameter D, looking at a wavelength λ, the fourquadrant phase mask coronagraph (4QPM, Rouan et al. 2000) and the vortex coronagraph (VC, Mawet et al. 2005) perfectly cancel an onaxis, unresolved point source, and detect companions as close as 0.9–1.75 λ/D (0.9 for the 4QPM and a vortex with a topological charge 2, and 1.75 for a VC with a topological charge 4).
The 4QPM coronagraph is a second order coronagraph: its offaxis transmission goes as θ^{2}, where θ is the angular distance to the star. This makes the 4QPM coronagraph quite sensitive to jitter and to the finite size of the star. Like the 4QPM coronagraph, a VC with a topological charge 2 is a second order coronagraph. A charge 4 VC is a fourth order coronagraph, however: its offaxis transmission goes as θ^{4} instead of θ^{2}. It is thus much less sensitive to the finite stellar size and to low order aberrations, which would otherwise limit the performance of the VC as they do with the 4QPM.
Another fourth order phase mask coronagraph is the eightoctant phase mask (8OPM, Murakami et al. 2008; Carlotti et al. 2009), which is to the 4QPM, as the charge 4 VC is to the charge 2 VC. The 4QPM and the 8OPM attenuate offaxis sources along their phase edges, however, and this limits the extent of the discovery space space around the star, especially in the case of the 8OPM.
Nonetheless, the performance of the 4QPM coronagraph and the VC is greatly reduced when the telescope is onaxis. For instance numerical simulations predict a 10^{4}−10^{5} contrast between 1 and 5 λ/D from the star, for a VC used with a 14% centrally obscured aperture such as one of the Very Large Telescope (VLT) 8mclass unit telescopes (UT). Onsky results have been obtained with a VC installed on VLT/NACO (Mawet et al. 2013a). Recently the VC has also been installed at two other telescopes: the Subaru telescope, and the Large Binocular Telescope.
A dualstage coronagraph can be used to cancel the diffraction effects of a circular central obscuration with a 4QPM (Galicher et al. 2011) or a VC (Mawet et al. 2011b), but it can only partially attenuate those of the secondary supports (or any other artifacts in the pupil plane). It also requires twice as many components, making the alignment more difficult and increasing the size of the coronagraph.
It is possible to apodize the aperture to avoid the diffraction effects that prevent the 4QPM or the VC to create highcontrast at a small IWA with an onaxis telescope.
As a matter of fact, a proof of concept for an apodized VC has already been successfully tested on the sky: a small, clear circular subaperture with a 1.5 m diameter has been used at the Palomar telescope with a VC (Serabyn et al. 2010; Mawet et al. 2011a). Unfortunately, this results in an effective resolution three times as small and a throughput eight times as small as what could have been obtained with the main telescope aperture.
It was shown in Carlotti (2013) that shaped pupils can be optimally designed for a given combination of an arbitrarilyshaped aperture and a phase mask. Examples of such designs were numerically optimized for one of the 8 m unit telescopes of the VLT, and for a 4QPM coronagraph, creating a few 10^{10} contrast at 1 λ/D (41 mas in the Hband) from an unresolved star with a system throughput of about 64%.
Mawet et al. (2013b) shows that a ringapodized vortex coronagraph (RAVC) can perfectly attenuate the onaxis light with a circular aperture with a circular central obscuration. Secondary supports are not taken into account, however. For a charge 2 VC this apodizer is composed of two rings: one is fully transmissive while the other is only partially transmissive. For a topological charge 4 vortex, a thin dark annulus separates the two previously described rings. The obscuration ratio of the aperture sets the radii of the rings and transmission of the outermost ring for which the throughput of the coronagraph is maximal.
The Lyot stop is directly constrained by the optimal rings radii. The same property characterizes the apodized 4QPM: the throughput of the coronagraph depends on the choice of the Lyot stop, namely the ratio of its central obscuration.
Fig. 1
Optical layout of an apodized vortex coronagraph. The apodizer, vortex phase mask, Lyot stop and camera are located in the successive pupil and image planes A, B, C, and D. 
Figure 1 displays the optical layout of an apodized vortex coronagraph: an apodizer is located in a pupil plane A. In a subsequent image plane B, the Fourier transform of the electric field that goes through the apodizer is multiplied by the complex amplitude of the phase mask. In the reimaged pupil plane C, i.e., the Lyot plane, the Lyot stop blocks the diffracted onaxis light. Finally, the science camera could be located in a reimaged focal plane D.
Shaped pupils can be numerically optimized for the VC as they have already been for the 4QPM coronagraph, but the complexity of the computation process is much higher than in the case of the 4QPM. This is due to the fact that the pupiltopupil transform described in Eq. (11) of Carlotti (2013) cannot be written as nested sums. It must be directly evaluated as a 2D sum instead of 2 nested 1D sums.
As a consequence the memory requirements scale as the fourth power of the number of points chosen to discretize the pupil, while they only scale as the square for the case of the 4QPM. While this does not fundamentally preclude any computation, it makes the optimization of pupil apodizers much more demanding in terms of numerical resources, than their 4QPM counterparts (Carlotti 2013).
In this paper, we build upon and devise new methods to solve the limitations of the early papers (Carlotti 2013; Mawet et al. 2013b) so as to design apodizers for 2D apertures and a VC.
A new computer dedicated to solving large scale optimization problems has made possible the computation of apodizers for the VC, with twofold symmetry masks computed over 512^{2} points in a few hours. This computer has a fourcore 3.6 Ghz processor and 64 GB of RAM. AMPL, a mathematical programming language is used to transcribe the optimization problem in a language that can be understood by one of the available solvers. In addition to the LOQO solver, the CPLEX and Gurobi solvers have been used. LOQO does not currently use more than one core of a processor, contrary to CPLEX and Gurobi, which can both use multiple cores. It should be noted that AMPL itself does not support parallel computing either.
To compute apodizers for the VC, we also had to modify the nature of the optimization problem. Indeed, there is a difference in the constraints set in the numerical optimization of apodizers for the VC and for the 4QPM coronagraph. Linear constraints were used for the 4QPM in Carlotti (2013). This was possible because, with an aperture with two axes of symmetry such as the VLT’s, the 4QPM creates a purely real amplitude in the Lyot plane – and applying the constraints on the amplitude or the intensity does not make a difference. This is not the case with the VC: constraints must be set on the intensity of the electric field, which makes them quadratic.
We present in this paper the properties of apodizers optimized for a charge 2 VC and charge 4 VC with a 32λ/D outer working angle, and designed to attenuate the onaxis light that goes through the Lyot stop by a 10^{6} factor, which result in a 10^{8} − 10^{10} contrast at 1–3λ/D from the star depending on the topological charge and the telescope aperture. An exhaustive number of apertures have been considered, with 10–30% central obscurations, and for each of them 0, 0.5 and 1% thick spiders. We compare the apodizers computed for spiderfree apertures to the ringapodizers, and we study the impact that an increasing spider thickness has on the performance metrics of the coronagraph.
Section 2 details the mathematical formalism used to derive the closed form expressions of the pupiltopupil transforms used to compute the electric field in the Lyot plane. Section 3 presents the optimization problem, and explains the methodology that is followed in the rest of the paper. Section 4 details the throughput and the inner working angle of the numerical optimizations. For spiderless apertures it compares them with the results presented in Mawet et al. (2013b). It also addresses the limitations of broadband observations. Section 5 tackles the manufacturing aspects of this coronagraph. Section 6 draws a conclusion to this paper.
2. Analytical expressions for the apodized VC
2.1. Formalism of apodized phase mask coronagraphs
Here we briefly recall the formalism introduced in Carlotti (2013) and extend it to the case of arbitrary charge vortex The electric field in the image plane before the phase mask is a function of the transmission of the apodized aperture A(x,y): (1)If M(u,v) describes the phase mask, then the electric field right after the mask is written as the product of E(u,v) and M(u,v). The expression of the electric field in the subsequent pupil plane is the Fourier transform of that product: (2)Equations (1) and (2) can be combined into a single expression: (3)The Fourier transform F(x,y) of the mask appears if Eq. (3) is written in a slightly different way:
(4)Equation (4) can thus be seen as being proportional to the convolution product of the aperture A(x,y) and the Fourier transform of the mask F(x,y): (5)If a closed form for the function F(x,y) can be derived, then Eq. (5) can be used to compute the electric field in the Lyot plane without explicitly computing the electric field in the intermediate image plane. The calculation of F(x,y) depends on the chosen phase mask M(u,v). It is more convenient to write the expression of F(x,y) as a double integral with respect to ρ and θ, where ρ and θ are the radial distance and the azymutal angle in the (u,v) plane. For vortex phase masks, the function M(ρ,θ) does not depend on ρ, and equals e^{iθl}, where l is the topological charge. This property makes it possible to proceed first with the integration with respect to ρ:
(6)Note that L represents the radius of the image plane mask (in units of λ/D). While it is usually assumed that the mask is infinitely large, it is necessary to specify its size in order to derive the closed form expression for Ψ(r) given in Eq. (6).
2.2. General expression for even topological charges
In order to derive this result, we start by rewriting Eq. (6), but instead of first integrating over ρ, we integrate first over θ:
(7)where we have first written x = rcosφ, y = rsinφ, and then proceeded with the change of variables . We recognize here an expression which is closely related to the definition of the Bessel function of the lth order. The expression of F(x,y) becomes: (8)and the final expression for F(x,y) is obtained after integrating over ρ: (9)where F is the generalized hypergeometric function.
2.3. Expression for the first few even topological charges
It is not difficult to derive for the first few even topological charge closed form expressions which only use Bessel functions. The expression of F(x,y) has previously been derived for the 4QPM and a VPM of topological charge l = 2: (10)In the case of a topological charge l = 4, the integration over θ gives the following result: (11)And for a topological charge l = 6, this expression becomes: (12)
3. Numerically optimized apodizers
Using Eqs. (10) and (11) derived in the previous section, and proceeding with the same discretization process as described in Carlotti (2013), apodizers have been optimized for different apertures. The mask is discretized over N points along each axis of the pupil plane. We assume a unit aperture diameter, and the distance between nearest neighbors is .
3.1. Set up of the optimization problem
The problem consists in maximizing the total transmission of the mask under the following constraints set on the extremum values of the transmission of the apodizer, and of the amplitude of the electric field in the Lyot plane: (13)where Δ_{A} is the region defined by the telescope’s aperture, Δ_{C} is the region defined by the Lyot stop, and c measures the attenuation of the intensity in the reimaged pupil plane in a logarithmic scale.
Since the vortex phase mask creates a complex amplitude in the Lyot plane, using linear constraints would force valid phasers in this plane to be located inside a square of size ×10^{−c/ 2}. To explore all valid solutions, one must let the valid phasers live inside a larger area: a circle of radius 10^{−c/ 2} (which circumscribes the previous square). This corresponds to the quadratic constraint that appears in Eq. (13).
Fig. 2
Example of a design for the Keck telescope, a charge 4 vortex coronagraph, and a circular Lyot stop: a) apodizer transmission; b) intensity in the Lyot plane (log scale); and c) comparison of the Keck aperture and the apodizer. In the last figure the apodizer is represented in white, while the Keck aperture is represented in gray (wherever the apodizer transmission is zero). 
3.2. Results in the case of noncircular apertures
Before studying exhaustively the properties of apodized vortex coronagraphs for various apertures, we first have to comment on the importance of the choice of the Lyot stop, and how the vortex phase mask constrains this choice. In particular, we want here to consider the case of noncircular apertures.
Many major groundbased and spacebased telescopes have circular apertures (for instance the VLT’s UTs, or the Hubble Space Telescope, among many others). The two Keck telescopes are segmented, however, and the shape of their outer edge is not circular (although their central obscuration is). The same can be said of the Gran Telescopio Canarias (GTC), which has a very similar pupil.
Other segmented telescopes are the future major goundbased telescopes such as the European Extremely Large Telescope (EELT), the ThirtyMeter Telescope (TMT), and a potentially large segmented space telescope. We should also mention the James Webb Space Telescope – whose primary mirror is also segmented – but the design of its coronagraphic instruments has long been finalized, and no modification can be expected.
The formalism that we have introduced allows us to design apodizers for any type of aperture, and we have done so by considering the aperture of the Keck telescopes. Both primary mirrors are composed of 36 hexagonal segments of 1.8 m (corner to corner). The size of this aperture is close to 10.9 m along its longest diameter, and 9 m along its smallest diameter, and it is commongrayly accepted that the aperture mean diameter is 10 m.
We have found that in the best case apodizers only have a 1–2% transmission if the Lyot stop transmissive area encompasses the outermost segments (those that lie in the 9−11 m diameter ring).
On the contrary, if the Lyot stop does not encompasses these outermost segments, i.e., if the Lyot stop outer edge is circular, then the transmission of the apodizers becomes significantly higher.
This is not too surprising: in its original design, the vortex coronagraph assumes a circular aperture noncircular artifacts, such as the spiders or a noncircular outer edge, make it difficult to retrieve a high contrast. While an apodizer can be used to help restoring a high contrast, it cannot entirely compensate for the fact that the vortex phase mask works preferentially with circular clear aperture.
As it is shown in Fig. 2, the transmission of an apodizer computed for the Keck telescope and a circular Lyot stop is almost zero everywhere outside the circle sets by the Lyot stop’s outer edge. No constrains was explicitly set on the apodizer, and the apodizer is only indirectly constrains by the Lyot stop properties.
With Keck, blocking light in the outermost segments means a maximum throughput 20% smaller than the throughput associated with the entire pupil, and a mild resolution loss of about 8% (assuming a 10 m diameter for the aperture).
The apertures of the EELT and the TMT are composed of a much higher number of segments, however, and making the outer edge circular for them does not come with the same throughput or resolution limitation. Besides, in the case of EELT, the secondary mirror will be circular. Its diameter will be such that it will only capture the inner, largest circular area possible, making the effective diameter of the telescope 37 m.
Since it appears necessary to restrict the Lyot stop’s outer edge to a circle, and that it indirectly constrains the apodizer transmission to be zero outside the same circle, we have chosen to only consider circular apertures and circular Lyot stop in the rest of this study.
3.3. Results in the case of circular apertures
Five linear obscuration ratios of 10, 15, 20, 25, and 30% have been considered. The secondary support structure is orthogonal and their thickness is either 0, 0.5, or 1%. The first case corresponds to a spiderfree centrally obscured aperture, which is useful to consider when comparing the properties of the numerically optimized apodizers with those of the ringapodizer which too are designed for an aperture with a central obscuration but no spiders. The 0.5 and 1% thicknesses have been chosen to match the thickness of current telescopes’ spiders. For instance the spiders of the Hale telescope at the Palomar observatory are 0.25% thick, and those of the 8 m unit telescopes at the VLT are 0.5% thick.
Overall, the Lyot stops’ geometry and the apertures’ geometry are similar. The diameter of the Lyot stops is 96% that of the aperture diameter. This is meant to prevent the diffraction effects of the finite size of the vortex phase mask. For the same reason, the secondary supports of the Lyot stops are twice as large as those of the apertures. The main difference between the Lyot stops and the apertures is the size of the central obscuration of the Lyot stops, which is larger for the Lyot stops than for the apertures, and is a key parameter that must be adjusted when looking for an optimal apodizer.
The two axes of symmetry displayed by the apertures and the Lyot stops are used to reduce the memory required during the optimizations. The transmission of the apodizers was optimized in one quadrant of the pupil plane over N = 256 points (along each axis). The amplitude was computed in one quadrant of the Lyot plane over M = 64 points. As the vortex mask was assumed to be 32λ/D in diameter, this number of points is enough to satisfy the NyquistShannon sampling theorem.
Note that for this number of points up to a third of the available memory was used. While further increasing it is possible, it was enough to assess the interest of the method by computing a large number of apodizers.
4. Performance
4.1. System throughput
As it was observed with the apodized 4QPM coronagraph and with the ringapodized VC, the throughput of the system depends on the size of the central obscuration of the Lyot stop: an optimal value exists for which the throughput is maximized.
4.1.1. Comparison with the ring apodizers
The ringapodized vortex coronagraph proposed in Mawet et al. (2013b) creates highcontrast in spite an aperture with a circular, central obscuration. The ring apodizers are formed of 2 or more rings with different transmissions. For instance, a charge 2 vortex requires a first ring with a t = 1 transmission, surrounded by a 0 <t< 1 second ring. In the case of a charge 4 vortex, a t = 1 first ring is surrounded by a t = 0 second ring, which itself is encircled by a 0 <t< 1 third ring. The transmissions and radii of the rings depend on the obscuration ratio of the aperture.
Fig. 3
Comparison of the throughputs estimates of the ring apodized VC (dotted lines) and the shaped pupil apodized VC (solid lines) for topological charges 2 and 4. Different measures of the throughput are illustrated. Pupil plane estimates are showed in red, while image plane estimates are showed in blue. The dotted black line represents the throughputs given in Mawet et al. (2013b). 
Ring apodizers offer an elegant and relatively simple approach to a complicated problem, but they do not take the diffraction effects of the secondary supports into account. Moreover both the central obscuration and the shape of the pupil are assumed to be circular, which is not the case of all telescopes, although we showed that numerically optimized apodizers designed for noncircular apertures were forced to be quasicircular by the VC.
The numerically optimized apodizers and the ring apodizers are designed to solve the same general problem, however, and it is interesting to try and compare how similar or not they can be.
A first element of comparison comes from looking at the optimal apodizers that have been computed for the spiderfree aperture. These apodizers look very much like the ring apodizers. The shaped pupils can also be divided into two or three areas, the only difference being that the transmission of the outer annulus is binary with the shaped pupils, whereas it has a gray transmission everywhere in the annulus of the ringapodizers.
Figure 3 shows how the system’s throughput of the ring apodizers change with the obscuration ratio of the aperture. The corresponding quantities for the optimal apodizers computed for the five 10–30% central obscuration spiderfree apertures are displayed for comparison.
Two types of throughput estimates are distinguished: pupilplane estimates and imageplane estimates. Pupilplane estimates are the ratio of the total planet intensity going through the Lyot stop, and the total planet intensity going through the plane located immediately before the apodizer. Imageplane estimates are the ratio of the total planet intensities measured in a small circular area centered on the planet with and without the apodized coronagraph. We chose this circular area to be 2λ/D in radius.
For comparison, the throughputs given in Mawet et al. (2013b) are computed in the pupilplane, although this computation was not performed by propagating light in a simulated system, but rather by relying on the analytic description of the apodized coronagraph.
A first observation is that the throughputs given in Mawet et al. (2013b) for charges l = 2 and l = 4 are almost identical to the pupilplane estimates we have computed by propagating a planetary companion in a simulated coronagraph. The very small differences that were observed are very likely due to the fact that the size of the VC had to be finite, and was set to 256 λ/D.
In the rest of this comparison analysis, we refer to the system throughput as the imageplane estimate of the throughput.
In the case of the charge l = 2 VC, the system throughput associated with the numerically optimized apodizers is always lower than what it is with the ring apodizers. The difference gets smaller with the central obscuration, however. The system throughput for the ringapodizers is 72% for a 10% central obscurations CO_{A}, and 30% for CO_{A} = 30%, and this decrease is quite linear. As for the system throughput for the numerically optimized apodizers, it is 62% for CO_{A} = 10%, and 27% for CO_{A} = 30%. It is thus 10% smaller for the smallest obscuration, and 3% smaller for the largest obscuration.
On the other hand, in the case of the charge l = 4 VC, the system throughput associated with the numerically optimized apodizers is smaller than that of the RAVC for very small obscurations, but it is significantly larger for large central obscurations. This difference goes from being − 4% for CO_{A} = 10%, to being 0 for CO_{A} = 13%, to being +10% for CO_{A} = 20%, and +17% for CO_{A} = 30%. For CO_{A} = 10% it is 56%. It is 41% for CO_{A} = 20%, and it is 30% for CO_{A} = 30%.
The fact that (a) the focal plane masks have a finite size (32λ/D), and (b) the remaining onaxis energy inside the Lyot stop is small but not zero may be the main reasons that explain why the throughput of the apodized charge 4 VC is higher than it could have been expected. This is outside the scope of this paper, but we intend to address this point in future work. One of the difficulty preventing us to present results on this topic in this paper comes from the current impossibility to compute apodizers for larger OWA than 32λ/D, as it requires a RAM quantity larger than what is available with the computers we have access to at the moment.
We already mentioned that, when looking for the values of CO_{LS}, we started with Lyot stops prescribed using the analytical expression of the ring apodizers and found a local optimum in that neighborhood. While not showed in Fig. 3, we have found that the CO_{LS} values are very close to their counterpart for the ring apodizers for l = 2, and close to but increasingly smaller for l = 4.
This confirms that the formalism used to design ring apodizers can and should be used to infer a first estimate of the obscuration ratio of the Lyot stop used in the numerical optimizations.
4.1.2. Additional impact of the spider thickness
To find the optimal obscuration ratio, we have computed apodizers for several Lyot stops, starting with the obscuration ratio predicted by the closed form expressions derived in Mawet et al. (2013b), and iterating until a maximum was reached.
Tables 1 and 2 list – for the topological charges l = 2 and l = 4 – the parameters of the optimal Lyot stops and of the apodizers that correspond to them, for the five apertures and the three spider thicknesses that were considered. To reach the maximum throughput, larger central obscurations CO_{A} require larger Lyot stop obscurations CO_{LS}, and this decreases both the transmission of the apodizer T_{A} and the maximum throughput of the coronagraph T_{max}.
The tables also display the transmission of the Lyot stop T_{LS}. It should be noted that the product of T_{A} and T_{LS} – which could be seen as the expected throughput of the coronagraph – is about 10% smaller than T_{max}: some of the offaxis light gets diffracted outside the Lyot stop, and is thus blocked by it.
Charge 2 VC: main parameters for the combination of apodizers and Lyot stops that give the highest system throughputs for the five aperture diameters CO_{A} and the three spider thicknesses t_{S}.
Increasing CO_{A} from 10 to 30% decreases T_{max} by about a third for the charge 4 AVC, no matter the spider thickness. For the charge 2 AVC, T_{max} decreases by 50% for the spiderless apertures, by 40% for t_{S} = 0.5%, and by 23% for t_{S} = 1%.
Given the same central obscuration, thicker spiders, however, greatly decrease the coronagraphic throughput. The impact of the spider thickness t_{S} on the throughput T_{max} increases with t_{S}: while there is only a 5–8% decrease of the throughput between coronagraphs designed for a spiderfree aperture and designs with 0.5% thick spiders, there is a 14–19% decrease of the throughput between the latter and 1% thick spiders. The smaller the central obscuration of the aperture, the larger the decrease.
Charge 4 VC: main parameters for the combination of apodizers and Lyot stops that give the highest system through puts for the five aperture diameters CO_{A} and the three spider thicknesses t_{S}.
Figures 4 and 5 show a total of twelve apodizers optimized for charge 2 and charge 4 VC, for apertures with 10 and 30% central obscurations, and 0, 0.5 and 1% thick spiders. The apodizers designed for the spiderfree apertures display concentric dark rings which thickness increases towards the center of the apodizer. The rings appear to be apodizing the central obscuration of the aperture.
For a charge 2 VC, one can noticed two different areas:

a first annulus that surrounds the central obscuration of the aperture, where the transmission is close to 1. The outer diameter of this annulus is set by the Lyot stop for which the mask is designed.

a second annulus than encircles the first, where, for small central obscurations, small blocking regions are located on a 2D lattice with a 1/OWA periodicity (in units of pupil diameter). Most of these regions are isolated, but some merge into larger structures.
For a charge 4 VC, three areas are noticeable:

a first annulus where the transmission is 1 almost everywhere.

a dark ring that surrounds the first region, which diameter is set by the Lyot stop for which the mask is designed.

a second annulus than encircles the two other regions. Here again, small blocking structures appear on the same 1/OWA periodic 2D lattice.
The presence of spiders creates a more complicated pattern, as the mask apodize both the central obscuration and the spiders. It is interesting to compare the patterns displayed by the apodizers designed for a charge 4 VC, and for apertures with 1% thick spiders and 10 and 30% central obscurations: in the case of the smaller obscuration, thick horizontal and vertical dark strips can be seen next to the spiders, i.e., the apodizer appears to be mainly apodizing the spiders. On the other hand, when looking at the other apodizer, thick dark concentric rings can be seen, and the apodizer appears to be mainly apodizing the central obscuration.
Fig. 4
Apodizers giving the highest system throughput for a charge 2 VC used with a 10% (top) and a 30% (bottom) central obscuration, and three spider thicknesses: 0, 0.5, and 1% (from left to right). 
Fig. 5
Apodizers giving the highest system throughput for a charge 4 VC used with a 10% (top) and a 30% (bottom) central obscuration, and three spider thicknesses: 0, 0.5, and 1% (from left to right). 
4.2. Inner working angle
The main advantage of phase mask coronagraphs is their very small IWA. As noted in previous papers, using an apodizer changes the usual offaxis transmission of the coronagraph.
The maximum throughput T_{max} is reached for sources at around 25–26 λ/D from the star, but T_{max}/ 2 is reached at much shorter distances (≈1–3λ/D depending on the topological charge, the central obscuration, and the spider thickness of the aperture). This halfthroughput is usually used to define the IWA of the coronagraph (as suggested in Guyon et al. 2006).
Fig. 6
Comparison of the relative offaxis transmission of the VC (black line) and an AVC (red line) designed for a CO_{A} = 20% obscured aperture with t_{S} = 0.5% spiders. The comparison is made for a) a charge 2 vortex and b) a charge 4 vortex. The transmissions are normalized with respect to the coronagraph throughput T_{max}. 
A first and important observation must be made: as illustrated in Fig. 6, and contrary to the offaxis transmission curve of a VC, the transmission of an AVC is not monotonic. Because of that, in some cases, using the IWA to fully describe the resolution properties of the AVC may be difficult.
For angular distances θ between 0 and 1–1.2λ/D, the transmission (normalized with respect to the maximum coronagraphic throughput T_{max}) follows a θ^{l} function, with l the topological charge of the VC. This confirms that, like the VC, the AVC is a second or fourthorder coronagraph.
In Fig. 6, the apodizer of the AVC is designed for a 20% obscured aperture with 0.5% thick spiders. The transmission increase rate slows down between 1 and 1.6λ/D for l = 2, and between 1.5 and 2.2λ/D for l = 4. This does not have much effect on the IWA of the charge 2 AVC: it is almost identical to what it is for the charge 2 VC (0.9λ/D). The same is not true for the charge 4 AVC, however: while 40% of the maximum throughput is reached for θ = 1.4λ/D, the IWA of this coronagraph is about 2.5λ/D (while it is 1.75λ/D with a charge 4 VC).
Figure 7 shows how the offaxis transmission for the five different central obscurations and the three different spider thicknesses that were considered. Note that these numbers are also reported in Tables 1 and 2. Again, the transmission slows down at about 1–1.2λ/D, before increasing again at about 1.5–2λ/D. The effect strengthens with the obscuration ratio and the topological charge. For l = 2, the transmission slows down for CO_{A} = 10%, and it stalls for CO_{A} = 30%. For l = 4, it stalls for CO_{A} = 10%, and it decreases for CO_{A} = 30%.
For l = 2 the IWA of the AVC is close to that of the VC for t_{S} = 0% and 0.5%, where the halftransmission is reached for θ = 0.9λ/D. The IWA only differs for t_{S} = 1% for which it increases with the obscuration, ranging between 1.1 and 1.5λ/D. For l = 4 the IWA of the AVC is larger than what it is with the VC. The halftransmission is reached for θ = 1.8–3.2λ/D for t_{S} = 0% and 0.5%, and for CO_{A} = 10–30%. For t_{S} = 1%, it ranges between 2.2 and 3.2λ/D.
Fig. 7
Transmission (in %) of an offaxis source as a function of the distance to the star (in units of λ/D), for a) and d) apertures with no spiders; b) and e) apertures with 0.5% thick spiders; and c) and f) apertures with 1% thick spiders. Figures a), b), and c) correspond to a charge 2 VC; and figures d), e), and f) correspond to a charge 4 VC. In each case, the transmission is given for the five central obscurations: 10, 15, 20, 25, and 30% (from top to bottom). The red dots indicate the IWA of each coronagraph. 
Associated with the apodizer, the vector vortex masks attenuate the onaxis light of the star in a given region of the Lyot plane. Figure 8 illustrates one of the 2D intensity distributions that have been computed. In the dark zones that appears in the Lyot plane, the intensity has been reduced to a few millionth of what it would have been without the apodizer.
The PSF of the coronagraph is computed by taking the Fourier transform of the product of the appropriate Lyot stop and the electric field observed in the Lyot plane. Normalizing this PSF is done by computing the propagation of an offaxis source (located as close to the OWA as it is possible without being attenuated). The intensity of the first PSF is divided by the maximum intensity of the second one. An example of such a PSF is showed in Fig. 9. One can see that a central structure extends up to 12–15λ/D from the star, where the intensity is lower than 3 ×10^{9}. Beyond this distance, the intensity stays close to 3 ×10^{10} (there is no point in looking further away than 32λ/D since that distance is also the radius of the vortex phase mask).
Tables 1 and 2 list the imaging properties of the 15 coronagraphs, and in particular the contrast measured at the IWA (it is an average value computed in an annulus with a mean radius that equals the IWA, and a width of 1λ/D). The contrast in the image plane is only indirectly constrained: the real constraints is set on the intensity of the electric field in the Lyot plane.
The contrast is not as low with the l = 2 AVC than it is with the l = 4 AVC. For l = 2 it is between 5 × 10^{9} and 8 × 10^{8}, and neither the obscuration ratio nor the spiders thickness appear to have a strong influence on the contrast (its mean value is 5 × 10^{8} for the spiderfree apertures, 1.5 × 10^{8} for the t_{S} = 0.5% apertures, and 5 × 10^{8} for the t_{S} = 0.5% apertures).
It can be noticed that the contrast increases with both CO_{A} and t_{S}. The contrast is close to 2–3 × 10^{10} for the spiderfree apertures, except for CO_{LS} = 30% where it is 6× 10^{10}. It becomes about 10 times higher for t_{S} = 0.5%, with a maximum of 7 × 10^{9}. Finally, for t_{S} = 1%, the maximum contrast becomes 2–3 × 10^{8}.
The contrast increase is partially explained by the fact that in the optimizations the intensity measured in the Lyot plane is not normalized by the transmission of the apodizer T_{A}. Nonetheless, T_{A} and T_{max}) do not change enough from one aperture to the other to explain the contrast increase, and the constraint set on the intensity should be adjusted to improve the contrast.
Fig. 8
Example of a 2D intensity distribution in the Lyot plane of the apodized VC. The intensity is displayed on a log scale for an aperture with a 15% central obscuration and 1% thick spiders. 
Fig. 9
Coronagraphic PSF (azymutal average) for the apodized vortex designed for a 15% centrally obscured aperture and 0.5% thick spiders. The intensity is displayed on a log scale. Angular distances are given in units of λ/D. 
4.3. Chromaticity
The most recent iterations of the vortex phase mask are partially achromatic: 10^{8} − 10^{9} contrast have been measured in a 10% band (Mawet et al. 2012). It is important to study the chromaticity of the apodized coronagraph as well.
The apodizer is designed for a focal plane mask with a diameter which is fixed when measured in units of λ/D. Hence, the physical size of the mask should scale with the wavelength, while in fact it stays the same. At shorter wavelengths, the mask will appear to cover a larger area that expected, and the opposite effect will occur at longer wavelengths.
As it can be seen in Fig. 10 a relatively bright ring surrounds the PSF of the optimal apodizer. This is quite understandable since the apodizer modifies the PSF to make it compatible with the vortex phase mask, which has a finite radius that in this case equals 32λ/D. Beyond that angular distance the amplitude of the electric field has no reason to match the requirements of the vortex phase mask. The impact on the performance of the coronagraph is thus expected to be stronger at shorter than at longer wavelengths.
It is also expected that chromatic effects may be weaker for larger vortex phase masks, as the intensity of the ring that surrounds the region for which the mask is designed will get weaker too. The current limitation of the mask radius is only due to the limited amount of RAM in the computer we used to perform the optimization (64GB). Much larger amounts are available in other machines, however, and larger masks could be computed. One should keep in mind that smaller mask radii may be responsible for higher coronagraphic throughputs, however.
Fig. 10
2D a) and azymutal average b) of the PSF of the optimal apodizer designed for a 10% central obscuration, 1% thick spiders, and a charge 4 vortex. In the azymutal average, the red dotdashed line denotes the highest spatial frequency captured by the mask at λ = 1.1λ_{0}, the blue dashed line does the same λ = 0.9λ_{0}, and the dotted green line does it at λ = λ_{0}. 
Figure 11 shows how the contrast changes with the bandwidth for an AVC designed for a CO_{A} = 15%, t_{S} = 0.5% aperture. The indicated contrast is the mean contrast inside an annulus centered on the star, and with a 2.3λ/D inner radius, and 26λ/D outer radius. At λ = λ_{0}, the physical size of the mask corresponds to the 32λ/D for which the mask is designed.
The contrast goes down with the bandwidth. A 10^{7} contrast is obtained with a 34% bandwidth, a 10^{8} contrast with an 18% bandwidth, and a 10^{9} contrast with a 2% bandwidth.
While a 1.5 × 10^{8} contrast in a 20% bandwidth might be enough for a groundbased instrument, additional constraints may have to be set on the Lyot plane intensity to increase the bandwidth and work with lower contrast. The electric field can be computed for two wavelengths instead of a single one, and to do that two mask sizes must be considered in the optimization problem. To make sure that highcontrast is obtained at the two extremities of a 20% bandwidth, one could for instance compute the electric field for 32 and 25.6λ/D focal plane masks. Additional constraints may have to be set inbetween these two extremum values.
Fig. 11
Mean contrast as a function of the bandwidth (in %) for a charge 4 VC, and an aperture with 15% central obscuration and 0.5% thick spiders. The contrast is measured between the IWA and the OWA of the coronagraph, inside an annulus. 
5. Manufacturing
The apodizers displayed in this paper are not absolutely binary, but it has been shown in previous papers (Carlotti et al. 2011; Carlotti 2013) that solutions to these optimization problems tend to be binary, i.e, that the mask transmission can be rounded without affecting the contrast as long as enough points, say 512 in one quadrant of the pupil plane for a 10^{9} contrast, are used to discretize the apodizer.
A binary apodizer can be manufactured as a quasiachromatic component. The most common way to manufacture these devices uses photolithography. Several microdots apodizers (Martinez et al. 2009; Sivaramakrishnan et al. 2010) have been manufactured with this technique for the coronagraphs of the gemini planet imager (GPI, Macintosh et al. 2007) project, and the spectropolarimetric highcontrast exoplanet research (SPHERE, Beuzit et al. 2008) project.
Laboratory experiments conducted for the development of the optics of the SPICA coronagraphic instrument (Enya et al. 2008) showed that transmissive binary apodizers could provide a 8 × 10^{8} contrast in the visible. The substrate that is used can introduce chromatic wavefront errors, as well as internal reflections. The smallest features in the binary apodizers developed for SPICA are similar in size to those of the apodizers presented in this paper.
Several 2D shaped pupils have been manufactured in the last few years. The size of their smallest feature is virtually identical to that of the binary apodizers developed for SPICA. In particular two masks have been designed for the Subaru telescope, and they are now part of the SCExAO instrument (see Martinache & Guyon 2009; Carlotti et al. 2012). Both have been manufactured using photolithography, with a resolution of about 10 μm (with about 1800 points across the apodizer). The high reliability of the etching process is illustrated in Fig. 12. The smallest features are reproduced with fidelity, including 1pixel large details.
Fig. 12
Comparison of the theoretical and actual transmissions of a portion of a 2D shaped pupil manufactured for the Subaru telescope. Left: microscope photography of the mask. Right: corresponding desired profile. Each pixel has a 10 μm size. 
Contrasts lower than 10^{8} will require other types of binary apodizers, however. They must either be substrateless (and thus freestanding) or reflective.
The possibility of using MEMS or MOEMS devices such as microshutter arrays or micromirror arrays is particularly compelling at it would allow a dynamic control of the pupil plane amplitude, in addition to replacing several apodizers with a single device. This is especially interesting when considering the absence each night of several randomly selected mirrors in the pupil of the EELT (and potentially the TMT and other future telescopes with a large number of segments). Although it is not clear whether this technology is currently mature enough to be used for highcontrast imaging, such devices have been extensively tested for spectrophotometry instruments. Located in an image plane they allow the dynamical selection of multiple targets.
Deep reactionion etching – the same technique used to manufacture JWSTNIRSpec’s microshutter arrays – can also be used to manufacture transmissive, freestanding masks. The ripple masks used in the highcontrast imaging laboratory at Princeton University, and in the highcontrast imaging testbed (HCIT) at the jet propulsion laboratory (JPL) were manufactured using this technique (Belikov et al. 2007). In this case small openings of a few tens of microns in size would be etched in a silicon wafer, with an orthonormal grid overlaying the mask to make it freestanding (and taken into account in the design process). Such masks would be achromatic, and would not introduce any wavefront aberration. Because of the nonzero width of the overlaying grid, the transmission of the mask would be slightly reduced, however.
Finally, reflective apodizers are currently in development at JPL (Balasubramanian et al. 2013) in the context of the shaped pupil coronagraph developed for WFIRSTAFTA. A black silicon coating is deposited on a flat, reflective substrate to mask the nonreflective area of the apodizer.
6. Discussion and conclusion
This paper builds upon and devises new methods to solve the limitations of two early papers:

Mawet et al. (2013b) presented 1D ring apodizers analytically optimized to create highcontrast with a VC and a 1D circular, centrally obscured aperture, without spiders.

Carlotti (2013) presented 2D apodizers numerically optimized to create highcontrast with a 4QPM coronagraph and a 2D centrally obscured, arbitrarily shaped aperture, with spiders.
We have showed that apodizers can be numerically optimized in 2D to help restoring the highcontrast imaging capabilities of vortex coronagraphs when used with obscured apertures with arbitrary shapes and spiders. While the ringapodizers have gray transmissions, these apodizers resemble the shaped pupils in that they have binary transmissions, and that they can be manufactured using the same processes.
Closed form expressions for even topological charges have been derived and used to directly compute the electric field in the Lyot plane as a function of the electric field in the first pupil plane. These pupiltopupil transforms have the advantage of not explicitly sampling the intermediate image plane and the vortex phase mask located there. One of the constraints of using these transforms is that the angular extent of the mask is limited. Given the computer currently used, it is possible to numerically optimize apodizers that are discretized over several hundred points on each axis of the pupil plane.
Because it is 2D, our formalism can be applied to the case of a segmented, noncircular aperture such as the pupil of one of the Keck telescopes. In practice the apodizer transmission is almost zero everywhere outside the inscribed circle defined by this noncircular aperture. Moreover, to maximize the transmission of the apodizer, the Lyot stop for which the apodizer is optimized must also be limited to that inscribed circle. This limitation most probably comes from the properties of the vortex mask, which is initially supposed to be used with a circular aperture. Apart from the apertures of the Keck telescopes, the GTC, or the JWST, only the apertures of the TMT and the EELT feature a noncircular outer edge, but their large number of segments make them almost circular, and this limitation does not affect them much.
Hence, we have considered circular obscured apertures with spiders. Apodizers have been computed for five different central obscuration ratios: 10, 15, 20, 25, and 30%. For each of these obscurations, we have considered three different spider thicknesses (0, 0.5 and 1%). The spiders form an orthogonal pattern. Apodizers have been optimized for charge l = 2 and charge l = 4 vortex phase masks. We were able to compute apodizers for 64λ/D OWA masks, but we have chosen to set the OWA to 32λ/D since the substantially shorter time taken by the optimizations, and the smaller required amount of RAM, make it possible to solve several optimization problems at the same time, and thus explore a large number of different cases. This was necessary to exhaustively characterize the performance of the coronagraph in terms of throughput, IWA, and contrast.
In the optimization problems that we solve, the constraints are set on the intensity of the electric field computed using the pupiltopupil transform. We have chosen to set the attenuation of the onaxis light that goes through the Lyot stop to 10^{6}. In practice this creates a 10^{8} − 10^{10} contrast in the image plane depending on the transmission of the apodizer. The Lyot stops for which the coronagraphs have been designed have also been optimized, although their optimization has not been performed together with the optimization of the apodizers. Like all the apertures that we have considered, all Lyot stops are circular.
In the case of the spiderfree apertures, the results of the numerical optimizations are similar to the ringapodizers computed using the closed form expressions derived in Mawet et al. (2013b). Both types of apodizers are characterized by multiple rings with different mean transmissions. A similar aspect was found in our masks: the radii of their rings are close, and so are their transmissions.
For topological charges l = 2, the throughputs are slightly smaller than those of the ringapodized vortex coronagraphs. Like the charge l = 2 vortex phase mask, the 4QPM is a θ^{2} coronagraph, and it is interesting to make an indicative comparison between the apodizer computed for the 4QPM in Carlotti (2013) for the aperture of the VLT. The throughput of this optimal apodized 4QPM is 65%, while the throughput of the optimal solution for an AVC computed for the closest look alike aperture in the present paper (with a 15% central obscuration, and 0.5% thick spiders) is 46%, which is substantially lower. This difference appears solely associated with the type of focal phase mask.
For topological charges of l = 4, the throughputs of the optimal solutions are close to or larger than those of the RAVC, especially for large central obscurations. For instance, for a 20% obscured aperture, the throughput of the numerically optimized apodized l = 4 VC is 42%, while it is about 32% with the RAVC. For a 30% obscuration, these throughput becomes 30% and 13%, respectively. We explain this difference by (a) the finite OWA of our vortex masks; and (b) the nonzero onaxis intensity that remains in the Lyot plane. The influence on the performance of the coronagraph of these two degrees of freedom has not been characterized yet, and doing so will be one of our next objectives.
Unlike the 1D apodizers showed in Mawet et al. (2013b), the 2D apodizers that we have presented take the aperture spiders into account. We chose to consider two nonzero thicknesses: 0.5 and 1% of the diameter, the former being the thickness of the secondary supports of the 8m unit telescopes at the VLT. The throughput changes significantly with the spider thickness. While the throughputs associated with the 0.5% thick spider apertures are – for both topological charges – about 3–10% and 5–6% smaller than those associated with spiderfree apertures, the throughputs associated with the 1% thick spider apertures are 9–33% and 16–28% smaller, respectively.
In addition to the horizontal and vertical alignments that are required by apodized coronagraphs, taking the spiders into account requires a clocking alignment as well. The spiders for which the apodizer is designed can be oversized with respect to the aperture spiders to accommodate clocking errors.
For small offaxis angles θ, the transmission of offaxis sources follows a θ^{l} function, which confirms the fact that, like the VC, the AVC is an lth order coronograph. The offaxis transmission of the AVC is not as monotonic as it is with the VC, however. The rate at which the transmission increases slows down, stalls or even decreases at around θ = 1 − 1.2λ/D, and it only resumes its maximum value at around θ = 1.5 − 2λ/D. The strength of this effect goes with the central obscuration.
For l = 2, this does not affect the IWA much: it remains close to 0.9λ/D, except for the 1% thick spiders for which the IWA becomes 1.1–1.5λ/D. The same is not true for l = 4. In this case the IWA of the AVC is larger than it is with the VC, and it is mostly affected by the obscuration ratio. It is close to 1.75λ/D for a 10% obscuration, and increases to 2.5λ/D for a 20% obscuration, and to 3.2λ/D for a 30% obscuration. It should be noted, however, that it is difficult to describe the resolution properties of the coronagraph by solely referring to the usual definition of the IWA, i.e., the distance at which 50% of the light of the offaxis light is transmitted. In particular the angle associated with a 40% transmission – about 1.4λ/D – is very similar to that of the VC for obscuration values up to 20%.
The active compensation of aperture discontinuities (ACAD) presented in Pueyo & Norman (2013) is a pupil mapping concept that uses two DMs to reduce the effective thickness of the spiders without loosing photons. ACAD could be used to implement the ringapodizers presented in Mawet et al. (2013b). The apodizers presented in this paper are a compelling counterpart to ACAD, which could be affected by the malfunction of some of the DM actuators. In this view, they offer a more conservative solution at only a moderate cost in throughput. Because the spider
thickness has an important influence on the throughput, apodizers designed for the VC could also take advantage of ACAD.
Although shaped pupils are binary apodizers – which are inherently achromatic – these apodizers are designed for a fixed mask radius measured in units of λ/D, and the fixed physical radius of the mask causes the PSF at different wavelengths to be attenuated in different ways. The chromatic effects may be small enough in some cases. For instance a 1.5 × 10^{8} maximum contrast is obtained in a 20% bandwidth for an aperture with a 15% central obscuration and 0.5% spiders (resembling that of the VLT’s UT). To make sure that lower contrasts are obtained in a large bandwidth, apodizers will have to be optimized for multiple mask sizes.
References
 Balasubramanian, K., Wilson, D., White, V., et al. 2013, in SPIE Conf. Ser., 8864 [Google Scholar]
 Belikov, R., Give’on, A., Kern, B., et al. 2007, in SPIE Conf. Ser., 6693 [Google Scholar]
 Beuzit, J.L., Feldt, M., Dohlen, K., et al. 2008, in SPIE Conf. Ser., 7014 [Google Scholar]
 Carlotti, A. 2013, A&A, 551, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Carlotti, A., Ricort, G., & Aime, C. 2009, A&A, 504, 663 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Carlotti, A., Vanderbei, R., & Kasdin, N. J. 2011, Opt. Express, 19, 26796 [NASA ADS] [CrossRef] [Google Scholar]
 Carlotti, A., Kasdin, N. J., Martinache, F., et al. 2012, in SPIE Conf. Ser., 8446 [Google Scholar]
 Carlotti, A., Kasdin, N. J., Vanderbei, R. J., & Riggs, A. J. E. 2013, in SPIE Conf. Ser. [Google Scholar]
 Enya, K., Abe, L., Tanaka, S., et al. 2008, A&A, 480, 899 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Galicher, R., Baudoz, P., & Baudrand, J. 2011, A&A, 530, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Guyon, O. 2003, A&A, 404, 379 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Guyon, O., Pluzhnik, E. A., Kuchner, M. J., Collins, B., & Ridgway, S. T. 2006, ApJS, 167, 81 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Guyon, O., Martinache, F., Belikov, R., & Soummer, R. 2010, ApJS, 190, 220 [Google Scholar]
 Kasdin, N. J., Vanderbei, R. J., & Belikov, R. 2007, Comptes Rendus Physique, 8, 312 [Google Scholar]
 Kasdin, N. J., Carlotti, A., Pueyo, L., Groff, T., & Vanderbei, R. 2011, in SPIE Conf. Ser., 8151 [Google Scholar]
 Macintosh, B., Graham, J., Palmer, D., et al. 2007, Comptes Rendus Physique, 8, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Malbet, F., Yu, J. W., & Shao, M. 1995, PASP, 107, 386 [NASA ADS] [CrossRef] [Google Scholar]
 Martinache, F. & Guyon, O. 2009, in SPIE Conf. Ser., 7440 [Google Scholar]
 Martinez, P., Dorrer, C., Kasper, M., Boccaletti, A., & Dohlen, K. 2009, A&A, 500, 1281 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mawet, D., Riaud, P., Absil, O., Baudrand, J., & Surdej, J. 2005, in SPIE Conf. Ser. 5905, ed. D. R. Coulter, 502 [Google Scholar]
 Mawet, D., Mennesson, B., Serabyn, E., Stapelfeldt, K., & Absil, O. 2011a, ApJ, 738, L12 [NASA ADS] [CrossRef] [Google Scholar]
 Mawet, D., Serabyn, E., Wallace, J. K., & Pueyo, L. 2011b, Optics Letters, 36, 1506 [NASA ADS] [CrossRef] [Google Scholar]
 Mawet, D., Pueyo, L., Lawson, P., et al. 2012, in SPIE Conf. Ser. 8442 [Google Scholar]
 Mawet, D., Absil, O., Delacroix, C., et al. 2013a, A&A, 552, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mawet, D., Pueyo, L., Carlotti, A., et al. 2013b, ApJS, 209, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Murakami, N., Uemura, R., Baba, N., et al. 2008, PASP, 120, 1112 [NASA ADS] [CrossRef] [Google Scholar]
 Newman, K., Belikov, R., & Guyon, O. 2013, in Am. Astron. Soc. Meet. Abstracts, 221, 149.31 [Google Scholar]
 Pueyo, L., & Kasdin, N. J. 2007, ApJ, 666, 609 [NASA ADS] [CrossRef] [Google Scholar]
 Pueyo, L., & Norman, C. 2013, ApJ, 769, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Pueyo, L., Kay, J., Kasdin, N. J., et al. 2009, Appl. Opt., 48, 6296 [NASA ADS] [CrossRef] [Google Scholar]
 Riggs, A. J. E., Groff, T. D., Carlotti, A., et al. 2013, in SPIE Conf. Ser., 8864 [Google Scholar]
 Rouan, D., Riaud, P., Boccaletti, A., Clénet, Y., & Labeyrie, A. 2000, PASP, 112, 1479 [NASA ADS] [CrossRef] [Google Scholar]
 Serabyn, E., Mawet, D., & Burruss, R. 2010, Nature, 464, 1018 [NASA ADS] [CrossRef] [Google Scholar]
 Sivaramakrishnan, A., Soummer, R., Oppenheimer, B. R., et al. 2010, in SPIE Conf. Ser., 7735 [Google Scholar]
 Soummer, R., Pueyo, L., Ferrari, A., et al. 2009, ApJ, 695, 695 [NASA ADS] [CrossRef] [Google Scholar]
 Soummer, R., Sivaramakrishnan, A., Pueyo, L., Macintosh, B., & Oppenheimer, B. R. 2011, ApJ, 729, 144 [NASA ADS] [CrossRef] [Google Scholar]
 Spergel, D. & Kasdin, J. 2001, in BAAS, 33, 1431 [Google Scholar]
 Vanderbei, R. J. 2012, Mathematical Programming Computation, 590 [Google Scholar]
 Vanderbei, R. J., Spergel, D. N., & Kasdin, N. J. 2003, ApJ, 590, 593 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Charge 2 VC: main parameters for the combination of apodizers and Lyot stops that give the highest system throughputs for the five aperture diameters CO_{A} and the three spider thicknesses t_{S}.
Charge 4 VC: main parameters for the combination of apodizers and Lyot stops that give the highest system through puts for the five aperture diameters CO_{A} and the three spider thicknesses t_{S}.
All Figures
Fig. 1
Optical layout of an apodized vortex coronagraph. The apodizer, vortex phase mask, Lyot stop and camera are located in the successive pupil and image planes A, B, C, and D. 

In the text 
Fig. 2
Example of a design for the Keck telescope, a charge 4 vortex coronagraph, and a circular Lyot stop: a) apodizer transmission; b) intensity in the Lyot plane (log scale); and c) comparison of the Keck aperture and the apodizer. In the last figure the apodizer is represented in white, while the Keck aperture is represented in gray (wherever the apodizer transmission is zero). 

In the text 
Fig. 3
Comparison of the throughputs estimates of the ring apodized VC (dotted lines) and the shaped pupil apodized VC (solid lines) for topological charges 2 and 4. Different measures of the throughput are illustrated. Pupil plane estimates are showed in red, while image plane estimates are showed in blue. The dotted black line represents the throughputs given in Mawet et al. (2013b). 

In the text 
Fig. 4
Apodizers giving the highest system throughput for a charge 2 VC used with a 10% (top) and a 30% (bottom) central obscuration, and three spider thicknesses: 0, 0.5, and 1% (from left to right). 

In the text 
Fig. 5
Apodizers giving the highest system throughput for a charge 4 VC used with a 10% (top) and a 30% (bottom) central obscuration, and three spider thicknesses: 0, 0.5, and 1% (from left to right). 

In the text 
Fig. 6
Comparison of the relative offaxis transmission of the VC (black line) and an AVC (red line) designed for a CO_{A} = 20% obscured aperture with t_{S} = 0.5% spiders. The comparison is made for a) a charge 2 vortex and b) a charge 4 vortex. The transmissions are normalized with respect to the coronagraph throughput T_{max}. 

In the text 
Fig. 7
Transmission (in %) of an offaxis source as a function of the distance to the star (in units of λ/D), for a) and d) apertures with no spiders; b) and e) apertures with 0.5% thick spiders; and c) and f) apertures with 1% thick spiders. Figures a), b), and c) correspond to a charge 2 VC; and figures d), e), and f) correspond to a charge 4 VC. In each case, the transmission is given for the five central obscurations: 10, 15, 20, 25, and 30% (from top to bottom). The red dots indicate the IWA of each coronagraph. 

In the text 
Fig. 8
Example of a 2D intensity distribution in the Lyot plane of the apodized VC. The intensity is displayed on a log scale for an aperture with a 15% central obscuration and 1% thick spiders. 

In the text 
Fig. 9
Coronagraphic PSF (azymutal average) for the apodized vortex designed for a 15% centrally obscured aperture and 0.5% thick spiders. The intensity is displayed on a log scale. Angular distances are given in units of λ/D. 

In the text 
Fig. 10
2D a) and azymutal average b) of the PSF of the optimal apodizer designed for a 10% central obscuration, 1% thick spiders, and a charge 4 vortex. In the azymutal average, the red dotdashed line denotes the highest spatial frequency captured by the mask at λ = 1.1λ_{0}, the blue dashed line does the same λ = 0.9λ_{0}, and the dotted green line does it at λ = λ_{0}. 

In the text 
Fig. 11
Mean contrast as a function of the bandwidth (in %) for a charge 4 VC, and an aperture with 15% central obscuration and 0.5% thick spiders. The contrast is measured between the IWA and the OWA of the coronagraph, inside an annulus. 

In the text 
Fig. 12
Comparison of the theoretical and actual transmissions of a portion of a 2D shaped pupil manufactured for the Subaru telescope. Left: microscope photography of the mask. Right: corresponding desired profile. Each pixel has a 10 μm size. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.