The Nancy Grace Roman Space Telescope will carry a coronagraph instrument (CGI) that will serve as a demonstrator for technologies needed for future high-contrast imaging missions in space, including deformable mirrors (DMs) to correct high-order wavefront errors that would otherwise limit the achievable contrast. The CGI has three baselined interchangeable observing configurations, one of which is a bowtie shaped pupil coronagraph for high-contrast spectroscopy. We present the flight designs for two closely related mask configurations of the bowtie shaped pupil coronagraph: a baseline 0-deg mask configuration for the technology demonstration and a 60-deg mask configuration contributed by the NASA Exoplanet Exploration Program. The shaped pupil mask and Lyot stop for each mask configuration result from an iterative process that maximizes the core throughput subject to constraints on other performance metrics, such as the contrast: a linear program optimizes the shaped pupil mask for a given Lyot stop, and the optimization repeats for various Lyot stops until the highest-throughput combination is identifiable. The flight designs for the baseline and rotated mask configurations have core throughputs of 4.50% and 3.89%, respectively, at 4 λD and are robust to conservative estimates of potential pupil errors such as misalignments and manufacturing errors. If these estimates are exceeded in flight, the DMs can be used to mitigate the effects of the excess error. |
1.IntroductionSince time immemorial, people have looked to the stars and wondered if there might be life elsewhere in the universe. This question has appeared in some form in every decadal survey for astronomy since 1972,1–6 and over the decades, multiple techniques have been developed to detect exoplanets, including those in the habitable zone of their parent star.7–13 Two of these techniques, transits and direct imaging, can be used to spectroscopically analyze the light from the exoplanet, thereby enabling a search for biomarkers that suggest the possibility of life. Currently, a space-based direct imaging mission provides one of the best paths to the ultimate goal of searching for life on an exoplanet similar to Earth, an exoplanet of Earth-like size in the habitable zone of a Sun-like star.13 Large space-based direct imaging mission concepts, such as the Habitable Exoplanet Observatory (HabEx)14 and the Large UV/Optical/IR Surveyor (LUVOIR),15 are being developed with this goal in mind, and a 6-m class infrared/optical/ultraviolet telescope for high-contrast imaging and spectroscopy is listed in the most recent decadal survey as the highest-priority recommendation for space frontier missions.6 However, the challenge of directly imaging and spectroscopically characterizing an Earth-like exoplanet is not insignificant; there are currently numerous technical gaps that will need to be addressed.16,17 To advance our general understanding of high-contrast imaging in space, NASA’s Nancy Grace Roman Space Telescope will carry a coronagraph instrument (CGI) that will serve as a demonstrator for some of the technology that will be needed for future missions. Notably, the CGI will be the first space-based coronagraph to use deformable mirrors (DMs) to correct high-order wavefront errors that would otherwise limit the achievable contrast, to incorporate low-order wavefront sensing and control, and to use precision-fabricated masks specifically designed to maximize the throughput.18,19 As a technology demonstrator, the CGI has three baselined interchangeable observing configurations, one of which is for high-contrast spectroscopy. (The other two observing configurations are for imaging: a narrow-field-of-view hybrid Lyot coronagraph designed for a 10% spectral band centered at 575 nm, and a wide-field-of-view shaped pupil coronagraph designed for a 10% spectral band centered at 825 nm; imaging over other spectral bands is possible but unsupported.20–22) The spectroscopy configuration uses a shaped pupil coronagraph to achieve a contrast of within a two-sided bowtie-shaped region spanning 3.0 to for wavelengths across a 15% spectral band centered at either 660 or 730 nm, and the coronagraph masks are followed by a slit-prism spectrograph with a spectral resolution of at the center of the selected spectral band.20,21,23 Due to the shape of the high-contrast region (also called the dark hole), this coronagraph is referred to as the bowtie shaped pupil coronagraph. The bowtie shaped pupil coronagraph has two closely related mask configurations: a baseline mask configuration that will be used for the technology demonstration, and a mask configuration that is rotated 60 deg relative to the baseline mask configuration. (To accommodate two different 15% spectral bands, there are two variations of each of these mask configurations. However, because the two variations of a given mask configuration differ only in their physical scaling, we will not discuss them individually in this paper. A full list of every designed configuration for the CGI may be found in a paper by Riggs et al.21) The rotated mask configuration is not part of the technology demonstration but contributed by the NASA Exoplanet Exploration Program (ExEP) so that it can be available for use if scheduling and funding constraints allow, and it is important to note that the rotated mask configuration is being accommodated on the grounds that it does no harm to the baseline design and causes no schedule delay.24 The rationale for contributing the rotated mask configuration is that it improves the annular field of view coverage and increases the probability of detecting our canonical target exoplanet, 47 Ursae Majoris c (47 UMa c); with the two bowties and a -deg roll, the detection probability during the nominal technology demonstration and beyond is nearly as high as it would be if the coronagraph had a 360-deg field of view due to the combination of the two bowtie openings, the bowties precessing, and orbital motion.25–27 (However, we emphasize that no observing time is currently allocated to the rotated mask configuration.) The two mask configurations each produce an appropriately rotated dark hole by using a set of three masks: a shaped pupil mask for the first pupil plane and bowtie-shaped masks for the first focal plane and the Lyot plane (the focal plane mask and Lyot stop, respectively), as shown in Fig. 1. The focal plane mask design is determined by a set of requirements on the field of view and rotated to match the specified orientation; it is identical for the two mask configurations except for the orientation. The shaped pupil mask and Lyot stop, by comparison, are products of a design survey that searches for the combination with the highest throughput. In this survey, a shaped pupil mask is optimized for a particular Lyot stop (Sec. 2), and the optimization is repeated for a range of Lyot stops (Sec. 3). The final designs for the baseline and rotated mask configurations have core throughputs of 4.50% and 3.89%, respectively, for an off-axis source at . These designs are robust to conservative (three-sigma) estimates of potential pupil errors (such as misalignments and manufacturing errors), but in the event that these estimates are exceeded, the DMs that correct for high-order wavefront errors may be used to mitigate the effects of the pupil errors as well (Sec. 4). 2.Shaped Pupil Mask Design ProcessDue to computational constraints, the shaped pupil mask and the Lyot stop are designed using an iterative optimization process that, subject to constraints on the contrast and other performance metrics, seeks to identify the combination with the highest core throughput , which is defined in the official CGI project requirements28 as the ratio of the power passing through the half-max core of the point spread function (PSF) in the final image plane, , to the power passing through the pupil plane in the absence of a shaped pupil mask, : Because the PSF for the bowtie shaped pupil coronagraph has multiple bright lobes, it is important to note that the half-max core of the PSF is defined as the set of all points in the image plane for which the intensity is at least half of the maximum value, not just the points in the main lobe (Fig. 2).The throughput value is dependent on the location of the point source used to calculate the PSF, and we have selected a point source at . This value was chosen because it corresponds to the location of 47 UMa c, our canonical target exoplanet that was chosen because it is a confirmed exoplanet around a star with a V magnitude of almost 5 and a small angular diameter.29 The emphasis on maximizing the core throughput was driven by the fact that it can be one of the primary factors limiting the ability to take spectroscopic measurements; however, in scenarios dominated by speckle instability, the effects of adjusting the throughput are more muted. This metric changed from a more signal-to-noise ratio (SNR) metric (that accounted for PSF sharpness) when the CGI switched from having an integral field spectrograph (IFS) to a slit spectrograph. For the IFS, the area of the PSF was a concern because observations were dark-current limited. With the switch to the slit spectrograph, the planet PSF is no longer broken up across several lenslets, thereby reducing the PSF area on the detector by a factor of a few. As a result, the challenge is to increase the throughput to reduce the photon noise of the planet itself (which is assumed to be faint, with a contrast on the order of ). Throughput is also an easy-to-use proxy for sharpness, while sharpness makes optimizations intractable. The first step of the design process is to construct an appropriate starting point for the shaped pupil mask optimization by creating a representation of the telescope pupil with the correct orientation (Sec. 2.2) and adding padding to compensate for potential pupil errors, such as clocking, rotation, translation, and magnification (Sec. 2.3). The resulting padded telescope pupil serves as the shaped pupil mask optimization starting point throughout an entire design survey, but the padding may be adjusted for different surveys. To complete a design survey, an optimization routine determines the highest-throughput shaped pupil mask for a given Lyot stop (Sec. 2.1), and the optimization is repeated for a range of potential Lyot stops until the highest-throughput combination can be identified (Sec. 3). 2.1.Optimizing a Shaped Pupil MaskThe bowtie shaped pupil masks are optimized using a linear program that maximizes the peak pixel value for an off-axis planet (a proxy for the throughput) while also constraining the contrast and the sensitivity to tip and tilt. This linear program (written in AMPL30) treats all of the equations for propagating light through the coronagraph as constraints, and it incorporates tip and tilt with the addition of separate propagations for each type of aberrated wavefront. A thorough description of how to model the propagation of light through the coronagraph and the specific techniques and equations used by the linear program may be found in papers by Zimmerman and Riggs and are not repeated here.31–33 Directly constraining the sensitivity to tip and tilt leads to shaped pupil masks that are more robust to residual pointing jitter and resolved stellar diameters, but it also reduces the throughput. It is therefore important to balance the benefits of reducing the tip and tilt sensitivity against the benefits of increasing the throughput. The linear program includes a coefficient in the aberrated wavefront propagations that allows the amount of tip and tilt in an optimization to be adjusted, with higher amounts producing shaped pupil masks with lower sensitivity. For the final flight shaped pupil masks, the value of this coefficient was selected based on the results of a limited design survey conducted by the CGI Integrated Modeling Team that evaluated the overall performance of masks optimized with various amounts of tip and tilt included.34 Because the linear program simulates propagating light through the coronagraph to evaluate each potential shaped pupil mask design, it requires knowledge of the coronagraph masks in each plane as well as the desired shape of the dark hole. As a result, two-dimensional arrays representing the shaped pupil mask starting point, the focal plane mask, and the Lyot stop are required inputs to the linear program, and the shape of the dark hole is specified by a set of constraints. To keep the optimization problem tractable, we require -axis symmetry for each of the coronagraph masks and the dark hole and exploit additional symmetry wherever possible.35 The -axis symmetry requirement significantly reduces the size of the optimization problem, but it also means that care must be taken when defining the arrays for the coronagraph masks. For the baseline mask configuration, the dark hole and almost all of the coronagraph masks are symmetric about the axis due to the choice of orientation; the bowties for the focal plane mask, Lyot stop, and dark hole each have two axes of symmetry, and these axes align with the and axes for the 0-deg orientation. The shaped pupil mask, by comparison, is symmetric only because the linear program assumes it to be and outputs only one half of the mask. Because the telescope pupil is slightly asymmetric, we use a symmetrized version as the shaped pupil mask optimization starting point to ensure that the shaped pupil mask completely covers the telescope pupil as assumed (Fig. 3). For the rotated (60-deg) mask configuration, neither the dark hole nor any of the coronagraph masks are symmetric about the axis (as illustrated in Fig. 1). To satisfy the symmetry requirement, we optimize a shaped pupil mask for the rotated mask configuration by rotating the telescope pupil 60 deg counterclockwise and keeping the baseline orientation for the remaining planes; the shaped pupil optimization starting point is therefore a symmetrized version of the rotated telescope pupil (Fig. 3). It is also worth noting that although we can optimize shaped pupil masks for a 90-deg mask configuration because the axes of symmetry for all three bowties again align with the and axes, the throughput is significantly lower due to the less-advantageous orientation of the dark hole with respect to the struts supporting the secondary mirror (Fig. 4). 2.2.Constructing an Appropriately Rotated Unpadded Telescope PupilThe first step to generating the starting point for the shaped pupil mask optimization is to construct a representation of the telescope pupil that has been rotated appropriately for the desired mask configuration. This representation is based on a map of the nominal CGI entrance pupil that was generated using the FRED raytracing software program.36–38 With 7980 and 8054 illuminated pixels along the and axes, respectively, this pupil map has a much higher resolution than the shaped pupil mask optimization process needs (1000 pixels along an axis is sufficient), and because this map is provided as a single array for the entire pupil, it does not enable the manipulation of individual elements, such as the struts, directly. However, the ability to arbitrarily modify any pupil element facilitates adding specialized pupil padding tailored to provide the desired degree of robustness to individual pupil errors (Sec. 2.3); although pupil padding can be added for global errors (translation, magnification, and clocking) by manipulating the provided array as a whole, adding specific pupil padding for component-level errors (misalignments and errors in the sizes of the struts, central obscuration, and tabs) requires working with the individual components directly. The ability to manipulate individual pupil elements also facilitates studying the impacts of both the pupil padding and any pupil errors (Sec. 4). It is therefore preferable to work with a modifiable, lower-resolution representation generated by a telescope pupil model that simplifies the true telescope pupil into a collection of well-defined lines and shapes that can be individually manipulated as desired. In the telescope pupil model, the pupil is divided into multiple components that are each modeled individually: the primary mirror, the secondary mirror, the six struts that support the secondary, and three tabs (which create small obscurations next to the secondary). The primary and secondary mirrors are each represented by an ellipse, and the struts are each represented (at least initially) by a rectangle. The tabs are then represented by sections of ellipses that fall between the appropriate struts. The sets of parameters that define the individual ellipses and rectangles are determined independently with one exception: the two tabs that are not located on the axis for a pupil orientation of 0 deg are assumed to share the same ellipse. A list of the parameter values and an explanation of the methods used to infer them are provided in the Appendix. Once defined, these component models are combined to form the full telescope pupil model (Fig. 5). To provide important cross-checks and additional insights into the behavior of the telescope pupil as it is rotated (for optimizing a 60-deg shaped pupil mask) and padded (modified by blocking additional regions of the pupil plane to enforce symmetry and account for potential pupil errors), we use two independent methods to construct the specific telescope pupil models for the baseline and rotated mask configurations. (The two mask configurations were also designed semi-independently, with the baseline mask configuration primarily designed at the Jet Propulsion Laboratory, JPL, and the rotated mask configuration primarily designed at Princeton, as described in Sec. 3). For the baseline mask configuration, the telescope pupil model is generated with PROPER,39 a library of optical propagation routines developed at JPL, as the superposition of elliptical and rectangular obscurations. For the rotated mask configuration, the telescope pupil model consists of a set of algebraic equations (explained below) that allow the individual edges of the obscurations to be manipulated directly, and the equations are written in a way that allows the pupil to be rotated as needed. Although not identical, the telescope pupil representations resulting from these two models agree very closely when generated for the same orientation, with only negligible differences (Fig. 6). In the algebraic model, we generate a binary representation of the telescope pupil by first defining sets delineated by ellipses and lines that describe the clear area of the unobscured pupil, , and the individual obscurations due to the secondary, struts, and tabs, and then applying these set definitions to a grid of points. By binary, we mean that a point in the pupil plane is either fully clear or fully obscured, with clear points assigned a value of 1 and obscured points (as well as points in the pupil plane that fall outside ) assigned a value of 0. The total obscuration is the union of all of the individual obscurations, and the clear area of the pupil is the set of points that fall within but outside , as illustrated in Fig. 5. To allow the pupil to be rotated arbitrarily, the perimeters of the primary and secondary are described by ellipses of the form where, in the absence of any pupil errors or padding, is the semiminor axis length, is the semimajor axis length, is the orientation of the desired mask configuration (0 deg for the baseline mask configuration and 60 deg for the rotated mask configuration), and and are the and coordinates, respectively, that specify the location of the ellipse center. Because any collected light must land on the primary mirror, consists of the set of points that fall within the primary’s perimeter: where the subscript indicates parameter values specific to the primary. For the secondary, the inequality follows the same pattern, but the meaning is reversed; the set names are chosen such that, for any component of the modeled telescope pupil, is the set of points within (or obscured by) . The set , therefore, describes the obscuration due to the secondary instead of the clear area outside the secondary: where the subscript sec indicates values specific to the secondary. The obscurations due to the tabs are defined similarly, but additional logic is required to restrict each of these obscurations to the subset of points within the ellipse that fall between the appropriate struts and outside .To determine the obscuration due to the struts, it is useful to start by modeling a generic strut as a rectangle of length and half width centered on the origin. A particular strut can then be modeled by rotating this generic rectangle by an angle , where is the nominal orientation of the specified strut, and translating the result so that its center is located at , where and are the and coordinates, respectively, that specify the location of the strut’s center in the absence of any pupil rotation. This rectangular strut model is useful for determining the locations of the strut’s four corners in the absence of any pupil padding and misalignments. Because the addition of pupil padding can modify the shape of the strut’s obscuration so that the long edges are no longer parallel, it is ultimately beneficial to model the strut by using the corners to determine the lines that define the strut’s longer sides and specifying the obscuration as the set of points that fall between these lines (truncated by the perimeter of the primary on one side and the secondary on the other). The specific definition for this set varies depending on the particular strut and the overall pupil orientation because these factors affect which corners define the upper of the two lines as well as how to bound the set to prevent it from poking through the opposite side of the secondary. 2.3.Padding the PupilOnce the telescope pupil model has been constructed, the next step in generating the starting point for the shaped pupil optimization is to add the pupil padding. By padding, we mean that additional regions of the pupil plane are blocked to enforce symmetry (symmetry padding) or to ensure that the shaped pupil mask completely covers the telescope pupil even if errors are present (error padding). The symmetry padding is required to satisfy an assumption made during the optimization process, and the amount that is needed depends on the pupil orientation (Sec. 2.3.2). The error padding, however, is more discretionary; although we can apply enough error padding to compensate for the maximum expected pupil error, we can also use less if there are other means of compensating for at least part of the potential pupil error (Sec. 2.3.1). Both types of pupil padding may be applied, but the application process varies depending on the type of padding and the method used to generate the telescope pupil model. 2.3.1.Padding for pupil errorsBecause the shaped pupil mask is optimized for a specific representation of the telescope pupil, problems can arise when there is a mismatch between the modeled and true telescope pupils. In particular, the contrast and throughput are no longer guaranteed to reach their design values if part of the true telescope pupil is not covered by the shaped pupil mask. This scenario can be caused by manufacturing uncertainties or misalignments, so a conservative approach to preventing the performance degradation is to add padding to the shaped pupil optimization starting point to account for these potential errors up to some specified amount. Although this error padding improves the robustness to pupil errors, it comes at the cost of some of the core throughput. The list of errors and their tolerances therefore plays an important role in the performance of the final design. For both the baseline and rotated mask configurations, we pad for a total of seven errors. Three of these are due to global misalignments (translations, clocking, and magnification), and three are due to component-level manufacturing uncertainties (for the struts; the secondary, also called the central obscuration; and the tabs). The final error, the rolloff, accounts for the fact that light reflected from the very outer edge of the primary is unusable due to the process by which the mirror was ground and polished. The conservative tolerance values assigned to each of these errors are identical for the baseline and rotated mask configurations with one exception: the baseline mask configuration uses a slightly higher value for the translation tolerance based on a more conservative interpretation of a lateral pupil misalignment requirement (Table 1). This choice was made to slightly increase the robustness of the baseline mask configuration and reduce its technical risk. For the rotated mask configuration, which is contributed by ExEP so that it can be available for use if scheduling and funding constraints allow and is not part of the technology demonstration, the smaller value was selected to slightly increase the core throughput (by two-tenths of a percentage point). Table 1Tolerances used for error padding. The diameter D is the full aperture diameter including the rolloff, not the diameter of the smaller aperture once the rolloff is masked.
The error padding is applied to the shaped pupil optimization starting point in one of two ways, both of which account for worst-case scenarios in which all of the errors are present simultaneously at the most extreme values specified by the tolerances. For the baseline mask configuration, the error padding is applied by compiling multiple versions of the array that represents the modeled telescope pupil. First, the component-level error padding and the rolloff padding are applied by modifying the telescope pupil model; the component-level error padding is applied by increasing the widths of the struts and the sizes of the secondary and tabs by the amounts specified in the tolerances, and the rolloff padding is applied by decreasing the radius of the primary. The global error padding is then applied by using PROPER to generate versions of the modified telescope pupil for a range of possible translation, clocking, and magnification combinations and taking the union of the padding cases. For the rotated mask configuration, the error padding is applied by modifying the sets of equations used to generate the individual elements of the telescope pupil model. Following the conventions established for the baseline mask configuration, we first account for the component-level errors and rolloff and then pad for magnification, clocking, and translation. With all of the types of error padding incorporated, the telescope pupil model is no longer a simple collection of ellipses and rectangles, but it is still relatively straightforward. The first modifications are made to the definition of with the goal of padding the primary’s perimeter so that the unusable area is always completely blocked by the shaped pupil mask even in the presence of pupil errors. If the primary had been modeled as the set of points within a circle centered at the origin, the modifications would have been as simple as reducing the circle’s radius to a value that accounts for the rolloff tolerance as well as the magnification and translation tolerances. Conservatively, this value is where is the radius of the unpadded circle, is the negative magnification tolerance, is the absolute value of the translation tolerance, and is the rolloff tolerance. The negative magnification tolerance is used because the unusable portion of the primary comes into view as the telescope pupil is demagnified, and the translation tolerance is multiplied by the square root of two to reflect a conservative interpretation that the telescope pupil can translate by the tolerance amount in and simultaneously.However, because the primary is modeled as the set of points within a general ellipse, we must also account for clocking and the projection effects that modify the shape of the unusable region. We pad for clocking errors by considering the two boundary-case ellipses in addition to the single ellipse that defined the primary. These boundary-case ellipses are rotated by the maximum and minimum values specified by the clocking tolerance, and the sets of points within these ellipses are defined as where the signs in the subscripts of , , and reflect the sign of the clocking value used to define the boundary ellipse; is the absolute value of the clocking tolerance; the factors of and scale the rolloff tolerance to the appropriate amount for each axis of the boundary ellipse; and and are the coordinates locating the center of the boundary ellipse: The clear area of the unobscured primary is then the intersection of the sets of points contained within each of the boundary ellipses and the original ellipse [Fig. 7(a)]: where represents the set of points within the original, unrotated ellipse [Eq. (3)].Although the secondary is also modeled as the set of points within a general ellipse, the modifications to are the opposite of the modifications to because is an obscuration instead of a clear area. Instead of reducing the size of the boundary ellipses and taking the intersection of the sets of points that fall within each ellipse, we now increase the size of the boundary ellipses and take the union of the sets of points that fall within these ellipses and all of the clocking cases that they bound [Fig. 7(b)]. Because the endpoints of the major axis of the secondary’s ellipse trace arcs as the ellipse is rotated from one boundary case to the other, we can represent the union of the sets of points that fall within all clocking cases as the union of the sets of points that fall within the boundary ellipses and the sets of points and that fall within the circular sectors swept out by the major axis: The sets are defined as where, as before, the signs of , , and reflect the sign of the clocking value used to define the boundary ellipse; is the positive magnification tolerance; is the absolute value of the alignment tolerance for the secondary; and and are the coordinates locating the center of the boundary ellipse:To define and , we first define a boundary line between the two circular sectors swept out by the major axis: Then, we consider each circular sector in turn. Although the specific equations and inequalties that define the circular sectors depend on , the overall method for defining and does not. For , we locate the endpoints of the traced arc, define the lines and that connect these endpoints to the origin , and identify the circle of which arc is a part. Then, we define as the set of points within this circle that fall between and and are on the correct side of . The set is defined similarly.The three sets , , and that represent the obscurations due to the tabs are modified in the same manner as . The additional logic that bounds , , and between the appropriate struts is unaffected by the padding, although it does require the boundaries of the padded struts to be defined first. The specific inequalities and equations are orientation-dependent and vary by tab. The six sets that represent the obscurations due to the individual struts are the most complicated to define because there are multiple details that depend on the specific strut and the pupil orientation. The first step is to follow the procedure outlined at the end of Sec. 2.2, with the strut half width increased by to account for the strut tolerance. For simplicity, let and be the boundary lines that define the longer sides of strut , and let be the set of points that fall between and , bounded by the perimeter of the primary on one end and the perimeter of the secondary on the other (Fig. 8). Then, we pad for each of the global errors in turn by identifying the boundary cases that move and the farthest outside and updating the definitions of the lines and accordingly. For example, to pad for magnification errors, we magnify one line by and the other by , choosing the case that increases the strut width. Then, to pad for clocking errors, we rotate one of the new lines by and the other by , again choosing the case that increases the size of the obscuration. Finally, we pad for translation errors by shifting the resulting lines by in and in , with the sign for each shift chosen based on which moves the specified boundary line away from . The signs chosen for one line will be the opposite of those for the other. 2.3.2.Padding for symmetryBecause the shaped pupil mask optimization algorithm requires that the pupil be symmetric about the axis (as explained in Sec. 2.1), the naturally asymmetric telescope pupil model is forced to be symmetric through the addition of pupil padding. This padding is added by manipulating the array that represents the modeled telescope pupil, so, for the analytically modeled telescope pupil, the symmetry padding is applied after any error padding. The set of pixels that correspond to the pupil asymmetries, , is identified by subtracting a mirror image of the left half of the pupil from the right half. Padding is then added to the right half of the pupil by setting the pixel value equal to zero wherever . Similarly, padding is added to the left half of the pupil by setting the pixel value equal to zero wherever the mirror image of is negative. The two padded halves are then recombined to form the symmetrized version of the telescope pupil. Although it is possible to achieve the same result by simply padding one side of the pupil and combining the result with its mirror image, we track the amount of padding added to each side of the pupil in order to transparently explain where each pixel of padding originates. Unlike the error padding, the symmetry padding is affected by the pupil orientation due to the lack of perfect rotational symmetry. (As noted in Sec. 2.1, we optimize shaped pupil masks with the pupil rotated by the specified amount when needed to make the problem tractable; in reality, the focal plane mask, Lyot stop, and dark hole are rotated.) The differences are noticeable even between the 0-deg, 60-deg, and -deg orientations, which are all nearly symmetric about the axis (Fig. 9). The -deg orientations have significantly more symmetry padding than the 0-deg orientation for three reasons, all of which are due to the CGI being mounted slightly off-axis: the ellipses that describe the outer diameters of the primary, the secondary, and the tabs are no longer symmetric about the axis; the two tabs that do not fall on the axis for the -deg orientations are no longer fit to the same ellipse; and the strut pairs are less symmetric. Ultimately, these changes lead to a throughput loss of about half a percentage point compared with the throughput of the equivalent 0-deg mask configuration. For the rotated mask configuration, the -deg orientation was selected over the -deg orientation due to its slightly higher symmetry and therefore slightly better throughput. 3.Flight Design SelectionWhile the shaped pupil mask optimization determines the best shaped pupil mask for a given Lyot stop, the design survey identifies the Lyot stop that leads to the highest-throughput design. The Lyot stop for the bowtie shaped pupil coronagraph is a bowtie-shaped mask rotated 90 deg relative to the focal plane mask; this orientation aligns the Lyot stop with the reflective regions of the shaped pupil mask, and the bowtie shape is well-matched to the regions illuminated by off-axis sources32 (Fig. 10). A Lyot stop is specified by three main parameters: the inner and outer diameters of the opening and the opening angle [Fig. 11(a)]. These parameters determine the basic perimeter of the bowtie, which consists of four line segments and four arcs. To improve manufacturability and reduce the risk of damage during launch, the basic bowtie shape is modified slightly: the eight sharp corners where a line segment and an arc intersect are replaced by fillets [Fig. 11(b)]. Each fillet is a section of a circle of radius that is tangent to the line segment at one end and tangent to the arc at the other, allowing the perimeter to smoothly transition from the line segment to the arc. Although is adjustable, it is held constant unless being optimized explicitly. For the design surveys, we considered Lyot stops with inner diameters of to , outer diameters of to , and opening angles of 80 deg to 100 deg. (The diameter is the full aperture diameter, not the smaller diameter once the rolloff is masked.) Due to time constraints, this parameter space was sampled coarsely at first and then successively finely until a maximum-throughput design could be identified in a parameter cube of resolution in the inner and outer diameters and 1 deg in the opening angle; at this resolution, the throughput changes very slowly around the maximum point, and additional refinement is unnecessary. Although the design surveys for the baseline and rotated mask configurations (Secs. 3.1 and 3.2, respectively) all followed this basic procedure, the specific process varied because the two mask configurations were designed semi-independently (with the baseline mask configuration primarily designed at JPL and the rotated mask configuration primarily designed at Princeton) and under fluctuating time constraints. Despite the differences in the design methodologies, the flight designs for the baseline and rotated mask configurations (Sec. 3.3) are strikingly similar, as would be expected if the telescope pupil were perfectly rotationally symmetric. 3.1.Design Surveys for the Baseline Mask ConfigurationTo quickly narrow the parameter space for the baseline mask configuration, the design process began with a preliminary design survey that coarsely sampled the set of possible Lyot stops while generating shaped pupil masks with a lower resolution than the full manufacturing resolution. Although the lower-resolution () shaped pupil masks are not purely binary as required, they offer comparable performance to full-resolution () shaped pupil masks in simulations while also being significantly faster to optimize. (The shaped pupil mask optimization routine does not guarantee that the resulting shaped pupil mask will be perfectly binary, but when the shaped pupil mask is optimized at the full manufacturing resolution, the contrast does not degrade when the pixel values are rounded to produce a purely binary mask.) The results of this survey suggest that the highest-throughput design corresponds to a Lyot stop with an inner diameter of to , an outer diameter of , and an opening angle of 87 deg to 88 deg (Figs. 12 and 13). In the interest of time, the parameter space for subsequent design surveys was therefore narrowed to inner diameters of to , outer diameters of to , and an opening angle of 88 deg; the opening angle was fixed due to time constraints and the minimal benefit of varying it near the maximum-throughput point. In response to concerns about launch stresses, we conducted an additional low-resolution design survey to determine the radius of curvature for the Lyot stop fillets. The results of this survey indicate that although a radius of leads to the highest throughput, the radius can be increased to without a significant throughput loss (Fig. 14). We therefore selected a radius of to reduce the risk of damage to the Lyot stop during launch. For this radius and an opening angle of 88 deg, a Lyot stop with an inner diameter of and an outer diameter of leads to the highest throughput for the baseline mask configuration, as confirmed by a final high-resolution design survey (Fig. 15). 3.2.Design Surveys for the Rotated Mask ConfigurationThe design surveys for the rotated mask configuration were performed concurrently with the design surveys for the baseline mask configuration and followed a similar process, incorporating information from the baseline design survey results as needed. For example, once the final value for the radius of curvature for the Lyot stop fillets, , was determined from the results of a design survey for the baseline mask configuration (Fig. 14), the value of used for the rotated mask configuration was updated to match. All of the design surveys for the rotated mask configuration generated full-resolution shaped pupil masks, and the differences from one survey to the next consisted of updates to the telescope pupil, tolerances, and other design parameters. Due to time constraints, the parameter space (consisting of possible values for the inner and outer diameters of the Lyot stop and the opening angle) was explored iteratively with successively targeted mini surveys until the maximum-throughput point could be identified, and the final design survey was initially centered on the maximum-throughput point of the penultimate survey (Figs. 16 and 17), which differed only by using a slightly more conservative translation tolerance. The final survey produced two designs with equivalent throughputs (Figs. 18 and 19), but the degeneracy was broken once was updated to its final value. Time constraints prohibited running a separate design survey to incorporate this last update, so the shaped pupil masks for the two highest-throughput point designs were reoptimized with updated Lyot stops; the resulting design with the higher throughput was selected for the rotated mask configuration. 3.3.Final Designs for the Baseline and Rotated Mask ConfigurationsThe flight designs for the baseline and rotated mask configurations are very similar because the orientation of the dark hole relative to the struts is nearly equivalent for the two cases. The shaped pupil masks look almost like rotated versions of each other (Fig. 20), and the Lyot stops for the two mask configurations use very similar parameters: an inner diameter, outer diameter, and opening angle of , , and 88 deg, respectively, for the baseline mask configuration and , , and 89 deg, respectively, for the rotated mask configuration. However, while the two designs are very similar, they have different throughputs; the baseline mask configuration has a core throughput of 4.50% at , whereas the rotated mask configuration has a core throughput of 3.89% at . [For reference, the throughput for the CGI input pupil alone is 34.4% when calculated using Eq. (1); note that this definition neglects all reflectivities and transmittivities and accounts only for the diffraction of the pupil features.] The difference in throughput is due to differences in the amount of pupil padding: the starting point for the shaped pupil optimization for the rotated mask configuration required more symmetry padding (Sec. 2.3.2), and the aperture loss due to the additional symmetry padding was not fully offset by using a slightly less conservative translation tolerance. A full-plane optimization that did not require symmetry likely would have produced baseline and rotated designs with more comparable throughputs because there would have been no need for the symmetry padding, but this type of optimization was not tractable due to the large design space. To briefly comment on the integration time required for spectroscopy with the baseline mask configuration, we first note that specific integration time estimates depend strongly on how conservatively the CGI performance is modeled. The most favorable reflected light planets observable by the CGI are very challenging targets; they have separations near the inner working angle and flux ratios of a few at best, and they are near the performance limit of the CGI, where contributions from random noise and speckle noise are both significant. There may be very few reflected light planets (such as HD 219134 h and 47 Uma c) for which spectroscopy is possible, with integration time estimates ranging from about one day to several weeks. For the interested reader, the Roman Coronagraph Exposure Time Calculator40 is a publicly available tool for estimating the integration times for different planetary system configurations. 4.Compensating for Excessive ErrorsThe flight designs for the baseline and rotated mask configurations are designed to mask the telescope pupil even with errors up to the amounts specified by the tolerances. However, these tolerances are three-sigma tolerances instead of not-to-exceed values, and it is possible (though not expected) for the pupil errors to be large enough that the telescope pupil is partially visible behind the shaped pupil mask, leading to a degradation in the contrast. Although using more pupil padding could have increased the robustness to pupil errors, it is preferable to compensate for any excessive errors by using the DMs that are part of the wavefront sensing and control system for the CGI. 4.1.Costs of Additional Pupil PaddingBecause the pupil padding is a permanent feature of the shaped pupil mask, the throughput losses due to the padding are present even if the pupil errors never approach the tolerance values. These losses can be significant, but they are not evenly distributed among the various types of pupil padding. To illustrate the general trends that describe how the different types of padding contribute to the total throughput loss, we generate two sets of telescope pupils for the rotated mask configuration with varying types and amounts of padding and calculate the amount of clear aperture for each pupil relative to that of the unpadded pupil (a quickly calculated proxy for the throughput loss). First, using the tolerance values, we pad for only one error type at a time to determine which types of padding most significantly affect the amount of clear aperture for the rotated mask configuration’s flight design. Then, we start with the fully padded telescope pupil that we used to optimize the flight design and increase one type of padding beyond its tolerance value to illustrate the cost associated with the additional padding. As shown in Table 2, the rolloff padding is by far the dominant contributor to the aperture loss, followed by the padding for clocking, the struts, magnification, and translation. The padding for the secondary and the tabs, by comparison, minimally affects the aperture. If we increase the padding for any error beyond its tolerance, the amount of clear aperture decreases. The loss is minimal if we increase the padding for the tabs or the secondary, but it is noticeable for the other cases (Fig. 21). The loss is particularly dramatic if we increase the padding for magnification, and this is due to the corresponding change in the rolloff padding. [To prevent the unusable portion of the primary from coming into view, the rolloff padding is increased based on the padding for magnification, translation, and clocking. The clocking tolerance affects the rolloff padding because the modeled pupil is elliptical instead of circular, as explained in Eqs. (6)–(8) and Fig. 7(a)]. With the rolloff padding held constant, additional padding for clocking errors becomes the most costly. Table 2Impact of pupil padding on the clear aperture.
4.2.Using the DMs to Compensate for Excessive ErrorsIn addition to being potentially costly, padding beyond the three-sigma tolerances is unnecessary. To demonstrate this point, we evaluate how excessive errors affect the initial contrast (the contrast calculated without the DMs present and with no wavefront errors) and examine how well the DMs are able to compensate for the excessive error in some of the worst cases. For these tests, we start with a worst-case representation of the telescope pupil that has every error present at its tolerance value. In the case of magnification, which has an asymmetric tolerance, we choose the positive value, and it is also worth noting that we treat magnification as a resizing of the binary telescope pupil; we do not incorporate any edge blurring. Although this representation of the telescope pupil is still completely covered by the shaped pupil mask, the simultaneous presence of all error types at the levels used for padding reduces the extent to which the padding for one type of error can mask a different error. From this starting point, we increase a single error beyond its tolerance and calculate the initial contrast for the rotated mask configuration. Finally, we identify the worst of these error cases and evaluate how well the DMs are able to achieve a target contrast of across a 15% band using an electric field conjugation (EFC) algorithm.41 This performance goal is related to the CGI requirements, and EFC is the baseline wavefront sensing and control algorithm for the CGI. The most dramatic degradation in the initial contrast occurs as the clocking is increased beyond its tolerance. As the clocking increases 3 mrad past its tolerance (to 150% of its tolerance value), the initial contrast degrades by an order of magnitude (Fig. 22). Achieving the same amount of degradation for any of the other error cases that we consider, by comparison, requires increasing the error to of its tolerance. However, this contrast degradation can be mitigated at least partially by the DMs. In our example, even as much as 5 mrad of excessive clocking can be corrected fully by the DMs (Table 4 and Fig. 23), which is fortuitous because clocking is one of the costliest errors to pad. To simulate the DM corrections, we used a suite of modular optical propagation and coronagraph modeling algorithms developed at Princeton. These algorithms implement Fresnel and Fraunhofer propagation using a matrix Fourier transform approach based on that of Soummer et al.,42 and the DM surfaces are modeled by superposing influence functions, which describe how a single actuator affects the DM surface when pistoned. For each actuator, we scale the influence function provided with PROPER43 by the product of the gain (which we assume is for all actuators) and the commanded voltage, and then we sum all of the individual influence functions (each centered over their respective actuators) to produce the modeled surface. To match the DMs used in the CGI, we represent the two DMs as actuator arrays separated by 1 m;20 selected technical information about the CGI DMs is provided for reference in Table 3. Table 3Selected technical details for the DMs. This information (and more) can be found in a paper by Poberezhskiy et al.20
The extent to which excessive errors can be corrected depends in part on the magnitude of the error and the choice of the dark hole. It is also worth noting that the ability to compensate for pupil errors is also dependent on the ratio of the square of the DM width to the DM separation, with smaller ratios preferred, but this ratio is fixed for the Roman Space Telescope. To illustrate the general patterns that describe how well the errors can be corrected, we consider two of our worst excessive error cases (11 and 14 mrad of total clocking) and use an EFC algorithm to improve the contrast over either a two-sided dark hole or a single-sided dark hole (Tables 4 and 5 and Figs. 23 and 24). Because the specific number of iterations required depends upon the size of the dark hole and the magnitude of the pupil error, the algorithm is set to run for 30 iterations, the maximum number of allowed wavefront sensing and control iterations planned for digging the initial dark holes on orbit for the baseline configurations of the hybrid Lyot coronagraph (HLC), the bowtie shaped pupil coronagraph, and the wide-field-of-view shaped pupil coronagraph. The postcorrection values shown in Tables 4 and 5 correspond to either the first iteration for which the contrast is beneath the target across the entire band or to the iteration for which the contrast has converged as close to the target as possible. Table 4Control results for 11 mrad of total clocking (5 mrad beyond the tolerance).
Table 5Control results for 14 mrad of total clocking (8 mrad beyond the tolerance).
As we would expect, smaller errors require fewer control iterations and are more likely to be successfully corrected over the full bandwidth for a two-sided dark hole (Table 4 and Fig. 23). As the magnitude of the error increases, there is eventually a point at which there is simply too much error for the DMs to correct completely across the full band for a two-sided dark hole. At this point, the edges of the band fail to reach the target contrast. As the error continues to increase, the contrast at the outer edges of the band continues to deteriorate, and the fully corrected portion of the band shrinks in toward the central wavelength. However, reducing the size of the dark hole by switching to a single-sided dark hole can enable the DMs to fully correct the error across the band once again (Table 5 and Fig. 24), and there are contributed field stops included in the CGI that allow for single-sided bowtie dark holes if that becomes necessary. It is also worth noting that, for each of the example cases, the actuator heights and commanded voltages (shown in Table 6) are well within their limits (Table 3). For the cases in which the DMs successfully corrected the error, the largest absolute change in actuator height was 9.4 nm, which is about 2% of the stroke limit. For the example in which the DMs were unable to completely correct the error, the largest absolute change in actuator height was 22.4 nm. The root-mean-square (RMS) stroke for the nominal HLC design is about 30 nm RMS before accounting for any optical aberrations, whereas for the shaped pupil coronagraphs it is nominally zero because the apodizer does most of the contrast generation. Therefore, a change in actuator height for the bowtie shaped pupil coronagraph in these cases with misalignments is well within the DM actuation capabilities. Overall, there is a range of excessive errors that can be completely corrected by the DMs across the full band, and beyond this range, the contrast degrades gracefully. Table 6Maximum and minimum voltage commands and actuator heights for the clocking error cases.
5.SummaryCurrently, a space-based direct imaging mission offers one of the best paths to searching for life on an Earth-like exoplanet, and the level of interest in such a mission is underscored by the most recent decadal survey’s recommendation of a 6-m class infrared/optical/ultraviolet telescope for high-contrast imaging and spectroscopy as the highest priority for space frontier missions. To help close some of the numerous technical gaps that need to be addressed, the CGI on the Nancy Grace Roman Space Telescope will serve as a demonstrator for high-contrast imaging technologies including DMs to correct high-order wavefront errors and precision-fabricated masks specifically designed to attenuate the relative intensity of the star off-axis while maximizing the core throughput. The CGI has three interchangeable observing configurations, one of which is a bowtie shaped pupil coronagraph for high-contrast spectroscopy. The bowtie shaped pupil coronagraph has a baseline 0-deg mask configuration for the technology demonstration and a 60-deg mask configuration contributed by ExEP. Each mask configuration consists of a shaped pupil mask, a bowtie-shaped focal plane mask, and a bowtie-shaped Lyot stop. The designs of the shaped pupil mask and the Lyot stop are determined through an iterative process: a linear program optimizes the shaped pupil mask for a given Lyot stop, and the optimization is repeated for a range of potential Lyot stops until the highest-throughput combination is identifiable in a parameter space cube of resolution in the inner and outer diameters of the Lyot stop and 1 deg in the opening angle. The resulting flight designs for the baseline and rotated mask configurations have core throughputs of 4.50% and 3.89%, respectively. To provide some robustness to errors, the shaped pupil mask optimization directly incorporates tip and tilt, and the optimization starting point is padded so that the telescope pupil is fully masked for conservative estimates of misalignments and manufacturing errors. This padding is specifically tailored to a list of error tolerances as opposed to being generated by rotating, translating, and scaling a provided high-resolution pupil array, and this new approach for generating robust shaped pupil designs is made tractable by a completely parametric model of the telescope pupil (which we also use to study the costs of adding various types of robustness and the impacts of pupil errors not masked by the padding). We have developed two versions of this model: one that generates a pupil by superposing elliptical and rectangular obscurations in PROPER, and one that consists of a set of algebraic equations that allow the individual edges of the obscurations to be manipulated directly. These equations are written in a way that allows the pupil to be rotated, facilitating designs for rotated mask configurations. There is also some padding to enforce -axis symmetry, which is required to keep the optimization problem tractable, and the amount of symmetry padding depends on the mask configuration orientation. However, the pupil padding costs throughput, and the throughput difference between the baseline and rotated mask configurations is largely due to the higher amount of symmetry padding required by the rotated mask configuration. The throughput loss due to pupil padding is also unevenly distributed among the various types of error padding, with the rolloff padding costing the most throughput and the padding for symmetry, the tabs, and the central obscuration costing the least. Because the throughput losses due to pupil padding are fixed, it is important to balance the need for pupil padding against the throughput cost. The baseline and rotated mask configurations are padded to conservative three-sigma estimates of various pupil errors, which is expected to provide sufficient robustness. However, if the pupil errors exceed the estimates, the DMs can be used to compensate at least partially. In our simulations, excessive clocking errors lead to the most dramatic contrast degradation by far, but the DMs are able to completely compensate for the error across a 15% band in even the worst cases (11 mrad and 14 mrad of total clocking, with padding for only 6 mrad). The one caveat is that a single-sided dark hole is required for the 14 mrad case, while a two-sided dark hole can be used for the 11 mrad case. In general, the DMs are able to compensate for excessive pupil errors up to a point, after which the contrast degrades gracefully; the fully corrected portion of the band shrinks inward toward the central wavelength, and the corrected contrast at the ends of the band degrades as the error increases. Future missions with obscured apertures may benefit from extending the use of the DMs a step further: selectively reducing the pupil padding to improve the throughput and relying on the DMs to compensate for the unpadded error as needed. The successful implementation of this approach and the extent to which the padding can be reduced depend strongly on the ratio of the square of the DM width to the DM separation, in addition to depending on details including the specific telescope pupil obscuration and the magnitudes of the expected errors. The throughput gains depend on the specific situation as well, but preliminary results suggest the gains can be as high as a percentage point or two in some cases. Continued investigations of the throughput gains achievable by reducing the pupil padding and relying on the DMs to compensate for pupil errors as needed, and the factors that affect the success and usefulness of this approach, are the subjects of future work. 6.Appendix: Pupil Model ParametersThe parameters for the two telescope pupil models were determined independently through separate trial-and-error approaches with different goals, and due to tight time constraints, much of the work was done by hand. For each model, we inferred a set of parameter values that reproduced the provided high-resolution array. Then, as a final step, we shifted the coordinates to center the modeled primary. (The pupil position in the provided array is arbitrary.) For the PROPER-based model (used to develop the baseline mask configuration), the overall goal was to find a set of parameters that led to the best fit, which minimized the summed squared difference between the provided high-resolution pupil array and the PROPER-generated model. The first step was to determine the approximate lateral position and size (width or radius) of each pupil component by eye. Then, this initial pupil model was refined through an iterative grid search in which one to three parameters were varied at a time for each individual obscuration. For the algebraic model (used to develop the rotated mask configuration), the overall goal was to find a conservative set of parameters that led to a pupil with the smallest obscuration that still completely covered the obscuration in the provided pupil. The first step was to determine a set of parameters that produced a pupil with an obscuration that covered every pixel of the obscuration in the provided pupil; this was easily checked by subtracting the provided pupil from the modeled pupil. Then, using a plot of the difference of the modeled and provided pupils as a guide, the parameters were refined by hand for one pupil component at a time. Changes that reduced the number of pixels inside the modeled obscuration were kept so long as the provided obscuration was completely covered, and changes that increased the number of obscured pixels were discarded. The process ended when any further adjustments either increased the number of obscured pixels or uncovered any part of the obscuration in the provided pupil. Although the parameters for the two pupil models were inferred through different methods with different goals, the final numerical values are very similar, as shown in Table 7. The differences in the angles of the last three struts reflect differences in the angle ranges used to define the nominal strut orientations. It is also worth noting that while the algebraic model briefly uses the strut length, the exact value is unimportant as long as it is positive. Table 7Pupil model parameters for the baseline and rotated mask configurations.
AcknowledgmentsPart of the work reported here was performed at Princeton University under NASA Contract NNG16PJ29C, and part was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004). ReferencesNational Research Council, Astronomy and Astrophysics for the 1970’s, Volume 1: Report of the Astronomy Survey Committee, National Academy Press, Washington, D.C
(1972). Google Scholar
National Research Council, Astronomy and Astrophysics for the 1980’s, Volume 1: Report of the Astronomy Survey Committee, National Academy Press, Washington, D.C
(1982). Google Scholar
National Research Council, The Decade of Discovery in Astronomy and Astrophysics, National Academy Press, Washington, D.C
(1991). Google Scholar
National Research Council, Astronomy and Astrophysics in the New Millennium, National Academy Press, Washington, D.C
(2001). Google Scholar
National Research Council, New Worlds, New Horizons in Astronomy and Astrophysics, The National Academies Press, Washington, D.C
(2010). Google Scholar
National Research Council, Pathways to Discovery in Astronomy and Astrophysics for the 2020s, National Academies Press, Washington, D.C
(2021). Google Scholar
A. Wolszczan and D. A. Frail,
“A planetary system around the millisecond pulsar PSR1257 + 12,”
Nature, 355 145
–147
(1992). https://doi.org/10.1038/355145a0 Google Scholar
M. Mayor and D. Queloz,
“A Jupiter-mass companion to a solar-type star,”
Nature, 378 355
–359
(1995). https://doi.org/10.1038/378355a0 Google Scholar
I. A. Bond et al.,
“OGLE 2003-BLG-235/MOA 2003-BLG-53: a planetary microlensing event,”
Astrophys. J., 606 L155
–L158
(2004). https://doi.org/10.1086/420928 ASJOAB 0004-637X Google Scholar
G. Chauvin et al.,
“A giant planet candidate near a young brown dwarf,”
Astron. Astrophy., 425 L29
–L32
(2004). https://doi.org/10.1051/0004-6361:200400056 Google Scholar
M. W. Muterspaugh et al.,
“The phases differential astrometry data archive. V. Candidate substellar companions to binary systems,”
Astron. J., 140 1657
–1671
(2010). https://doi.org/10.1088/0004-6256/140/6/1657 ANJOAA 0004-6256 Google Scholar
W. J. Borucki,
“KEPLER Mission: development and overview,”
Rep. Prog. Phys., 79 036901
–50
(2016). https://doi.org/10.1088/0034-4885/79/3/036901 RPPHAG 0034-4885 Google Scholar
National Academies of Sciences, Engineering, and Medicine, Exoplanet Science Strategy, National Academies Press, Washington, D.C
(2018). Google Scholar
, “Exploring new worlds – understanding our universe,”
(2019) https://www.jpl.nasa.gov/habex/pdf/HabEx-Final-Report-Public-Release-LINKED-0924.pdf Google Scholar
, “The LUVOIR mission concept study final report,”
(2019) https://asd.gsfc.nasa.gov/luvoir/reports/ Google Scholar
B. Crill and N. Siegler,
“Space technology for directly imaging and characterizing exo-Earths,”
Proc. SPIE, 10398 103980H
(2017). https://doi.org/10.1117/12.2275697 PSISDG 0277-786X Google Scholar
B. Crill and N. Siegler,
“Exoplanet Exploration Program 2019 technology plan appendix,”
(2019). https://exoplanets.nasa.gov/internal_resources/1123/ Google Scholar
H. Tang et al.,
“The WFIRST coronagraph instrument phase B optical design,”
Proc. SPIE, 11117 111170C
(2019). https://doi.org/10.1117/12.2530506 PSISDG 0277-786X Google Scholar
N. J. Kasdin et al.,
“The Nancy Grace Roman Space Telescope Coronagraph Instrument (CGI) technology demonstration,”
Proc. SPIE, 11443 114431U
(2020). https://doi.org/10.1117/12.2562997 PSISDG 0277-786X Google Scholar
I. Y. Poberezhskiy et al.,
“Roman Space Telescope Coronagraph: engineering design and operating concept,”
Proc. SPIE, 11443 114431V
(2021). https://doi.org/10.1117/12.2563480 PSISDG 0277-786X Google Scholar
A. J. E. Riggs et al.,
“Flight mask designs of the Roman Space Telescope Coronagraph Instrument,”
Proc. SPIE, 11823 118231Y
(2021). https://doi.org/10.1117/12.2598599 PSISDG 0277-786X Google Scholar
J. Trauger et al.,
“Hybrid Lyot coronagraph for WFIRST-AFTA: coronagraph design and performance metrics,”
J. Astron. Telesc. Instrum. Syst., 2
(1), 011013
(2016). https://doi.org/10.1117/1.JATIS.2.1.011013 Google Scholar
T. Groff et al.,
“Roman Space Telescope CGI: prism and polarizer characterization modes,”
Proc. SPIE, 11443 114433D
(2021). https://doi.org/10.1117/12.2562925 PSISDG 0277-786X Google Scholar
V. Bailey,
“Contributed mask summary,”
(2020). Google Scholar
V. Bailey,
“Possible mask options,”
(2020). Google Scholar
E. Nielsen,
“WFIRST: 47 UMa c and the Bowtie,”
(2020). Google Scholar
E. Nielsen,
“Update of additional bowtie orientations 47 uma detectability,”
(2020). Google Scholar
R. Boucarut,
“RST Coronagraph Instrument Requirements Document (CGIRD),”
Goddard Space Flight Center, Greenbelt, Maryland
(2021). Google Scholar
V. Bailey,
(2020). Google Scholar
R. Fourer, D. M. Gay and B. W. Kernighan, AMPL: A Modeling Language for Mathematical Programming, 2nd ed.Duxbury Thomson(2003). Google Scholar
N. T. Zimmerman et al.,
“Shaped pupil Lyot coronagraphs: high-contrast solutions for restricted focal planes,”
J. Astron. Telesc. Instrum. Syst., 2 011012
(2016). https://doi.org/10.1117/1.JATIS.2.1.011012 Google Scholar
A. Riggs et al.,
“Shaped pupil coronagraph design improvements for the WFIRST coronagraph instrument,”
Proc. SPIE, 10400 104000O
(2017). https://doi.org/10.1117/12.2274437 PSISDG 0277-786X Google Scholar
A. J. E. Riggs, G. Ruane and B. D. Kern,
“Directly constraining low-order aberration sensitivities in the WFIRST coronagraph design,”
Proc. SPIE, 11117 111170F
(2019). https://doi.org/10.1117/12.2529588 PSISDG 0277-786X Google Scholar
B. Nemati,
(2020). Google Scholar
R. J. Vanderbei,
“Fast Fourier optimization,”
Math. Prog. Comput., 4 53
–69
(2012). https://doi.org/10.1007/s12532-011-0034-8 1867-2949 Google Scholar
A. Riggs,
“Summary of decisions on CGI pupil usage,”
(2020). Google Scholar
“2020-05-13 CGI entrance pupil 8k as-built uncoated optics.fits,”
(2020) http://github.com/ajeldorado/falco-matlab/blob/master/lib/masks/WFIRST/PhaseC/pupil_CGI-20200513_8k_binary_noF.png Google Scholar
, “FRED overview,”
(2020) https://photonengr.com/wp-content/uploads/2020/04/FREDOverview19.40.2.pdf Google Scholar
J. E. Krist,
“PROPER: an optical propagation library for IDL,”
Proc. SPIE, 6675 66750P
(2007). https://doi.org/10.1117/12.731179 PSISDG 0277-786X Google Scholar
S. R. Hildebrandt,
“Roman coronagraph exposure time calculator,”
https://roman.ipac.caltech.edu/sims/ETC.html 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
(1), 011009
(2016). https://doi.org/10.1117/1.JATIS.2.1.011009 Google Scholar
R. Soummer et al.,
“Fast computation of Lyot-style coronagraph propagation,”
Opt. Express, 15
(24), 15935
–15951
(2007). https://doi.org/10.1364/OE.15.015935 OPEXFF 1094-4087 Google Scholar
J. Krist et al.,
“PROPER: an optical propagation library for IDL, Python, and Matlab,”
(2021). https://proper-library.sourceforge.net Google Scholar
BiographyJessica Gersh-Range is an associate research scholar at Princeton’s High Contrast Imaging Laboratory, where she optimized shaped pupil Lyot coronagraph masks to support technology development efforts for the Nancy Grace Roman Space Telescope. She was the recipient of a NASA Graduate Student Researchers Program fellowship at Marshall Space Flight Center, and she also worked at the Space Telescope Science Institute as a graduate student. She received her BA degree in physics from Swarthmore College and her PhD in mechanical engineering from Cornell University. A J Eldorado Riggs is an optical engineer at the Jet Propulsion Laboratory. His research focuses on the high-contrast imaging of exoplanets, in particular, mask optimization and wavefront sensing and control for the Coronagraph Instrument on the Nancy Grace Roman Space Telescope. He received his BS degree in physics and mechanical engineering from Yale University and his PhD in mechanical and aerospace engineering from Princeton University in 2011 and 2016, respectively. N. Jeremy Kasdin is the Eugene Higgins professor of mechanical and aerospace engineering, emeritus, at Princeton University and the former assistant dean for engineering at the University of San Francisco. He received his PhD from Stanford University in 1991. After being the chief systems engineer for NASA’s Gravity Probe B spacecraft, he joined the Princeton faculty in 1999, where he researched high-contrast imaging technology for exoplanet imaging. From 2014 to 2016, he was the vice dean of the School of Engineering and Applied Science at Princeton. From 2014 to 2021, he was the adjutant scientist for the coronagraph instrument on NASA’s Nancy Grace Roman Space Telescope. |