This paper introduces an analytical method to calculate segment-level wavefront error (WFE) tolerances to enable the detection of faint extra-solar planets using segmented-aperture telescopes in space. This study provides a full treatment of the case of spatially uncorrelated segment phasing errors for segmented telescope coronagraphy, which has so far only been approached using ad-hoc Monte Carlo (MC) simulations. Instead of describing the wavefront tolerance globally for all segments, our method produces spatially dependent requirement maps. We relate the statistical mean contrast in the coronagraph dark hole to the standard deviation of the WFE of each individual segment on the primary mirror. This statistical framework for segment-level tolerancing extends the Pair-based Analytical model for Segmented Telescope Imaging from Space (PASTIS), which is based uniquely on a matrix multiplication for the optical propagation. We confirm our analytical results with MC simulations of end-to-end optical propagations through a coronagraph. Comparing our results for the Apodized Pupil Lyot Coronagraph designs for the Large Ultraviolet Optical Infrared telescope to previous studies, we show general agreement but we provide a relaxation of the requirements for a significant subset of segments in the pupil. These requirement maps are unique to any given telescope geometry and coronagraph design. The spatially uncorrelated segment tolerances we calculate are a key element of a complete error budget that will also need to include allocations for correlated segment contributions. We discuss how the PASTIS formalism can be extended to the spatially correlated case by deriving the statistical mean contrast and its variance for a non-diagonal aberration covariance matrix. The PASTIS tolerancing framework therefore brings a new capability that is necessary for the global tolerancing of future segmented space observatories. |
1.IntroductionImaging Earth-like exoplanets and searching for biomarkers is one of the key science objectives in space astronomy for the next decade. The close proximity of such planets to their host star, and a flux ratio on the order of at visible wavelengths makes this a challenging task. These two goals can be achieved using large-aperture telescopes for large light collecting areas and high angular resolution, in combination with static and dynamic starlight suppression techniques with coronagraphs and wavefront sensing and control (WFS&C) methods. The invention of the coronagraph1 synthesized with early ideas for the direct imaging of planets2 have led to several space mission concepts being developed toward this goal today. The Habitable Exoplanet Observatory3 and the Large UV Optical Infrared Surveyor4,5 (LUVOIR) are two space-based concepts recently studied by NASA as possible future flagship missions. Their primary science objective is the direct detection and spectral characterization of habitable Earth-like planets6 and the search for life; they require primary mirror diameters of 4 to 15 m. Meanwhile, the ground-based community is preparing for the era of Extremely Large Telescopes (ELTs) where 30- to 40-m class telescopes such as the Thirty Meter Telescope7 (TMT), the Giant Magellan Telescope8 (GMT), and the European Extremely Large Telescope9 (E-ELT) will be equipped with coronagraphs and extreme adaptive optics systems to search for and characterize exoplanets.10 What unites all of these observatories is that they have significantly larger primary mirrors than their respective space-based and ground-based predecessors. This poses a number of problems that need to be solved, including considerations about overall mass, cost, and plausible launch vehicles for space-based missions. The logical consequence to this is that most of these observatories will have segmented primary mirrors, much like the Keck telescope11 or the James Webb Space Telescope (JWST).12,13 This will allow for lighter-weight backplanes and foldable structures for launch purposes, or even in-situ space assembly.14,15 Telescope segmentation introduces additional diffraction effects in the focal plane,16,17,18,19 and sources for wavefront errors (WFEs) due to segment misalignments and lighter mirror structure deformation in the form of localized segment-level aberration modes. All WFEs degrade the imaging performance in a high-contrast system20,21 as they generate light residuals all over the focal plane. Such WFE will directly impact the performance of the coronagraph instrument. There are a number of coronagraph designs that were developed specifically to maximize performance on telescopes with arbitrary apertures, which includes secondary obscurations, spiders and segmentation gaps.22–27 All high-contrast instruments that aim at very high-contrast such as what is necessary to detect Earth-like planets will deploy strategies that combine static coronagraph masks in pupil and focal planes with active control of the electric field28–31 to create a zone of deep contrast in the final image plane, the dark hole (DH). To enable such wavefront control techniques, several methods for focal plane wavefront sensing have been developed32,33 to feed into a whole system of sensors and control loops that constitute the high-contrast instrument. Even after careful cophasing of the segmented aperture and implementation of WFS&C techniques that reach the required star attenuation level, there will always be some residual errors due to drifts in the system (e.g., from thermal instabilities). These changes to the mechanical structure and in the optical train will have a direct effect on the observability of a faint point source, as a sufficient signal-to-noise ratio is needed for detection within confidence limits.34 As a consequence, these high-contrast goals with segmented apertures impose severe requirements not only on static wavefront quality, but also stability requirements on the WFE and the overall mechanical structures of the telescope. There are various works that have tried to quantify these wavefront stability requirements for high-contrast imaging, both with and without segmented apertures in mind. The Nancy Grace Roman Space Telescope (formerly known as the Wide-Field Infra-Red Survey Telescope, WFIRST) is a 2.4-m monolithic space telescope with a large central obscuration and six thick, non-radial support struts35 that render high-contrast imaging particularly challenging.36 Bound to launch in 2025, it will provide technology demonstrations for stellar coronagraphy at contrast levels37 with the Roman Space Telescope Coronagraphic Instrument (CGI). Going to a segmented telescope introduces an increased number of degrees of freedom that will influence the final contrast. While there are solutions that aim to maintain a good contrast in the DH across integration times by means of continuous WFS&C,38 the problem must also be approached from an overall engineering perspective.39–41 In particular, the direct effects of segmentation on the final coronagraphic contrast42,43 are of interest in the context of high-contrast imaging, and there is an ongoing effort to characterize and quantify the requirements for such ultra-stable telescopes.44–47 Studies performing Monte Carlo (MC) end-to-end (E2E) simulations48,49 have confirmed the strict WFE requirements of a couple of tens of picometers over tens of minutes to enable the search for faint extra-solar planets, and analytical methods for the derivation of coronagraphic performance specifications have been proposed.50,51 One thing that all of these studies have in common is that they define global WFE tolerances over the entire telescope pupil, where the segments have a random contribution to the overall aberrations. In this paper, we focus on analytically defining requirements on a segment-to-segment basis instead, using the Pair-based Analytical model for Segmented Telescope Imaging from Space (PASTIS),52–54 which models the DH average contrast of a coronagraph on a segmented telescope as a function of the segment aberrations. We first introduce a new semi-analytical (SA) calculation method for the PASTIS matrix.55 Then, we show how to compute the statistical mean of the contrast using the PASTIS modes and extend the model inversion to calculate segment-level WFE requirements for a given target contrast. A full error budget that aims at maintaining a particular DH contrast will contain WFE contributions both from spatially correlated and uncorrelated segments in the telescope pupil. The impact of aberrations made of correlated segments on coronagraph contrast has been studied in various cases, e.g., low-order Zernike modes and high-frequency checkerboard-like patterns in the pupil.48,56,57 Aberrations made of spatially uncorrelated segments on the other hand have so far mostly been addressed in E2E simulations where the segments’ amplitudes had equal standard deviations.43,49,50 In this paper, we tackle the uncorrelated contribution, and establish analytically how to allocate WFE contributions to all segments individually. This addresses an essential component in the overall error budget, which had not been formally established yet. We then use the PASTIS approach to also generalize this to the correlated case. In Sec. 2, we recall the development of the analytical propagation model and how the underlying PASTIS matrix was initially built through an analytical calculation. We then generalize the matrix calculation to all coronagraphs and segmented apertures with an extension to the SA matrix calculation, which eliminates the post-calibration step that used to be performed on a perfect coronagraph model. We show that the average contrast is always a quadratic function of the aberrations and drop the requirement of having a symmetrical DH. In Sec. 3, we perform the model inversion and validate the SA matrix. Further, we show that the PASTIS modes can be used to define a statistical framework for the analysis, additionally to their deterministic relation to the DH contrast. In Sec. 4, we derive the statistical mean contrast and its variance from two separate components that describe the imaging properties of the coronagraph instrument on the one hand, and the thermo-mechanical segment statistical correlations on the other hand. We proceed with the calculation of independent segment-based WFE requirements and how to validate them in a statistical sense, and we show how to apply this to correlated segments. All simulations in Sec. 2 to Sec. 4 are done with a narrow-angle Apodized Pupil Lyot Coronagraph (APLC) on the primary pupil of LUVOIR-A (see Fig. 1 and Sec. 5) at a wavelength of 500 nm, however, these methods can be applied to any combination of coronagraph and segmented telescope. A full demonstration of the PASTIS analysis is given in Sec. 5, where we calculate these segment tolerances for the case of three different APLC designs for LUVOIR-A and highlight some consequences of this approach. In Sec. 6, we discuss our results and compare them to previously derived requirements and how they can be used in observatory error-budgets, and in Sec. 7, we report our conclusions. Fig. 1(a) LUVOIR-A design aperture with a diameter of 15 m. (b) Narrow-angle apodizer for the LUVOIR-A APLC, intended for exoplanet characterization. It uses a FPM with a radius of (with the wavelength and the telescope diameter). (c) Resulting coronagraphic image, with a DH from 3.4 to and an average normalized intensity of , which is the coronagraph floor in the absence of optical aberrations. ![]() It is noteworthy that the main metric of the PASTIS model is the spatial average raw contrast in the DH (normalized coronagraphic intensity to peak of direct image), which is what we refer to as “contrast” throughout this paper, as opposed to a spatially dependent quantity. We also want to point out how we differentiate between this spatially averaged DH intensity, the “average DH contrast” on the one side, and a statistical mean (expectation value) of this averaged contrast over many optical propagations on the other side, the statistical “mean contrast.” 2.PASTIS Model of Telescope Segment-Level Aberrations in High-Contrast CoronagraphyThe PASTIS model was initially established for a perfect coronagraph using an analytical propagation model for aberrated pairs of segments;54 the application to real coronagraphs required a second-step numerical calibration. Here, we generalize the model to any coronagraph on any segmented aperture geometry using an SA derivation of the PASTIS matrix. We also show that the validity of the PASTIS results is not limited to symmetrical DHs but extends to non-symmetrical ones as well. Independently of the way the PASTIS matrix is calculated (analytically or semi-analytically), the derivations and conclusions that we build on the PASTIS approach retain their analytical power and potential. 2.1.Matrix Formalism to Calculate the Average Dark Hole ContrastThe goal of PASTIS is to model coronagraphic images in the presence of optical aberrations on a segmented primary mirror, which can be represented for example using localized Zernike polynomials. This basis is an obvious possible choice since segment-level piston, tip/tilt, focus, and astigmatism are naturally occurring aberrations from segment misalignments, for example in three-mirror anastigmat (TMA) designs such as JWST12,58 or LUVOIR.4,49 Although beyond the scope of this paper, the PASTIS approach can also be applied directly to any other function basis, for example to represent mirror WFEs induced by thermo-mechanical effects (mounting and backplane deformations).59–61 In this paper, we simply expand the phase aberration in the segmented pupil as a sum of local (segment-level) Zernike polynomials [Eq. (9) in Ref. 54]: where is the pupil plane coordinate, is the phase from the segmented primary, and is the total number of segments, indexed by . The is the Zernike coefficient with Noll index62 up to the maximum Zernike , and is the ’th Zernike. In this paper, we limit the study to a single Zernike mode (piston; index ) as shown in Fig. 2. Hence, we drop the index in all consecutive equations, but the PASTIS methodology is applicable to any Zernike mode, combination thereof, or other types of segment-level modes.Fig. 2Piston pair aberrations on a segmented pupil (top) and the resulting image plane intensity distributions in the DH (bottom), using the narrow-angle APLC for LUVOIR-A in an E2E propagation model. The left three panels show different interference pairs with corresponding Young-like interference fringes, while the right panel shows a random distribution of local piston on all segments of the pupil and the resulting image plane intensity. All plots appear on the same scale. ![]() In high-contrast coronagraphy, the best contrast is not typically obtained for the perfect aperture without any aberration, but more commonly in the presence of a wavefront control solution using deformable mirrors.30,46 Therefore, we are studying the response of the coronagraphic system to a perturbation around that solution. Defining as the phase solution for best DH contrast and as the segmented perturbation, the phase can be divided into High-contrast coronagraphy requires exquisite wavefront quality around the DH solution, and therefore, we assume the small aberration regime for , where the electric field is well approximated as an affine function of the phase: . The phase is zero where the pupil aperture is zero, and is a complex pupil that includes the wavefront solution to produce the static DH (with both phase and amplitude contributions, and including static errors). It is noteworthy that the phase is not necessarily small.30,31 Using Fourier optics for a scalar description of the electric field and of its propagation, the coronagraph propagation can be represented by a linear operator . This is a valid assumption for Lyot-style coronagraphs, for example an APLC56,63,64 such as the one illustrated in Fig. 1 for the LUVOIR-A coronagraph design, or a vortex coronagraph.65,66 High-order vortex designs would need a special treatment for the specific low-order modes (e.g., defocus) they reject perfectly,27 but the tolerancing of such global modes is not the main purpose of PASTIS anyway. We can hence express the intensity distribution in the final image plane as with the image plane coordinate. This intensity is therefore the sum of three terms [Eq. (16) in Ref. 54]:In most cases of interest, we will be working in a symmetrical DH. It can be shown that the spatial average of the linear cross-term in Eq. (4) over a symmetrical DH is zero [Appendix A in Ref. 54]. This simplifies Eq. (4) to a quadratic function of the phase: The main metric used in this paper is the spatial average contrast over the extent of the DH, , so using , we can express the average DH intensity as Using the expression for the phase decomposition from Eq. (1) in Eq. (6), we can derive the intensity as a function of all aberrated segment pair combinations: and therefore:This double sum combines all pairs of segments where segment has an aberration amplitude of the localized phase aberration . These cross-terms from each aberrated pair of segments are very similar to Young interference fringes, and this forms the basic idea behind the PASTIS model.54 The orientation and periodicity of these fringes depend on the separation and orientation of the according aberrated pair, as shown in Fig. 2. It is important to note that the pair-wise model is not an ad-hoc idea to build the model by pairs. It derives from the fact that we expand the primary mirror phase on a discrete number of segments. Since we build a propagation model for the intensity, the “pairs” simply appear in Eq. (8) from all the cross-terms when calculating the square modulus of the electric field in Eq. (7). Equation (8) can be readily re-written as a matrix multiplication: where is the average contrast in the DH, the coronagraph floor (i.e., the average contrast in the DH at best contrast with , in the absence of phase perturbations), is the PASTIS matrix with elements , is the aberration vector of the local Zernike coefficients on all discrete segments and its transpose. The elements of the PASTIS matrix in Eq. (9) therefore directly identify asWhile this derivation is always true in the most common case of a symmetrical DH, there are coronagraph designs that produce half-sided DHs.67 We can show that the quadratic dependency of the contrast on the phase perturbations remains true in this most general case. We rewrite Eq. (4) in a similar matrix form as Eq. (9), but preserving the linear term: where is a vector that does not need to be expressed explicitly here. If we take the derivative of this equation and solve for the aberration vector that provides the minimum contrast , we can identify and , and therefore eliminate the linear term by performing a simple change of variable:This quadratic expression is similar to Eq. (9), but with a segmented mirror perturbation solution that improves contrast compared to the case without aberrations. As discussed above, we also assume a wavefront control solution with deformable mirrors to be included in the term (and hence ), with both amplitude and phase contributions. Therefore, this guarantees that the best contrast in the presence of that wavefront control solution and DH is obtained for , which in turn means that any arbitrary segment aberration vector will always degrade the contrast. Note that this does not preclude to have a non-zero static segmented correction included as part of the term . This is equally true in broadband light: When summing over wavelengths, the quadratic nature of Eq. (12) remains true, albeit with different coefficients , and . We have shown that the average DH contrast is always a quadratic function of a segmented phase perturbation , which can be discretized into a per-segment aberration amplitude vector , coefficients on a modal basis. We can calculate this average DH contrast for any aberration vector directly, using the PASTIS matrix expression [Eq. (9)]. This is particularly interesting and efficient since it does not require E2E optical simulations and only involves simple linear algebra. Furthermore, this analytical expression can be inverted to establish a segment-level WFE budget that meets a given level of contrast. This will be detailed in Sec. 4. 2.2.Semi-Analytical Calculation of the PASTIS MatrixThe PASTIS matrix can be calculated using the original analytical approach for a perfect coronagraph, then calibrated numerically for a real coronagraph and to include pupil features (e.g., support structures) [Eq. (20) in Ref. 54]. This approach was validated against an E2E model for the 36-segment ATLAST telescope pupil with an APLC [Fig. 7 in Ref. 54] to within an error of 3%. Here, we introduce another way to calculate the PASTIS matrix using an E2E simulation68 of the average DH contrast for all individually aberrated segment pairs, from which we can identify semi-analytically the matrix elements in Eq. (10). This presents the advantage of enabling a direct calculation of the matrix for any telescope geometry, any coronagraph, and any choice of segment-level aberrations (including fully numerical ones such as segment figures induced by thermo-mechanical effects). The phase for each segment pair is expressed as a Zernike aberration: We denote by the average DH contrast, on the pair of segments , that can be calculated numerically for a small wavefront aberration like and compared to the quadratic expression of Eq. (6) under the linear expansion of this phase term: The elements of the PASTIS matrix [Eq. (10)] can then be identified directly in Eq. (14) as where the diagonal terms of the PASTIS matrix are and the off-diagonal elements (This equation corrects a sign error in a previous conference proceeding.[Eq. (12) in Ref. 55]):For simplicity, we choose the same calibration aberration amplitude for both segments. Throughout our analytical development above, is in units of radians, as is a phase. Since the PASTIS matrix can be normalized to any units though, the units of the aberration amplitude can be chosen freely in the computation of Eqs. (16) and (17). The units of the PASTIS matrix are therefore in contrast per square of units of (contrast having no physical dimension), which is consistent with Eq. (9). It is noteworthy that in the presented case in Fig. 3, the units of the aberration amplitude is waves. The aberration amplitude has to be chosen such that the global pupil aberration it results in yields an average DH contrast higher than the contrast floor, but small enough to remain in the small phase aberrations linear regime. This will be discussed further in Sec. 2.3. Fig. 3(a) SA PASTIS matrix of the 120 segment LUVOIR-A design with the narrow-angle APLC. The matrix is symmetric by construction and the blue streaks correspond to negative values. The diagonal elements show which segments have more impact on the contrast than others. The outermost segments (60 to 120) have lower matrix values because of the darker apodization for these segments. (b) This is clearly visible in the superimposed image of the apodizer on the segments. ![]() The matrix is symmetric by definition since . Off-diagonal elements of the PASTIS matrix [Eq. (17)] can be negative, which is not an issue since the only constraint for the matrix is to be positive semi-definite to ensure positive eigenvalues, since they correspond to each eigenmode’s contrast. This will be discussed in detail in Sec. 3. We could potentially calculate the matrix elements [Eq. (10)] directly by calculating those complex electric field quantities. Usually though, full E2E simulators that calculate the image plane intensity are readily available and necessary for multiple other reasons. This means that choosing to calculate the PASTIS matrix through image plane intensities makes it more flexible and portable to other simulators. More importantly, working with intensities allows us to measure an empirical PASTIS matrix without the estimation errors and computational overheads of using an electric field estimator, allowing this theory to be experimentally tested. In summary, the PASTIS matrix is constructed in two steps: (1) calculate aberrated images and their corresponding DH average contrast for each pair of aberrated segments , and (2) use these contrast values to identify analytically the elements of the PASTIS matrix based on Eqs. (16) and (17). Here, the numerical calculation of these aberrated images for pairs of segments using an E2E simulator (see Appendix A) replaces the analytical expression of Young fringes between pairs of segments.54 This approach provides more accuracy, flexibility, and generality for use with any coronagraph and telescope geometry, since the analytical approach has to be calibrated using a numerical simulation anyway. 2.3.Validating the Semi-Analytical PASTIS MatrixThe SA PASTIS matrix for the narrow-angle LUVOIR APLC is calculated following Sec. 2.2 and shown in Fig. 3. The PASTIS matrix shows how some segments have a higher impact on the final contrast than others. This is visible along the diagonal, which records the contrast contribution from each individual segment alone. For example, segments 60 to 120 have a lower contrast contribution, as they correspond to the darker areas of the apodizer on the outer two rings of the aperture [see Fig. 3(b)]. This effect is also visible on the innermost ring of hexagons. We can also notice streaks of negative values in the matrix in the off-axis areas, as discussed in Sec. 2.2. We validate the SA PASTIS matrix by comparing the PASTIS contrast obtained with the matrix formalism of Eq. (9) to the contrast from an E2E simulator using the same inputs. We show the comparison in Fig. 4. The coronagraph floor for this particular APLC design in the absence of aberrations is . The PASTIS model starts to diverge from the E2E calculation at large WFE root-mean-square (RMS) where the linear approximation of the phase breaks down. It is noteworthy that the choice of of the wavelength (used in the presented example, and is on the order of 1 nm in the visible range) on a single segment yields a global pupil WFE RMS of waves, which translates into an average DH contrast just above the coronagraph floor, but keeps it in the small aberration regime. Here, the accuracy of the SA matrix approach is significantly higher than that of the fully analytical matrix because the construction of the PASTIS matrix is based on the actual E2E simulation as opposed to a post-calibrated analytical fringe model. Fig. 4Average DH contrast as a function of WFE using both an E2E simulator (dashed red) and the PASTIS matrix propagation (solid blue). In a hockey stick graph behavior, the contrast is limited by the coronagraph itself at low WFEs corresponding to the flattened-out curve to the left (at ). From about waves to waves of WFE RMS, the contrast is limited by segment phasing aberrations. In this range the estimation error of PASTIS is 0.06% compared to the reference E2E model. A calibration aberration per segment of 1/500 wave, on a 120 segment pupil, translates to a global WFE of waves when calibrating the PASTIS matrix diagonal, or with two simultaneously aberrated segments, which is in the small aberration regime of the model, just above the coronagraph floor. The curves shown are obtained as the mean of the same 20 random realizations for each RMS value, both for the E2E simulator and the PASTIS propagation. At large WFEs (close to 0.1 waves RMS) the linear approximation breaks down and the two curves no longer match perfectly. ![]() 3.Model Inversion and Statistical Mean Contrast DerivationOnce the PASTIS matrix has been calculated, Eq. (9) gives a fully analytical expression of the DH average contrast for any random segment-level aberration . This makes PASTIS particularly well suited for error budgeting analyses, compared to otherwise computation-intensive MC analyses. More interestingly, this analytical model can be inverted to determine the pupil plane aberration vector that meets a specific average contrast target , using an eigendecomposition of the PASTIS matrix. We also show that the model inversion to obtain the target contrast as a function of eigenmodes is achieved both in a deterministic and statistical sense. 3.1.Eigendecomposition of the PASTIS Model and Mode-Segment RelationshipThe PASTIS matrix is square and symmetric by construction, and therefore diagonalizable. We perform the eigendecomposition: where U is unitary, hence invertible and . The columns of are the eigenmodes of the PASTIS matrix , which can be written as column vectors of and is the total number of eigenmodes (which is equal to the total number of segments ), indexed by . is a diagonal matrix whose diagonal elements are the eigenvalues of the matrix ; it is the diagonalized PASTIS matrix . The analysis of the eigenmodes provides information about the critical modes of the system that can be used to place tolerances on segment cophasing and stability. The full set of modes of the LUVOIR-A primary with the narrow-angle APLC is shown in Fig. 5, and a selection of modes in Figs. 7, 8, and 9.Fig. 5All PASTIS modes for the LUVOIR-A narrow-angle APLC, for local piston aberrations, sorted from highest to lowest eigenvalue. The modes are unitless, showcasing the relative scaling of the segments to each other, and between all modes. They gain physical meaning when multiplied by a mode aberration amplitude in units of WFE or phase. Their respective eigenvalues and hence relative impact on final contrast is shown in Fig. 6. ![]() Fig. 6Eigenvalues, or sensitivity of contrast to mode index , for the piston PASTIS modes of the LUVOIR-A telescope with the small FPM coronagraph design, shown in Fig. 5. Note how the PASTIS matrix and modes do not depend on the target contrast, but they do on the choice of telescope geometry and coronagraph, making them the proper modes of the optical system. ![]() Fig. 7Low-impact modes with high tolerances for the narrow-angle APLC on the LUVOIR-A telescope, for local piston aberrations. These modes have little impact on the final contrast—they are similar, but not equal, to discretized Zernike modes and the coronagraph rejects them very well by design. ![]() Fig. 8Mid-impact modes with medium tolerances for the narrow-angle APLC on the LUVOIR-A telescope, for local piston aberrations. These modes have medium impact on the final contrast, relatively speaking. These modes show mostly low spatial frequency features except for high spatial frequency components in the parts of the pupil where the apodizer covers most of the segments. ![]() Fig. 9High-impact modes with low tolerances for the narrow-angle APLC on the LUVOIR-A telescope, for local piston aberrations. These modes have the highest impact on the final contrast. They consist entirely of high spatial frequency components in the parts of the pupil where the apodizer (and other pupil plane optics) are the most transmissive. ![]() The eigenvalues shown in Fig. 6 indicate how much each mode contributes to the final image contrast if applied to the pupil in their natural normalization, without any imposed weighting. This figure shows that the high-spatial frequency modes (to the left) have a much higher impact than the lower-spatial frequency modes (to the right). The PASTIS modes form an orthonormal basis set that allows us to express any arbitrary, segment-based pupil plane aberration as a linear combination of the modes with mode weighting factors : This can also be written as indicating a basis transformation between the mode-basis and the segment basis. The inverse basis transformation is thus given by . This relationship demonstrates the physical equivalence of working in the mode-basis or in the segment basis, as we can transform any expression in one space into an expression of equivalent meaning in the other space. We further explore the physical meaning of the PASTIS modes in Sec. 3.2 and 3.3.3.2.Contrast as a Function of the EigenmodesThe mode weights will depend on how much each individual mode contributes to the final contrast, and their associated eigenvalue. Inserting Eq. (20) into Eq. (9) allows us to define this relationship: and finallyThe final contrast is therefore the sum of all squared mode weights, multiplied by their respective eigenvalue. Since the modes contribute independently to the final contrast (they are orthonormal by construction), we can define a per-mode contrast as and obtain that the total contrast is the sum of all individual contrast contributions:We can then find the ’th mode weight that gives the allocated contrast contribution as Equation (25) gives the weighting factor for each PASTIS mode when it has a particular contrast contribution . We can illustrate this expression by calculating the mode weights corresponding specifically to a uniform contrast contribution of the overall target contrast over all modes, . Then we calculate them as where the is the particular set of mode weights in the case of a uniform contrast allocation across all modes. The resulting mode weights for a total contrast allocation of are shown in Fig. 10.Fig. 10PASTIS mode weights for the uniform contrast allocation across all modes. The low-index modes to the left, which correspond to high spatial frequencies, have a lower WFE tolerance than the low spatial frequency modes with high index to the right. These mode amplitudes are inversely proportional to the eigenvalues associated with each mode [Eq. (26)], and they scale the modes such that each of them contributes the same contrast to the overall target contrast. The cumulative contrast response of the modes multiplied by these weights is shown in Fig. 11. ![]() In Fig. 11, we confirm the validity of the mode weights by showing the average DH contrast from an E2E propagation of the cumulative WFE for all modes. The linearity of the plot, as well as the end value at the target contrast validates the uniform contrast allocation to each PASTIS mode from Eq. (26). Fig. 11Cumulative contrast from all PASTIS modes when allocating the total contrast uniformly across all modes (see Sec. 3.2). For instance, the measured contrast corresponding to the first 60 accumulated weighted modes is about . We multiply all modes by their respective mode amplitude and propagate them cumulatively to the image plane, both with the E2E simulator (dashed red) and the PASTIS propagation (solid blue). Without any of the modes applied, we get the contrast floor from the coronagraph , while application of all modes together yields the requested target contrast, here . Each mode is allocated an equal contrast contribution to the final contrast, which results in a linear cumulative contrast curve. Note how neither line starts at the coronagraph floor because the lowest-index mode already adds a contrast contribution on top of the baseline contrast. The corresponding PASTIS mode weights to obtain this uniform allocation of contrast per mode, is shown in Fig. 10. ![]() 3.3.Statistical Mean of the Contrast from Mode AmplitudesIn this section, we analyze the properties of the model in the statistical sense to prepare a framework for the segment-level error budget in Sec. 4. We extend the formalism from purely deterministic mode weights to random variables. We obtain the statistical mean contrast by substituting Eq. (23) into Eq. (24) and taking the mean: Assuming zero-mean normal distributions of the PASTIS modes, we can readily identify their variance as We verify this for the uniform contrast allocation per mode [Eq. (26)] in an E2E MC simulation where we draw random samples of the PASTIS mode coefficients , following zero-mean normal distributions with standard deviations , according to the uniform contrast allocation from Eq. (26): Figure 12 shows the result of such an MC simulation where we sum each randomly weighted, individual set of modes to a unique wavefront map, propagate it to the image plane with the E2E simulator and measure the spatial average contrast in the DH. We validate that the statistical mean of all these contrast values is the target contrast for which we calculated the vector of mode standard deviations in the first place, in this case . Fig. 12Validation of the uniform contrast allocation across PASTIS modes with an E2E MC simulation, drawing random sets of mode weights from a normal distribution with standard deviations [Eq. (28)], corresponding to an equal contrast allocation per mode [Eq. (26)] with weights . The overall WFE of one realization is the sum of all weighted modes for each set, and is propagated in an E2E simulation. The histogram represents 100,000 realizations of the average DH contrast for a target contrast of . ![]() In summary, PASTIS provides an analytical model to go from a set of segment aberrations to the DH average contrast in the coronagraphic image. The calculation of the PASTIS matrix eigenmodes allows to invert this model: we can set a target contrast and allocate WFE amplitudes to each of the system’s eigenmodes to reach that target contrast. Since they form an orthonormal set of modes, their individual contrast contributions add independently, each with its own sensitivity, as a fraction of the average contrast in Sec. 3.2, on top of their statistical description in Sec. 3.3. These optical modes contain the full information of the image formation system, including the apodizer, coronagraph, primary geometry, and other optical components in the system, and are distinct from and unaware of the mechanical behavior of the telescope. By this nature, we can use them to understand the fundamental limitations for high contrast imaging with a segmented aperture, which will in sequence be further constrained by thermo-mechanical properties of the telescope as we describe in Sec. 4. 4.Segment-Level Tolerance StatisticsSection 3 details how we can construct mode-based requirements that satisfy a target on the spatial average contrast in the image plane. This relationship is not only valid in a deterministic sense (relating a mode amplitude directly to its contrast contribution), but also in the statistical sense (relating the standard deviation of random modes to the overall DH contrast). The PASTIS modes form an orthonormal basis where each mode contributes independently to the DH average contrast: the contrast of a sum of weighted modes is equal to the sum of the contrast contributions for all weighted modes. The PASTIS matrix and its eigenmodes fully describe the optical propagation through the system in terms of WFE effects from a segmented aperture on the average contrast in the DH with a coronagraph. However, they do not contain any information about thermal or mechanical effects necessary to describe the final performance of a given segmented observatory. A sound framework to develop error budgets on segmented apertures therefore requires the combination of both the optical response of the telescope and coronagraph (encapsulated in the PASTIS modes), and the thermo-mechanical response of the telescope and the observatory (encapsulated in the segment aberration covariance matrix). The goal of this section is to combine the information about the imaging formation through the coronagraph on a segmented mirror with the thermo-mechanical properties of the observatory, to establish requirements to reach a target DH contrast. 4.1.Statistical Mean Contrast and Its Variance in Segmented CoronagraphyWe can calculate the statistical mean contrast of the DH spatial average directly from Eq. (9), exploiting the fact that the trace of a scalar is the scalar itself, and that : and finally where is the segment covariance matrix, containing the as-built thermo-mechanical correlations between segments. Equation (31) allows us to calculate the statistical mean of the average DH contrast directly from the knowledge of the segment covariance matrix, no matter if there is correlation between the segments or not, combining the imaging properties of the high contrast imaging system, contained in , with the thermo-mechanical behavior of the instrument contained in .Similarly, we can derive an analytical expression for the variance of the DH contrast. Assuming that follows a zero-mean Gaussian distribution, the variance for Eq. (9) takes the very simple form69 (Theorem 5.2c): These two equations provide an unambiguous closed form derivation of the mean contrast and its variance from the optical model of the imaging system (encapsulated in ), and from the thermo-mechanical properties of the telescope (captured by ). The PASTIS matrix knows nothing of the thermo-mechanical effects of the observatory and is obtained by diffractive modelling of the coronagraph, while the segment covariance matrix comes from thermal and mechanical modeling of the observatory and is completely detached from the image formation system of the telescope. The two matrices together ( and ) fully describe the statistical response of the coronagraph system to a particular WFE allocation on segments and therefore allow to establish a set of top-level requirements on segment tolerances for an observatory. The enabling aspect of Eqs. (31) and (32) for segment-level tolerancing is that the trace is invariant under a basis transformation. It follows that if either one of the two matrices, the PASTIS matrix or the thermo-mechanical covariance matrix, is expressed in its diagonal basis, the expressions for the contrast mean and variance simplify greatly, as we show in the following sections. The segment tolerancing can thus be achieved either by diagonalizing , or by doing so with . We have treated the case of diagonalizing the PASTIS matrix in Sec. 3, where we describe the analytical framework for segmented telescope tolerancing in the diagonal basis that most naturally describes the optical sensitivity of the system to the DH contrast. In the following two sections, we turn to a basis that diagonalizes the segment covariance matrix instead, permitting us to perform the tolerancing on appropriate system modes. 4.2.Uncorrelated Segment-Level RequirementsIn finding a diagonal basis for the thermo-mechanical matrix, the easiest case is when is already diagonal, which physically corresponds to independent segments on the primary mirror. In this case, the diagonal elements of , namely, the segment variances , fully describe the effect of the primary mirror segments on the DH contrast, and the statistical mean of the contrast in Eq. (31), , finds a simple expression similar to Eq. (27): Similarly to the mode-based error budget presented in Sec. 3, we now want to find a segment-based error budget to formulate the WFE limits on each segment that reach a specific statistical mean target contrast . Turning to a statistical mean contrast allows us to define a similar allocation of contrast contributions to all segments as we did statistically (and deterministically) in the PASTIS mode basis [Eq. (24)]. The most straightforward way of doing this is to allocate the target contrast equally to all segments: If we define as the standard deviation of the WFE on the ’th segment, comparably to Eq. (28): then by combining the three previous equations we obtain the per-segment WFE requirement for this particular contrast allocation:The expression in Eq. (36) lets us calculate a per-segment requirement for all individual segments in the pupil of a coronagraph instrument, given a statistical mean target contrast. The main assumption for this is that the segments are independent from each other, and that we have access to the statistical mean value of the contrast. While the mean contrast is easily measurable on images through averaging, this might not be the case for an ultra-stable facility such as LUVOIR. However, the statistical mean contrast is an important quantity to perform segment-level WFE tolerancing, especially with regards to mirror manufacturing. In a more physical sense, we know that the intensity or contrast is proportional to the variance of the WFE. Therefore, the total final contrast over the full pupil is proportional to the sum of the segment variances, which is also proportional to the sum of the contrast for all segments, conforming with Parseval’s theorem. We validate this independent segment-level error budget for three different LUVOIR coronagraphs in Sec. 5. 4.3.Case of Correlated SegmentsWhile the assumption of statistically independent segments brings insights into WFE tolerances of segmented mirrors for coronagraphic imaging, it is not general enough to encompass all possible modes for such telescopes where large-scale thermo-mechanical drifts occur (for example, backplane “flapping” mode around the folding motion of the primary mirror). Here, we discuss extensions of the PASTIS approach to the case of correlated segments. In the case of independent segments, the covariance matrix of the aberration vector is a simple diagonal matrix holding the segment variances . However, the covariance matrix is no longer diagonal for correlated segments, because of mechanical coupling for example due to large-scale backplane deformations. In this case, the statistical mean contrast and its variance remain analytically computable with Eqs. (31) and (32). The tolerancing can be done by performing an eigendecomposition on , which will diagonalize it and provide an orthonormal set of eigenmodes that describe the mechanical perturbations of the telescope system, which is also known as the Karhunen–Loève basis. By writing , we obtain the PASTIS matrix in this new basis, , through the transformation matrix , where () still describes the optical properties of the system. In this basis, Eq. (31) takes the form: Since the thermo-mechanical covariance matrix is diagonal, we identify its diagonal elements as the thermo-mechanical mode variances . Like for Eqs. (27) and (33), this simplifies Eq. (37) yet again to which allows us to make a reasonable allocation of contrast contributions across all thermo-mechanical eigenmodes. In the same way Eq. (36) calculates a per-segment variance in the basis of independent segments, we can use Eq. (38) to tolerance the per-mode variances to any given target contrast, albeit this time for individual mechanical eigenmodes. Similarly, a transformed expression can be found for the contrast variance in Eq. (32):In this most general case of correlated segments, the knowledge of the segment-level covariance matrix and its diagonal eigenbasis supersedes the simpler description in terms of segment-level variances that is only possible in the uncorrelated case (Sec. 4.2). It allows us to express the mean contrast explicitly as a function of variances of thermo-mechanical modes that can be toleranced in a similar fashion to what was done in Sec. 3. We have presented a quantitative, fully analytical method to calculate segment-level tolerances for a high-contrast instrument on a segmented aperture telescope. These follow directly from the PASTIS matrix for which we provided a new, SA way for its calculation that exploits a numerical simulator to compute the effects of the segments on the intensity in the image plane. We encode the optical and thermo-mechanical properties of the observatory separately, with the PASTIS matrix and the segment covariance matrix , which when put together allow for the analytical calculation of the expected mean contrast and its variance. These equations are invariant under a basis transformation, which permits us to find an appropriate diagonalized basis to derive individual WFE tolerances. This can either be done by diagonalizing , as we showed in Sec. 3, or by finding a diagonal basis for . A special case is given if is naturally diagonal due to independent segments on the segmented mirror; in this case, we derive a per-segment requirement map by following the analytical framework set forth in Sec. 3. In the more general case of correlations between the segments, due to thermo-mechanical properties of the telescope, we diagonalize and use the same analysis principles in the Karhunen–Loève basis of , which allows us to calculate per-mode WFE tolerances. While the deformation matrix will be acquired through thermo-mechanical modeling and can include thermal, vibrational, or gravitational perturbations, the PASTIS matrix and the PASTIS modes give insight into the purely optical properties of the observatory, and the sensitivity of the optical system to contrast. In the next section, we validate the segment-level error budget in the case of independent segments (Sec. 4.2) for three different APLC designs for LUVOIR. 5.Application to LUVOIR WFE TolerancingThe LUVOIR study4 has two point-design cases (LUVOIR-A and LUVOIR-B), 15 and 8 m in diameters, respectively, each designed as a TMA and containing a suite of scientific instruments that include coronagraphs. The LUVOIR-A CGI70 includes a suite of three numerically optimized APLC coronagraphs71 with focal plane mask (FPM) diameters that maximize the exo-Earth yield in both detection and characterization.72,73 The optical train of an APLC [Fig. 5 in Ref. 54] contains an apodizer in the pupil plane that modulates the optical beam in amplitude, a FPM occulting the on-axis point-spread-function (PSF) core, and a Lyot stop in the subsequent pupil plane that blocks the light diffracted at the FPM.56,63,64 Of the three LUVOIR APLCs, we used the smallest FPM coronagraph, or narrow-angle coronagraph, for theory validation in the previous sections. It is typically used for spectroscopic characterization in the wavelength band where molecular oxygen and water can be detected (0.76 and ). Planet detection can however be performed at shorter wavelengths (e.g., around ) where a given angular size corresponds to a larger inner working angle in diffraction resolution units (). This larger inner working angle corresponds to a larger FPM, and a larger FPM allows for apodizer designs with a higher throughput and a more robust coronagraph design, which is where the trade-off between the three designs (narrow-, medium-, and wide-angle) is made. The three LUVOIR-A APLC designs are shown in the top row of Fig. 13. The corresponding FPM have radii of 3.50, 6.82, and 13.38 , respectively, followed by a hard edge annular Lyot stop, whose inner and outer diameters are 12.0% and 98.2% of the circumscribed diameter of the apodizers. The resulting coronagraphic DH sizes are 3.4–12, 6.7–23.7, and 13.3–46.9 respectively, with a coronagraph floor of , , and for the three designs. Fig. 13Top: The three baseline apodizer designs for LUVOIR-A, a narrow-angle (left), medium-angle (middle) and wide-angle (right) mask (details see Sec. 5). Bottom: Segment tolerance maps for narrow-angle (left), medium-angle (middle) and wide-angle (right) APLC designs on LUVOIR-A for a target contrast of , at a wavelength of 500 nm. All three tolerance maps are shown on the same scale. It is noteworthy how each segment value denotes the standard deviation of a zero-mean normal distribution from which the segment aberrations in WFE RMS are drawn. The minimum and maximum values of these maps are, from left to right: 7 and 116 pm, 25 and 93 pm, and 92 and 181 pm. ![]() In this section, we present a full analysis to obtain segment requirements for these three LUVOIR APLCs, and validate the results by performing MC simulations with an E2E simulator. We also take a deeper look into the narrow-angle APLC by analyzing the PASTIS mode-based decomposition of the individual-segment requirements. This monochromatic analysis was performed at a wavelength of 500 nm, which is the lower limit wavelength for the LUVOIR coronagraphs and where we expect to detect planets. 5.1.Segment Requirements and Monte Carlo Simulations for Three APLC DesignsWe first calculate the PASTIS matrix for each of these three APLC designs, according to the methodology described in Sec. 2.2. We can then establish a segment-level error budget in the assumption of uncorrelated segments, according to Eq. (36). In Fig. 13, we show the resulting segment requirement maps for a target contrast of and for all three APLC designs. It is important to note that these requirement maps do not represent the WFE over the segmented pupil, but instead show the standard deviations on the tolerable WFE for each segment to retrieve, statistically, the desired mean target contrast. In this sense, the maps in Fig. 13 are a prescription for the drawing of random segment WFE realizations such as the examples shown in Fig. 14. These random maps are then propagated with the E2E simulator and their average contrast values build the MC histograms in Fig. 15. One big takeaway point from Fig. 13 is that the segment requirements are not uniform across the pupil, but clearly follow the apodization of the coronagraph mask. The PASTIS matrix holds knowledge of the optical effect of not only the segments but also the coronagraph instrument on the final contrast, so by including that knowledge into the derivation of the segment constraints we obtain a requirement map optimized for that particular instrument. Moreover, we can observe a direct trade-off between the coronagraph apodization and the per-segment requirements—the more aggressive the apodization and the lower the throughput, the more we can relax the requirements on the more concealed segments within one coronagraph. However, more aggressive apodization usually comes with smaller FPM coronagraphs that filter low-order modes less, which will lead to more stringent overall requirements. This leads to a direct trade-off between FPM size, throughput and WFE requirements (see also Sec. 6). Fig. 14Four random segment-based WFE maps drawn from a zero-mean normal distribution and the per-segment standard deviations from the left prescription map in Fig. 13, for the narrow-angle APLC design and a target contrast of . After each random map is created, we propagate it through the E2E simulator and record its average contrast to build the left MC simulation in Fig. 15. ![]() Fig. 15Validation of the independent segment tolerancing with E2E MC simulations, for different target contrasts, using the narrow-angle APLC design. Each segment in one of the 100,000 WFE realizations is drawn from a zero-mean normal distribution with standard deviation . The dashed-dotted lines mark the target contrast of each case, which are successfully recovered by the mean values of the histograms, in accordance with their analytical calculation in Eq. (31). The dotted lines mark the 1-sigma confidence limits of this contrast distribution, which are , and for the three target contrasts , and respectively, and they accord with the numbers calculated by Eq. (32). ![]() These requirement maps can be calculated for any target contrast in the range of validity of the PASTIS model, which we discussed in Sec. 2.3. We can verify Eq. (36) by running MC simulations with the E2E simulator across a grid of different coronagraph instruments and target contrasts. Using a range of target contrasts , , and on the narrow-angle baseline LUVOIR APLC design, we first calculate the segment constraints (the requirement maps for and are not shown here, but they show the same spatial distribution over the segments as in Fig. 13, only different by a proportionality factor). We draw the WFE amplitude for each individual segment from a zero-mean normal distribution and its standard deviation : We use these random aberration amplitudes on all segments to compose a WFE map on the segmented pupil and then propagate this WFE map through the E2E simulator to measure the resulting spatial average contrast in the DH. Doing this 100,000 times for each target contrast case, we obtain the histograms shown in Fig. 15. The mean of the resulting MC simulations clearly recovers the target contrast for which the segment requirements have been calculated, which is indicated by the dashed-dotted line. Both these mean values, and the standard deviations, indicated with the dotted lines in Fig. 15, agree with the theoretical values calculated analytically from the segment covariance matrix [Eqs. (31) and (32)]. Also, we have verified the correct recovery of the same range of target contrasts by means of MC simulations for the other two APLC designs shown in Fig. 13 (resulting histograms not shown in this paper). 5.2.Modal Analysis of the Segment-Based RequirementsThe segment requirement maps were obtained assuming a uniform contrast allocation across all segments [Eq. (34)]. We also assumed statistically independent segments, so that their correlation matrix was diagonal. Here, we further explore this uniform error budget in the segment basis by analyzing the corresponding distribution of proper system modes of the optical system, the PASTIS modes. Using the transformation matrix from the eigendecomposition of the PASTIS matrix , we can calculate the corresponding covariance matrix in the PASTIS mode basis with . Given this linear transformation, if the covariance matrix is diagonal in one space, we do not expect it to be diagonal in the other space. The matrix , obtained from the diagonal segment covariance matrix assembled from the requirement map, is shown in Fig. 16(a). Figure 16(b) compares the extracted standard deviations for PASTIS modes along the diagonal of , with the PASTIS mode weights previously calculated in Sec. 3. Although is not diagonal, this is a legitimate comparison. The PASTIS matrix is always diagonal in its own eigenbasis, expressed as matrix in Sec. 3.1. This is why the average contrast expression from the statistical mode weights in Eq. (27) only requires the diagonal elements of the covariance matrix , no matter whether it is diagonal or not, i.e., whether the mode weights show some correlations or not. The difference with respect to the mode weights [Eq. (26)] obtained under the assumption of a uniform contrast allocation per mode (Fig. 10) is notable: the mode weights of low mode index have increased tolerances, which is very interesting from a system design point of view, while large index modes (above ) are strongly attenuated in the PASTIS mode basis error budget for independent segments. This is also clearly visible in Fig. 17 where the contrast contribution per mode is relatively flat at a low mode index but drops to negligible contributions at high-index modes. For the case of the flat contrast allocation across modes, this figure shows a flat line at for comparison (dashed line). This effect can be well understood by looking back at Fig. 5, where high-index modes appear to be very similar to low-order Zernike modes, therefore, having highly correlated segments, and the low-index modes appear as high-spatial frequencies, i.e., with more uncorrelated segments. Therefore, the construction of a segment-level error budget for uncorrelated segments creates a modal distribution with extremely low weights on the PASTIS modes that have highly correlated segments (high-index modes), as seen in Fig. 16. Also, since the mode contrast contribution is directly related to the mode weight [Eq. (25)], the same effect is visible in the allocated contrast per mode (Fig. 17). Fig. 16(a) Covariance matrix , calculated from the diagonal covariance matrix in segment-space, , with . Although there are clearly some correlations present between the high-index PASTIS modes in the top right corner (low spatial frequencies), this does not matter as long as we are in the PASTIS segment basis, where the PASTIS matrix is diagonal. When this is the case, the mean contrast only depends on the diagonal elements of [Eq. (27)]. (b) PASTIS mode amplitudes for the case of independent segments in WFE RMS (solid red). They are extracted from the mode covariance matrix , after constructing an error budget assuming independent segments that contribute equally to the total contrast, at 500 nm. Overlapping (dashed gray), we can see the mode weights from the uniform contrast allocation to all PASTIS modes from Fig. 10. We can clearly see how compared to that flat allocation, the independent-segment error budget increases the tolerances of low-index modes (left) that have less segment correlation, and dampens the tolerances of high-index modes (right) that are highly correlated, low-spatial frequency modes. ![]() Fig. 17Contrast per individual PASTIS mode when derived from the error budget in which all segments contribute independently and equally to the final contrast (solid blue). High-index modes to the right, which correspond to low-spatial frequencies and therefore highly correlated segments (see Fig. 5), are highly attenuated and contribute negligible amounts to the contrast. The uniform contrast allocation across all modes at is indicated with the dashed gray line (the contrast floor has been removed in both curves). ![]() To illustrate this further, we calculate a cumulative contrast plot similar to Fig. 11, which was initially obtained for a uniform contrast allocation per mode. The new result is shown in Fig. 18, where we can see that it is no longer linear, i.e., the modes no longer contribute equally to the total mean contrast. The slope of the blue curve is indicative of the allocated tolerances for each mode: the low-spatial frequency PASTIS modes on the righthand side now contribute significantly less to the final contrast, while the first modes contribute more, while still resulting in the exact cumulative target contrast. This is consistent with the behavior discussed in Figs. 16 and 17. Fig. 18Cumulative contrast of the PASTIS modes for the error budget in which all segments contribute independently and equally to the final contrast (solid blue), compared to the case in Fig. 11 where all modes contribute the same contrast (dashed gray). The high-index modes have negligible contrast impact (see also Fig. 17) as they correspond to low-spatial frequency, highly correlated segments. This plot also confirms the assumption that the mode covariance matrix is nearly diagonal. ![]() 6.DiscussionThe results we obtain for the three LUVOIR APLC designs in Sec. 5, under the assumption of statistically independent segments, span more than one order of magnitude from the most constrained segment on the small APLC to the most relaxed segment on the large APLC (7 to 181 pm). These results show not only a dependency on the coronagraph design but also a wide range of segment requirements within the pupil for one single coronagraph. The segment with the most stringent requirement can tolerate a standard deviation of 7 pm local piston error on the narrow-angle APLC for a target contrast of , at a wavelength of 500 nm, which is comparable to previous results for segment-based piston errors on LUVOIR-A (10 pm at a wavelength of 575 nm49), while studies on other apertures and with other segment numbers quote similar numbers.43,50 However, we show that segments in other parts of the pupil have a much higher local piston tolerance, a standard deviation up to 116 pm on the small APLC design, which suggests that not all segments need to be held to the same tolerance level. Instead, we can relax the segment-level requirements on those segments that do not influence the average DH contrast as much, while still obtaining the same statistical mean contrast. Local relaxation of the WFE limits on certain parts of the pupil can be exploited for example for the backplane mechanical design and observatory-level control strategy. In particular, the tolerances will also depend on the total number of segments in the pupil. While this has not been studied systematically in this paper, PASTIS can enable such work. For example, with fewer segments in the aperture, the spatial frequencies corresponding to segment misalignments will be lower. Therefore, WFE from these misalignments will be more filtered by the coronagraph, which might lead to increased tolerances. Inversely, with more segments in the pupil, the highest spatial frequencies from segment misalignments will not be filtered by the FPM and might thus result in lower tolerances, as is already the case for the most sensitive PASTIS modes. Wavefront aberrations from a mirror with fewer segments will be typically more filtered by the coronagraph. However, more and smaller segments will introduce high spatial frequencies that will diffract light into the image beyond the outer working angle, decreasing the impact of each segment misalignment on the DH. Overall, the number of segments in the pupil will raise competing effects that will influence the overall tolerances and will be highly influenced by the type of coronagraph.74 In Sec. 5.1, we briefly mentioned the observed trade-off between coronagraph throughput, FPM size, and per-segment requirements. For a given coronagraph design, we observe relaxed requirements for those segments in the pupil that are more concealed by the apodizer (i.e., more black in the apodizer image). The overall tolerance over the entire pupil also increases with the size of the FPM. The larger the FPM, the higher the rejection in particular of low-order spatial modes, which correspond to high mode indices where the modes are similar to low-order, Zernike-like global modes. This results in higher mode weights, which becomes obvious in Fig. 10, where these high-index modes on the right side show higher WFE tolerances. When moving to the independent-segment error budget in Fig. 16(b), this effect becomes less obvious as the low-order mode tolerances get dampened due to our assumption of uncorrelated segments contributing equally to the mean contrast. The projection of these uncorrelated segments onto the mode basis favors high-spatial frequency modes. This leaves the low-order modes statistically weak, as they would otherwise contribute to inter-segment correlation. However, this is only true for these fully uncorrelated segment-level WFE contributions. This uncorrelated error budget is over-constraining the low-order (high-index) modes (see Fig. 16(b)), therefore not taking advantage of the coronagraphic rejection of these spatially correlated modes. A complete error budget will need to allocate contrast contributions separately between the correlated and the uncorrelated components of the WFE. The final result will therefore have a modal weight distribution in-between the solid (fully uncorrelated) and dashed lines (uniform contrast across modes) in Fig. 16(b). Additionally, large-FPM APLCs have higher throughput apodizers. This means that their larger coronagraphic rejection (associated with the larger FPM) contributes more to the WFE tolerance relaxation than the apodizer throughput itself. Therefore, a true optimization of the WFE tolerances will be a trade-off between the FPM size and the fraction of apodization in the pupil. Further, this tolerancing work introduces new design considerations for high-contrast instruments, which is the optimization of the coronagraphic component with respect to segment phasing tolerances. Such an optimization will aim to release the segment tolerances while keeping a reasonable contrast goal, with the ultimate goal to maximize exoplanet yield, which should be explored in future work. In a realistic telescope of course, the segments are typically correlated due to the deformations of the backplane structure. The development of technologies that support increased wavefront stability of segmented telescopes is actively being worked on today.45 This includes precise methods for thermo-mechanical modeling and measuring of such deformation effects on the segmented primary mirror. The underlying segment correlations can be provided either as a segment covariance matrix, or in the form of thermo-mechanical system eigenmodes, in which case we can directly use Eq. (37) for a tolerancing analysis, after expressing the PASTIS matrix in this new basis. While local piston errors have been shown to have the largest impact on the contrast,49,50 we can generate a PASTIS matrix for other local Zernike modes as well (e.g., tip/tilt, focus, astigmatism). The feasibility of this has already been shown in the analytical approach54 and should hence be regarded as a mere functional addition. Moreover, if we have knowledge of telescope-design and hardware dependent local aberration modes (e.g., from effects such as adhesive shrinkage, bulk temperature, coatings, or gravity), these can be used as well to build the PASTIS matrix and derive their corresponding segment-level requirements. Instead of evaluating the tolerance levels mode by mode, we can also calculate a multi-mode PASTIS matrix that incorporates combinations of local modes (e.g., piston–tip–tilt, or combinations of custom modes) to derive segment tolerances that will take into account that more than one distinct local aberration mode is contributing to the overall WFE. Another application of extended PASTIS matrices is the generalization to high-spatial frequency effects (e.g. from polishing). Instead of building a PASTIS matrix with pair-wise Zernike aberrations of segments, we can use sine waves locally on the segments in lieu of Zernikes. Each spatial frequency and orientation would then be a new local mode, and we can use many of them to make a multi-mode PASTIS matrix in the same way as with any other local modes. While a generalization to a continuous distribution of frequencies to build tolerances in terms of a power spectral density might be possible, it is beyond the scope of the present discussion. Nevertheless, the generalization to a few sine-wave frequencies (e.g., with their corresponding speckles localized at the inner and outer working angle, or in the middle of the DH) is a direct and straightforward extension of the present illustration and would provide meaningful input for tolerancing purposes of polishing errors. The presented tolerancing model provides WFE limits on the segments, but it does not define how these limits are to be maintained. Relying purely on the mechanical stability of the telescope will not be enough to stay within these requirements and therefore an active optics system will be needed to measure WFE deviations and compensate for the residuals. Such an active optics system will include WFS&C and signal-to-noise considerations, as the wavefront sensor needs enough photons to provide an accurate wavefront estimate and correction.75 Finally, the application of PASTIS to ground-based observatories is possible but will have to take additional effects in account. While future large segmented telescopes (TMT, E-ELT, and GMT) will reach contrast levels of to within the next decade,76 which is sufficient to enter the high-contrast regime that the PASTIS model can be applied to, future work will have to include the effects of residual turbulence to truthfully represent coronagraphic observations on those observatories. As a first common approximation for this purpose, the coronagraphic PSF can be expressed as the sum of a static and a dynamic contribution,77 so the PASTIS analysis can be used for a characterization of the static part. 7.ConclusionsThe goal of PASTIS, as established by Leboulleux et al.,54 is an analytical propagation model to calculate the average DH contrast in a coronagraphic system, in the presence of segment-level aberrations. This is achieved with a closed-form expression in Eq. (9) that depends exclusively on the PASTIS matrix , acting on the aberration amplitudes on all segments. In this paper, we extended the calculation of the matrix to an SA approach, where the optical propagation of segment aberrations is performed numerically before assembling the matrix analytically [Eqs. (16) and (17)]. This makes the model more accurate as it includes all details of the optical system as provided with the E2E simulator. We also show that the model holds even for a non-symmetrical DH. The SA PASTIS approach is therefore a flexible tolerancing tool that can be adapted readily to any telescope geometry or coronagraph, as shown in Sec. 5. It can be used to study trade-offs between coronagraph designs that will provide certain tolerance distributions over the segments on the primary mirror, and telescope-level engineering constraints implemented in other parts of the observatory. We used the model to derive analytical expressions for the statistical mean and variance of the average DH contrast [Eqs. (31) and (32)]. This opens the possibility for WFE tolerancing of a segmented observatory. In addition to the optical properties modeled by the PASTIS matrix , these expressions involve the segment-level covariance matrix that describes the thermo-mechanical properties of the telescope. Indeed, deformations of the backplane structure typically lead to correlated segment poses (e.g. backplane flapping modes around the folding motion of the primary mirror). It is this combination of the optical with the thermo-mechanical characteristics that lays the foundation to a complete and analytical method for the tolerancing of segmented aberrations. The key to calculate WFE requirements with this framework is to find the diagonal basis of either of the two matrices, or . The PASTIS matrix can easily be diagonalized by means of an eigendecomposition, in which case the tolerancing can be performed on the PASTIS eigenmodes. They form an orthonormal set of modes, representing the proper optical system modes of the given observatory and coronagraph instrument, and allow to analyze its fundamental limitations in terms of WFE propagation in segmented aperture coronagraphy. These eigenmodes contribute additively to the image plane average contrast according to their mode-level tolerances, which correspond to the standard deviations per mode , associated with the statistical mean of the average DH contrast. Additionally, if we have information on the thermo-mechanical behavior of the instrument or of the whole observatory embodied in the segment aberration covariance matrix (e.g., through finite-element simulations), we can choose to work in the eigenbasis of these thermo-mechanical perturbations, a.k.a. the Karhunen–Loève basis, to perform the tolerancing. This allows us to put requirements on structural deformations that impact the segmented primary mirror (e.g., backplane and mirror support structures) and can be modeled as rigid-body motions at the segment level. A better approximation of these thermo-mechanical modes will be obtained with a multi-mode PASTIS matrix by combining multiple local Zernikes and/or ad-hoc local aberration modes, as described in Sec. 6. In the simplified case of independent segments, the covariance matrix contains only diagonal elements, representing the individual segment WFE variances . We have built a segment-based error budget by allocating equal contribution to contrast from all segments. This allowed us to calculate WFE requirements for all segments individually, building segment requirement maps as we showed in Sec. 5 for the three APLC designs of the LUVOIR-A telescope. The advantage of this method is two-fold: First, rather than calculating WFE tolerances globally over the entire pupil, we can obtain a WFE standard deviation per segment, which can locally lead to a relaxation in requirements. Second, we do not need to perform full MC simulations that evaluate different realizations of WFE maps for that purpose, instead, we can calculate these requirements analytically in one single step. The analysis presented in this paper is statistical but static, i.e., without temporal evolution; the extension to dynamical drift rates depends on the observing scenario and wavefront control strategy, which will put the PASTIS propagation model on different time scales.45,46 Future work will address such dynamic analysis methods for continuous WFS&C cases on ultra-stable telescopes. 8.Appendix A. LUVOIR-A End-to-End SimulatorThe LUVOIR-A E2E simulator used in this paper is fully written in Python and uses the HCIPy78 package for optical propagations and the Lyot coronagraph implementation. The simulator inherits from the HCIPy SegmentedMirror class that provides the capability to individually actuate the segments in a telescope pupil. The core of the simulator is an abstract class that combines this segmented actuation with an APLC, which gets inherited by the the main LUVOIR-A simulator that defines the specific aperture geometry, apodizers, FPMs and Lyot stops of the LUVOIR-A APLCs. The source files (aperture, apodizers, and Lyot stop) of this APLC suite are credited to the Segmented-aperture Coronagraph Design and Analysis study, and the user can switch between the three designs with a single keyword parameter. The simulator calculates the electric field at each pupil and focal plane of the optical system, including the final image plane of the coronagraph. These electric fields are then returned by the simulator. Correct optical propagation of various coronagraphs, including the APLC, has been implemented by various other simulators, as have been segmented mirrors. The problem that this particular simulator solves and makes available is that of proper sampling of the pupil for the segmented actuation with respect to the pupil plane coronagraph optics: the apodizer and the Lyot stop. It is very important that the pixels of the individual segments overlap perfectly with the apodizer, otherwise the coronagraph will not perform well. Both the aperture as well as the apodizer design and Lyot stop are input parameters to the simulator, hence the correct pixel-to-pixel mapping between the aperture and the apodizer are guaranteed. Currently, the segmented mirror control includes local piston, tip and tilt, but can be easily extended to other modes and modal bases by implementing newer versions of HCIPy, which is ongoing refactoring. The simulator version used for this paper is published within the PASTIS Python package68 and available on GitHub. AcknowledgmentsThis work was co-authored by employees of BALL AEROSPACE as part of the the Ultra-Stable Telescope Research and Analysis (ULTRA) Program under Contract No. 80MSFC20C0018 with the National Aeronautics and Space Administration (PI: L. Coyle), and by STScI employees under corresponding subcontracts No.18KMB00077 and No. 19KMB00102 with Ball Aerospace (PI: R. Soummer, Sci-PI: L. Pueyo). This work was also co-authored by employees of the French National Aerospace Research Center ONERA (Office National d’Études et de Recherches Aérospatiales), and benefited from the support of the WOLF project ANR-18-CE31-0018 of the French National Research Agency (ANR), as well as the internal research project VASCO. This work was also supported in part by the National Aeronautics and Space Administration under Grant 80NSSC19K0120 issued through the Strategic Astrophysics Technology/Technology Demonstration for Exoplanet Missions Program (SAT-TDEM; PI: R. Soummer), and by the Segmented-aperture Coronagraph Design and Analysis funded by ExEP, under JPL subcontract No.1539872, and by the STScI Director’s Discretionary Research Funds. The authors acknowledge the contributions of Kathryn St. Laurent, James Noss and Emiel Por to the design and optimization of the LUVOIR-A apodizers. I. Laginja is thankful to Théo Jolivet and Anand Sivaramakrishnan for helpful discussions. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to reproduce, prepare derivative works, distribute copies to the public, and perform publicly and display publicly, or allow others to do so, for United States Government purposes. All other rights are reserved by the copyright owner. This research was developed in Python,79 an open-source programming language, and made use of the Numpy,80,81 Matplotlib,82,83 Astropy84–86 and HCIPy78 packages. The software for all analysis was published in PASTIS,68 a modular, open-source Python package for segment-level error budgeting of segmented telescopes. ReferencesB. Lyot,
“The study of the solar corona and prominences without eclipses (George Darwin Lecture, 1939),”
MNRAS, 99 580
(1939). https://doi.org/10.1093/mnras/99.8.580 Google Scholar
N. G. Roman,
“Planets of other suns,”
Astron. J., 64 344
–345
(1959). https://doi.org/10.1086/108038 ANJOAA 0004-6256 Google Scholar
B. S. Gaudi et al.,
“The Habitable Exoplanet Observatory (HabEx) mission concept study final report,”
(2020). Google Scholar
The LUVOIR Team,
“The LUVOIR mission concept study final report,”
(2019). Google Scholar
M. R. Bolcar,
“The large UV/optical/infrared (LUVOIR) surveyor: engineering design and technology overview,”
Proc. SPIE, 11115 111150P
(2019). https://doi.org/10.1117/12.2530456 PSISDG 0277-786X Google Scholar
A. Roberge, M. R. Bolcar and K. C. France,
“Telling the story of life in the cosmos: the LUVOIR telescope concepts,”
Proc. SPIE, 11115 111150O
(2019). https://doi.org/10.1117/12.2530475 PSISDG 0277-786X Google Scholar
L. Simard et al.,
“Thirty meter telescope science instruments: a status report,”
Proc. SPIE, 9908 99081V
(2016). https://doi.org/10.1117/12.2234058 PSISDG 0277-786X Google Scholar
J. Fanson et al.,
“Overview and status of the Giant Magellan Telescope project,”
Proc. SPIE, 10700 1070012
(2018). https://doi.org/10.1117/12.2561852 PSISDG 0277-786X Google Scholar
S. Ramsay et al.,
“The ESO Extremely Large Telescope instrumentation programme,”
Proc. SPIE, 11203 1120303
(2020). https://doi.org/10.1117/12.2541400 PSISDG 0277-786X Google Scholar
O. Guyon et al.,
“How ELTs will acquire the first spectra of rocky habitable planets,”
Proc. SPIE, 8447 84471X
(2012). https://doi.org/10.1117/12.927181 Google Scholar
T. S. Mast, J. E. Nelson and W. J. Welch,
“Effects of primary mirror segmentation on telescope image quality,”
Proc. SPIE, 0332 123
–133
(1982). https://doi.org/10.1117/12.933512 PSISDG 0277-786X Google Scholar
D. S. Acton et al.,
“Wavefront sensing and controls for the James Webb Space Telescope,”
Proc. SPIE, 8442 84422H
(2012). https://doi.org/10.1117/12.925015 PSISDG 0277-786X Google Scholar
M. D. Perrin et al.,
“Updated optical modeling of JWST coronagraph performance contrast, stability, and strategies,”
Proc. SPIE, 10698 1069809
(2018). https://doi.org/10.1117/12.2313552 PSISDG 0277-786X Google Scholar
R. S. Polidan et al.,
“Servicing and assembly: enabling the most ambitious future space observatories,”
Proc. SPIE, 10698 1069825
(2018). https://doi.org/10.1117/12.2313013 PSISDG 0277-786X Google Scholar
L. M. Bowman et al.,
“In-space assembly application and technology for NASA’s future science observatory and platform missions,”
Proc. SPIE, 10698 1069826
(2018). https://doi.org/10.1117/12.2311413 PSISDG 0277-786X Google Scholar
P. A. Lightsey and M. Chrisp,
“Image quality for large segmented space telescopes,”
Proc. SPIE, 4850 453
–460
(2003). https://doi.org/10.1117/12.461220 PSISDG 0277-786X Google Scholar
N. Yaitskova, K. Dohlen and P. Dierickx,
“Analytical study of diffraction effects in extremely large segmented telescopes,”
J. Opt. Soc. Am. A, 20 1563
–1575
(2003). https://doi.org/10.1364/JOSAA.20.001563 JOAOD6 0740-3232 Google Scholar
M. Troy and G. A. Chanan,
“Diffraction effects from giant segmented mirror telescopes,”
Proc. SPIE, 4840 81
–92
(2003). https://doi.org/10.1117/12.457985 PSISDG 0277-786X Google Scholar
S. Itoh et al.,
“Point spread function of hexagonally segmented telescopes by new symmetrical formulation,”
MNRAS, 483 119
–131
(2019). https://doi.org/10.1093/mnras/sty3052 Google Scholar
I. J. Crossfield and M. Troy,
“Segment aberration effects on contrast,”
Appl. Opt., 46 4533
–4540
(2007). https://doi.org/10.1364/AO.46.004533 APOPAI 0003-6935 Google Scholar
N. Yaitskova and M. Troy,
“Rolled edges and phasing of segmented telescopes,”
Appl. Opt., 50 542
(2011). https://doi.org/10.1364/AO.50.000542 APOPAI 0003-6935 Google Scholar
A. Sivaramakrishnan and N. Yaitskova,
“Lyot coronagraphy on giant segmented-mirror telescopes,”
Astrophys. J. Lett., 626 L65
–L68
(2005). https://doi.org/10.1086/431460 AJLEEY 0004-637X Google Scholar
P. Martinez et al.,
“Comparison of coronagraphs for high-contrast imaging in the context of extremely large telescopes,”
Astron. Astrophys., 492 289
–300
(2008). https://doi.org/10.1051/0004-6361:200810650 AAEJAF 0004-6361 Google Scholar
R. Soummer et al.,
“Apodized Pupil Lyot coronagraphs for arbitrary apertures. II. Theoretical properties and application to extremely large telescopes,”
Astrophys. J., 695 695
–706
(2009). https://doi.org/10.1088/0004-637X/695/1/695 ASJOAB 0004-637X Google Scholar
O. Guyon et al.,
“High performance Lyot and PIAA coronagraphy for arbitrarily shaped telescope apertures,”
Astrophys. J., 780 171
(2014). https://doi.org/10.1088/0004-637X/780/2/171 ASJOAB 0004-637X Google Scholar
N. T. Zimmerman et al.,
“Lyot coronagraph design study for large, segmented space telescope apertures,”
Proc. SPIE, 9904 99041Y
(2016). https://doi.org/10.1117/12.2233205 PSISDG 0277-786X Google Scholar
G. Ruane et al.,
“Performance and sensitivity of vortex coronagraphs on segmented space telescopes,”
Proc. SPIE, 10400 104000J
(2017). https://doi.org/10.1117/12.2274508 PSISDG 0277-786X Google Scholar
P. J. Bordé and W. A. Traub,
“High-contrast imaging from space: speckle nulling in a low-aberration regime,”
Astrophys. J., 638 488
–498
(2006). https://doi.org/10.1086/498669 ASJOAB 0004-637X Google Scholar
L. Pueyo et al.,
“Optimal dark hole generation via two deformable mirrors with stroke minimization,”
Appl. Opt., 48 6296
(2009). https://doi.org/10.1364/AO.48.006296 APOPAI 0003-6935 Google Scholar
J. Mazoyer et al.,
“Active correction of aperture discontinuities-optimized stroke minimization. I. A new adaptive interaction matrix algorithm,”
Astron. J., 155 7
(2018). https://doi.org/10.3847/1538-3881/aa91cf ANJOAA 0004-6256 Google Scholar
J. Mazoyer et al.,
“Active correction of aperture discontinuities-optimized stroke minimization. II. Optimization for future missions,”
Astron. J., 155 8
(2018). https://doi.org/10.3847/1538-3881/aa91d7 ANJOAA 0004-6256 Google Scholar
T. D. Groff et al.,
“Methods and limitations of focal plane sensing, estimation, and control in high-contrast imaging,”
J. Astron. Telesc. Instrum. Syst., 2 011009
(2016). https://doi.org/10.1117/1.JATIS.2.1.011009 Google Scholar
N. Jovanovic et al.,
“Review of high-contrast imaging systems for current and future ground-based and space-based telescopes: Part II. Common path wavefront sensing/control and coherent differential imaging,”
Proc. SPIE, 10703 107031U
(2018). https://doi.org/10.1117/12.2314260 PSISDG 0277-786X Google Scholar
R. G. Lyon and M. Clampin,
“Space telescope sensitivity and controls for exoplanet imaging,”
Opt. Eng., 51 011002
(2012). https://doi.org/10.1117/1.OE.51.1.011002 Google Scholar
J. Krist et al.,
“An overview of WFIRST/AFTA coronagraph optical modeling,”
Proc. SPIE, 9605 960505
(2015). https://doi.org/10.1117/12.2188361 PSISDG 0277-786X Google Scholar
B. Nemati, J. E. Krist and B. Mennesson,
“Sensitivity of the WFIRST coronagraph performance to key instrument parameters,”
Proc. SPIE, 10400 1040007
(2017). https://doi.org/10.1117/12.2274396 PSISDG 0277-786X Google Scholar
D. Savransky and D. Garrett,
“WFIRST-AFTA coronagraph science yield modeling with EXOSIMS,”
J. Astron. Telesc. Instrum. Syst., 2 011006
(2016). https://doi.org/10.1117/1.JATIS.2.1.011006 Google Scholar
L. Pogorelyuk and N. J. Kasdin,
“Dark hole maintenance and a posteriori intensity estimation in the presence of speckle drift in a high-contrast space coronagraph,”
Astrophys. J., 873 95
(2019). https://doi.org/10.3847/1538-4357/ab0461 ASJOAB 0004-637X Google Scholar
L. Feinberg et al.,
“Ultra-stable segmented telescope sensing and control architecture,”
Proc. SPIE, 10398 103980E
(2017). https://doi.org/10.1117/12.2272810 PSISDG 0277-786X Google Scholar
H. P. Stahl,
“Advanced mirror technology development (AMTD): year five status,”
Proc. SPIE, 10401 104010O
(2017). https://doi.org/10.1117/12.2275190 PSISDG 0277-786X Google Scholar
M. East et al.,
“Picometer level stability of a mounted mirror assembly,”
Proc. SPIE, 10698 106980Y
(2018). https://doi.org/10.1117/12.2314332 PSISDG 0277-786X Google Scholar
N. Yaitskova and K. Dohlen,
“Tip–tilt error for extremely large segmented telescopes: detailed theoretical point-spread-function analysis and numerical simulation results,”
J. Opt. Soc. Am. A, 19 1274
–1285
(2002). https://doi.org/10.1364/JOSAA.19.001274 JOAOD6 0740-3232 Google Scholar
M. T. Stahl, S. B. Shaklan and H. P. Stahl,
“Preliminary analysis of effect of random segment errors on coronagraph performance,”
Proc. SPIE, 9605 96050P
(2015). https://doi.org/10.1117/12.2190160 PSISDG 0277-786X Google Scholar
L. E. Coyle, J. S. Knight and M. Adkins,
“Edge sensor concept for segment stabilization,”
Proc. SPIE, 10698 1069869
(2018). https://doi.org/10.1117/12.2312224 PSISDG 0277-786X Google Scholar
L. E. Coyle et al.,
“Large ultra-stable telescope system study,”
Proc. SPIE, 11115 111150R
(2019). https://doi.org/10.1117/12.2525396 PSISDG 0277-786X Google Scholar
L. Pueyo et al.,
“The LUVOIR extreme coronagraph for living planetary systems (ECLIPS) I: searching and characterizing exoplanetary gems,”
Proc. SPIE, 11117 1111703
(2019). https://doi.org/ 10.1117/12.2530722 PSISDG 0277-786X Google Scholar
P. Hallibert et al.,
“Technologies for large ultra-stable optical missions: current perspectives and developments at ESA,”
Proc. SPIE, 11116 1111610
(2019). https://doi.org/10.1117/12.2529320 PSISDG 0277-786X Google Scholar
D. B. Moore and D. C. Redding,
“Picometer differential wavefront metrology by nonlinear Zernike wavefront sensing for LUVOIR,”
Proc. SPIE, 10698 1069841
(2018). https://doi.org/10.1117/12.2312600 PSISDG 0277-786X Google Scholar
R. Juanola-Parramon et al.,
“The LUVOIR extreme coronagraph for living planetary systems (ECLIPS) II. Performance evaluation, aberration sensitivity analysis and exoplanet detection simulations,”
Proc. SPIE, 11117 1111702
(2019). https://doi.org/10.1117/12.2530356 PSISDG 0277-786X Google Scholar
B. Nemati et al.,
“The effects of space telescope primary mirror segment errors on coronagraph instrument performance,”
Proc. SPIE, 10398 103980G
(2017). https://doi.org/10.1117/12.2273072 PSISDG 0277-786X Google Scholar
B. Nemati et al.,
“Method for deriving optical telescope performance specifications for Earth-detecting coronagraphs,”
J. Astron. Telesc. Instrum. Syst., 6 039002
(2020). https://doi.org/10.1117/1.JATIS.6.3.039002 Google Scholar
L. Leboulleux et al.,
“Sensitivity analysis for high-contrast missions with segmented telescopes,”
Proc. SPIE, 10400 104000M
(2017). https://doi.org/10.1117/12.2274347 PSISDG 0277-786X Google Scholar
L. Leboulleux et al.,
“Sensitivity analysis for high-contrast imaging with segmented space telescopes,”
Proc. SPIE, 10698 106986H
(2018). https://doi.org/10.1117/12.2313904 PSISDG 0277-786X Google Scholar
L. Leboulleux et al.,
“Pair-based Analytical model for Segmented Telescopes Imaging from Space for sensitivity analysis,”
J. Astron. Telesc. Instrum. Syst., 4 035002
(2018). https://doi.org/10.1117/1.JATIS.4.3.035002 Google Scholar
I. Laginja et al.,
“Wavefront error tolerancing for direct imaging of exo-Earths with a large segmented telescope in space,”
Proc. SPIE, 11117 1111717
(2019). https://doi.org/10.1117/12.2530300 PSISDG 0277-786X Google Scholar
M. N’Diaye et al.,
“Apodized Pupil Lyot coronagraphs for arbitrary apertures. V. Hybrid shaped pupil designs for imaging earth-like planets with future space observatories,”
Astrophys. J., 818 163
(2016). https://doi.org/10.3847/0004-637X/818/2/163 ASJOAB 0004-637X Google Scholar
E. S. Douglas et al.,
“WFIRST coronagraph technology requirements: status update and systems engineering approach,”
Proc. SPIE, 10705 1070526
(2018). https://doi.org/10.1117/12.2314221 PSISDG 0277-786X Google Scholar
D. S. Acton et al.,
“James Webb Space Telescope wavefront sensing and control algorithms,”
Proc. SPIE, 5487 887
–896
(2004). https://doi.org/10.1117/12.551846 PSISDG 0277-786X Google Scholar
L. E. Coyle et al.,
“Ultra-stable telescope research and analysis (ULTRA) program, phase 1 report,”
(2019). http://astrostrategictech.us/pdf/projectfiles/Reports/17-SLSTD17-0003_FR_Ball_2019_04_Public.pdf Google Scholar
M. East, C. Sullivan and C. Wells,
“Ultrastable mirror assembly design (Conference presentation),”
Proc. SPIE, 11116 111160Z
(2019). https://doi.org/10.1117/12.2529575 PSISDG 0277-786X Google Scholar
C. Wells and M. East,
“Mirror segment sensitivity analysis and performance budgets for LUVOIR ULTRA (Conference presentation),”
Proc. SPIE, 11116 1111603
(2019). https://doi.org/10.1117/12.2529405 PSISDG 0277-786X Google Scholar
R. J. Noll,
“Zernike polynomials and atmospheric turbulence,”
J. Opt. Soc. Am. (1917-1983), 66 207
–211
(1976). https://doi.org/10.1364/JOSA.66.000207 Google Scholar
R. Soummer, C. Aime and P. E. Falloon,
“Stellar coronagraphy with prolate apodized circular apertures,”
Astron. Astrophys., 397 1161
–1172
(2003). https://doi.org/10.1051/0004-6361:20021573 AAEJAF 0004-6361 Google Scholar
M. N’Diaye, L. Pueyo and R. Soummer,
“Apodized Pupil Lyot coronagraphs for arbitrary apertures. IV. Reduced inner working angle and increased robustness to low-order aberrations,”
Astrophys. J., 799 225
(2015). https://doi.org/10.1088/0004-637X/799/2/225 ASJOAB 0004-637X Google Scholar
G. Foo, D. M. Palacios and G. A. Swartzlander,
“Optical vortex coronagraph,”
Opt. Lett., 30 3308
–3310
(2005). https://doi.org/10.1364/OL.30.003308 OPLEDP 0146-9592 Google Scholar
D. Mawet et al.,
“Ring-apodized vortex coronagraphs for obscured telescopes. I. transmissive ring apodizers,”
Astrophys. J., 209 7
(2013). https://doi.org/10.1088/0067-0049/209/1/7 Google Scholar
E. H. Por,
“Phase-apodized-pupil Lyot coronagraphs for arbitrary telescope pupils,”
Astrophys. J., 888 127
(2020). https://doi.org/10.3847/1538-4357/ab3857 ASJOAB 0004-637X Google Scholar
I. Laginja,
“PASTIS: v2.0.0 semi-analytical PASTIS for LUVOIR,”
(2020). Google Scholar
A. C. Rencher and G. B. Schaalje, Linear Models in Statistics, 2nd ed.Wiley-Interscience(2008). Google Scholar
L. Pueyo et al.,
“The LUVOIR architecture “A” coronagraph instrument,”
Proc. SPIE, 10398 103980F
(2017). https://doi.org/10.1117/12.2274654 PSISDG 0277-786X Google Scholar
E. H. Por et al.,
“Exploiting symmetries and progressive refinement for apodized pupil Lyot coronagraph design,”
Proc. SPIE, 11443 114433P
(2020). https://doi.org/10.1117/12.2562571 PSISDG 0277-786X Google Scholar
C. C. Stark et al.,
“Lower limits on aperture size for an ExoEarth detecting coronagraphic mission,”
Astrophys. J., 808 149
(2015). https://doi.org/10.1088/0004-637X/808/2/149 ASJOAB 0004-637X Google Scholar
C. C. Stark et al.,
“ExoEarth yield landscape for future direct imaging space telescopes,”
J. Astron. Telesc. Instrum. Syst., 5 024009
(2019). https://doi.org/10.1117/1.JATIS.5.2.024009 Google Scholar
I. Laginja et al.,
“Predicting contrast sensitivity to segmented aperture misalignment modes for the HiCAT testbed,”
Proc. SPIE, 11443 114433J
(2020). https://doi.org/10.1117/12.2560113 PSISDG 0277-786X Google Scholar
L. Pueyo et al.,
“Stability error budget for exo-earth imaging with a large segmented telescope in space,”
in American Astronomical Society Meeting Abstracts #233, 158.24,
(2019). https://aas233-aas.ipostersessions.com/default.aspx?s=F3-CD-DE-E7-D7-CD-24-BC-8B-57-9C-0B-DA-73-A1-00&guestview=true Google Scholar
M. E. Kasper et al.,
“EPICS: the exoplanet imager for the E-ELT,”
Proc. SPIE, 7015 70151S
(2008). https://doi.org/10.1117/12.789047 PSISDG 0277-786X Google Scholar
M. Ygouf et al.,
“Simultaneous exoplanet detection and instrument aberration retrieval in multispectral coronagraphic imaging,”
Astron. Astrophys., 551 A138
(2013). https://doi.org/10.1051/0004-6361/201220318 AAEJAF 0004-6361 Google Scholar
E. H. Por et al.,
“High Contrast Imaging for Python (HCIPy): an open-source adaptive optics and coronagraph simulator,”
Proc. SPIE, 10703 1070342
(2018). https://doi.org/10.1117/12.2314407 PSISDG 0277-786X Google Scholar
T. E. Oliphant, A Guide to NumPy, 1 Trelgol Publishing(2006). Google Scholar
S. van der Walt, S. C. Colbert and G. Varoquaux,
“The NumPy array: a structure for efficient numerical computation,”
Comput. Sci. Eng., 13 22
–30
(2011). https://doi.org/10.1109/MCSE.2011.37 Google Scholar
J. D. Hunter,
“Matplotlib: a 2D graphics environment,”
Comput. Sci. Eng., 9 90
–95
(2007). https://doi.org/10.1109/MCSE.2007.55 Google Scholar
T. A. Caswell et al.,
“matplotlib/matplotlib: Rel: v3.1.1,”
(2019). Google Scholar
T. P. Robitaille et al.,
“Astropy: a community Python package for astronomy,”
Astron. Astrophys., 558 A33
(2013). https://doi.org/10.1051/0004-6361/201322068 AAEJAF 0004-6361 Google Scholar
A. M. Price-Whelan et al.,
“The astropy project: building an open-science project and status of the v2.0 core package,”
Astron. J., 156 123
(2018). https://doi.org/10.3847/1538-3881/aabc4f ANJOAA 0004-6256 Google Scholar
The Astropy Collaboration,
“Astropy: a core python package for astronomy,”
(2018). Google Scholar
BiographyIva Laginja is a PhD candidate with ONERA, the French aerospace lab and previous member, now collaborator, of the Russell B. Makidon Optics Laboratory at the Space Telescope Science Institute in Baltimore, USA. She received BSc degrees in astronomy and physics at the University of Vienna, Austria, and a MSc degree in astronomy and instrumentation at Leiden University, Netherlands. Her research interests include WFS&C, high contrast imaging and coronagraphy for exoplanet detection. Rémi Soummer is an associate astronomer at STScI. He received his doctorate from the University of Nice in 2002 and has been working in the field of high contrast imaging and instrumentation for the detection and characterization of exoplanets ever since. He is currently the head of the Russell B. Makidon Optics Laboratory and working on a coronagraph demonstration for future large-segmented aperture space telescopes. Laurent M. Mugnier is a research director at ONERA, where he is in charge of data processing activities of the High Angular Resolution team. His research interests include astronomical and satellite observation from the ground with adaptive optics, exoplanet detection, Earth observation, retinal imaging for ophthalmology, free-space optical telecommunications. He has been the (co)advisor of more than 30 PhD students. He has authored 8 contributions to textbooks and more than 75 papers in peer-reviewed international journals. Laurent Pueyo is an associate astronomer at STScI. He received his doctorate from Princeton University in 2008 and previously worked as a NASA Fellow at JPL, California and as a Sagan Fellow at JHU. His research focuses on imaging faint planets around nearby stars. He has pioneered optical technologies that allow astronomers to take images of other planetary systems and has developed data analysis methods now standardly used to study extrasolar planets. Jean-François Sauvage is a researcher specialized in adaptive optics since almost 15 years at ONERA, and has worked in collaboration with LAM (Laboratoire d’Astrophysique de Marseille) for 5 years. He specialized in the development of eXtreme AO systems for astronomy (in the SPHERE instrument), and in particular he developed COFFEE, a new focal-plane wave-front sensor dedicated to high-contrast imagers. Lucie Leboulleux is a postdoctoral researcher working at LESIA/Observatoire de Paris in the field of direct imaging of exoplanets. Her work aims at developing analytical models to replace E2E simulations of instruments that have an enormous numerical burden and are very time consuming. She also worked on the HiCAT testbed at the STScI on the experimental validation of a wavefront sensor. She recently received the National Rising Talents L’Oréal-UNESCO For Women in Science Award. Laura Coyle is a senior optical engineer at Ball Aerospace. She received her BA degree in physics from Colgate University in 2009 and a PhD in optical sciences from the University of Arizona in 2014. She has over ten years of experience in the field of precision optical alignment and metrology. Her recent work has focused on developing technologies for segment active sensing and control and advancing coronagraph performance-based stability budgeting for optical systems. J. Scott Knight is a senior staff consultant in Missions Systems Engineering at Ball Aerospace. His current focus is payload architecture for large space systems. He has been involved with segmented space telescopes for over 15+ years concentrating on system engineering, modeling, test design for ground-based verification and flight commissioning. |