Space Telescope Imaging Spectrograph (STIS) is a second generation instrument for Hubble Space Telescope (HST) that was installed in early 1997 (Ref. 1) during Servicing Mission Two. It possesses a comprehensive collection of imaging and spectroscopic modes that span from the far ultraviolet to the near-infrared (IR). It is currently the oldest operating instrument on board HST that primarily obtains science observations. It retains a significant proportion of science time on HST, with 19% of all general observer (GO) orbits during cycles 23 to 25.
The high-contrast imaging mode on STIS includes a focal plane mask designated 50CORON. The focal plane mask consists of two wedges (A and B) joined perpendicularly, along with a wide occulting bar at the top of the field of view (FOV), as projected on-sky in HST’s science aperture instrument frame,2 and a narrow occulting finger on the right edge of the FOV (see Fig. 1). On the one hand, positions on the wedges are designated according to the width of the wedge at that location: for example, WEDGEA1.0 corresponds to the location on the vertical wedge with a width of 1.0″. On the other hand, the other apertures have designed lengths in the names: a wide occulting bar named BAR10 is 10″ long, and an occulting finger denoted BAR5 has a length of 5″. BAR5 was bent during the assembly of STIS and was only recently commissioned for high-contrast imaging via the outsourced calibration Program 12923 (PI: A. Gáspár).3 It is now a fully supported aperture location for STIS. The occulters are combined with a hard-edged Lyot stop in the pupil plane that is present for all imaging modes of STIS. The combination of the occulters and Lyot stop creates a simple coronagraph that marginally suppresses the diffraction pattern of the telescope.4
In this paper, we investigate the limits to STIS’s performance with respect to total-intensity, visible-light, high-contrast imaging. We determine what sets its limits on contrast performance at small inner working angles; in addition, we can assess its implications for what might be possible with second-generation high-contrast imagers in space, e.g., the coronagraphic instrument (CGI) proposed for the Wide Field Infrared Survey Telescope (WFIRST) mission.6,7 In Sec. 2, we talk about general strategies to obtain high contrast with STIS, building upon previous results;8,9,10 therein we also include the development of a coronagraphic noise model, which can be used to predict the performance of a given observing strategy; together with an analysis of the impact that charge transfer inefficiency has on high-contrast imaging with CCDs. In Sec. 3, we present a calibration program recently obtained with BAR5 to demonstrate high-contrast capabilities with STIS that rival current ground-based near-IR-optimized coronagraphs, which are equipped with extreme adaptive optics (AO) systems, as well as a discussion of STIS’s sensitivity to point sources and circumstellar disks. In Sec. 4, we describe a handful of common science cases that may require STIS’s coronagraphic aperture and the considerations required. In Sec. 5, we apply our new STIS coronagraphic noise model to investigate the potential of WFIRST to conduct high-contrast imaging beyond its dark hole, and in Sec. 6, we list the areas in which STIS is unique compared to other high-contrast imagers. We finish with our conclusions in Sec. 7.
High-Contrast Imaging with STIS
In their simplest form, high-contrast imaging techniques seek to minimize the diffracted light from a brighter central astrophysical object, allowing fainter point sources or extended sources to be detected at small angular separations from the host object.
High-contrast imaging with STIS requires two major steps to accomplish this task. The first is to obtain the image of a target, which is achieved by placing the target behind the occulters—this prevents light from the core of the stellar point spread function (PSF) from saturating and bleeding on the CCD detector. The raw contrast in such images is typically not sufficient for most observations since a significant halo of the PSF (also called the wings of the PSF) is present. The second step is to subtract off an astrometrically coregistered and intensity scaled reference PSF to enhance contrast. A final optional step is to obtain additional images of the target at different spacecraft orientations to obtain information about the sky scene that is occulted by the wedges, bars, and diffraction spikes of the target.
High-contrast imaging programs with STIS in the first few years of its operation primarily focused on a strategy of obtaining images of a target at two separate spacecraft orientations and with a single observation of a reference star.4 This strategy minimizes the total number of HST orbits required to obtain a final sky image (3), but at the expense of angular coverage, inner working angle, and exposure depth at large angles. In this observing scheme, the reference star for the target must also be very closely matched in spectral type (or spectral energy distribution, SED): the ideal target and PSF reference star (a.k.a. PSF-template star) should have a difference in or colors . These early programs found that the contrast achieved with a given aperture location is the same–thus paving the way for combining aperture positions to maximize angular coverage or sensitivity. The 50CORON aperture mode is operated without any filters and the effective bandpass of the aperture is defined by the red and blue wing sensitivity cut-offs of the CCD detector. The broad effective bandpass thus creates different PSF wings as a function of the host star’s SED. Other high-contrast techniques were also investigated, such as coronagraphic spectroscopy with the STIS long slit gratings,11 and placing bright sources outside of the FoV to attempt deep imaging of nearby stars, such as Eri.12
With the refurbishment of STIS in 2009 during the fourth HST servicing mission (SM4),13,14 STIS became the sole space-based coronagraph in operation on HST since ACS/HRC was not repaired during the mission. Since SM4, multiple programs have used the STIS 50CORON aperture to obtain high contrast at large separations to directly image and characterize the orbit of Fomalhaut b, to directly image at high signal-to-noise debris disks previously unobserved in the visible, and to find very faint disks.10,15,16 Many of these programs adopted observations at multiple spacecraft orientations (known as multiroll differential imaging, MRDI), and some used the central star itself as a reference, akin to ground-based angular differential imaging techniques.17,18
Concurrently, postprocessing techniques that suppress speckle noise, such as the locally optimized combination of images (LOCI) algorithm,19 the Karheunen–Loève image projection (KLIP) algorithm,20 and non-negative matrix factorization21 have been developed and shown to be successful at pushing contrasts closer to the photon noise limit of a given coronagraph. These techniques have been shown to be effective with archival NICMOS images22–25 from its calibrated archive26 and STIS27 coronagraphic images.
Strategies that Maximize Angular Coverage
There exist several strategies for obtaining full 360-deg circum-azimuthal field coverage of an object with STIS, but given the multiple options, we provide a subset of possibilities to give observers a flavor of what is possible. We detail some common examples in Sec. 4. Wedge and BAR position selection is primarily dictated by the desired inner working angle and the brightness of the star. However, due to the occulting mask obscurations and the presence of unapodized stellar diffraction spikes in the image scene, a significant portion of the stellocentric field is unimaged in a single-spacecraft orientation. This can be mitigated using multiple spacecraft orientations and/or by choosing a combination of aperture positions; see Ref. 28 for the observation of the Pic disk with three roll angles and two aperture positions on each wedge.
In another example, Fig. 2(a) shows the coverage in an image with a hypothetical HST program that utilizes a sequence of six relative spacecraft orientations (, 0 deg, , , , ) in six separate orbits to maximize angular coverage with an inner working angle of 0.45″ from a combination of shallow images at WEDGEA0.6 and deeper images at WEDGEA1.0. This strategy, including two additional orbits for a reference PSF star, has been used to characterize debris disks around a variety of stars.10 One scheduling consideration for this approach is the fact that allowed spacecraft orientations are different through the year and restricted—thus limiting opportunities to realize six significantly different orientations within a single cycle or year. Often, a gap in observations is required between two sets of orientations. This gap can be 3 to 6 months, depending on the location of the target on the sky. Figure 3 demonstrates a roll angle report of HD 38393 generated with version 26.1 of the Astronomer’s Proposal Tool (APT), which can aid in scheduling specific orientations on the sky for a given scene and target since it lists all possible scheduling opportunities without regard to additional constraints caused by other programs or the availability of HST guide stars.
An approach that obtains a similar inner working angle with fewer orbits is to combine the BAR5 and WEDGEA1.0 aperture locations as demonstrated in recent images of TW Hya as part of Program 13753. In this case, nearly 360-deg coverage down to 0.5″ was achieved with three spacecraft orientations (0, ), and partial angular coverage down to an inner working angle of 0.25″.29 The coverage map for three visits of that program is shown in Fig. 2(c). For edge-on circumstellar disks, a combination of BAR5 and WEDGEB1.0 may be a preferable choice to obtain shallow images of circumstellar material close to the star and deeper images further from the star with a minimum of spacecraft orientations since these two aperture locations are parallel to each other.
Another approach that achieves nearly 360-deg coverage with only two orientations at small inner working angles with short exposure times is to combine BAR5 with lower left and lower right locations of BAR10 (Fig. 4). Although BAR5 is a fully supported aperture position with STIS, which can be selected in APT, the lower left and lower right positions of BAR10 can be utilized by executing a POS-TARG, or small angle maneuver from the BAR10 central position. The BAR10 positions have been demonstrated to be moderately effective3 but have not yet been combined with BAR5 for any existing programs. However, by combining the three positions, one could obtain in principle 360-deg coverage to an inner working angle of 0.2″, as demonstrated in Fig. 2. One potential risk to this approach is that flatfield images of the 50CORON aperture indicate up to a 3% decrease in sensitivity close to all wedge positions, but it is particularly noticeable near BAR10, where such sensitivity decreases can extend up to 10 to 20 pixels from the BAR10 edges. Currently, these sensitivity deviations are not corrected within the calSTIS pipeline.30 Users are encouraged to download coronagraphic flatfield images from the archive if their scientific program requires a high degree of photometric repeatability. In most cases, PSF wing intensity nonrepeatability likely will dominate over these effects.
At least one program has also experimented with placing targets via a POS-TARG at the edge of the BAR5 aperture (i.e., Program 13725, PI: P. Kalas). This will trade angular coverage in a single exposure for currently unknown performance of the edge of BAR5—since the occulting finger is bent, it is currently unclear what impact that may have on photometric accuracy or repeatability of contrast, compared to the nominal BAR5 position that has been investigated.
Strategies that Reach Highest Contrast
Several strategies can help reach higher contrast achieved with STIS coronagraphy. These include strategies for selecting the best reference stars and careful selection of exposure times to maximize the total counts in the PSF without driving the region of interest into saturation. Additionally, we review here the sources of noise that are present in coronagraphic observations, including noise terms that are generally not accounted for within the traditional signal-to-noise ratio (SNR) calculations that are performed for the direct imaging of sources.
If the science program requires a reference star, matching the color of the reference star across the STIS bandpass is critical for ensuring the best contrast results. Details of recommended color matches are given in Ref. 9, but we review here—typically colors must be within 0.05, with closer color matches and spectral type matches working the best. Great care should be taken in considering the SED of the target and reference over the full bandpass of STIS, which runs from to (Sec. 2), suggesting that better stellar matches can be made if one considers the full SED of an object in that range, or at the very least between and if more complete photometry is not available. This is especially true with very blue or very red stars, where the effective bandpass of the PSF is skewed to wavelengths shorter or longer than . Interestingly, the broad sensitivity also means that STIS is more sensitive to substellar objects than a typical visible light imager since it retains sensitivity to photons from 0.8 to , where and dwarfs still retain significant flux.
In the case of reference star selection, an ideal candidate is reasonably close in the sky to the target, usually within 10 deg. This minimizes the thermal changes in the observatory between the target star and the reference when observed in contiguous orbits without other intervening telescope repointings. If multiple orbits on the target are executed, it is best to have the reference star in either the second or third orbit in the sequence since most of the thermal settling of the telescope can occur within the first orbit after a slew. Large thermal swings can drive defocus from differential heating across the support structure, which is the cause of telescope “breathing.”
Table 1 lists observed count rates in the PSF wings of a selection of stellar spectral types in units of and using the BAR5 aperture location. The count rates were determined by taking azimuthally averaged surface brightness profiles of observed STIS PSF wings while avoiding diffraction spikes, measuring the median count rate at a given radius, and multiplying by the detector gain. Calculations of the count rate in the peak pixel come from the STIS exposure time calculator (ETC).31 Note that the maximum count rate close to the mask edge can vary from this median value by as much as 20%, users should be aware of this if they do not wish to saturate directly at the mask edge, especially if they are executing small dithers behind a mask aperture position. Table 1 provides a handy rule of thumb for the maximally allowed exposure time for a given BAR or WEDGE position assuming a target of 80% full well and a saturation limit of for gain . Previous observations with STIS show that the actual full well varies somewhat across the detector32 and is thus a consideration when doing a more careful calculation for how long to expose a target. The 20% margin provides safety for subpixel dithers behind the aperture position and/or miscentering of the target.
The observed count rates as a function of radius for the stellar halo of different spectral types, avoiding diffraction spikes and occulters. Sx are the count rates in the PSF wings at a x″.
Another consideration is the selection of an exposure time for individual read outs of the CCD. The upper limit to the photon limited contrast possible in a single orbit is dictated by the total exposure time, and the efficiency with which a readout can be obtained. The typical readout time of the STIS CCD is roughly 1 min but can be as small as 20 s if subarrays with a small number of rows are used. Overheads thus limit the total counts in the PSF per orbit under the restriction of not saturating for a given mask position. If we assume that we design a program that has counts close to the full well at the edge of the BAR5 mask position in each exposure (), we can calculate the total integration time available for observing a star within a typical 50-min HST orbital visibility window (), accounting for readout and data management overheads and assuming a subarray with a readout time of 20 s. We select this particular subarray since it represents a compromise between quick read-out time and spatial coverage. The total exposure times and counts in the PSF per orbit are presented in Table 2. Total exposure times were calculated using version 26.1 of the APT in the phase II proposal mode, and the count rates were scaled from Table 1. Users should adopt these values as a rough rule of thumb and do their own more detailed calculations for their specific program requirements.
Effective number of counts as a function of radius within a typical Hubble orbital visibilty of 50 min and assuming typical overheads resulting in a total integration time Tint. We assume observations are taken in gain G=4 mode. Observers should determine the specific exposure times obtained for a given target and its visibility window.
Table 2 demonstrates that STIS will obtain similar number of total counts within the PSF wings (to within 50%) within 1 orbit and with the gain () of the CCD set to 4 for stars with . is in general preferred for high-contrast imaging due to the fact that more counts can be collected within a pixel before it saturates. The range of stellar magnitudes is not dissimilar to the range of magnitudes that are effective with high-order AO coronagraphs such as VLT/SPHERE33 or the Gemini Planet Imager.34 It should be noted that for very bright stars high-contrast imaging is incredibly inefficient due to the overhead associated with subarray readout times and buffer dumps—for a 50-min visibility window only 19.2 s in total is spent on source for a . One significant difference to the ground is that STIS’s Strehl ratio does not degrade for fainter sources. Thus for objects with , STIS’s contrast performance is dictated by the total exposure time available to expose the wings of the PSF relative to the background limited sensitivity of the exposures.
Calculating STIS Raw Coronagraphic Performance (without PSF-Subtraction)
Currently, the STIS ETC does not include a straightforward way to estimate the raw performance of STIS coronagraphic modes for a given exposure time and target star. For STIS imaging modes, the ETC does include noise sources such as solar system zodiacal light, CCD dark current (which has slowly been growing with time due to its exposure to the harsh radiation of space), and CCD read noise. Light from the occulted star’s PSF wings, however, is by far the largest contributor in particular at small stellocentric angles, but it is currently not considered in the ETC. Further, stochastic noise from quasistatic aberrations and jitter are currently not considered. A simple analytical case for noise sources can be assumed following the methodology of a typical background limited observation. We also use information provided in the Hubble Space Telescope ETC User Manual,35 we first define the noise coming from the STIS CCD detector :
We then define the important background noise sources for space observatories, including light from the Solar System zodiacal dust disk () and scattered light from Earthshine (), also in units of :
We assume the expected detector noise properties of STIS projected to mid-Cycle 26 for the middle of the detector and with gain . In practice, Earthshine and zodiacal light are much smaller than other noise sources.
For high-contrast imaging, one also needs to take into account the photon noise from the PSF wing intensity (), as well as the “speckles” from both static aberrations as well as slowly varying wavefront errors (). Static speckles can typically be removed via subtraction of a reference star. Slowly varying speckles correlate with the thermal and focus state of the telescope, the so-called telescope breathing. For STIS, because of its wide bandpass, the speckles are actually spoke-like in shape and extend radially (see Fig. 8 of Ref. 4). The average observed intensity is comprised of the static, aberrated coronagraphic PSF intensity as well as speckles from high-frequency telescope jitter and speckles from quasistatic errors .8,36 The photon noise from these three contributions is given by
Note that can be determined from a log–log interpolation of Table 2, or by referring to the publicly available, azimuthally averaged PSF intensities provided on the STIS Instrument website37 for a more detailed profile. These intensities are calculated relative to the peak pixel of the STIS PSF and are per pixel intensities. The other intensities are assumed to be much less than and are ignored.
We finally construct the total noise and put it into units of analog to digital units (i.e., the counts within a STIS calibrated exposure obtained from the MAST archive) via multiplication of the gain () of the CCD, which is typically either 1 or 4. We also account for reference star subtraction with the intensity scale factor , which corresponds to () where is the ratio of the flux from the target and reference star observations, and , which corresponds to :
Estimating Systematic Speckle Behavior in STIS High-Contrast Images
It is important to note that has not formally been studied with STIS. Ground-based AO systems typically have speckle intensities that follow a modified Rician probability distribution function (PDF) or a Weibull PDF, rather than a Gaussian PDF of the mean intensity.17,36,38 Recent work has also shown that with extreme-AO fed systems, atmospheric decorrelations, and internal optics decorrelations grow exponentially within the first tens of seconds and then grow linearly on the timescale of minutes to hours.39 Atmospheric driven speckles have intensity distributions as a function of angular separation that mimic a Moffat profile halo with a full-width half-max that is roughly proportional to the characteristic length of the atmospheric turbulence.40 We expect that the majority of speckles in STIS will be generated by a combination of jitter, which operates on short timescales, and focus changes, which operate on orbital visibility timescales ().
In this case,8,36 the noise associated with the speckles can be directly related to the speckle intensities of each type as well as their cross terms over the total number of independent-speckle realizations, :
Directly quanitfying each term in Eq. (5) requires careful analysis and assumes that these terms have well-defined timescales. Since most coronagraphic exposures are taken with of a few seconds to a few hundred seconds, one can estimate by calculating noise in each pixel and measuring how much it exceeds that expected from the PSF wings alone. Such an analysis by definition ignores the possible spatial correlations due to the spoke-like shape of STIS speckles. Accounting for spatial correlations is beyond the scope of this paper, but we assume that one cannot average over pixels in a small photometric aperture due to this correlation.
We investigated the measured variance of speckles in several STIS coronagraphic images using a combination of the BAR5, WEDGEA0.6, and WEDGEA1.0 aperture locations. In particular, we studied the intensity root-mean-square (RMS) in the HD 38393 () observations with BAR5 from Program 14426 (PI: J. Debes), and HD 30447 (), HD 141943 (), HD 35841 (), and HD 191089 () from Program 13381 (PI: M. Perrin); therefore, our observations sample the state of the telescope between 2013 and 2015. For the Program 13381 targets, we investigated the noise in both WEDGEA1.0 and WEDGEA0.6 aperture locations (see Table 3). Since each of these targets had multiple exposures of the CCD within a single orbit in at least two orbits, we first calculated an RMS noise per pixel as a function of time for the images of one star. We azimuthally averaged the resulting 2-D RMS over time to construct a radial distribution of RMS. To determine the contribution from speckle noise, we made a few assumptions. First, we assumed that the majority of the noise observed is due to deterministic photon statistics as estimated from the STIS CCD detector parameters, the brightness of the star, and the observed average PSF wing intensity with static aberrations as determined by Eq. (3) and assuming that the speckle noise is 0. At large angular distances for short-exposure times, we see excellent agreement with these predictions. We also assumed that the speckle noise in one detector pixel is given by , where is the flux from the central astrophysical object in the peak pixel in units of , and is the speckle RMS dependence as a function of the angular separation to the central object. This functional form for the speckle variance as function of angular distance from the star is roughly equivalent to the normalized speckle variance discussed in Eq. (36) of Ref. 36. This means that the variance we measured should scale both with exposure time, and with the brightness of the central object.
Stars used for calculation systematic speckle noise at BAR5, WEDGEA0.6, and WEDGEA1.0.
Taking our RMS measures per star, we calculated the residual noise above the predicted photon and detector noise, especially in regions where the noise was much larger than the predicted variance. We then scaled the residual noise by both the exposure time (which ranged from 0.2 to 540 s) and the expected flux in electrons for each star as calculated by Psynphot for the relevant apparent magnitude and spectral type of each target (the stellar spectral types were F3–G2). Figure 5 demonstrates our resulting empirical measure of for our targeted stars, with the exception of HD 141943. In the case of HD 141943, the observed noise closely matched the predicted noise from just the PSF wings and the detector, precluding an accurate fit of additional systematic residuals and implying, for those observations, the speckle noise was a factor of two lower than that compared to the other observations. We thus expect our results for the majority of the stars are “typical” and that, while some variation exists due to the state of the telescope, the variability does not exceed a factor of two for nominal telescope operations for that time period. More detailed estimates on the distribution of speckle noise levels and what they might depend on as a function of the STIS lifetime is left for future work. We recommend that anyone using this approach to predict performance to take coronagraphic observations from the MAST archive as close in time as possible to verify the latest speckle behavior of STIS and to build some margin, if possible, into their observational design. Figure 5 also demonstrates that our assumption on the dependence on stellar flux and exposure time were correct over a fairly typical range of stellar brightnesses and exposure times relevant for typical STIS observations, and that the primary dependence on angular separation, can be fit empirically. By averaging together all the existing residual noise curves, we fit a function to the angular dependence of the speckle RMS with a form of , where , , and is in units of milliarcseconds. The constant is equal to the platescale of STIS, 50.7 mas. The uncertainty of each noise measure was typically between 20% and 50% across the four stars. The dependence on angular radius is similar to the dependence of sensitivity limits derived in deep STIS coronagraphic images of HD 207129.41 Putting this together, for our noise model would become:
This encapsulates our assumptions about how the speckle variance scales with exposure time and implicitly assumes that speckles are likely spatially correlated on small angular scales. This is appropriate given the fact that a single STIS CCD pixel is close to , where is the central wavelength of the CCD bandpass and multiple pixels will encompass multiple speckles. Additionally, since the sensitivity of the CCD spans from 2000 to 10,000 Å, individual speckles can span a large number of STIS pixels in the radial direction. As an example, if an aberration caused a speckle at a radius of , it would span from 68 to 343 mas, or the equivalent of 1.3–5 pixels in radial extent.
We now include Eq. (6) in Eq. (3) to obtain the final estimate of the noise. In the next section, we compare this to a suite of observations from Program 14426.
The physical origin of this variance profile, which roughly tracks the power-law decrease of the STIS PSF wings is likely caused by speckle pinning.36 Since the STIS wedges and bars do not optimally suppress the diffraction pattern, STIS coronagraphic images approximate direct imaging at a high Strehl ratio.
We note that a similar procedure could be implemented during the James Webb Space Telescope (JWST) commissioning in order to quantify speckle statistics for that observatory, as well as for other future high-contrast imaging missions. A dedicated calibration program would image sources with increasing exposure times to characterize the raw uncorrected speckle intensity field under differing wavefront control algorithms as well as with reference star subtraction. The performance could then be monitored over the mission lifetime to aid in providing planning information to GOs. It is often a significant stumbling block to successful high-contrast imaging, if the noise statistics of speckles are not properly accounted for.
An additional subtlety of the speckle noise is in optimizing an observing strategy when one is attempting to use a reference star to obtain high contrast. Using a reference star with a different brightness and exposure time can drastically increase residual speckle noise because it will then be dominated by the speckles for one of the stars. The residual halo described in STIS observations of HR 879942 is due in part to the very different exposure times and apparent magnitudes of HR 8799 and its selected reference star.
Validation of the Coronagraphic Noise Model
We validate our noise model by comparing the per pixel RMS of 24 exposures of HR 8799 taken as part of Program 12281 (PI: M. Clampin) to predictions. We calculate the RMS over time per pixel as described in Sec. 2.4. The exposure time for each CR-SPLIT was 20 s, with 8 CR-SPLITs per exposure. We compared the azimuthally averaged RMS as a function of time by excluding regions containing the WEDGEA1.0 mask and HST diffraction spikes to the analytical noise formulation given in Eqs. (1)–(6). Figure 6 shows the results, along with the individual components of the noise. At this exposure time, the dark and sky backgrounds are negligible compared to the noise in the PSF wings and due to speckles. We selected this dataset because it was taken close in time to the stars for which we measured the speckle noise directly. Looking at the ratio of the radial noise profile relative to the predicted profile, we determine that the median observed noise was 8% higher than predicted between 0.5″ and 3.0″. Close to the star, we overpredicted the expected speckle noise by 40%. It is important to note that this noise per pixel is not sufficient to predict the contrast sensitivity—one will need a defined photometric aperture as well as an accounting of the degradation of contrast at the inner edge due to incomplete azimuthal coverage due to the occulters and the residual diffraction spikes of the target.
For observation planning, one will also need an estimate of the target source flux and its accompanying photon noise to conduct a full calculation to determine an optimum exposure time. We investigate these complexities further for point sources in Secs. 3.4 and 3.6 for extended sources such as circumstellar disks.
Effects of Charge Transfer Efficiency Degradation and Potential Impacts to Science
Due to its 21-year lifetime in low Earth orbit since 1997, the STIS CCD has suffered from cumulative radiation damage. The damage primarily takes the form of increasing numbers of hot pixels on the detector, increasing dark rate, and increasing numbers of defects in the silicon substrate that “trap” electrons generated from a photon striking the detector. These trapped electrons are then released later than expected on the detector during readout, resulting in an artificial trailing of light away from the readout direction and a degradation in charge transfer efficiency (CTE).43 In the case of a point source on STIS, this will appear as a trail of light that points to the bottom of the image. For an extended source, charge traps will create a low-level halo that would be asymmetrically stronger FOR a few pixels below the source.
CTE degradation in high-contrast imaging is a phenomenon that a PSF halo observed with different exposure times will have a subtly different two-dimensional surface brightness distribution. The impact of degrading CTE on high-contrast imaging is of particular interest to the design and operation of the CGI within NASA’s proposed WFIRST mission.44 Although CGI will likely utilize a different type of detector than STIS, the requirements on contrast will be much more stringent.
We investigate the upper limits to CTE degradation impacts with STIS by investigating the behavior of images of the same source taken within the same orbit but with differing exposure times. One example of this is observations of HR 8799 taken in Program 12281. Both 20- and 60-s exposures were obtained with the WEDGEA1.0 aperture position of HR 8799, resulting in a factor of 3 difference in the total counts on the CCD from the PSF wings. Additionally, the program selected a PSF reference star, HIP 117990, which is a factor of 2 fainter than HR 8799. For the PSF exposures, each readout had an exposure time of 318, resulting in a factor of 2.6 and 7.8 more counts on the detector than for the 60- and 20-s HR8799 exposures, respectively.
We constructed images of four cases: a mean combination of 20-s exposures subtracted by 20-s exposures from an adjacent orbit (case A), a combination of the PSF images subtracted from the 60-s exposures (case B), a combination of the 60-s exposures subtracted from the 20-s exposures (case C), and a combination of the PSF images subtracted from the 20-s exposures (case D). In this way, we investigate the impact of differing total counts on the detector between reference PSFs with increasing scale factors of 1, 2.6, 3, and 7.8 relative to the target for the four cases. From case A to case D, one would expect that the PSF reference star would be less impacted by CTE degradation and one would expect positive residuals preferentially pointing toward the bottom of the detector.
Figure 7 shows the resulting subtraction images for each case, demonstrating the presence of increasing residual surface brightness at small stellocentric angular radii as the multiplicative factor of the counts on the detector for the PSF relative to the target increases. Although we note an increasing residual halo in the images, this is not an entirely clean experiment since there could be color and focus differences between HR 8799 and its reference [although in this case the color difference between the two objects is relatively small ]. There appear to be vertical asymmetries in the residuals that are most noticeable along the diffraction spike subtraction residuals and near the bottom of the PSF. At an angular radius 1″, the median surface brightness of the residuals is 2, 6, and 9 counts for cases B, C, and D, respectively compared to estimated to be in the peak pixel of HR 8799. These numbers correspond to residual emission per pixel from CTE degradation to be no more than at a contrast of between 2 and . In the worst cases, where a specifically sharp PSF structure exists the emission can be a factor of two to four brighter (such as at the feature located at the 10 o’clock position in Fig. 7.
Based on these results, we recommend that users select PSF reference stars that are close in brightness or slightly brighter and to choose exposure times that ensure the same number of counts on the detector within the PSF wings to better than a factor of 2. This should limit the impact of CTE degradation to a surface brightness that is lower than relative to the peak flux of the star at 1″ as well as minimize noise from differing levels of speckle noise. These observations will also have the advantage of sampling the breathing within an orbit in a similar way to better match PSF images to the target and minimize subtraction residuals.
The impact of CTE for high-contrast imaging has been an area of focus for the CGI of the WFIRST mission. The CGI will be a technology demonstration of extreme high-contrast imaging at small inner working angles that will likely use a photon-counting EMCCD, which will still be vulnerable to CTE degradation due to radiation damage. More work will need to be done to understand whether slight differences in counts on an EMCCD detector could result in similar behavior to the upper limits we place above. Since we do not significantly detect any degradation after more than 20 years in space, it is likely that the largest impact to CGI will be in the loss of signal from faint sources rather than contrast degradation.
High-Contrast Imaging at Small Inner Working Angles with BAR5
The BAR5 aperture location on the 50CORON aperture has a physical half width of 0.15″, and it has been demonstrated to have an effective inner working angle of 0.2″. As the aperture with the smallest inner working angle on 50CORON, BAR5 has achieved reliable high-contrast results for bright circumstellar disks.3 Program 14426 was designed to achieve a working contrast at the BAR5 inner working angle with a goal of obtaining contrasts of or better, comparable in visible-light to near-IR contrasts achieved on the ground.
HD 38393 ( Lep A) is a nearby 1.3 Gyr old45 F6V46 star at 47 in the Ursa Major moving group,48 with a common proper motion K dwarf companion ( Lep B) at a projected separation of 95″ (855 au).49,50 It was chosen due to its availability for scheduling, its brightness, and its relative dearth of known companions or circumstellar material. HD 38393 was previously directly imaged with the MMT/CLIO Instrument.51,52 The tabulated contrast sensitivity curve53 for HD 38393 shows that substellar objects more massive than are ruled out for this star beyond au. HD 38393 has not been investigated for radial velocity variations.54 Due to its proximity to Earth, HD 38393 has been the target of Spitzer photometric and Keck nulling interferometer searches for dust, showing a lack of any dust around the system with limits to relative IR luminosities for cold dust and for warmer dust, respectively.55,56
Program 14426 was designed to observe the star HD 38393 over the course of nine orbits.57 Each orbit was at a different spacecraft orientation, and the orbits were arranged in three sets of three continuous orbits separated by roughly one or two months. The first set of observations were executed on October 21, 2015, with the final set of visits executing on February 26, 2016. The final visit was also used to perform a test of the updated Science Instrument Aperture File2 position of BAR5 to validate the supported aperture position in APT starting with version 24.1.
HD 38393 was imaged with subpixel dithers in a grid behind the BAR5 aperture: on the one hand, small miscenterings on the order of 0.25 pixel () in BAR5 can cause significant differences in the amount of light leaking around the edge of the BAR5 mask, and subpixel dithering provides margins against mask and target misregistration with imperfect initial acquisition, with then improved PSF rejection;3 on the other hand, it allows for an operational test of a similar dithering strategy proposed for the JWST coronagraphs.58 For classical PSF subtraction, subpixel dithers can allow better matching of the images between a target star and reference star; for the statistical postprocessing techniques (e.g., LOCI, KLIP, and NMF), it provides additional diversity to the PSF library. This strategy has been suggested to mitigate the tendency of coronagraphic masks on the MIRI and NIRCam instruments to degrade in performance with the expected target acquisition accuracy of JWST.58,59 Simulations of the MIRI four quadrant phase mask contrasts with subpixel grid dithers showed up to a factor of improvement over raw contrast58 and a factor of improvement for NIRCam coronagraphic modes.59 Being able to test this observing strategy on-sky before the launch of JWST provides a useful test of this approach.
Subpixel dithering relies highly on the accurate centroiding of all the images. For STIS, this is straightforward. Because the masks and Lyot stop in the instrument are not optimized, significant diffraction spikes from the telescope support structure remain in coronagraphic images. Although this is not ideal for imaging purposes, it does serve as a unique external probe of the stellar centroid. The diffraction spikes are nearly at 45-deg angles from the vertical direction on the detector. Therefore, it is possible to measure the vertical location of the spike on the detector as a function of detector column and linearly fit the spikes to infer a stellar centroid. If one subtracts off the PSF halo and uses a mask that encompasses only the diffraction spikes, accurate centroids can be obtained with this method, typically with uncertainties of better than 0.05 pixel or 2.5 mas. Another centroiding method is the Radon transform-based line integral.60 For STIS, one can apply an angular separation-dependent correction map to the spikes, which can weigh the contribution from different pixels, then perform the line integral along 45 deg and 135 deg only27 and obtain similar results as the previous method.
High-Contrast Image Postprocessing for Point-Source Detection
A large number of STIS high-contrast observations have been obtained with classical high-contrast imaging methods that rely on the stability and repeatability of the stellar PSF over periods of several hours. Other postprocessing methods exist as well, and it is useful to consider the benefits of each type of technique in the context of STIS observing. In this section, we investigate constructing high-contrast images of HD 39383 using both classical PSF subtraction and KLIP postprocessing.20
Classical MRDI subtraction
For all the nine visits57 of HD 38393 with STIS, in each visit, the HD 38393 images from the same dither position were aligned with the “X marks the spot” method10 and registered to a common centroid, then median combined into a single image for the dither position. PSF subtraction for each dither position within each visit was done using adjacent orbits: for example, visit 4 used either the visit 5 or visit 6 image at the same dither point as a PSF reference. The better subtraction was determined by eye, and a mask for each dither point was constructed. All nine separate spacecraft orientations were derotated and combined. This procedure is the most reminiscent of MRDI and is particularly well suited for STIS. In the case of Program 14426, our use of the host star as the reference guarantees that PSF subtraction residuals due to potential chromatic mismatches in otherwise used template stars are fully eliminated. Otherwise, as discussed in Sec. 2, if one uses a reference star, contrast performance may be degraded depending on how much the reference star departs from the target star’s SED across the STIS bandpass.
Karheunen–Loève image projection subtraction
The KLIP algorithm is used to reach the highest contrast to detect point sources.20 To perform it onto the 810 total images (9 orientations dither positions CCD readouts), they were aligned by performing Radon transform-based line integral using centerRadon during a systematic postprocessing for the entire 50CORON archive in Ref. 21. To minimize the influence from the BAR5 occulter and residuals from the diffraction spikes, masks were created on subarrays of the data and applied. For each subarray image, we took the centerRadon-determined positions of the star in the detector coordinates, then interpolated and cut a mask of the same size from the whole FoV STIS mask created in Ref. 29.
A reference library of PSFs was generated for each spacecraft orientation (90 images) from the other spacecraft orientations (720 images) to determine the contrast achieved within 24 pixel (): exterior to this angular separation, point-sources can be detected to the photon limits. Following the KLIP algorithm, we generated model PSFs based on KLIP components of the STIS quasi-static PSF and subtracted these from the target image (Sec. 3.4).
Sensitivity to Point Sources
To test the STIS contrast sensitivity performance of our HD 38393 observations with classical MRDI and KLIP, we inserted artificial point sources using calculated TinyTim PSFs61,62 using a G0 spectral type SED. To quantify the contrast limit for point source at different positions in the STIS images and simulate real observations, only one point source was injected then reduced for SNR calculation. This process is then sequentially performed for point sources at different physical locations on a grid in the polar coordinates: an angular separation, i.e., radius, ranging from 0.2″ to 1.2″ (4 pixel to 24 pixel) with a step of 0.05″ (1 pixel), and an azimuthal angle from 0 deg to 360 deg with a step of 30 deg. We scaled the point sources relative to a source with the equivalent of a total aperture corrected band magnitude of 26.05 ( with STIS). To determine if a point source was recovered, we obtained aperture photometry with a square centered at its location and determined whether the source was detected at an through a binary search. For each on-sky location, the iterative binary search started with two contrast limits and : in each step, we first inject a planet signal at that location with a trial contrast that equals the average of current upper and lower limits, then we reduce the data and measure the SNR for that planet. If its SNR is , then the new upper limit is set to that trial contrast; otherwise, the new lower limit equals the trial contrast. We repeatedly update the upper or lower limits until when they converge, i.e., when the SNR for the trial contrast is within 0.01 from 5. To calculate the final contrast for a specific radius, we took the median of the contrasts obtained for all the azimuthal samples.
For KLIP, we also injected the interpolated TinyTim point source PSFs to simulate the dithered observations of one point source in all nine orientations as described at the beginning of this section. For the images from each orientation, the images from the other eight orientations were treated as references to construct the KLIP components; we selected the first 400 components to reach the flat plateau of the KLIP residual variance20 for a trade-off between maximizing quasistatic noise suppression and minimizing the over-fitting of random noise. For each injected point source, we used the 400 components, we modeled the 90 point-source-injected images in the other orientations and subtracted the model from the images. We then median-combined all the 810 residual images and calculated the SNR for the point source. For a planet, we calculate its SNR assuming the small sample statistics:63 we first obtain the signal by measuring the total counts of a 9 pixel box () centered at the planet location and subtract it by 9 times the median of the pixels that are within 2 to 4 pixel from the planet. We then obtain the standard deviation for a ring of pixels that have stellocentric separation that are within of that of the planet (excluding the region for the planet, consequently, a total of pixel) and multiply that number by to obtain the noise for the planet. The SNR is thus the signal divided by the noise.
To calculate the contrast sensitivity curve with KLIP, we adopt the above injection and reduction procedure for one point source and calculated the SNRs for the sources injected on the above-mentioned polar grid. We performed a binary search until an SNR of 5 is reached. For each radius, the contrasts at different azimuthal angles were median combined as in the MRDI scenario.
The contrast sensitivity curves of KLIP and MRDI compared to predicted noise limits are presented in Fig. 8. This figure shows the resulting point source sensitivities for the MRDI and KLIP reductions compared to the raw PSF wing intensity and the predicted noise model limited contrast :Figure 8 shows that the contrast sensitivity beyond 0.7″ matches that expected for noise from the PSF wings alone and that KLIP is very close to the theoretical limit under our assumptions of speckle statistics with the exception of right near the inner working angle.
Effect of Subpixel Dithering and Postprocessing on Contrast Performance
HD 38393 was observed in Program 14426 with a dithering pattern in one pixel during a visit, with a dithering step of (). To explore the improvement of point-source contrast with subpixel dithering, we generated contrast curves for four groups; for each telescope orientation:
For each of the four groups, there are a total of 81 images (), and the images from different dithering scenarios are matched to their dithering setups. We performed the identical KLIP point source contrast search as in Sec. 3.4 but only using 40 components to reflect the shrinking of the number of images. We first calculated the contrast curves for the four groups, then obtained the improvement factors over the nondithering case by dividing the contrasts of nondithered contrast by the dithered contrasts and present them in Fig. 9.
From the improvement curves, we can see that both vertical dither and horizontal dithering improve the contrast limit, and a (with a step size of , ) (This is the recommended pattern and step size in Ref. 3.) dithering strategy improves the performance by a factor of at small inner working angles. These results are similar to those expected for the same strategy with JWST/NIRCam coronagraphic modes.59
Sensitivity to Circumstellar Material
The observations of HD 38393 can also be used to gauge BAR5 performance in the detection of starlight scattered by circumstellar material.
As in Sec. 3.4, we will determine STIS’s sensitivity to spatially resolved disk-scattered starlight using artificial sources injected into the images and processed via KLIP. We assume a secure detection of disk-scattered starlight when the SNR in a given pixel is . This is acceptable when, by definition, a disk-signal is spread over at least several resolution elements (and typically much larger), so the SNR of the total disk-scattered starlight is (For the extended structure in this paper, we have ignored the spatially correlated noise. Positively correlated noise will increase the noise level, our disk surface brightness limit is, therefore, optimistic. To reduce such noise, we used 1-pixel wide ring models. See Ref. 64 for a proper treatment of correlated noise using forward modeling.). Thus we would require that any detectable disk-scattered starlight must have an in STIS pixels or an angular area of . In practice, we calculate the signal to be the element-wise median of the derotated reduced images and the noise to be the element-wise standard deviation of the derotated reduced images.
We injected artificial disk-scattered starlight in the form of a circular face-on ring with a width of 1 pixel (thus unresolved in the radial direction but avoids spatially correlated noise that impacts wider rings) and convolved it with a TinyTim point source STIS PSF, injected it to the 90 images of one orientation at HD 38393, then used the 720 images from the other eight orientations as the references and performed KLIP subtraction; this identical process was performed for all of the nine orientations, then derotated and median combined, which was used to calculate the average SNR in the ring. The actual observational design of the 14426 Program is not well suited for the detection of disks at low-inclination angles (i.e., nearly face-on), due to the lack of a PSF template and the use of MRDI. However, our approach efficiently determines STIS’s azimuthally averaged sensitivity to spatially resolved disk-scattered starlight assuming a more advantageous experimental design is used, such as the use of a large reference library of PSFs from stars that do not possess extended disks.
Assuming the scattering particles in the ring are spectrally neutral (i.e., gray scattering), we infer the scattered light sensitivities in units of -band surface brightness in Fig. 10. We show two possible sensitivity lines because the sensitivity is dependent on the convolution of the STIS PSF and the spatially resolved disk-scattered starlight. For radially unresolved disk-scattered starlight in a narrow ring-like disk (such as what we have modeled), the PSF washes the disk signal out by the enclosed energy fraction in the peak pixel. As the disk-scattered starlight approaches several resolution elements in total radial extent, the signal approaches the true surface brightness (the lower curve in Fig. 10). This brackets the range of possible disk morphologies. We also have plotted the expected surface brightness of hypothetical exozodi around Ceti (47) and Fomalhaut (65). The exozodi disks are assumed to have constant optical depth, with their radii proportional to . The hypothetical models are 300 times brighter than the Solar System zodical cloud at 1 au (i.e., ).66 As in Ref. 66, we equate this surface brightness to 1 “Zodi,” which is not to be confused with similar definitions that are based in the amount of IR radiation from the Solar System’s zodiacal cloud. We note that this result is for “predicted” performance with BAR5 and does not represent the actual sensitivity to scattered light for HD 38393 due to self-subtraction. That is beyond the scope of this paper.
Recommendations for Specific Scientific Cases
Navigating the myriad aperture locations and formulating an optimal combination of them can be a fraught and confusing process. We thus try to detail some recommendations for common scientific cases with respect to high-contrast imaging with the 50CORON aperture. Users are encouraged to consult with the HST helpdesk if their scientific case does not match those listed here. In general, if the scientific object of interest resides at , it is recommended to use the BAR5 location. If the target resides beyond , it is recommended to use the WEDGEA1.0 location as that has the most archival data available for postprocessing libraries. For very low-surface brightness objects beyond , it is recommended to select the longest possible exposure times while still not saturating the detector at and to select the WEDGEB2.5 position,27 which has the most archival observations of the wider wedge positions. We point out particular HST GO programs that have utilized similar approaches to those detailed here—all phase II programs are publicly available in MAST for users to study.
Point Source at Small Inner Working Angles
Let us assume that a young substellar object is the target of a STIS program and it has an SED similar to GL 229B (a T dwarf, which is available as a nonstellar object option in the HST ETC). Its separation is 0.3″, and it is in orbit around an M5V star. From the STIS imaging ETC, such an object would have a contrast of relative to its host in the 50CCD bandpass and thus would be detectable with STIS and the BAR5 position using an approach similar to Program 14426: ADI + postprocessing. The user would need to use the noise model to assess the total exposure time needed and would likely require at least three separate spacecraft orientations, taking care to select the ORIENT constraint in APT to place the companion perpendicular to BAR5’s long axis. Assuming the PA of the companion is 30 deg, ORIENTS of 75 deg and 255 deg would be optimal, based on the rule of thumb for the WEDGEB and BAR5 locations that one selects ORIENTS relative to a specific PA by adding 45 deg or 225 deg. At this angular distance, one would select relative orientations of to the nominal ORIENT to ensure that the companion was not significantly self-subtracted.
Disk at Moderate Inner Working Angles with a Known Position Angle
Let us assume that the science target is an edge-on disk with a constant surface brightness of beyond 1″ and its major axis has a PA of 125 deg. In this case, WEDGEA0.6 or WEDGEA1.0 would be appropriate depending on the brightness of the source and how quickly the PSF saturates at the inner working angle of each mask. In this case, one would select ORIENTS rotated 90 deg from the WEDGEB and BAR5 recommendation since the WEDGEA occulter is perpendicular to the other occulters’ long axis. For users seeking to understand bright material close to the inner working angle yet also to detect faint nebulosity far from the target, a combination of WEDGEA0.6 and WEDGEA1.0 would provide larger dynamic range. This is similar to the strategy laid out in Program 12228. In addition to a selection of aperture location and orientation, it is paramount to obtain observations of a reference star with a similar magnitude, similar SED, and relatively nearby in the sky to ensure a similar thermal evolution of the observatory. If three orientations were planned, a fourth orbit would be needed to take observations of the reference star concurrent with the science target.
Faint Point Source at Large Angular Distance
For this case, we assume a very faint object with a contrast of relative to a magnitude star at a separation of 9″, which is similar to Fomalhaut b, as studied by Program 12576 (PI: P. Kalas). In this case, sensitivity and contrast are paramount over inner working angle and so the bright star is placed behind the WEDGEB2.5 aperture and allowed to saturate well beyond the full well. Even so, exposure times are short and of order 67 s per CR-SPLIT. Since the object is far from the star, small changes in spacecraft orientation for MRDI ensure that self-subtraction is not an issue for the point source target. MRDI and postprocessing can be used to achieve the desired contrast level over 12 orbits.
In this case, we assume observations of a near face-on disk that ranges in surface brightness from at 0.3″ to at 4″, where full angular coverage and photometric repeatability is desired. This is similar to GO 13753 (PI: Debes), which observed TW Hya with a combination of BAR5 and WEDGEA1.0. In this case, ORIENT constraints are not needed, but a careful scheduling of multiple orbits with ORIENT offsets is essential. In this case, the combination of two aperture positions that are perpendicular to each other affords nearly full azimuthal coverage with just three spacecraft orientations close to the inner working angle (see Fig. 2). In the case of face-on observations of fainter disks near the inner working angle, the reference star should be very close in SED and in apparent magnitude to ensure minimal speckle residuals.
Implications for WFIRST/CGI Full Field of View Imaging
One unknown feature of future space-based coronagraphs, such as the CGI on WFIRST, is the exact speckle behavior of these telescopes once they are launched. In this work, we have demonstrated that space-based observations with HST/STIS have residual speckle noise that is well characterized by a steeply varying function with stellocentric angle. To that end, we can use HST observations to investigate the potential performance of future space-based missions in the absence of active wavefront control. Such an investigation is important for two reasons. First, typical dark holes created by active wavefront control usually are limited by the size (actuator count) of their deformable mirrors, thus restricted to a range of angular separations that can be typically smaller than the full angular scale of an extended object, including most spatially resolved protoplanetary and debris disks known today. Second, there exists a nonzero probability that the deformable mirrors used in such an instrument in space may become severely degraded or stop working without the opportunity for repair. Finally, investigating the speckle behavior of operating space-based observatories can better inform input models to wavefront control algorithms, particularly if a physical basis for the observed speckle evolution is determined. Understanding the potential performance at large angular distances or in the absence of a dark hole is useful for defining new observing modes and science cases for future high-contrast imaging, particularly for resolved circumstellar material at moderate angular separations from a star.
We investigate the possible performance of a CGI-like instrument without a dark hole. To do this, we rely on parameters from the noise and detector model of CGI developed to estimate its sensitivity to exoplanets,67 discussed below and including our STIS-derived speckle variance. WFIRST/CGI will possess a similar primary mirror size, thermal properties that should be more stable than HST, and lower pointing jitter due to a low-order wavefront sensor and a fast steering tip/tilt mirror. It is thus a reasonable and conservative approach to extrapolate performance from STIS. We note that Ref. 67 assumed a systematic speckle noise source that depends on the total exposure time while we have demonstrated that uncorrected speckles combine proportional to the speckles within a single readout. This is primarily due to the fact that for STIS the time scale for speckles to change is fairly short, compared to the assumption that systematic speckle variations in CGI with wavefront correction will be as long or longer than the total exposure time on a source. This assumption will be straightforward to test during the CGI commissioning phase.
To do our calculation, we modify our existing noise model in Eq. (1) to include the impact of the EMCCD excess noise factor, the gain of the EMCCD (), and a clock induced charge noise term ():67
Next, we recalculate the expected count rate from the star on the detector, if unocculted, since we need to relate this to the estimated speckle noise and the predicted PSF wing intensity (). We will relate this to the flux of the peak pixel on the detector, just as we did in Eq. (3):Table 4 lists relevant values for assumed CGI parameters used in our calculations (taken from Ref. 68).
Assumed values for CGI performance.
The exact design for the hybrid Lyot coronagraph (HLC)69 and its performance without a dark hole is not currently known, as the design has not been finalized and few models of the performance without deformable mirrors have been generated at large angular separations. In principle, the detector FOV diameter could be as large as 9″, and CGI will likely take observations both in total intensity and polarized intensity light. We, therefore, assume that there will be an observing mode in total intensity light, the PSF intensity at angles beyond 1″ will approximate those of STIS while conserving flux with a different plate scale (, ), and the CGI PSF will be depressed by a factor due to the presence of a focal plane mask and Lyot stop (i.e., the HLC or the shaped pupil coronagraph, ). Thus we assume .
We estimate the suppression factor based on HLC designs that approximate the range in far-field PSF wing intensity in the presence of no wavefront correction when compared with the azimuthally averaged STIS PSF wings (E. Cady, private communication). At 1″, the PSF intensity is roughly a factor of 3 smaller than STIS, whereas at larger radii the intensity is a factor of 10 fainter after accounting for differing plate scales. We thus investigate the STIS PSF suppressed by a factor of 3 and 10 to bracket the likely on-orbit performance. A significant uncertainty is to what degree the mirror roughness power spectrum density might have on the far-field wings of the PSF.8 Finally, we add our empirical prescription for speckle intensity, accounting for the fact that in this case is equal to .
Figure 11 shows the resulting surface brightness sensitivity in as a function of radius between STIS’s 50CCD bandpass and our hypothetical HLC performances for a solar type star. We find that CGI performance is expected to be better than STIS, despite the latter instrument’s larger bandpass and platescale. The CGI performance at large angular radii can be further improved by binning pixels: the CGI spatial resolution is slightly oversampled by the detector in the nominal 575-nm HLC band, and noise from binning should be negligible. Similarly, STIS observing efficiency with short exposure times is quite low due to spacecraft buffer management, Earth occultations, and detector readout, so effective additional performance can also be improved by obtaining efficiencies much higher than HST. At an L2 orbit and with faster detector readout times, WFIRST/CGI is likely to achieve higher efficiencies.
Our initial investigation here has two takeaways. First, WFIRST CGI can reach or exceed HST/STIS visible light high-contrast imaging beyond the nominal HLC dark hole provided the HLC suppresses the PSF in excess of the STIS coronagraphic wedges. Second, this is particularly important for the study of previously resolved large disks since only a handful of objects have measurements of disk photometry in multiple visible wavelengths, particularly interior to 1″. CGI thus would have the ability to characterize a large number of bright circumstellar disks between 0.1″ and 3″ that have previously only been observed with STIS and/or ACS.
Advantages of STIS versus Ground-Based High-Contrast Imaging
STIS retains unique capabilities despite its age relative to newer ground-based coronagraphs. First, the broad bandpass of the STIS CCD ensures that for solar-type and early-type stars, the effective wavelength of the observations is shorter than 5800 Å, complementing near-IR observations from the ground. Especially for circumstellar disks, most ground-based observations are currently executed by observing in a differential polarization mode, or with heavy postprocessing that rarely retains low-spatial frequency information. STIS in combination with multiple spacecraft orientations and roll-differential imaging with a reference star PSF retains the full information of an image in total intensity light. STIS also has an unrivaled isoplanatic FoV compared to most ground-based coronagraphs, extending to distances of up to 25 to 50″ from a star depending on the aperture location.
STIS is well-suited for very deep imaging in the visible for circumstellar material or faint companions and compares favorably with ground-based observatories. An example is SPHERE/ZIMPOL images of TW Hya’s disk compared to those obtained with STIS, for similar total exposure times. Disk surface brightnesses fainter by a factor of a few hundred out to several arcseconds were detected in the STIS images compared to ZIMPOL.29,70 Other examples include the detection of very low-surface brightness debris disks such as HD 202628 and HD 207129.41 The -band surface brightnesses of these disks are and , respectively (see Fig. 12).
The stability of space also provides for very precise absolute flux measurements of high-contrast images, which can be difficult to do from the ground due to constantly changing atmospheric conditions. An example can again be drawn from STIS images of TW Hya, which showed differences in azimuthally averaged disk surface brightness over a 15-year time period,29 or with Pictoris, which showed variation over 18 years.28
Finally, STIS is not limited to stars brighter than a certain apparent magnitude, as is the case for AO-fed coronagraphs. There are practical limits to the faintness of a primary target, which depend somewhat on the desired contrast one wants to achieve and the exposure time for its acquisition. The limiting detectable magnitude of an object in 50 min of STIS integration time is roughly .72 For an object with , the limiting contrast for a point source at most angular radii will be , due to the background sensitivity of the image. Observers should keep this in mind if they are observing faint primary targets.
High-contrast imaging with STIS remains an active part of the unique observing modes of the HST. The recent commissioning of the BAR5 aperture location coupled with modern observing strategies and postprocessing provides for high contrasts, small inner working angles, stable space-based PSFs, and deep sensitivities over a large FOV. This combination of capabilities over a total intensity visible bandpass allows complementary observations of objects with ground-based extreme AO-fed coronagraphs, JWST, and archival NICMOS data.
The authors would like to acknowledge the careful and painstaking work of A. Gáspár in their commissioning of BAR5, and the leadership of the late B. Woodgate in the design and construction of STIS without whom this work would not be possible. We also acknowledge useful conversations with B. Nemati, B. Macintosh, E. Choquet, E. Douglas, L. Pueyo, and J. Lomax as well as careful feedback from the two reviewers that aided in the clarity of this paper. G. S. acknowledges support from HST GO program Nos. 13786 and 15219. This research has made use of the Direct Imaging Virtual Archive (DIVA), operated at CeSAM/LAM, Marseille, France. B. R. acknowledges the computational resources from the Maryland Advanced Research Computing Center (MARCC), which is funded by a State of Maryland grant to Johns Hopkins University through the Institute for Data Intensive Engineering and Science (IDIES).
http://www.stsci.edu/hst/observatory/apertures/siaf.html June 2019). Google Scholar
http://www.stsci.edu/hst/stis/documents/isrs/2017_03.pdf Google Scholar
https://archive.stsci.edu/prepds/laplace/index.html September 2019). Google Scholar
http://etc.stsci.edu/etcstatic/users_guide/1_2_1_snr.html January 2019). Google Scholar
http://www.stsci.edu/~STIS/coronagraphic_bars/GO12923raw.dat April 2019). Google Scholar
https://gea.esac.esa.int/archive/documentation/GDR2/ Google Scholar
http://www.hopewriter.com/Astronomyfiles/Data/SurveyPaper/ July 2019). Google Scholar
http://www.stsci.edu/cgi-bin/get-proposal-info?id=14426 February 2019). Google Scholar
https://wfirst.ipac.caltech.edu/sims/Param_db.html May 2019). Google Scholar
https://doi.org/10.1117/1.JATIS.2.1.011019 Google Scholar
John H. Debes received his BA degree in physics in 1999 from Johns Hopkins University and his PhD in astronomy and astrophysics from Pennsylvania State University in 2005. He is an ESA/AURA astronomer and STIS branch manager at the Space Telescope Science Institute. He has authored or co-authored more than 70 refereed papers and one book chapter. His interests include high-contrast imaging of exoplanets and debris disks, as well as the study of planetary systems that survive postmain sequence evolution.
Bin Ren received his BS degree in physics and his BSc degree in mathematical statistics from Xiamen University in 2013 and his MA degree in physics and astronomy from Johns Hopkins University in 2016. He is a PhD candidate in the Department of Physics and Astronomy and an MSEng candidate in the Department of Applied Mathematics and Statistics at Johns Hopkins University. His interests include postprocessing methods for high-contrast imaging and characterization of circumstellar systems.
Glenn Schneider is an astronomer at Steward Observatory and was the project instrument scientist for HST’s near-infrared camera and multiobject spectrometer. His research interests are centered on the formation, evolution, and characterization of extrasolar planetary systems and have focused on the direct detection of substellar and planetary mass companions to young and nearby stars and the circumstellar environments in which such systems arise and interact. He has led several STIS high-contrast imaging programs, including the commissioning of the BAR5 aperture location.