Overview and Reassessment of Noise Budget of Starshade Exoplanet Imaging

High-contrast imaging enabled by a starshade in formation flight with a space telescope can provide a near-term pathway to search for and characterize temperate and small planets of nearby stars. NASA's Starshade Technology Development Activity to TRL5 (S5) is rapidly maturing the required technologies to the point at which starshades could be integrated into potential future missions. Here we reappraise the noise budget of starshade-enabled exoplanet imaging to incorporate the experimentally demonstrated optical performance of the starshade and its optical edge. Our analyses of stray light sources - including the leakage through micrometeoroid damage and the reflection of bright celestial bodies - indicate that sunlight scattered by the optical edge (i.e., the solar glint) is by far the dominant stray light. With telescope and observation parameters that approximately correspond to Starshade Rendezvous with Roman and HabEx, we find that the dominating noise source would be exozodiacal light for characterizing a temperate and Earth-sized planet around Sun-like and earlier stars and the solar glint for later-type stars. Further reducing the brightness of solar glint by a factor of 10 with a coating would prevent it from becoming the dominant noise for both Roman and HabEx. With an instrument contrast of 1E-10, the residual starlight is not a dominant noise; and increasing the contrast level by a factor 10 would not lead to any appreciable change in the expected science performance. If unbiased calibration of the background to the photon-noise limit can be achieved, Starshade Rendezvous with Roman could provide nearly photon-limited spectroscopy of temperate and Earth-sized planets of F, G, and K stars<4 parsecs away, and HabEx could extend this capability to many more stars<8 parsecs. (Abridged)


Introduction
Direct imaging of exoplanets from space holds promise to write a new chapter in astronomy and planetary science. With most of the exoplanets discovered to date in tightly bound orbits of their host stars, and thus uninhabitable unless the stars are much fainter than the Sun, direct imaging would detect planets in the habitable zones 1 of more Sun-like stars. If some of these planets are small and predominantly rocky in composition, they may have environments hospitable for life.
One of the primary goals of exoplanet direct imaging is to search for temperate and small planets of nearby stars and study the chemical composition of their atmospheres with spectroscopy.
A starshade working in tandem with a space telescope provides one of the best near-term opportunities to achieve this goal. Starshade is an external occulter flown along the line of sight from a telescope to a target star. With the shape designed to mimic the optical effects of an optimally apodized screen, a starshade can create a "deep shadow" where the starlight is suppressed by 10 10 . 2,3 The telescope kept in this shadow would be able to detect planets and disks around the star at very high contrast. A starshade that prevents the starlight from entering the telescope would allow many simplifications of the telescope optics. For example, precise wavefront control would not be necessary, which also reduces the number of reflections before feeding the light to a detector, and thus increases the optical throughput of the instrument. The costs of these benefits are the added complexity of formation flying and the complications involved in launching and deploying the large and optically precise starshade.
Two advanced mission concepts to discover Earth-like planets in the habitable zones of Sunlike stars being considered at NASA would use starshades as one of the starlight suppression techniques. The Roman Space Telescope 4, 5 will be capable of collecting starlight reflected by large exoplanets with its coronagraph instrument. The Starshade Rendezvous Probe, 6 an advanced mission concept, would further enable Roman to search for Earth-sized planets in the habitable zones of ∼ 10 nearby stars, with the possibility to obtain their limited-bandwidth spectra at a moderate resolution (∼ 70). The Habitable Exoplanet Observatory (HabEx, 7 ), a concept of a 4-m space telescope with a starshade, has the main objective to search for Earth-sized planets in a larger stellar sample and obtain their spectra in a wider band and at a higher resolution. Starshade Rendezvous with Roman and HabEx, both with a starshade, would have the spectral characterization of small planets in the habitable zones of nearby stars as the key and probably limiting science objective.
To enable these potential exoplanet science missions, NASA's Exoplanet Exploration Program (ExEP) is executing a directed and focused activity, the Starshade Technology Development Activity to TRL5 (S5). For S5, TRL5 is defined as demonstrating critical performance in relevant environments at the subsystem level with medium fidelity prototypes. The technology development plan of S5 1 adopts formation flying and observation scenarios of the Rendezvous and HabEx mission concepts as the baseline. Completion of S5 would bring starshade technologies to TRL5 for both concepts. Specifically, S5 includes experiments and analyses to demonstrate small-scale starshade masks that could reach 10 −10 instrument contrast at the inner working angle (IWA) at a flight-like Fresnel number, to develop an optical edge for the starshade petals that would limit scattered sunlight (i.e., solar glint) to acceptable levels, to demonstrate the ability to sense the lateral offset between the starshade and the telescope to an accuracy of 30 cm, and to demonstrate the ability to design and manufacture the starshade mechanical elements that could meet the contrast requirement. Together with S5, ExEP has chartered a Science and Industry Partnership to engage the broader science and technology communities during the execution of the S5 activity.
Given the completion of most of S5's technology milestones on instrument contrast, 8, 9 solar glint, 10 and formation flying, 11 we are motivated to revisit the noise budget of starshades' application in exoplanet imaging. While estimates of exoplanet yields from starshade-assisted imaging have been published, 12-14 these works have not mapped the technology progress to science performance with the level of detail of the present paper. Particularly, previous works have not explicitly included solar glint in their noise budget. The purpose of this paper is to update the expected performance of starshade-enabled exoplanet imaging in light of new constraints on starlight and stray 1 https://exoplanets.nasa.gov/internal resources/1033/ light suppression resulted from S5 work and to assess the noise budget of the spectral characterization of temperate and small planets of nearby stars. We will focus on the performance parameters that are directly related to exoplanet imaging, and defer the assessment on mechanical precision tolerance and stability -which controls the instrument contrast -to an error budget analysis 15 and the S5 technology development plan. We will focus on revealing the dominating noise term under a wide range of realistic planet scenarios, to guide the priorities of future development. We will first evaluate the science performance on the assumption that the background could be calibrated to the photon-noise limit, and then discuss the impact of temporal and spatial variability on background calibration.
The paper is organized as follows. We first describe the model used to derive the S/N of a starshade-enabled exoplanet observation based on performance parameters in § 2. In § 3, we discuss a range of stray light sources that may enter the telescope and are not included as a potentially dominant term in the noise budget. § 4 presents the expected S/N for observing nearby planetary systems with the current performance demonstrated by S5 and with potential future development.
We discuss imperfect background calibration due to temporal and spatial variability and the sensitivity of exozodi levels in § 5, and conclude with future prospects in § 6.

STARSHADE PERFORMANCE MODEL
where N P , N S , N G , N E , N Z , and N D are counts from the planet, the star, the solar glint, the exozodiacal dust, the local zodiacal dust, and detector noise, and C is the instrument contrast. The parameters α and β in Eq. (1) result from background subtraction, which will be discussed in Section 2.1. The stray light sources shown in Fig. 1 but not included in Eq. (1) will be discussed in Section 3. The counts are defined by a photometric aperture, which in turns corresponds to the point spread function (PSF) of the telescope. We assume a photometric aperture diameter of λ/D, which is consistent with the S5 milestone reports. [8][9][10] Our choice of photometric aperture would encircle 46% of the flux from a point source, and is smaller than Ref., 14,16 which is 1.4λ/D.
The S/N yielded from the λ/D photometric diameter is approximately 7% less than the theoretical maximum achieved at 1.4λ/D in the background-limited regime. The aperture size in reality may eventually be controlled by the fixed pixel scale of the detector. Parameters used to estimate the contribution of the background sources and the expected S/N are summarized in Table 1.

Background Removal
The quality of background subtraction is characterized by the parameters α and β in Eq. (1).
When β = 0, the background can be subtracted to the photon-noise limit, where α describes the fractional increase from the photon noise that would manifest in the result. For example, if image processing can be approximated by subtracting two adjacent pixels of equal background, one with the planet and the other without, α = 2. For another example, if many more pixels can be used to characterize the background, α would approach unity. One may thus reasonably expect that with largely smooth and static background, β → 0 and α ∼ N P is the count from the planet and N B is the count from the dominant background or noise term.
Eq. (1) can be converted to the following form: where S/N 0 is the S/N when β = 0. Eq. (2) expresses how much S/N would degrade due to imperfect background calibration. Several insights can be observed from Eq. (2). First, when K 1 (i.e., the planet dominates over the background), the S/N is not prone to degradation due to imperfect background calibration. Second, when K 1 (i.e., the background dominates over the planet), or at large S/N 0 (i.e., large ∆T ), the asymptotic S/N would be S/N → K/β = N P /βN B .
Third, at the critical case of β = K S/N 0 , the S/N would become S/N 0 / √ 2, i.e., degraded by a factor of √ 2. Note that the value of β can be different for each source of background; as we will show later, one source would usually dominate for each observation. Here we neglect potential interference between the noise sources and its potential contribution to the β term.
The analysis here tells us that the fundamental limit of planet detection is determined by the flux ratio between the planet and the dominant background source (K in Eq. 2), and how well image processing can subtract the background to the photon-noise limit (β in Eq. 2). Causes for not achieving the photon-noise limit include detector systematics (as is the case for transit observations with Hubble and Spitzer 17 ) and temporal variability of speckles in coronagraphic direct imaging. 18 These causes do not apply to future starshade direct imaging because of the use of EMCCD and the decoupling between starlight suppression from telescope optics. The capability of deep imaging provided by starshade may however render other causes to be the limiting factor, for example, the temporal variability of residual starlight and solar glint and the spatial distribution of exozodiacal dust. We adopt α = 2 and β = 0 in the analyses that follow, and come back to discuss this point in Section 5.

Residual Starlight
Based on the S5 milestone reports, 8,9 we adopt a contrast ratio of 10 −10 at the IWA resulting from an imperfect starshade. As the habitable zones of nearby stars are often substantially larger than the separation as approximately θ −3.4 , where θ is the off-axis angle. This is slightly steeper than the expected Airy pattern drop of θ −3 for a filled aperture, due to the distribution of the residual starlight (i.e., some near the center of the starshade and some localized near the petals). Lastly, an off-axis companion beyond the IWA does not experience a change in transmission for starshade observations. 20

Solar Glint
While the sun will not be on the telescope-facing side of the starshade, the starshade's edges will scatter sunlight towards the telescope via a combination of diffraction, diffuse reflection, and specular reflection. 21 When these scattering mechanisms are considered together, the telescope will see the scattered sunlight coming mainly from localized regions on a few petals where the  Table 1 Parameters adopted in this work. *The magnitude is defined such that a 25-magnitude glint would have the same energy flux in the wavelength band in parenthesis as a solar-spectrum point source that has a magnitude of 25 relative to an A0V star, analog to the Johnson magnitude system. optical edge is aligned for specular reflection (Fig. 1). These will appear as two broad lobes due to the telescope's finite spatial resolution, i.e., the solar glint. The solar glint is a unique effect in direct imaging using a starshade.
We adopt the brightness of the solar glint measured by S5 with a razor-sharp, amorphous metal edge, 10 and use SISTER to calculate its angular dependency (Fig. 2). The expected magnitudes of the solar glint are calculated by combining scattering measurements of prototype optical edge segments and the optical models of the starshade. 10 We adopt a "worst case" scenario, i.e., the maximum solar angle (∼ 83 • ) and 95% confidence upper limit of the brightness at the IWA. 10 The expected brightness of the solar glint for HabEx is approximately one-order-of-magnitude less than that for Roman. This is because the starshade of HabEx would be much more separated from the telescope than the starshade from Roman ( Table 1). The solar glint features a hump that peaks at the IWA and quickly drops as the angular separation increases, and thus impacts planet search near the IWA most significantly.
Recognizing the potential impact of the solar glint on science performance, S5 and its Science and Industry Partnership have been actively seeking improvement on the optical edge technology. A multi-layer, thin-film coating has recently emerged as a highly promising design, and experiments have indicated that the coating would result in a solar glint brightness lower than the uncoated design by approximately one order of magnitude. 22 We will also study the science performance with the solar glint brightness suppressed by the coating.

Exozodiacal Light
Results from the most sensitive exozodiacal dust survey indicate that the majority of nearby Sunlike stars have relatively low habitable-zone dust levels, with the best-fit median to be 3 times the Solar-System level, while some stars (e.g., Eridani) are significantly more dusty. 23 In this work, we assume "1 zodi" corresponds to 22 mag arcsec −2 in the V band at the habitable zone 16 and typically assume "3 zodis" in the analysis. This assumption was also adopted by mission concept studies and exoplanet yield analyses. 6,7,14 We discuss the sensitivity of the exozodi levels in Section 5.3. We also assume that the brightness of the exozodiacal light is independent of the distance or the stellar type if evaluated at the habitable zone 2 , besides a factor from the spectral shape of the star. 16 Scaling to the habitable zone, we approximate the brightness of the exozodiacal light as dependent on the semi-major axis as a −2.44 . This scaling is based on example "Zodipic" simulations with a Solar-System dust density profile. 19

Detector Noise
We model the detector noise as the combination of dark current and clock-induced charge. 14 The EMCCD would have effectively zero read noise. We assume the frame rate to be 6.73 times the count rate of the brightest pixel. 14 This is to ensure that the probability of two photons arriving at the brightest pixel to be less than 1%, at the expense of increasing the clock-induced charge.
The number of pixels for each spectral element is assumed to be 42 for Roman and 56 for HabEx.
These are estimated for spectral characterization at 700 nm with the assumptions of a PSF core of 4 pixels, dispersed into 6 pixels per spectral elements, and the detector providing Nyquist sampling at 400 nm for Roman and 300 nm for HabEx. 14

Overview of Stray Lights
Here we overview other noise terms not accounted for in the primary noise budget (Eq. 1).

Astrophysical Background
The very dark shadow created by the starshade will reveal not only planets but also faint stars and distant galaxies. The Exo-S final report 24 provided a detailed analysis of background star and galaxy confusion and suggested mitigation strategies. If planets within 5 AU from the parent stars are potentially detectable (farther planets would often be too faint in reflected light), the largest area for planet search would be ∼ 6 arcsec 2 . One can expect ∼ 2 distant galaxies down to V∼ 31 25 and ∼ 0.2 stars in this area depending on the galactic latitude. 24 Stars with known companions that have small angular separation and would impact direct imaging have been excluded from the target lists of Starshade Rendezvous 6 with Roman and HabEx. 7 Spectra may provide clues to tell the planets apart from these background sources, and revisits and the detection of common proper motions with the parent star would be required to confirm the planets.

Micrometeoroid Holes
Micrometeoroids can penetrate the starshade's multi-layer opaque optical shield (OS). The baseline OS design consists of three evenly spaced layers of Black Kapton. Some large or high-velocity micrometeoroids can produce through-holes on the starshade. During a science integration, onaxis starlight can pass directly through the fraction of particle trajectories aligned to the starshade normal to disperse only via diffraction towards the telescope (Fig. 1). Off-axis sunlight instead requires multiple reflections within the OS before exiting to disperse in the Lambertian fashion towards the telescope (Fig. 1). It is thus necessary to consider micrometeoroid holes and the transmission of starlight and sunlight through them.
We have developed a model to estimate micrometeoroid holes and the resulting brightness levels for Starshade Rendezvous with Roman. The mass-flux distribution of micrometeoroids at L2 is estimated with the Grun model. 26 The Grun model does not account for seasonal meteor showers that bring elevated fluxes by 1 -2 orders of magnitude. It may be necessary to orient the starshade to a near edge-on direction during one or two showers per year, for a total of approximately one month per year. The lost observation time can be mitigated to some extent by the timing of retargeting maneuvers with long coast periods.
The minimum particle mass required to enter the OS is computed with a single layer ballistic equation from NASA's Preferred Reliability Practices document for micrometeorite protection (Standard PD-EC-1107). The minimum particle mass to pass the middle layer and then exit the third layer are computed with a two-layer ballistic equation. 27 The exit computation conserva-tively neglects the benefit of the middle layer. The incoming flux is assumed isotropic in direction considering variable starshade pointing throughout the mission. The particle specific density is conservatively bounded by a constant of 2.5 g cm −3 for considering the ability to penetrate the OS and 1 g cm −3 for considering the size of holes produced. Entry holes diameters approximately match the particle diameter for large particles and can be greater than the particle diameter by up to a factor of five for small particles. 28 When a high-velocity particle penetrates the first layer, the particle and the shield material will vaporize, creating an expanding gas cloud. A subset of the gas cloud can then penetrate the middle and exit layers. We estimate that the exit and middle-layer hole diameters can grow by up to a factor of thirty from the original particle diameter to account for the gas cloud expansion. 29 Our model indicates that the hole areas after three years on-orbit would be 0.1 parts-per-million (ppm) by surface area entry holes on both sides, 50 ppm exit holes on both sides, and 150 ppm middle-layer holes. The expected number of entry holes is about 5 × 10 5 , while the expected number of exit holes is < 400. The exit holes are produced only by large particles.
We estimate that ∼ 10% through-holes would provide a direct path for starlight to the telescope.
The leakage is limited by the entry hole area on both sides to yield 0.02 ppm of effective area.
Feeding this area to our optical performance model of starshades, 30 we estimate that the starlight leakage due to micrometeoroid holes to correspond to a residual starlight contrast of 10 −13 , lower than the allowed residual starlight contrast by three orders of magnitude.
We then estimate the upper bound of solar transmission to be the product of the porosity factors of the layers and obtain 4 × 10 −13 . The solar leakage is also proportional to micrometeoroid hole area, but not necessarily the number of holes. This transmission corresponds to a brightness of solar leakage at the telescope that has a visual magnitude of 39.6 at the IWA at 700 nm, with the brightness falling off with the off-axis angle at approximately the same rate as for the solar glint.
We thus expect the solar leakage due to micrometeoroid holes to be dimmer than the solar glint by > 10 magnitudes.
While the estimates apply for Starshade Rendezvous with Roman, they can be scaled for HabEx considering a longer (5-year) mission lifetime and a higher angular resolution (defined by λ/D).
The distance from the starshade to the telescope does not appear in this scaling because the solid angle per resolution element is independent of the distance. We estimate that the brightness from micrometeoroid holes on HabEx's starshade would have approximately the same magnitude as the Starshade Rendezvous, and remain much dimmer than the solar glint.

Bright-body Reflections
Bright celestial bodies positioned on the telescope side of the starshade may have a portion of their light reflected towards the telescope, with the brightness falling off with the off-axis angle at approximately the same rate as for solar glint. The starshade presents mostly Black Kapton to the telescope and we evaluate bright-body reflections with its BRDF data. To estimate the magnitude of bright-body reflections, we start from a Solar-System planet's absolute magnitude, and consider the worst-case orbit phasing to derive the brightness incident at the starshade. We then take into account the angular size of the starshade and its reflectivity in the direction of the telescope. Lastly, we distribute the brightness into each resolution element on the starshade, as the starshade would have nearly uniform brightness in this problem. Table 2 provides a detailed walk-through of our estimates of bright-body reflections for Solar-System bodies.
Light from Earth, Moon, and Venus can be reflected towards the telescope at grazing incidence angles, but only at the extremes of starshade's orbit position and pointing ( Fig. 1 Table 2 Estimates of Earth, Venus, and Jupiter brightness reflected by the starshade to the telescope at their orbital positions and the starshade's position and pointing that maximize these stray light sources. The estimates use the parameters of Starshade Rendezvous with Roman, and the result for HabEx would scale by the solid angle per resolution element. *The absolute magnitude of a Solar-System object is defined as the apparent magnitude that the object would have if it were 1 AU from both the Sun and the observer, and in conditions of ideal solar opposition. 31 maximally can appear at about the same magnitudes of 31.7 mags at the IWA (Table 2). To have this magnitude, Venus must be at the quadrature orbital phase viewed from the starshade when the starshade is at one extreme of its L2 orbit, a combination that would rarely occur. The Moon will be much dimmer than Earth and Venus for starshade reflection to the telescope.
Light from Mars and Jupiter can be reflected towards the telescope at nearly normal incidence angles ( Fig. 1). We find that Jupiter will never exceed a brightness of 31.9 mags at the IWA (Table   2), and Mars will be much dimmer.
In addition to the Solar-System bodies, the integrated light from the center of the Milky Way will likely be the brightest object that can be reflected towards the telescope by the starshade.

Other Stray Light
Other stray light may include secondary solar reflections, fluorescence, thruster exhaust solar scatter, and stray light produced by the telescope. The first three potential sources are subjects of starshade technology and our preliminary analyses indicate that these light sources likely contribute to the background by no more than a fraction of the solar glint. S5 is conducting extensive stray light analyses to evaluate possible stray light paths and, if necessary, adjust the starshade design and the observation constraints to mitigate their impacts.
The current starshade design is intended to preclude any sunlight reaching the telescope after only a single reflection, except for the optical edges (i.e., the solar glint). However, secondary solar reflections, where sunlight can reach the telescope after two reflections, are possible with certain particular light paths and they are being studied by S5. Preliminary analysis suggests an acceptable brightness level. For example, some stray light paths involve an out-of-plane deformed petal and stray light analyses are being developed to establish a tighter, yet still readily achievable requirement of how far the petal can deform out-of-plane versus the current specification based on diffraction performance alone. For another example, some solar leakage may occur through small gaps between the optical shield at the inner disk to petal interface. This interface is dynamic at the end of disk deployment as a hard tie-down point rotates into node plates at the ends of each inner disk facet. A labyrinth seal is designed to attenuate the sunlight with a preliminary verification of performance, and a more thorough stray light analysis is planned to confirm this design.
Black Kapton used in the starshade's OS is known to glow in response to high energy solar electrons due to material fluorescence. Informed by past studies for JWST, we estimate a relatively dim magnitude of 33.3 at the IWA, i.e., dimmer than the solar glint (with coating) by two orders of magnitude. Efforts are underway to better understand the input electron energy levels and evaluate the fluorescence of Black Kapton and other material refinements.
Periodic thruster firings to maintain formation flight are expected about once every 10 minutes. 11 Exhaust particles will scatter sunlight at a level that can saturate detectors. An operating concept to prevent detector saturation is for the starshade to notify the telescope of an imminent thruster firing and for the telescope to switch to a fast detector read-out mode for some conservative duration (e.g., 10 seconds). This is more than enough time for the bulk of the exhaust to leave the field of view and only causes a modest loss of observing time. Lastly, large exhaust particles might "loiter" around the starshade or stick to the exhaust nozzle to be pushed off at relatively low velocity at the next thruster firing. The optical density of these loitering particles would be much lower than the high-velocity exhaust particles, and detailed simulations are required to assess their brightness.

RESULTS AND ANALYSES
We choose stars that represent the target lists of Starshade Rendezvous 6 and HabEx 7 in this study.
As seen in Fig. 3 For each representative star, we evaluate the detectability of planets of varied size and planetstar separation. Rather than focusing on 1-R ⊕ planets that receive Earth-like insolation, we include larger planets at colder or hotter orbits as potential search targets. This is motivated by the diversity of planets revealed by current exoplanet searches and improved understanding of the habitable zones. Kepler planet surveys and planetary mass measurements available for a subset of the detected planets have indicated two populations of small planets. 33 The planets in the < 1.7 R ⊕ population are mostly rocky, 34 and the planets in the 1.7 ∼ 3.0 R ⊕ population could be planets with H 2 /He gas envelopes 35,36 and/or substantial water layers. 37,38 The larger planets, if they have moderate-size atmospheres, can host liquid-water oceans and thus be potentially habitable. 39,40 The "habitable zones" for planets with H 2 -dominated atmospheres can be substantially more separated from the parent stars than those for planets with N 2 -and CO 2 -dominated atmospheres. 39,41 We focus on the spectral characterization of exoplanets in this paper. We adopt a requirement of S/N=20 per spectral element, as this is a conservative estimate of the precision needed to measure atmospheric abundances from the reflected-light spectra and potentially distinguish the types of planets. [42][43][44][45][46][47] The expected integration time to detect planets in the broadband will be substantially less than what is shown in this section. Also, all results shown in this section assume unbiased calibration of the background to the photon-noise limit. We will come back to discuss background calibration in Section 5.  Table 1). We make the following observations from Fig. 4. First, exozodiacal light is the dominant noise term for planet searches around nearby F and A stars. These early-type stars have widely separated habitable zones, and thus the habitable-zone planets would have small planet-to-star contrast. Also, because we assume the exozodiacal brightness scales to the habitable zone (the 1-AU equivalent in this context), the expected exozodiacal light brightness near the IWA is high for these early-type stars.

Current Performance of Starshade Rendezvous
Second, for the search of small (< 2.5R ⊕ ) planets near the habitable zones of the nearest late  Table 1.
G and K stars (i.e., τ Ceti and Indi A), the planet itself is likely the dominant noise term. In other words, the current performance expected for Starshade Rendezvous with Roman could provide a photon-limited detection, rather than background-limited detection, for temperate and small planets around these most favorable targets. This is because the stars are less luminous than the Sun, which makes the habitable zones closer to the stars (0.4 -0.7 AU). The planet-star contrast at their habitable zones would be 7 -20×10 −10 for 1-R ⊕ planets. Coupled with the brightness of these stars (visual magnitude of 3.5 -4.7), the brightness of 1-R ⊕ planets in their habitable zones would have a visual magnitude of ∼ 26.3, i.e., much brighter than a "mean" scenario of 30. As such, the common belief that the detection would be mainly limited by exozodiacal light does not apply for these most favorable targets. They are also the prime targets from an integration time standpoint: to spectroscopically characterize (with S/N=20) a 1-R ⊕ planet near the habitable zone only requires integration of a few hundred hours, and the integration time reduces to ∼ 100 hours for larger rocky planets with 1.5-R ⊕ radius. Among these stars, Eridani has an especially bright dust disk that would probably prevent the search of small planets. 23 τ Ceti have two super-Earth-mass planets (M sin i ∼ 3.9M ⊕ ) near the habitable zone detected by radial-velocity measurements, 48,49 and an outer dust disk (> 6 AU) detected by far-infrared and radio observations. [50][51][52] The constraints on the orbital elements and masses of these planets would aid future characterization.
Third, for slightly farther stars (5 -6 parsecs, e.g., 40 Eridani, δ Pavonis, 82 Eridani, σ Draconis), the search for small planets near the habitable zones would be limited by exozodiacal light and solar glint. For the later-type stars whose habitable zones are close to the IWA, the search for habitable-zone planets are particularly affected by solar glint. With the current performance 10 and assuming calibration of the solar glint to the photon-noise limit, to characterize a 1-R ⊕ planet near the habitable zone requires integration of a few thousand hours, i.e., likely infeasible for a realistic space mission. Characterizing a 1.5 ∼ 2.5-R ⊕ planet, however, only requires an integration time of 100 ∼ 1000 hours. Large rocky planets and more volatile-rich planets could thus be studied with spectroscopy, and thus these "super-Earths" represent the near-term opportunity for the search of habitable worlds on these stars. 53 Finally, for stars ∼ 8 parsecs away, planet characterization near the habitable zones would limit to planets not smaller than ∼ 2 R ⊕ , and the integration time required for smaller planets would be well longer than 1000 hours. The solar glint continues to be the dominant noise term, except for the F and A stars. As discussed in Section 2.3, one may reasonably expect the actual solar glint brightness to be approximately one order of magnitude lower than the nominal value, due to coating technologies. 22 Also, we have used a "worst" scenario in the analysis presented in Fig. 4, in that the Sun was at the maximum angle from the line of sight (83 • ) and the planet was close to the position where the solar glint was maximized. An "average" scenario where the solar angle is less (∼ 53 • ) and the planet is away from the glint maximum would result in a reduction in the solar glint brightness by a factor of ∼ 3. 10 Figure 5 shows the impact on the science performance from a 10-fold reduction on the solar glint brightness. We find that this improvement would drastically enhance the prospect of the search for small planets around later-type nearby stars. The solar glint would no longer be the dominant noise term for the planet observations around these stars. The exozodiacal light becomes the dominant noise term in place. Consequently, a decrease in the required integration time by  Fig. 4, except that the brightness of solar glint is lower than the nominal value by a factor of 10 and the residual starlight contrast is 3 × 10 −9 (i.e., higher than the nominal value by a factor of 30). a factor 2 ∼ 3 for planet search near the IWA would be expected (Fig. 6). This improvement is particularly substantial for exoplanet science, because this would reduce the time needed to characterize 1-R ⊕ planets in the habitable zones of stars 5 ∼ 6 parsecs away from a few thousand to approximately 1000 hours (i.e., from infeasible to marginally feasible), and would reduce the time needed to characterize 1.5-R ⊕ planets (i.e., large rocky planets) from ∼ 1000 hours to a few hundred hours (i.e., from marginally feasible to highly feasible). This improvement would also extend the search for planets around 8-parsec, Sun-like stars from 2-R ⊕ planets to 1.5-R ⊕ planets,

Impact of Future Development
i.e., into the regime of rocky planets. Together, these impacts would enlarge the stellar sample that would be amenable for searching rocky planets in the habitable zones by a factor of at least three, and thus a substantial improvement of the science prospect.
For starlight suppression, the nominal performance shown in Fig. 4 does not have residual starlight as the dominant noise term anywhere in the explored parameter space. We are thus motivated to assess the tolerance of exoplanet characterization on a degraded residual starlight contrast.
We have repeated the analysis with a contrast level that is 3, 10, 30, and 100 times higher than the   Table 1 in these simulations. With the additional optical edge coating, 22 the impact of the solar glint would be eliminated.  (Table 1). Based on nominal performance parameters, solar glint would be the dominant noise term for observing Earth-sized and smaller planets in the habitable zones of the latest-type stars in each distance group. This is because the habitable zones of these stars approach the IWA, where the solar glint is the brightest.

Expected Performance of HabEx
If the optical edge coating is applied and the solar glint brightness is reduced by a factor of 10, this impact of the solar glint would be eliminated.
We see that photon-limited detection and characterization of planets can be expected in the majority of the parameter space explored. Particularly, HabEx would be able to perform photonlimited planet characterization for Earth-sized planets in the habitable zones of nearby stars as far as ∼ 6 parsecs. The integration time needed to achieve an S/N=20 spectrum at the resolution of R = 140 is ∼ 10 hours for the closest (3 -4 parsecs) stars and 50 ∼ 100 hours for stars 5 -6 parsecs away. These estimated integration times indicate more than one-order-of-magnitude improvement over the Starshade Rendezvous with Roman. The improvement comes from not only the larger telescope, but also higher throughput, better angular resolution to reduce the exozodiacal light interference, and lower solar glint brightness.
Exozodiacal light starts to affect the observations of Earth-sized planets around stars ≥ 8 parsecs away. Fig. 8 shows that the exozodiacal light is the dominant noise term for planets < 1.5 R ⊕ in the habitable zones of ∼ 8-parsec stars. If the exozodiacal light can be subtracted to the photonnoise limit, however, the planets may still be observed and characterized with S/N=20 spectra within a few hundred hours. If the parameter α in Eq. 1 is further reduced from two to unity, the integration time in the exozodiacal-dominant regime would be smaller than what is shown in Fig.   8 by a factor ∼ 2.
The potential targets for HabEx include more distant stars up to ∼ 15 parsecs (Fig. 3). While not explicitly shown in Fig. 8, we estimate the impact of solar glint and exozodiacal light on the performance of planet observations around these farther stars. For the photometric aperture does not have solar glint as the main noise term and will be increasingly impacted by exozodiacal light for stars 8 -15 parsecs away.

Requirements for Background Calibration
To this point we work on the assumption that the background can be removed by imaging processing to the photon-noise limit (i.e., β = 0 in Eq. 1). Since the solar glint and the exozodiacal light dominate over the planetary light in the noise budget in substantial fractions of Roman's and HabEx's search spaces for small planets (Figs. 4, 5, and 8), the ability to calibrate these noise terms is important for delivering the science capability described in Section 4. Here we first quantify the precision needed for background calibration. Fig. 9 shows the ratio between the planetary flux and the background flux from solar glint and exozodiacal light (i.e., N P /βN B , or K in Eq. 2). As discussed in Section 2.1, the maximum S/N achievable is N P /βN B when β is nonzero. Therefore, the ratio shown in Fig. 9, divided by the desired S/N (i.e., 20 in this work), is the maximum tolerable residual fraction (β) of the background after calibration.
For the closest stars (< 4 parsecs), the flux from a 1-R ⊕ planet in the habitable zone is similar to the flux of the exozodiacal light and is greater than the flux of solar glint by more than one order of magnitude (Fig. 9). This means that the exozodiacal light must be calibrated to 5% to allow an S/N of 20. For farther stellar systems (5 ∼ 6 parsecs), the requirement for the calibration precision of exozodiacal light becomes 2%, and that for the calibration of solar glint is 5%. For stars ∼ 8 parsecs away, the limiting science objective for a reasonable integration time with Roman would be to characterize 1.5-R ⊕ planets in the habitable zones (Fig. 5). For this, the requirement for the calibration precision of exozodiacal light and solar glint remains 2% and 5%, respectively.
Taken together, the residual exozodiacal light should be less than 2% (i.e., β exozodi ≤ 0.02) and the residual solar glint should be less than 5% (i.e., β glint ≤ 0.05) to avoid significant adverse impact on the spectral characterization of small and temperate planets. At these critical precision levels of background calibration, the final S/N would be degraded from the S/N with photon-noise-limit calibration by a factor of distribution. In-situ production of the dust by collisions and giant impacts would lead to bright and localized dust "clumps", 56 but these events are rare in mature planetary systems, which is the case for most stars in the target lists of Roman and HabEx. Dust particles may also be transported to the habitable zones by Poynting-Roberson drag, 57 or delivered as comets scattered by outer planets and then sublimated. 58 The latter two processes would eventually produce a smoothly varying dust density profile, 59 and may enable high-precision background removal. However, dust particles may be trapped in mean-motion resonance with the planets 60, 61 and become "clumps." In addition to spatial structures, exozodiacal dust particles should have very different spectral shapes than the planets, which may provide another way to distinguish them. We stress that unlike other background terms, the exozodiacal signal and the density profile it implies is a science objective in its own right. Its spatial distribution and origin would be studied together with the planets.

Variability of the Background
Another potential challenging aspect of background calibration is that the background may vary during the long integration typically required for planet characterization. As demonstrate in Section 4, the solar glint is the dominant background term compared to other stray light sources and residual starlight. The solar glint brightness is expected to vary, as it depends on the lateral position of the starshade with respect to the line of sight. The next-in-the-line stray light source, reflection of earthshine and the Milky Way (Section 3.3), can be variable as well. 62 The solar glint and other stray light to starshade exoplanet imaging is the analog of "speckles" to coronagraph exoplanet imaging, and the variability of speckles has driven the design of the coronagraph instrument on Roman. 63 Here we provide a high-level estimate on the impact of background variability on starshade exoplanet imaging. We approximate the background variability as a random variation, because the variation of the starshade's lateral position would have a shorter timescale (< 10 minutes 11 ) compared to typical integration times. A change in the starshade's position is equivalent to a change in the off-axis angle for this problem, and the brightness of the residual starlight and the solar glint would have a fractional change as a function of the off-axis angle (Fig. 2). We thus approximate the background flux with a log-normal distribution with the parameter 0 and σ (Lognormal(0, σ)) with respect to the median value. We perform Monte Carlo simulations of photons arriving onto two photometric apertures on the detector. Both apertures receive photons from the varying background and one of the pixel additionally receive photons from the planet. The count difference between the two apertures -and the statistical distribution of it -thus indicates how well the planet is detected. Fig.   10 shows the resulting distributions for different levels of variability.
First of all, we see that the mean count in the subtracted aperture remains the same with and without the variability in Fig. 10. This indicates that the random variability would not bias the measurement of the planetary flux. Rendezvous with Roman. The demonstrated capability of formation flying should be able to control the starshade to stay well within this range during science operation. 11 The allowed range of variability (i.e., half an order of magnitude) appears to be also greater than the variability of earthshine. 62 Therefore, the variation in the solar glint and residual starlight caused by the starshade's motion in formation flying, as well as the variation in the brightness of the earthshine, should not result in substantial degradation of the S/N of planet measurements.

Sensitivity of Exozodi Levels
We have assumed the exozodi level to be 3 zodis in this work, guided by the best-fit median value from the hitherto most sensitive exozodiacal dust survey. 23  Of particular interest is where nearly photon-limited planetary spectroscopy may be expected with good background calibration. If the exozodi level increases from 3 zodis to 10 zodis, the boundary between the planet-dominated regime and the exozodi-dominated regime would increase from ∼ 1 R ⊕ to ∼ 1.8 R ⊕ for Roman observing the nearest (< 4 parsecs) stars (Fig. 5). The boundary would increase to > 3 R ⊕ for farther stars. This means that an increase of the exozodi level could prevent Roman Rendezvous to perform photon-limited planetary spectroscopy for any target stars. HabEx on the other hand is less susceptible to high exozodi levels. If the exozodi level increases from 3 zodis to 10 zodis, HabEx would still be able to obtain photon-limited planetary spectroscopy for Earth-sized or larger rocky planets (< 1.7 R ⊕ ) around stars < 6 parsecs away (Fig. 8). The density and structure of the exozodiacal disks of nearby stars will continue to be an important area of exploration.

Implications on Mission Designs
The analyses presented in this paper assume that the observations of the planets would take place at the orbital phase angle of π/3. In reality, we do not know planets and their orbital elements we suggest that mission designs should set a high priority to maximize the search completeness of planets in the habitable zones of these four stars, and particularly Indi A, and encourage precursor efforts such as radial-velocity measurements that prioritize these stars.
We have used the reference wavelength of 700 nm for the analyses presented in this work, and it provides a representation of the "green" band of Starshade Rendezvous with Roman (615-800 nm 64 ) and the UV-Visible band of HabEx (0.3-1.0 µm 7 ). Roman may also have imaging and spectroscopy capabilities in the "blue" band (425-552 nm 6 ). The blue band offers a few advantages over the green band, with less bright solar glint by 2 magnitudes 10 and a smaller PSF and thus less exozodiacal light by a factor of ∼ 2. These advantages would reduce the needed integration time to achieve S/N and loosen the requirements on background calibration. The blue band however requires a larger separation between the starshade and the telescope (37.2 Mm compared to 25.7 Mm for the green band), leading to longer distances and higher fuel and time cost to maneuver the starshade from target to target. Also, the blue band likely has fewer spectral features of interest from exoplanet atmospheres than the green band. [42][43][44][45][46][47] In addition, HabEx would perform nearinfrared spectroscopy (1.0-1.8 µm) on selected "high-interest" targets. 7 The magnitude of the solar glint has not been evaluated for the HabEx starshade in the near-infrared band, nor has the optical edge coating been designed to cover the near-infrared band. Future studies are required to quantify the solar glint and its impact on the science performance of HabEx in the 1.0-1.8 µm band.
Lastly, the planet search space that would be fully accessible for Starshade Rendezvous with Roman would be large rocky planets (∼ 1.5 R ⊕ ) and potentially water worlds (∼ 2.5 R ⊕ ). For all stars in the target lists, these planets would likely be within the reach from the perspectives of integration time and background calibration. The search completeness has not been assessed for these "super-Earths", but we suspect that they would have reasonable completeness. Therefore, detailed designs of Starshade Rendezvous with Roman would need to prioritize searching and characterizing large rocky planets and water worlds in the habitable zones, and relevant science investigations would need to further quantify the requirements for their orbital determination and atmospheric abundance retrieval.

CONCLUSION AND PROSPECTS
We provide an overview and reassessment of the noise budget of exoplanet imaging and spec-   Table 1.

5
The same as Fig. 4, except that the brightness of solar glint is lower than the nominal value by a factor of 10.
6 Improvement in the required integration time to characterize (S/N=20 in a 10-nm wavelength channel) planets around a nearby star, from a 10-fold reduction of the solar glint brightness from the nominal value. This is essentially a zoom-in view of Fig. 5 for sigma Draconis. The red dashed line corresponds to the IWA, and the blue dashed line corresponds to the planet-star separation for receiving Earth's insolation, i.e., the 1-AU equivalent.

7
The same as Fig. 4, except that the brightness of solar glint is lower than the nominal value by a factor of 10 and the residual starlight contrast is 3 × 10 −9 (i.e., higher than the nominal value by a factor of 30).