Translator Disclaimer
29 October 2019 Astrometry with the Wide-Field Infrared Space Telescope
Author Affiliations +

The Wide-Field Infrared Space Telescope (WFIRST) will be capable of delivering precise astrometry for faint sources over the enormous field of view of its main camera, the Wide-Field Imager (WFI). This unprecedented combination will be transformative for the many scientific questions that require precise positions, distances, and velocities of stars. We describe the expectations for the astrometric precision of the WFIRST WFI in different scenarios, illustrate how a broad range of science cases will see significant advances with such data, and identify aspects of WFIRST’s design where small adjustments could greatly improve its power as an astrometric instrument.



The wide field of view (FoV) and stable, sharp images delivered by the Wide-Field Imager (WFI) planned for the Wide-Field Infrared Space Telescope1 (WFIRST; Ref. 2) make it an excellent instrument for astrometry, one of the five major discovery areas identified in the 2010 Decadal Survey. WFIRST has two main advantages over other spacecraft missions: it can precisely measure very faint stars (complementary to Gaia); and it has a very wide FoV, complementary to the Hubble Space Telescope (HST) and the James Webb Space Telescope (JWST). Compared to HST, WFIRST’s wider FoV with similar image quality will provide many more astrometric targets in a single image but also hundreds more anchors to the astrometric reference frame in any field, including both background galaxies and stars with precise positions in the Gaia catalog.3,4 In addition, WFIRST will operate in the infrared (IR), a wavelength regime where the most precise relative astrometry has so far been achieved with adaptive optics images from large ground-based telescopes (e.g., 150  μas from Keck).5 WFIRST will provide at least a factor of three improvement in astrometry over the current state of the art in this wavelength range while spanning an FoV thousands of times larger. WFIRST is thus poised to make major contributions to multiple science topics in which astrometry plays an important role. In most cases, these contributions can be achieved without major alterations to the planned mission or instrument. In this paper, we summarize a few of the many compelling science cases where WFIRST astrometry could prove transformational, and then outline the areas where a small investment of attention now will ensure that WFIRST’s impact on this science is significant.


Expected Astrometric Performance of the WFIRST Wide-Field Imager

The astrometric performance of the Wide-Field Channel (WFC) in the WFI on WFIRST will depend on numerous elements, including hardware characteristics, the stability of the platform, characterization of the optics and of the detector (down to the individual pixel), the ability to design observations with the required properties for reference stars, and calibration of both the point-spread function (PSF) and of the geometric distortion (GD) of the focal plane. Here, we summarize the assumptions used in this work; Table 1 collects some performance estimates for different types of observations.

Table 1

Approximate expected astrometric performance of the WFIRST WFI for different types of observations. All estimates are for well-exposed point sources (refer to Ref. 2 and Table 4 for depths of the core survey programs). See the referred sections for details about the assumptions leading to each number. These estimates are order-of-magnitude only.

ContextEstimated performanceSections
Single-exposure precision0.01 px; 1.1 mas1.1
Typical GO program (100 exposures of one field)0.1 mas1.1
Absolute astrometry accuracy0.1 mas3
Relative PMs derived from HLS25  μasyr15.1
Relative astrometry, EML survey (per image)1 mas5.2
Relative astrometry, EML survey (full survey)310  μas5.2
Spatial scanning, single scan10  μas2.4
Spatial scanning, multiple exposures1  μas2.4
Centering on diffraction spikes10  μas2.4

Reference 6 discusses astrometric science in the WFIRST exoplanet microlensing (EML) survey field but makes very optimistic sensitivity estimates and ignores potential sources of systematic errors. For the purpose of this document, we assume that the single-exposure precision for well-exposed point sources is 0.01 pixel or about 1.1 mas; this is consistent with current experience on space-based platforms such as HST, as long as a comparable level of calibration activities are carried out.7 Depending on the platform stability and the quality and frequency of calibrations, this precision can be substantially improved by repeated, dithered observations as

Eq. (1)

where Δ is the single-exposure point-source precision in pixels, η is the plate scale, and N is the number of observations. We assume that a gain by a factor 10 (to 0.1 mas; i.e., 100 exposures) can be achieved for a typical guest-observer (GO) program. The EML survey, which will obtain many thousands of observations of each source, has more stringent requirements (Sec. 5.2). For proper motion (PM) measurements, the achievable accuracy depends on the single-image precision and the time baseline T between the first and last images. In the case of N evenly spaced images

Eq. (2)


Improvements in astrometric measurements can also be obtained by special techniques, such as spatial scanning and centering on diffraction spikes, described in Sec. 2.4, and by improvements in the pixel-level calibration, as discussed in Sec. 4.2.

Thanks to its large FoV and the availability of accurate reference sources from Gaia, each WFC exposure can achieve an absolute positional accuracy of 0.1 mas or better (see Sec. 3 for details). Although WFIRST can directly measure only relative parallaxes and PMs within its FoV, the ability to use reference stars in common with Gaia will allow parallaxes and PMs to be converted to an absolute reference frame to a worst-case accuracy of 10  μas in parallax and 10  μasyr1 in PM.


Science with WFIRST Astrometry

The science enabled by astrometry with WFIRST spans size scales from the local group (LG) to exoplanetary systems and provides important contributions to all three astrophysics goals in the NASA Science Plan. In this section, we survey the range of science topics to which this instrument can make important contributions. The astrometric precision needed for each of the following science cases is listed in Table 2.

Table 2

Required astrometric precision (in units of both WFI pixels and μas) for the different science cases discussed in Sec. 2.

SectionsScience caseAstrometric precision
2.1Motions of dwarf satellites in the LG2.2×104  pixelyr125  μasyr1
2.2Motion of stars in the distant MW stellar halo2×104  pixelyr125  μasyr1
2.3Low-mass end of the subhalo mass function1.8×104  pixelyr120  μasyr1
2.4Detection and characterization of exoplanets9×105  pixel10  μas
2.5Structure of the MW bulge9×105  pixel10  μas
2.6Star formation in the MW4.5×104  pixelyr150  μasyr1
2.7Isolated black holes and neutron stars4.5×104  pixel50  μas
2.8Internal kinematics in GCs1.8×104  pixelyr120  μasyr1


Motions of Local Group Galaxies

The range and reach of WFIRST astrometry complement and extend Gaia and Large Synoptic Survey Telescope (LSST) astrometry. Figure 1 compares the reach of current and planned PM surveys to the PMs corresponding to known velocities and distances of LG objects. Since the known orbital and internal velocities refer in almost every case to the radial component, these are intended only to represent the order of magnitude one might expect for the PMs (indeed, as in the case of M31, the orbital PMs may be significantly smaller than that inferred by radial velocity measurements alone). From this figure, it is clear that to measure PMs of satellites beyond the Milky Way’s (MW) virial radius will require better precision than LSST can achieve, at larger distances (and thus fainter magnitudes). This is the window of opportunity for WFIRST.

Fig. 1

PMs accessible to various current and planned surveys and measurements, compared to the PMs corresponding to characteristic velocities and distances for LG objects. Shaded regions show the distances and PMs for single stars accessible to the Gaia (magenta) and LSST (purple) surveys, compared to the approximate reach of the WFIRST HLS field for bright K giants assuming 15 exposures over 5 years (yellow) and the additional reach for cross-matches with HST imaging (orange). The diagonal lines show the PM associated with several characteristic transverse velocities as a function of distance: the typical range of orbital velocities in the Galactic halo (magenta and cyan) and the typical internal velocity dispersion of a dSph galaxy (green). Thick vertical lines mark heliocentric distances to: the large (blue dashed) and small (blue dot-dashed) Magellanic Clouds, the edge of the MW halo (gray), M31 (red), and the approximate edge of the LG (green). Gray dotted vertical lines mark heliocentric distances to other dwarf galaxies in the LG, including satellites of the MW and M31.8 Current PM measurements by HST912 and Gaia13 for MW GCs and satellite galaxies and for individual stars in the Sgr tidal stream14 are plotted as black points/symbols.


WFIRST astrometry is a crucial component of constraints on the nature of dark matter (DM) from the orbital and internal PMs of the satellite galaxies of the MW. The orbits of dwarf satellites can be used to help map the MW’s own DM halo out to its virial radius and beyond, placing our galaxy into a cosmological context. A complete knowledge of the MW halo’s properties enables tests of DM models through comparisons with predictions from simulations for its mass and shape, accretion history, and the mass and orbit distributions of its satellite galaxies. HST PMs of dwarf satellite galaxies, and Gaia astrometry in the inner Galaxy, will both make great strides toward this goal; however, Gaia has insufficient depth and HST insufficient FoV to reach to the edge of the MW halo (where the total mass is uncertain to a factor of 4; see Ref. 15) or to obtain internal PM measurements for many dwarf galaxies, which are crucial to break velocity degeneracies and understand the small-scale distribution of DM (see also work by the Gaia Challenge group,16 summarized in Fig 2.2 of Ref. 17).

The HST proper motion (HSTPROMO) campaign has used the HST, which has image quality similar to that expected for WFIRST, to measure PMs of both bound objects and stream stars in the MW (e.g., Refs. 10 and 18). In the coming years, HST will set a PM baseline for many more distant satellites:19 by the time WFIRST is ready, this baseline will be roughly 8 to 10 years for these satellites and far longer for dwarf galaxies with earlier observations. With its larger FoV and more sensitive detectors, WFIRST should be able to expand these measurements to more distant galaxies and achieve better accuracy due to the larger number of calibration objects available (and to the establishment of Gaia’s astrometric frame, see Sec. 3). Figure 2 shows the estimated number of stars in each of the MW’s satellites that are accessible to WFIRST at the depth of the planned high-latitude survey (HLS, J<26.7, blue) core program, as well as for a possible spatial-scanning mode on WFIRST (H<16, red). The upper panel of the figure shows the expected tangential velocity error for each dwarf assuming PM precisions of 25  μasyr1. Internal PM dispersions, as was recently done for the Large Magellanic Cloud with HST,20 are reachable with WFIRST for the galaxies shown in cyan. These include three ultrafaint galaxies (Segue I, Draco, and Ursa Minor) that are sufficiently DM-dominated to distinguish between cold DM models (which predict a cuspy inner mass profile) and warm or “fuzzy” DM models (which predict cored inner profiles). Draco and Ursa Minor each have 80 stars at H<16 and their internal velocity dispersions are marginally resolved even at 25  μasyr1, making them particularly good cases for spatial scanning to improve the internal PM accuracies by an order of magnitude.

Fig. 2

(a) Estimated observational error in tangential velocity assuming PM precision of 25  μasyr1. Galaxies in cyan have estimated velocity errors (δμ) comparable to or less than their intrinsic velocity dispersions (σv).(b) Number of stars in LG dwarf satellite galaxies brighter than the limiting apparent magnitude of the WFIRST HLS (J<26.7, blue) and the limit for spatial scanning (H<16, red). Plot courtesy Matthew Walker.



Motions of Stars in the Distant MW Halo

Besides dwarf galaxies, the MW’s halo contains the tidally disrupted remains of previously accreted galaxies, known as tidal streams. We expect that tidal debris should extend to at least the virial radius of the MW,21 but currently the most distant MW halo star known is an M giant at around 250 kpc.22 The most distant known populations with statistical samples of stars (BHB, RR Lyr, and M giant stars) extend to around 150 kpc (half the virial radius) at the magnitude limit of current surveys,2325 but the WFIRST HLS fiducial depth will reach to the MW’s virial radius down to the main-sequence turnoff. The orbits of distant stars probe the extent and total mass of the MW dark halo; they also represent a unique population of recently accreted small galaxies. With the HLS’s projected depth, the transition between the MW’s and M31’s spheres of influence, and perhaps the splashback radius of the MW (e.g., Ref. 26), could also be detected. PMs from WFIRST are crucial to these endeavors since complete phase-space information for these stars is the best way to confirm that stars associated in position at large distances are from the same progenitor and to connect groups on opposite sides of the galaxy through their orbits, leading to constraints on the mass profile and flattening of the Galactic dark halo at large distances (e.g., Ref. 27).

High-velocity stars are another interesting target, whether for GO observations or as serendipitous objects in repeated survey fields. These stars’ orbits, which have an extremely wide radial range, can potentially also be used to constrain the overall shape and mass of the MW DM halo.28

At distances of 100 to 300 kpc, preliminary work shows that PM precision of 25  μasyr1 or better is required to eliminate outliers in groups and connect structures on opposite sides of the galaxy using their phase-space positions (Fig. 3; Secunda et al. in preparation). A similar precision would be needed to identify high-velocity stars at or near the Galactic escape velocity.

Fig. 3

Identifying stellar structures in the distant Galactic halo. (a) Groups of stars identified in a mock stellar halo29 in the range 100 to 300 kpc, using sky positions (shown) and distances only. (b) Same stars colored by progenitor galaxy. Green arrows highlight the contribution of interlopers to group 5 (dark orange) in (a). (c)–(e) View of the same groups in energy-angular momentum projection, which requires six-dimensional phase space information including PMs. (d) Some outliers are already identifiable at 25  μasyr1 and (e) structures are clearly distinguishable at 5  μasyr1 (green arrows). Groups 2 and 4 (dark and light blue, respectively) in (a) are from the same tidally disrupted progenitor galaxy but are found on opposite sides of the sky; with 25  μasyr1 precision they can be associated through orbit integration, reflected in (d) and (e) by their similar values of angular momentum (jz).



Constraining the Low-Mass End of the Subhalo Mass Function

Astrometry from pointed guest observer (GO) projects could provide several routes to understanding the distribution of low-mass substructure in the MW. The abundance (or lack) of low-mass structure is a key to differentiating between cold and warm DM models since cold DM predicts abundant substructure at small scales while in warm DM the mass function is cut off at a characteristic scale related to the intrinsic temperature of the DM particle (and hence to its mass in the case of a thermal relic).

One route is to search for perturbations to tidal tails from distant globular clusters (GCs) or dwarf galaxies. Current searches are focused on quantifying substructure in the spatial distribution of tidal debris but are limited by knowledge of the MW background/foreground in the region of the stream as well as by Poisson fluctuations in the star counts (e.g., Ref. 30) and by the uncertain dynamical ages of streams, which depend on modeling the orbit (e.g., Ref. 31). WFIRST’s capability to reach deep into the stellar main sequence (MS) at these distances will help mitigate the shot-noise issue. More importantly, obtaining astrometry of fields including tidal streams would allow superior selection of stream stars relative to the background/foreground, improving the sensitivity to density fluctuations, and allowing better constraints on the time when material was first tidally stripped. Streams commonly stretch tens of degrees over the sky, so WFIRST’s large FoV is uniquely well suited to this application. A thin stream usually has a velocity dispersion of 1 to 10  kms1, so to provide a useful PM selection for a stream at 50 kpc would require relative PMs to a precision of about 20  μasyr1.

Another possibility is to search for deviations in the apparent positions of quasars due to strong lensing by dark substructures. For a distant quasar lensed by a 108M subhalo at 50 kpc, the Einstein radius of the lens is roughly 20  μas (presuming a singular isothermal sphere). One route would be to look for the time-dependence of the lensing around a single quasar: for a subhalo moving at 200  kms1, the time to cross the Einstein radius is about 10 days. Alternatively, a wide field containing many quasars could be examined for statistical deviations between exposures taken at different times (separated by longer than 10 days). Either approach would require absolute astrometry (i.e., consistent between exposures) accurate to 20  μas.


Detection and Characterization of Exoplanets

Very accurate PM estimates will make it feasible to search for the astrometric signature of exoplanets around nearby stars. For competitive constraints on exoplanet masses and orbital parameters, an instantaneous precision of better than 10  μas is required. The best constraints can be achieved for the most nearby stars (d10  pc). A dedicated GO program that specifically observes those most promising targets with a flexible schedule is therefore complementary to the EML.

Because of their close proximity, the target stars are generally very bright. This makes high-precision astrometry possible by using one of two different methods: spatial scanning and diffraction spike modeling.


Spatial scanning

Spatial scanning involves intentionally slewing the spacecraft during integration to create extended tracks from bright target and reference stars in the field of interest. This spreads out the signal from each star over hundreds or thousands of pixels, thereby avoiding saturation while integrating orders of magnitude more photons, and averaging over pixel-level artifacts that may significantly affect pointed observations.32,33 Scans in different sky directions can be combined to yield high precision for both coordinates. HST has attained precisions of 20 40  μas with this technique, limited in part by the small number of available reference stars and the variation of the focal plane geometry on the orbital time scale of the telescope (1  h). Because of its larger FoV and more stable orbit, we expect that the WFIRST WFC will be able to achieve precisions closer to the noise limit, about 10  μas per exposure. By combining multiple exposures, it will then be possible to achieve a final relative astrometric accuracy of 1  μas.

Desirable slew rates for spatial scanning are 0.510  s1, roughly corresponding to 12 to 250 pixels per read; this is the length of the region over which the light from each star will be spread within one readout frame. The fast, nondestructive reads of the WFIRST WFC will allow a clean separation of the signal accumulated within each pixel from different sources at different times, greatly reducing the confusion due to overlapping trails that has affected applications of this technique using the Wide-Field Camera 3/Ultraviolet-VISible (WFC3/UVIS) detector on HST. Within the desired range of scanning speeds, it will be possible to observe unsaturated sources 7 mag brighter than the pointed-observation saturation limit, or HAB4 mag. The fastest available scan speed primarily affects the maximum brightness of the source that can be accomodated; slower scan speeds in the range 0.5 to 2  s1 can achieve essentially the same benefits but with a fainter saturation limit.

Both confusion effects and the signal-to-noise ratio for spatial scanning observations would benefit more than pointed observations from obtaining all independent reads for each exposure: unlike pointed observations, signal does not build up linearly over time in each pixel but is deposited there during the narrow time interval in which a star passes over that pixel. Extending the interval between available reads increases both the background accumulated in each pixel without a corresponding increase in the signal, and the time interval over which signal from different stars in the same pixel cannot be cleanly separated. The availability of intermediate reads for download is of course subject to mission-level limits on science telemetry, so the number of reads to download may need to be determined on a scene- and project-dependent basis. Finally, spatial scanning observations will most likely need to be obtained under gyro control, as the required motion of the spacecraft will quickly exceed the size of the guiding window. More details, including error budgets, will be included in an upcoming white paper (Casertano et al., in preparation).


Centering on diffraction spikes

A second potential strategy for obtaining highly accurate astrometry of very bright stars involves centering on diffraction spikes. The approach is facilitated by the properties of the WFIRST H4RG detectors, which, unlike charge-coupled devices (CCDs), do not show “bleeding” of excess charges from saturated pixels to their neighbors (see Sec. 4.2). Astrometric precisions of 10  μas or better are achievable with this technique with integrations of 100 s for stars with J=5 or, making use of the recently added optical R062 filter, R=6.34 (Diffraction-spike measurements are superior in the short-wavelength range because the diffraction spike is sharper. Unlike the core of the WFIRST PSF, the diffraction spike is well sampled even in the bluest WFIRST filter given a pixel scale of 0.11 mas.)

As for spatial scanning, measurement accuracy will likely be limited by systematic uncertainties, in particular the fidelity of corrections for optical distortion and pixel-level artifacts (cf. Sec. 4). Thus, performing several exposures per visit is beneficial and should be able to yield precisions of 10  μas or better even in the presence of residual systematics.


Detection of Earth-mass exoplanets

These estimates indicate that a dedicated GO program with visits to target fields separated by months and spread out over the lifetime of WFIRST could detect Earth-mass exoplanets astrometrically around the most nearby stars, in some cases even in their respective habitable zones. In addition, it can probe Neptune-class planets around more distant stars and, by adding earlier measurements from Gaia, rocky planets with periods of >10  yr. Such measurements would be strongly synergistic with radial velocity campaigns, improving the mass constraints and breaking degeneracies in several orbital parameters,35 and enabling mass estimates of the direct-imaging exoplanets of the WFIRST coronagraph and possible starshade occulter programs.34


Detailed Structure of the Inner Milky Way

Gaia will revolutionize our understanding of MW structure in the outer parts of the MW, including the halo. However, Gaia has a very limited view of the inner MW due to the significant extinction in the Galactic plane at optical wavelengths (Fig. 4 shows an illustrative example using a simulated galaxy) as well as crowding (not accounted for in Fig. 4). WFIRST will probe significantly deeper into the inner MW and be less seriously affected by crowding than Gaia, allowing us to map the structure and kinematics of this region and complement the Gaia view. As an example, the EML survey will obtain precise parallaxes and ultraprecise PMs for over 50 million stars in a small area of the Galactic bulge, enabling a detailed analysis of their kinematics and density distribution. Currently, studies of bulge stellar populations are limited by the quality of the PM and the need to remove foreground disk stars, typically achieved via kinematic or photometric filters (see, e.g., Ref. 43). Both are statistical in nature and do not provide a direct determination of the distance to individual stars. According to the current requirements, a mission-long astrometric accuracy of 10  μas or better (with a stretch goal of 3  μas) should be achieved at HAB=21.6 (EML 20; see Sec. 5.2).

Fig. 4

Simulated completeness of the distribution of red clump stars with distances |z|<500  pc from the Galactic plane, based on the MW-like simulated galaxy in Ref. 36. (a) Stars that Gaia would detect in the optical with G<20; (b) those seen by WFIRST in the IR with K<26. Gray contours in both panels show the density of the complete distribution on a logarithmic scale; the black cross marks the location of the Sun. The synthetic red-clump catalog was constructed by drawing stars in the range 0.48<MI<0.08, 0.8<VI<1.437 from the model isochrones in Ref. 38, distributed according to the age and stellar mass density of the simulated star particles.39 The three-dimensional extinction map in Ref. 40 was interpolated to determine apparent magnitudes and reddening, and Gaia G magnitudes were calculated using pygaia.41 This simulated view ignores the effects of crowding (which significantly affect Gaia in the plane but are anticipated to be relevant for WFIRST only within 0.5  deg of the Galactic Center; see Sec. 2.6) and does not include a prominent Galactic bar (see Ref. 42). Gaia will largely be limited to heliocentric distances <4  kpc in the plane, and WFIRST can measure parallaxes and velocities of stars well beyond the Galactic Center.


At comparable accuracy in relative parallaxes, distances to individual stars can be measured to 9% at the bulge (3% if the stretch goal is achieved), and useful distance discrimination should be obtained to significantly fainter magnitudes. The two tangential components of the space velocity can be recovered to the same accuracy (in this regime, the distance error is dominant over the PM error in deriving the space velocity). This information will enable a much cleaner determination of the kinematics of the bulk of bulge stars in the EML survey field and readily identify subgroups of stars—disk or bulge—with anomalous kinematics. If depth effects can be accounted for, the end-of-mission PM accuracy translates to a velocity precision of 1  kms1; together with the very large number of stars measured, this will permit a clear component separation of the spatially overlapping bulge and halo populations (see, e.g., Ref. 44) and potentially identifying complex structures such as the anomalous motions found in the X-shaped regions of the bulge (see, e.g., Ref. 45). In principle, Gaia will achieve comparable precision over all of the bulge but only for the bright red giants at G15 or brighter; the uncertainties will be considerably larger (2 orders of magnitude) at Gaia’s faint limit.

Regions within 0.5 deg of the Galactic Center will likely suffer from crowding at H<21, as has been seen with HST-WFC3IR studies of this region,46 which would limit the astrometric precision to >0.5  masyr1 for stars fainter than this. Beyond this region, the stellar density is not typically high enough to impact astrometric precision down to H<24. However, the exact determination of how the WFIRST astrometric precision will scale with stellar density, SNR, PSF knowledge, and survey depth will require image-level simulations in the future.


Star Formation in the Milky Way

With the advent of large IR surveys of the Galactic Plane, many new young star clusters have been identified. The most massive of these young clusters are ideal laboratories for studies of star and cluster formation, stellar evolution, and cluster dynamics, but detailed studies of these regions are hampered by high and spatially variable extinction, high stellar densities, and confusion with foreground and background stars. Many of these limitations can be overcome with the addition of PMs observed in the IR, to separate out the comoving cluster members from the contaminating field population.7,47 Furthermore, measurements of the internal velocity structure of star clusters provide constraints on the unseen stellar population from dynamical mass measurements, thereby informing cluster evolution models.48

WFIRST is ideally suited for studies of massive young clusters and star forming regions in the MW, given its wide FoV at IR wavelengths, high spatial resolution, and potential for precise photometry and astrometry. For rapidly moving populations in the center of the Galaxy, a PM precision of 0.5  masyr1 per star is needed to separate cluster members from field stars. To obtain internal velocities or separate clusters in the disk, a PM precision of 0.05  masyr1 or better is desired. Even higher astrometric precisions would enable searches for binaries and higher-order multiples.

An important factor to consider for this science case is that cluster members span a large range in brightness. The brightest and most massive cluster members in clusters beyond 4 kpc often have J=9 or brighter. Careful calibration of persistence, shorter integration times, or possibly narrow-band filters will be needed to reach both bright, massive members and faint, low-mass members of clusters.


Isolated Black Holes and Neutron Stars

Our Galaxy likely contains 107108 stellar mass black holes and orders of magnitude more neutron stars.49 Measuring the number and mass statistics of these stellar remnants will provide important constraints on the initial stellar mass function, the fate of massive stars and the initial-final mass relation, the star formation history of our Galaxy, and the fundamental physics of compact objects. WFIRST has the ability to find such objects in large numbers through gravitational microlensing when a background star passes behind the compact object and is magnified photometrically. However, only the addition of WFIRST astrometry will enable us to measure the precise masses of these objects through astrometric microlensing. The apparent astrometric shift of the background star due to microlensing, which is proportional to M1/2, is 1 milliarcsecond for a 10M black hole at 4 kpc lensing a background star at 8 kpc (Fig. 5). Thus, the necessary astrometric precision to detect isolated black holes is <150  μas; a factor of 2 to 3 better precision would also allow the detection of neutron stars.

Fig. 5

Astrometric shift of a background bulge star (source, d=8  kpc) lensed by a foreground compact object such as a black hole, neutron star, or white dwarf (lens, d=4  kpc). The astrometric shift changes as a function of the projected source-lens separation on the sky, u, in units of the Einstein radius. For the 10M case, the Einstein radius is 4 mas and the time for the source to cross the Einstein ring is typically >100  days.


While Gaia or ground-based adaptive optics systems may detect one or a few isolated black holes,50 WFIRST’s IR capabilities and monitoring of the Galactic Center and bulge fields will yield the much larger samples needed to precisely measure the black hole and neutron star mass function and multiplicity. Microlensing by massive objects typically has long timescales, with Einstein crossing times >100  days for black holes, so WFIRST astrometry should be stable on these timescales; i.e., routinely calibrated on sky if possible.


Globular Clusters

In the last decade, a wealth of revolutionary studies has dramatically changed the traditional view of GCs as the best examples of “simple stellar populations:” stars with the same age and chemical composition. The presence of multiple stellar populations (MPs) in GCs has been widely established along all the stellar evolutionary phases (e.g., Ref. 51, 52, and references therein): spectroscopic studies have found significant star-to-star variation in light elements (e.g., Ref. 53 and references therein), whereas high-precision photometry, mostly from HST data, has clearly revealed the presence of distinct sequences in color-magnitude diagrams (CMDs) at all wavelengths (e.g., Refs. 52 and 56; see also Fig. 6 of Ref. 57). Several GCs have also shown the presence of significantly He-enhanced subpopulations (e.g., Refs. 58 and 59) and even subpopulations with distinct iron content in a few cases such as ω Cen, M22, Terzan 5, M54, NGC 5824, and M2.6063 These observational findings present formidable challenges for theories of the formation and evolution of GCs and have inaugurated a new era in GC research in which understanding how multiple stellar systems form and evolve is not just the curious study of an anomaly but a fundamental key to understanding star formation.

Fig. 6

(a) The mF275W versus mF275WmF336W CMD of ω Cen, showing several subpopulations of stars in all evolutionary sequences (from Ref. 54). MPs reveal themselves in high-precision photometry from the UV to the near-IR. (b) mF160W versus mF110WmF160W CMD of an outer field of ω Cen, corrected for differential reddening. Field stars (orange dots) are identified using PMs. From Ref. 55.


Measuring the PMs of stars in GCs is the most effective way to constrain the structure, formation, and dynamical evolution of these ancient stellar systems and, in turn, that of the MW itself. High-precision HST astrometry of GCs is now becoming available for a large number of objects (e.g., Ref. 64), but current PM catalogs are limited by the small FoV of HST, either to the innermost few arcminutes or to pencil-beam locations in the outskirts. While most dynamical interactions do happen in the center of GCs, answering many outstanding questions will require high-precision PMs of faint cluster stars over wide fields, for which WFIRST is by far the best tool. Here, we discuss a few examples of such investigations.


Multiple-population internal kinematics

The PM-based kinematic properties of MPs have so far been characterized for only three GCs: 47 Tuc,65 ω Cen,66 and NGC 2808.67 The short two-body relaxation timescale in the inner regions of these clusters, where most observations have so far been focused, implies that any initial differences in the kinematic properties of different stellar populations have likely been erased. The cluster outskirts, however, have much longer relaxation timescales and could still retain fossil kinematic information about the early stages of cluster evolution (e.g., Ref. 68). The outer regions can thus provide a wealth of information and constraints on the formation and early dynamics of MP clusters, on the subsequent long-term dynamical evolution driven by two-body relaxation, and on the role of the Galactic tidal field in the outskirts of clusters. For instance, it has been shown that second-generation stars in 47 Tuc,65 NGC 2808,67 and ω Cen69 are characterized by an increasing radial anisotropy in the outer regions with respect to first-generation stars (see Fig. 7). Even further out, at distances approaching the tidal radius, the effects of the external tidal field are expected to lead to a more isotropic velocity distribution.

Fig. 7

(a) Radial (blue) and tangential (green) PM dispersions as a function of color for the bluest (left) and reddest (right) parts of the MS of 47 Tuc (top) and of the Small Magellanic Cloud (bottom) (from Ref. 65). (b) Deviation from tangential-to-radial isotropy (horizontal line) for the five subpopulations in NGC 2808. Vertical lines mark the locations of rh, 1.5×rh, and 2×rh (from Ref. 67).


Both WFIRST’s wide FoV and its improved sensitivity will make revolutionary steps forward in understanding the initial differences in MPs if sufficient PM accuracy can be achieved. Due to mass segregation, the most abundant stars in the outskirts of GCs are low-mass, faint main-sequence objects. Gaia can only measure stars as faint as the turn-off region in most clusters and, therefore, will not be able to provide enough statistics to properly characterize the kinematics of the outer cluster regions. The expected internal velocity dispersion of cluster stars near the tidal radius is of the order of 1 to 3  kms1, even for the most massive clusters. The PM error adds in quadrature, so it should be less than half the intrinsic velocity dispersion, i.e., 1  kms1, in order to measure dispersions in cluster outskirts. At the typical distance of Galactic GCs, 10  kpc, this translates into PM errors of the order of 20  μasyr1.


Energy equipartition

It is widely assumed that GCs evolve toward a state of energy equipartition over many two-body relaxation times so that the velocity dispersion of an evolved cluster should scale with stellar mass as σmη, with η=0.5. Recently,70 used direct N-body simulations with a variety of realistic initial mass functions and initial conditions to show that this scenario is not correct (see also Ref. 71). None of these simulated systems reaches a state close to equipartition: instead, over sufficiently long timescales the mass-velocity dispersion relation converges to the value η0.08 as a consequence of the Spritzer instability (see Fig. 8).

Fig. 8

Time evolution of the energy equipartition indicator η for single main-sequence stars in N-body simulations, from Ref. 70. The time along the abscissa is expressed in units of the initial half-mass relaxation time trh (0). Complete energy equipartition (η=0.5; dotted line) is never attained, confirming previous investigations based on stability analysis.


These intriguing results have just started to be observationally tested (e.g., Refs. 69 and 72). To measure η, a wide range of stellar masses must be probed. Again, this task is out of reach for Gaia because of its relatively bright magnitude limit, but WFIRST will easily measure high-precision PMs down to the hydrogen-burning limit (HBL; 0.08M) and out to the tidal radius, thus constraining both the current state of energy equipartition in a cluster and its past dynamical evolution. As in the previous case, PM errors of the order of 20  μasyr1 are needed.


Hydrogen-burning limit and the brown-dwarf regime

WFIRST will also make it possible to study the luminosity functions of GCs beyond the HBL and into the brown-dwarf regime. Close to the HBL, old stars show a huge difference in luminosity for a small difference in mass, resulting in a plunge of the luminosity function toward zero for stars with masses just above this limit. Stars in GCs are homogeneous in age, distance, and chemical composition (within the same subpopulation), so at the typical GC age of >10  Gyr, stars with masses below the HBL will have faded by several magnitudes relative to those above it, creating a virtual cutoff in the luminosity function (e.g., Refs. 73 and 74; see also Fig. 9).

Fig. 9

Deep near-IR CMD of the GC M4, from Ref. 74. The white-dwarf and brown-dwarf regions are labeled, and low-mass stellar models are over-plotted in green and red. The expected end of the H-burning sequence is marked with red dashed lines and a shaded area.


The best place to observe the properties of stars approaching the HBL is once again outside a cluster’s core region, where contamination by light from much brighter red-giant-branch stars is negligible. The brown-dwarf regime in GCs is unexplored ground, so many new and intriguing discoveries may be waiting for WFIRST. Due to the relatively low number density of low-mass MS and brown-dwarf stars in the cluster outskirts, WFIRST is the perfect astronomical tool for these investigations as well. Proper-motion-based cluster-field separation is needed to create clean samples of cluster members; PM errors of the order of a few tenths of masyr1 might be sufficient to separate cluster stars from field stars, whereas errors one order of magnitude smaller would also enable studies of the internal kinematics of cluster stars in these lowest-mass regimes.


WFIRST Absolute Astrometric Performance

High-precision absolute astrometric measurements with HST have typically been based on the positions of well-measured background galaxies within the FoV, but a new major improvement in absolute astrometry measurements is imminent. When WFIRST begins observing the universe, the European Space Agency (ESA) mission Gaia will already be complete, providing absolute astrometric positions of unprecedented accuracy everywhere on the sky. WFIRST will be able to make use of Gaia’s absolute astrometric reference frame to convert its relative astrometry to absolute astrometry.

In broad terms, there are two methods to determine the absolute astrometric accuracy of a scientific observation. Initial astrometric positions can be obtained from the telescope pointing information using guide star (GS) data (a priori). This information includes the celestial coordinates of the GSs and the locations of the scientific instruments relative to the GSs in the focal plane of the telescope. These positions can be refined based on information available after an observation is made, namely the positions of all sources with accurate coordinates in external catalogs (a posteriori). Although there will be a minimum of four GSs, perhaps as many as 18, used to guide WFIRST observations, many more fainter stars within each exposure can be used to improve the absolute astrometric precision with the a posteriori method.75

According to WFIRST’s design and operations concept, GSs can be placed on any of the 18 WFI detectors, but at most one GS will be assigned to a given detector. Assuming the nominal 10<H<15.6 bright and faint GS magnitude limits, analysis has shown (Ref. 76) that the Two-Micron All Sky Survey (2MASS) can provide at least one GS candidate for each of the 18 WFI detectors with a probability close to 100% (and at least 10 detectors will have a GS candidate brighter than H=12.6). While perhaps as few as four bright GS (H<13) will suffice for the attitude control system, the fine guidance system design supports the use of up to 18 GSs since each of the detectors will read out a “guide window” to keep the readout pattern for all 18 detectors synchronized, even if the guide window does not contain a GS. For the grism mode, the whole sky will be available if the faint limit for GSs is HAB<14, but if the faint limit is pushed to HAB<13 or even HAB<12, then WFIRST will not be able to perform spectroscopy in a few regions (0.001% and 3.5% of the sky, respectively) around the Galactic poles.

The WFIRST mission will begin operation in the second half of the 2020s, several years after Gaia has completed its nominal 5-year mission in July 2019. However, the Gaia mission has already been extended to the end of 202277 and could in principle be extended up to a total of 5 years beyond its nominal mission based on the depletion rate of consumables and the degradation rate of the main CCD camera.

Near the faint Gaia limit (19<G<20), PMs in the Gaia catalog will have an end-of-mission error (assuming the nominal 5-year baseline) of about 0.20.3  masyr1 (see Ref. 78). This translates into a position uncertainty of about 1.6 to 2.4 mas (or 0.015 to 0.02 WFIRST WFI pixels) at the start of WFIRST operations and about twice as much by year 5. Position errors of this size will have a significant impact if the goal is to achieve high-precision (to better than 0.01 pixels) absolute astrometric measurements. In the following, we provide expected estimates based only on catalog errors. All other sources of errors (geometric-distortion residuals, centroiding errors, etc.) are ignored (more in Ref. 75). Gaia’s extension for another 5 years (the estimated maximum possible) will improve on the following analysis significantly, not only with improved PMs (lowering uncertainties by a factor 22 compared to the 5-year baseline) but also by reducing the timespan over which positions must be extrapolated (another factor of 2 improvement, assuming WFIRST begins operation in 2025).

Single-epoch GO and Guest-Investigator (GI) observations of a random location on the sky may have to rely solely on the information contained in prior astrometric catalogs (in particular Gaia’s catalog) to determine the absolute position of their sources. Assuming an average per-star positional error of 2 mas (corresponding Gaia’s expected end-of-mission astrometric error for GGaia=19, extrapolated to the late 2020s) and ignoring all other sources of errors (e.g., geometric-distortion or source centroiding errors), it will be possible to obtain absolute a posteriori positions to better than 0.05 mas (or about 5×104 WFI pixels) over half the sky. In regions with the lowest stellar densities, the expected absolute position error increases to 0.1 mas (103 WFI pixels).

For the planned WFIRST mission surveys (the HLS and EML surveys), repeated WFIRST observations spanning several years can be used to improve Gaia’s PMs, especially at the faint end, and to derive absolute positions and PMs for many fainter sources. The astrometric precision for the planned surveys is expected to be significantly better than what can be done with Gaia alone but is difficult to quantify at this time.

For all stars suitable for the a priori method, Gaia’s expected end-of-mission astrometric error ranges between 10 and 80  μasyr1,78 but WFIRST will likely choose the four GSs among the brightest available sources. We estimate that at least 7 to 8 GSs in the range 10.0<H2MASS<10.2 will be available over half the sky. If these stars land on at least four different WFI chips (a near 100% chance; Ref. 76), assuming their average magnitude is H2MASS10.1 (corresponding to GGaia13.7), and assuming they have a Gaia-extrapolated position error of 0.15  mas in late 2020s, then the a priori method is expected to offer absolute position measurements at the 0.075 mas level or better (7×104 WFI pixels) for half the sky. For the entire sky, on the other hand, we always expect at least 7 to 8 GSs within any given WFI FoV if the faint limit is relaxed to H2MASS=12.4 (or roughly GGaia=15.8). This translates into an upper limit for the expected a priori astrometric error of 0.2 mas (or about 2×103 pixels).

The proposed 5-year extension to the Gaia mission would improve WFIRST astrometry substantially. In this case, each Gaia source will have twice as many measurements over twice the time baseline, providing an increase in precision by a factor of 22. Moreover, the time between missions, and hence the interval over which WFIRST would need to extrapolate Gaia’s positions, would be reduced by 2.5 years, resulting in an additional factor of 1.4 improvement at WFIRST’s first light.

Finally, the ability to guide the telescope using more than four GSs, optimally one in each detector, provides the means to monitor the stability of WFIRST focal plane solution for all such observations by comparing the observed relative positions of the GSs to their catalogued positions. The ground system can impose the requirement that these GSs have Gaia positions (and parallaxes and PMs). The GS positions are reported 6 times per second, much more frequently than the WFI full frame images, even more so considering that not all frames will be saved to the recorder. Moreover, the Gaia field stars in the full frame images will likely have saturated PSFs after the first few reads. Therefore, the GSs provide a unique opportunity to monitor the focal plane solution, which is critical for high accuracy astrometry. If significantly fewer than 18 detectors routinely host a GSs, then the stability of the focal plane solution may need to be accessed using dedicated calibration observations, with the results interpolated to the intervening science visits.

The uncertainty of the conversion of relative to absolute parallaxes and PMs depends on the number of reference sources and their individual Gaia measurement errors. In regions of low stellar density—e.g., near the South Galactic Pole but away from NGC 288—the Gaia DR2 catalog contains about 150 stars per WFIRST FoV with Gaia magnitudes 17<G<19. These stars have median Gaia DR2 uncertainties of 170  μas in parallax and 300  μasyr1 per component in PM. At the end of the Gaia mission, these uncertainties are expect to drop below 100  μas and 100  μasyr1, respectively, allowing a conversion to absolute parallax and PM with a worst-case error better than 10  μas in parallax and 10  μasyr1 in PM. Typical performance in areas with higher stellar density will likely be much better. While systematic issues still exist with Gaia parallaxes and PMs at the level of a few tens of μas,79,80 improved calibration and processing will likely reduce these substantially in future releases.



Here we consider what is most likely to have an effect on the astrometric performance of the WFI. We highlight areas where astrometry-specific considerations are especially important and can add significant extrascience capability with little to no extracost. Our recommendations are summarized in Table 3.

Table 3

Summary of main recommendations for astrometry.

4.1GDThis significant systematic error for astrometry, not currently covered by core science requirements, should be considered in calibration plans for the WFI.
4.2Pixel-level effectsGround-based calibration should be considered, based on results of current tests by several labs. A spatial scanning mode would mitigate these effects for bright stars.
4.3Filters and color dependenceLikely straightforward to calibrate, but should be aware of systematic effects.
4.4Readout hysteresisStraightforward to minimize based on experience with current generation of HxRG detector/amplifier combinations.
HLSOptimal to evenly space observations over full time of survey, to extent permitted by other requirements. Current example schedules vary in PM outcome by factor of 2.
EML surveyProgramming an occasional larger dither will significantly help calibrate for general astrometry. Largest possible time-spacing between first and last exposures is optimal; regular intermediate observations will increase understanding of long-term PSF variations.
GOThe TAC process should allow for multiyear GO proposals to optimize PM baselines. For proposals covering large sky areas, time between field revisits should be maximized.
4.6JitterThis may be an issue for WFIRST where it was not for HST, given large requirement (14 mas). Requirements of the HLS for galaxy shape determination should help.
4.7Data managementDownloading every read with no coadds for at least part of the FoV is highly desirable for spatial scanning observations (see Sec. 2.4). Downloads of GS postage stamps are crucial and inexpensive for PSF jitter correction.
4.8High-level data products and archiveAstrometry (linked to the Gaia frame) and astrometric uncertainties (including PSF centroiding error estimates) should be part of the high-level products. A requirement should be set on the astrometric uncertainty. The archive should allow for multiple upgradable astrometric solutions.


Geometric Distortion

GD is the most significant systematic contributor for astrometry that is not currently covered by explicit requirements for either the HLS or EML survey. A dedicated set of observations to autocalibrate the GD of the WFI is currently being considered. There are two main ways to solve for the GD: using previous knowledge of the stellar positions in the field from existing astrometric catalogs (the “catalog” method), or via autocalibration, in which stellar positions themselves are iteratively solved for together with the GD. Each of these approaches has different advantages and disadvantages, but both depend strongly on the precision with which the position of stars can be measured using appropriate PSFs (e.g., Refs. 8182.

The catalog method is less demanding of telescope time since it requires fewer images to calibrate the GD and monitor temporal variations, but it strongly depends on the quality of the astrometric catalog used as a reference since the residuals and the systematic errors present in the reference catalog can propagate in the GD solution. In addition, unless the reference catalog is based itself on images taken very close in time to the WFIRST calibration images, PMs (and their errors) can introduce significant residuals in the GD solution. A recent technical report on the catalog method (Ref. 92) also highlights the importance of using accurate PSFs that take into account for jitter and interpixel capacitance effects.

The autocalibration approach requires more images, and therefore more telescope time, but offers a self-consistent calibration solution and can be designed to be formally insensitive to proper-motion-related errors.

On-sky GD calibration has historically been performed using large dithered exposures (as wide as the FoV in some cases) of a homogeneously distributed, moderately dense stellar field. Stars in the EML survey fields are homogeneously distributed, and while their overall stellar density can be too high, this can be mitigated by using only the brightest stars in each field to calibrate the GD. If the exposure time has been carefully chosen, the bright, unsaturated stars will still be reasonably far apart from each other, and their surrounding neighbors typically a few magnitudes fainter, so that the bright stars can still be considered fairly isolated. The Baade window was successfully used by Refs. 83and 88 to calibrate the GD of two different ground-based, wide-field IR detectors (ESO WFI@MPG and HAWK-I@VLT, respectively). The GC ω Cen, another possible target field for WFIRST calibration, was used by Ref. 89 to calibrate the GD of the IR WFI VISTA InfraRed CAMera (VIRCAM) at the Visible and Infrared Survey Telescope for Astronomy (VISTA). The same technique of using only the brightest stars could be applied to the crowded regions in the core of ω Cen. On the other hand, the stellar density near the tidal radius (48, Ref. 93, 2010 edition) may be too low. Other GCs could also be used for calibration, but because of its overall high number of members and its large tidal radius, ω Cen is the best target.


Example of an autocalibration strategy

A possible autocalibration strategy for the WFIRST WFI could be modeled on the calibration described in Ref. 89 for the VIRCAM@VISTA detectors. The VIRCAM WFI comprises 16 2k×2k VIRGO detectors, for a total FoV of about 1.3×1.0 sq. deg, but the very large gaps between the chips bring down the effective FoV to 0.59 sq. deg. The calibration program (ESO proposal 488.L-0500(A), PI: Bellini) used a combination of small and large 5×5 dithers. Large dithers were used to cover the gaps between chips, monitor low-frequency distortions, and construct a single common reference system for all observations; the small dithers were included to enable independent modeling of the high-frequency residuals of the GD within each chip (more in Ref. 89). The choice of a 5×5 dither pattern was a compromise to obtain a sufficiently high-precision GD correction in a reasonable amount of telescope time. Extremely small dithers (from a few subpixels to a few pixels apart) are not strictly needed to characterize the PSF in well-populated star fields since nature distributes stars randomly with respect to the pixel boundaries (see also Ref. 94).

Figure 10 shows an example of a possible dither strategy for WFIRST following this plan. Figure 10(a) is a plan for a small 5×5 dither pattern that covers each WFI detector from corner to corner. In the top left panel, the black dots mark the center of each of the 25 dithers, with the detector layout of the central dither shown in red and other dithers in gray. The bottom left panel shows the resulting depth-of-coverage map on a logarithmic scale, with a maximum of 25 different images covering the same patch of sky. Most of the map is covered by at least 22 images, but only 12 to 15 of these come from the same chip so that the same star is typically imaged in 12 to 15 different chip locations. The 5×5 pattern never repeats the same shift along the X or Y direction, thus guaranteeing that the same stars will never fall on the same column or row, to minimize the impact of possible detector defects or degeneracies in the distortion solution.

Fig. 10

(a) Example of a small, 5×5 dither pattern that covers each WFI detector from corner to corner. (b) Example of a large, 9×5 dither pattern that covers the entire FoV from corner to corner. Top row: Dither pattern layout on-sky, with the center of each dither marked by a black dot, the WFI outline of the central dither shown in red, and outlines of the other dithers in gaey). Bottom row: Depth-of-coverage map (number of repeat observations as a function of position) for the assumed dither strategies, on a logarithmic scale. See Sec. 4.1 for details.


Figure 10(b) shows an example plan for a 9×5 pattern of large dithers that covers the entire WFI FoV from corner to corner. Because of the rectangular shape of the WFI FoV, nine dithers on the X axis are needed to cover the FoV with similar spacing to the five dithers along the Y axis. As for the small dithers, the layouts and centers of each pointing are shown on top and the resulting depth-of-coverage map is shown on the bottom. In this case, the layout results in a maximum of 39 different images at the center of the pattern.

The dither patterns shown in Fig. 10 all have the same telescope rotation angle, but in order to properly calibrate the skew terms of the distortion, a few observations of the same field at different roll angles would be highly beneficial. It is not obvious to suggest exactly how many of these rotated exposures should be taken, but sampling the full circle every 45 deg 60 deg should suffice. The total FoV covered by the large dither pattern in Fig. 10 allows for the central pointing to be rotated by any angle and still be fully within the covered region.

The proposed dither strategy makes use of 25 small dithers and 45 large dithers for a given filter, plus six or eight additional pointings (60 deg or 45 deg sampling, respectively) to constrain the skew terms, for a total of 76 to 78 distinct pointings. Experience calibrating the HST GD shows that convergence in the GD solution is achieved when stellar positions transformed from one image to another taken with a different pointing have rms residuals comparable to the stellar centroiding errors. Simulations to assess the precision of the GD correction as a function of the adopted dither strategy are ongoing to determine, among other things, whether fewer pointings than the example shown here could be sufficient.

An additional complication introduced by WFIRST’s large FoV is due to the use of tangent-plane projections: adopting the same projection point for images taken more than a few arcminutes apart results in significant positional transformation residuals. The GD calibration therefore has to be carried out on the celestial sphere rather than on any given tangent plane, adding an extralayer of complexity (see also Ref. 89).


Long-term monitoring of the GD solution

The EML survey is intended to characterize the PSF and fine-tune the GD solution of the WFIRST WFI, possibly including other sources of systematic effects such as intrapixel sensitivity variations. The current design of the EML survey employs small dithers and fixed rotation angle, suggesting that a satisfactory autocalibration of the GD using only EML survey images will be very hard to achieve. The catalog method could instead be used to fine-tune and monitor the GD solution, probably using Gaia as the reference catalog, but the lack of different telescope rotation angles may result in poorly constrained skew-term variations if these are present, as is the case for the WFC of the Advanced Camera for Surveys on HST. A preliminary investigation into the possibility of using EML-like simulated images of the Bulge, Gaia absolute stellar positions, and WebbPSF-based WFIRST PSF models (Bellini, in preparation) showed that: (i) stellar positions measured by PSF models that do not take into account jitter variations are significantly affected by pixel-phase-like errors (of the order of a few to several hundredths of a pixel); and (ii) the density of Gaia stars in the simulated EML survey field (about 5000 stars per chip) is adequate to solve for the GD.

Improved PSF models, either derived independently for each individual exposure or as a function of jitter rms, will address the pixel-phase issues and allow time-monitoring and fine-tuning of the GD solution. Jitter will vary with the reaction wheel speed, particularly at speeds that excite a structural vibration mode. Thus, the level of jitter is expected to change with time even on short time scales, but should have an rms well below 14 mas most of the time. Regardless, because of the time variability and because excitation of the telescope structure can cause line-of-sight jitter not sensed by the gyroscopes, the guide-star data will be extremely valuable for characterizing the jitter. It would therefore be very useful to downlink the reads of the guide windows together with each exposure, especially given the small amount of additional data involved. Furthermore, the jitter-dependent PSF models obtained for the filters employed in the EML survey will not necessarily apply equally well to other WFIRST filters. If this turns out to be the case, it would be helpful to include settling criteria that allow more stringent jitter rms constraints when images are taken for the purpose of calibrating and/or monitoring the GD. The current settling criteria include constraints on both position and angular rates, but additional criteria aimed at achieving better stability prior to calibrations would be helpful and should be investigated, since a smaller jitter rms implies smaller pixel-phase errors, helps in removing the degeneracy between centroiding accuracy and GD residuals, and would thereby make calibration more efficient.

Given that the EML survey will make use of only two of the WFI filters, it is important to note that filter elements are known to add significant contributions to the GD (e.g., Refs. 8687.88). Exposures taken with the other filters must be used to monitor the time dependency of the GD solution in those filters. In principle, all WFI exposures can, and probably will, be used for this purpose. This would typically be done with the catalog method, but autocalibration can be applied when properly dithered exposures are available so that any variation of the GD solution in all filters can be mitigated to the fullest extent possible.

The possibility of ground calibration of both the distortion and intrapixel sensitivity variations (see Sec. 4.2) should also be considered. As was found with HST, a successful astrometric calibration could reasonably be expected to improve WFIRST’s point-source localization, and therefore all astrometry-related measurements, by an order of magnitude. Such an improvement would multiply WFIRST’s reach in distance or velocity sensitivity for astrometry, thereby unlocking an entirely new space for discoveries.


Pixel-Level Effects


Quantum efficiency variations

Variations in the quantum efficiency (QE) within a single pixel can affect the accuracy of localization and therefore the astrometric precision. Reference 95 measured the intrapixel response function for the H2RG detectors to be flown on JWST, which are direct precursors of those planned for WFIRST. They found that the variations in the response per pixel [shown in Fig. 11(b)] appear to be mainly caused by redistribution of charges from pixel-to-pixel rather than by variations in pixel sensitivity. The most important effect in redistributing charge between pixels was the diffusion of charge to neighboring pixels, followed by interpixel capacitance (measured between 0% and 4%). They also find that the type of small defect visible in Fig. 11 occurs in roughly 10% of pixels. Additional testing of next-generation detectors more closely resembling those to be used for WFIRST is ongoing, but we expect that they will exhibit lower levels of variation.

Fig. 11

(a) Figure 4 of Ref. 95, showing the pixel response at 650 nm for an 8-pixel-square region of an H2RG detector. (b) Pixel offsets in a 128×128 region of an H2RG detector. (Figure courtesy Michael Shao.)



Placement error

To translate pixel-level effects into predictions for the precision of localization, Michael Shao’s group has made some preliminary measurements of the pixel placement error in H2RG detectors. The “effective” pixel position, which is defined as the location of the centroid of the QE within each pixel, was measured in these tests relative to an ideal coordinate system. Pixel offsets can have multiple causes, including the QE variations within a pixel discussed in Ref. 95 and in Sec. 4.2.1. These tests considered the pixel offset of a 128×128  pixel section of a H2RG detector and measured an rms 0.02  pixel offset error for a source measured in a single exposure, twice the value assumed in this work for single-exposure precision (see Sec. 1.1). Figure 11(a) shows the pixel offset in the X direction for the portion of the H2RG detector that was tested. By eye, the pixel placement errors appear to be random, in which case relative astrometry for two stars falling within the 128×128 region could be improved to better than 0.02 pixel by centroiding using the average of neighboring pixels. However, it is also possible that the neighboring 128×128 group of pixels are systematically offset from the group tested (this type of systematic shift of a group of pixels has been seen in CCDs), in which case the relative astrometric error cannot be improved by averaging. The tested region is so far not large enough to detect this type of larger-scale systematic error. Currently, the accuracy of pixel position measurement is estimated at 0.5%, but the group is in the process of more thoroughly verifying this as well as pushing toward 0.1% accuracy level. The main sources of error now are spurious fringes due to ghost reflections from, e.g., the dewar window, which have been minimized by tuning the laser over a broad range of wavelengths.

A 0.02-pixel position error corresponds to a single-image astrometric error of about 2 mas, so calibration of this effect will be important if WFIRST wants to deliver astrometry at the level of Gaia (10 to 100 times better) or even LSST precision (5 to 10 times better). Multiple dithered images and spatial scanning can be used to improve accuracy over the “raw” single-image error, depending on the brightness of the targets.


Ground-versus space-based calibration of subpixel effects

A limitation of calibrating subpixel effects once the telescope is in space is the issue of telescope jitter, which can change the PSF on a time scale of hours. The use of images of crowded fields to solve for subpixel errors in the detector relies on a stable PSF over a period of time long enough to collect sufficient photons to calibrate subpixel effects. The presence of time-variable telescope jitter prevents this from happening by many orders of magnitude. It is almost certain that the combination of jitter and photon noise will not allow on orbit calibration better than just assuming a perfect detector, given the measured 0.02 pixel errors in H2RG detectors.

It may not be as time-consuming to scan and calibrate this type of variation on the ground as indicated in Ref. 95. They measured intrapixel QE variations by scanning a spot image across 100  pixels using an extremely time-consuming process. For larger regions containing 104  pixels, the accuracy of measurements of pixel spacing using this approach will be limited by the positional accuracy of the translation stage used to perform the scan, which will likely be less accurate than the micrometer stage used by Ref. 95. Therefore, while this is a good approach to measure intrapixel QE, it is not sufficient for calibrating the dimensional accuracy of a large focal plane array for astrometry. The approach of scanning a spot across a pixel individually would be prohibitively slow for the WFIRST focal plane, which will contain 300 million pixels. Instead, the tests described in Sec. 4.2.2, which consider all pixels simultaneously, could potentially be scaled up to calibrate all the detectors before launch. Current estimates are that this scanning process can calibrate the focal plane array roughly 104105 times faster than the technique in Ref. 95. For a 300 megapixel camera such as the WFIRST WFI, such a calibration is estimated to take about 1 to 2 weeks, not including setup time.



Persistence of brightly illuminated regions is known to affect H2RG devices, especially in areas that have been saturated beyond the full-well depth (see the Wide Field Camera 3 Instrument Handbook,96 section 7.9.4). Characterizations of the persistence for both H2RGs and H4RGs are currently ongoing in several laboratories, which should establish a model for the persistence amplitude and decay time. Such a model can then be implemented to test the biases arising from persistence, which are relevant not only for precision astrometry but also for weak-lensing measurements in the HLS. It remains to be confirmed whether the model from ground-based testing is consistent with the persistence experienced in flight, for which several exposures of suitably bright stars should be sufficient.

If persistence is found to be problematic, a dark filter could be employed to block the incoming light during slews, or slew trajectories could be chosen to avoid bright stars. Experience from HST indicates that without a dark filter, persistence during slewing/tracking will be significant for stars brighter than fourth magnitude at the maximum slew rate. Given that there is 1 star brighter than sixth magnitude per 10 square degrees of sky, avoiding these sparsely distributed sources during slewing should be fairly straightforward.


Brighter-fatter effect

The brighter-fatter phenomenon is a well-known detector characteristic in CCDs whereby objects that are brighter have a larger PSF (i.e., are fatter; see e.g., Ref. 97). This complicates many PSF-dependent measurements and characterizations, which generally operate under the assumption that the PSF size is invariant due to changes in flux or exposure time. More recently, this effect has been observed in an H2RG detector similar to the H4RG detectors planned for WFIRST.98,99 There are currently laboratory efforts to understand and quantify this effect and corresponding efforts to software-based mitigation strategies. These efforts include experiment emulation for WFIRST at the Precision Projector Laboratory, a detector emulation and characterization facility designed to understand detector systematics at the demanding level required by weak lensing experiments.100 While we are at the early stages of quantitative studies of the brighter fatter effect on HxRG detectors, we expect there to be an induced astrometric bias, especially for undersampled point sources. Depending on the flux and position of the source relative to the pixel grid, this effect could redistribute flux from the PSF center to neighboring pixels asymmetrically, causing the observed centroid to shift by up to 1% of a pixel width, requiring mitigation procedures for sub-1% astrometric measurements.


Mitigation strategies

Spatial-scanning and diffraction-spike measurements (Sec. 2.4) distribute the photons over hundreds or thousands of pixels and are therefore more robust against pixel-level effects. Spatial scans are also robust against jitter (Sec. 4.6). A complete summary of requirements for an astrometric spatial scanning mode will be presented in an accompanying report (Casertano et al., in preparation).


Filters and Color-Dependent Systematics

Color-dependent systematic errors in the GD may have a significant impact on WFIRST. The bluest filters of the WFC3/UVIS detector on HST have color-dependent residuals of a few hundredths of a pixel (see Fig. 6 of Ref. 87). In that case, it is likely that the problem is due to a chromatic effect induced by fused-silica CCD windows within the optical system, which refract blue and red photons differently and have a sharp increase of the refractive index below 4000 Å. As a consequence, the F225W, F275W, and F336W filters are the most affected. WFIRST detectors are not based on silica CCDs, but the possibility of similar color-dependent systematic effects should be taken into account.

In addition, especially for wide-band filters, the GD affecting redder and bluer photons is likely to be slightly different due to diffraction and optical distortion. A test using WFC3/UVIS observations in the F606W filter of blue-horizontal-branch and red-giant-branch stars in ω Cen showed that the measured positions of blue and red stars are off by 0.002 UVIS pixels on average, with respect to their true positions. Filter-dependent residuals could therefore introduce small but still significant color-dependent systematic effects in wide-band WFIRST filters, including global (position-independent) effects. These filter-dependent systematics are expected to be stable over time, so their calibration should be straightforward.


Hysteresis in Readout Electronics

WFIRST should be aware that, on at least one existing instrument with HxRG detectors, the readout electronics are affected by significant hysteresis.88 The effect was first observed by J. Anderson in data taken by the High Acuity Wide field K-band Imager (HAWK-I) on the Very Large Telescope (VLT) while observing a standard field to help calibrate JWST and the HAWK-I detector itself, and later by Ref. 88 in all the fields observed with HAWK-I during the commissioning of the detector (see their Table 2). Figure 12 shows the results of Anderson’s investigation on the top and of Ref. 88 on the bottom. Plotted are the positional residuals of stars in the “distortion-free” master frame transformed into the raw reference system of the detector and the raw positions of the same stars in the raw reference system of the detector. The master frame was constructed iteratively using the autocalibration method to solve for the GD. In both cases, a periodic signal in the astrometric residuals is observed in the calibration data, believed to be caused by the alternating readout directions for the different amplifiers. A detailed description of the effects of hysteresis on astrometry and how to minimize them is presented in Sec. 5 of Ref. 88. For the H2RG chips of the HAWK-I detectors, the hysteresis effect produces ±0.035  pixels (or ±3.7  mas) of systematic error, which can be easily modeled and corrected for. The WFIRST detectors will have 32 amplifiers, and scientific full-frame images will make use of all of them. If similar hysteresis effects to those of the H2RG of the HAWK-I camera are present as well in the WFIRST detectors, they will likely be easily minimized as was done for HAWK-I.

Fig. 12

Hysteresis in the HAWK-I detectors on the VLT, precursors to those planned for WFIRST. (a) Results from an investigation by Anderson. Four horizontal strips of the detector were analyzed, as shown on the far left. The central columns show the distortion as a function of x position in each strip; the vertical scale is in pixels (1 pix = 100 mas). The raw distortion is shown to the left, whereas to the right are shown the residuals after subtraction of a smooth global polynomial, revealing a high-frequency periodic signal. On the far right are close-up views of the astrometric effect. The top row shows one polynomial-subtracted stripe; while the central row shows the residuals phase-wrapped by 128 columns. The resulting step function (left center row) has a half-amplitude of 0.035 pixel; the right center row shows residuals after subtraction. The unwrapped, corrected residuals are shown in the bottom row. (b) Results from Ref. 88. The left-hand panel shows δx residuals for each of the four HAWK-I detectors after the polynomial correction is applied. The right-hand panel shows a periodogram, with a period of 128 pixels, containing all the points plotted on the left, with the median shown as a solid red line. The dashed red lines are at 0 and ±0.05 HAWK-I pixels (about ±5.3  mas).




A small amount of extraattention to scheduling can give great payout for astrometric measurements. As suggested by Eq. (2), the ideal scheme for measuring PMs is to space revisits to the same field as evenly as possible over the longest possible time baseline, in the interest of minimizing both random and systematic errors. Maximizing the time baseline reduces the random error while increasing the number of epochs protects against systematics, which are often time-dependent with unknown correlations. Here, we discuss a few considerations for the core science and guest observing components of the WFIRST mission.


High-latitude survey

The observation of HLS fields is planned to take place over a 5-year baseline, but the detailed schedule is not yet finalized. Starting the HLS with an initial exposure of each field in year 1 and reobserving at least once in each field as late as year 5 would produce an astrometric survey that extends Gaia’s PM precision (25 to 50  μasyr1) to stars six magnitudes fainter than Gaia can reach. However, the current range of schedules considered for the HLS can affect this projected precision by factors of a few. We recommend breaking ties between otherwise equivalent programs by considering the time-distribution of revisits.


Exoplanet microlensing survey

The EML survey cycles through its 10 target fields at 52-s intervals using the wide filter, with one set in the blue filter every 12 h. Currently, there are six total seasons planned, half in spring and half in fall, in order to measure relative parallaxes with a similar level of accuracy to that of PMs. Some observations will be front-weighted at the start of the survey, but at least one season should be planned for the end of the mission to obtain the longest possible time baselines. Pointing and solar-panel orientation requirements are likely to separate seasons by almost exactly 6 months, which may cause some complications in calibrating time-dependent PSF effects (see Sec. 4.1).


Guest observing

Much of the science described in Sec. 2 will likely be carried out through GO and GI programs. It is therefore crucial that the benefit of the astrometric requirements for the two core science programs listed above be made available to GO/GI programs as well (by providing the calibration information for these programs to GOs/GIs and allowing for multiple astrometric calibrations within the archive) in order to achieve the promised precision. This need is discussed in depth in Sec. 4.8.

Since WFIRST is an IR telescope, it can provide astrometry for regions of the Galaxy, such as the disk plane and bulge, that are completely inaccessible to the current generation of optical astrometric instruments. The EML survey will cover one such region, but to take full advantage of WFIRST’s astrometric capabilities, it will be crucial to allow multi-year proposals from GOs in order to optimize for PM baselines, as it is currently done for HST.



The current WFIRST requirements impose a maximum jitter rms of 14  mas, far larger than HST’s jitter rms of 2– to 5 mas. Jitter of the size allowed by WFIRST’s requirements could have a significant impact on the shape of the PSF. Tightening this limit would clearly result in better astrometry, but would also translate into a higher cost for some of the telescope components and is thus beyond the mandate of our working group. We have thus considered how to mitigate the effect of jitter at the level set by existing requirements.

Preliminary simulations of WFIRST’s geometric-distortion corrections based on EML-like Bulge images,92 which make use of time-constant, spatially variable library PSF models adapted from WebbPSF for WFIRST,101 show that there is significant degeneracy (at the 0.02 to 0.05 pixel level) between the achievable geometric-distortion correction and pixel-phase errors in stellar positions.

One way to break the degeneracy is to spatially perturb the library PSF in each individual exposure, so as to tailor the library PSF to the particular jitter status of each image. Because of the GD, jitter-induced PSF variations are expected to affect the PSF of different WFI detectors in a different way. Using this workaround to calibrate the GD would require images with roughly 20 to 40 thousand bright and isolated sources, homogeneously distributed across the WFI FoV, in order to map local PSF variations on scales of 500 pixels or so.

Another possibility would be to exploit the enormous number of images in the EML survey to map PSF variations at different jitter rms values, thereby creating a jitter-sampled set of spatially variable model PSFs. This technique will be limited by the use of only two filters for the EML survey since it is unclear whether the models generated by these two filters will apply equally well to the rest of the filterset.


Data Management

Both pointed and spatially scanned astrometric observations will be constrained by the data downlink rate. For the WFIRST reference mission at L2, the downlink rate is estimated to be about 1.3 TB day1, barring addition of extra ground stations.102 Particularly for pointed astrometric observations, this limits the number of reads per exposure that can be downloaded. The current plan is to allow configuration of the options for averaging and saving reads, similar to JWST. For astrometry, two options are particularly important:

  • Downloads of GS postage stamps are crucial and inexpensive (0.1% overhead!) for PSF jitter correction.

  • The ability to download every read with no coaddition for at least part of the field is especially desirable for spatial scanning observations, as discussed in Sec. 2.4.


High-Level Data Products and Archive

The combination of a wide FoV, the plentiful availability of reference stars with extremely accurate astrometry from Gaia, and the resolution and stability of a space-based platform makes it possible for WFIRST to achieve extremely accurate absolute astrometry: better than 100  μas for essentially all WFC imaging products (see Sec. 3 and Ref. 75). Any field observed more than once by the WFI is thus a potential astrometric field, providing the community with a wealth of opportunities for high-precision astrometry studies in many different domains, including many of the science cases described in Sec. 2 of this report. However, for this potential to be realized for all users throughout the mission, we strongly recommend that the data processing pipeline and the archive incorporate from the outset the necessary elements to obtain, propagate, and maintain in practice the astrometric accuracy that the mission characteristics make possible in principle.

Many of the features that make WFIRST an excellent astrometric instrument, notably the requirements for excellent PSF modeling and thermal stability, are already addressed by the core science programs (the HLS and EML surveys). These requirements are summarized in Appendix A. Given their use of repeated visits to the same fields, both of these programs also offer an opportunity to produce excellent astrometry for all observed stars in their footprints with little extra analysis. We recommend that derived PMs should be provided as part of the object catalogs for both surveys once multiple epochs have been observed. In the HLS particularly, a further cross-match of stars to the LSST catalog would extend the LSST survey into the IR regime for the region covered by HLS, identify variable stars that can be used as standard candles (particularly for RR Lyrae stars, as the period-luminosity relation is much tighter in the IR) and allow for cross-validation of PMs.

Specifically, we recommend that the mission considers the possibility of achieving the following goals:

  • 1. The initial (a priori) astrometric information for each image should be based on GS positions known to Gaia accuracy, together with an accurate WFC GD model and the analysis of guide-star window data to extract accurate instantaneous GS positions.

  • 2. The WFC GD model should be verified, and updated if needed, with sufficient frequency to maintain no worse than 100  μas precision, on the basis of on-orbit data on the geometric stability of the WFC focal plane.

  • 3. The pipeline to generate level 2 products should include an a posteriori alignment step based on cross-matching sources found in each image with the Gaia catalog; this information should be incorporated as an alternate (preferred) astrometric solution in the image metadata.

  • 4. The accurate astrometry thus determined should be propagated to level 3 and level 4 data products.

  • 5. The Archive should have the ability to retain and distribute multiple astrometric solutions for each data product, together with their pedigree and uncertainty. The Archive could also incorporate community-provided astrometric solutions if deemed useful.

  • 6. The data management system should have the ability to update the astrometric information for higher-level products when the astrometric solution for the contributing data products is updated.

These recommendations are based on part on our experience with HST. The astrometric accuracy available for HST data early in the mission was originally limited by the quality of the GS positions; therefore, modest effort was placed into improving other components of the astrometric fidelity, such as the knowledge and time evolution of the relative positions of instruments and guiders. Now that substantially better positional accuracy is available for guide and reference stars, retrofitting the HST pipeline and archive to improve the final absolute astrometric accuracy of HST-processed data has proven complex and resource-intensive. We recognize that several of these recommendations go beyond the current science requirements and may exceed the baseline capabilities of the mission, as currently planned. However, incorporating these considerations, to the extent possible, into the design of the WFIRST data processing and archive systems will greatly improve the quality and accessibility of mission data for astrometry, improving science outcomes in this area and ultimately reducing total development costs when compared with adding similar capabilities at a later time.


Appendix A: Astrometry-Relevant Requirements for Core Science Programs

Here, we summarize the current state of relevant requirements for the two core-science surveys WFIRST will carry out, broken down into three categories: basic science requirements (BSRs), requirements for the high-latitude imaging survey (HLIS), and requirements for the EML survey.


A.1 Astrometry with the High-Latitude Survey

The current requirements related to astrometry, as of July 2017, include the following:

  • BSR 2: WFIRST WFI shall measure shapes of galaxies at z=0 to 2 in at least two bands and fluxes in at least four bands for photometric redshifts, at a depth equivalent to a 5-sigma point source detection at AB magnitude J<26.9 or H<26.7, with photometric accuracy of 1% and with rms uncertainties (in the shape measurement filters only) below 103 in the PSF second moment and below 5×104 in the PSF ellipticity, in the HLS imaging survey.

  • HLIS 7: Obtain photometry, position, and shape measurements of galaxies in three filters (J, H, and F184), and photometry and position measurements in one additional color filter (Y; only for photo-z).

  • HLIS 8: Obtain S/N18 (matched filter detection significance) per shape/color filter for galaxy effective radius reff=180 mas and AB mag = 24.7/24.6/24.1 (J/H/F184).

  • HLIS 9: Determine PSF second moment to a relative error of 9.3×104  rms (shape/color filters only).

  • HLIS 10: Determine PSF ellipticity to 4.7×104  rms (shape/color filters only).

  • HLIS 11: The 50% encircled energy (EE50) radius of the PSF 0.12 (Y band), 0.12 (J), 0.14 (H), or 0.13 (F184) arcsec.

The reference dither pattern for the HLS is a set of 3 to 4 dithers of a size intended to cover the chip gaps, repeated to tile the field and in each of four filters: Y, J, H, and F184. A second pass over each field follows 6 months later at a different roll angle.

The projected, approximate bright and faint point-source limits of the HLS are summarized in Table 4. The faint limits are as stated in the requirements above. The bright limits were estimated based on the statement in the reference mission that pixels in the HLS will be read nondestructively every 5.4 s, and assuming that any pixels that saturate before the fourth such read will be hard-saturated (see the discussion in Ref. 75). Using the GS ETC, assuming 25% of light in the central pixel for all filters and 65,000 electrons as the saturation level, the values listed in the table give the approximate bright limit (probably accurate to within 0.3 to 0.4 magnitudes).

Table 4

Bright and faint point-source limits for the WFIRST HLS (AB magnitudes), based on the requirements described in Sec. 5.1 A.1.

Bright limit15.515.615.715.3
Faint limit26.926.7


A.2 Astrometry with the Exoplanet MicroLensing Survey

The current astrometry-related requirements being discussed (as of June 29, 2017) for the EML survey are:

  • EML 8: Relative photometric measurements in the primary microlensing filter that have a statistical S/N of 100 per exposure for a HAB=21.6 star.

  • EML 14: The EE50 radius of the PSF in the wide filter shall be <0.15.

  • EML 19: The relative astrometric measurements shall have a statistical precision of 1 mas per measurement for HAB=21.6.

  • EML 20: Relative astrometric measurements will have systematic precision of 10  μas over the full survey (stretch goal of 3  μas).

Currently, both narrow (two-pixel wide) and large (10-in. wide) possibilities for dithering are being explored for this core project. From the perspective of astrometric calibration, large dithers are crucial to accurately measure skew and kurtosis in the wings of the PSF.

Parallaxes and PMs over the bulge field are part of the mission of this program to characterize the masses of the star-planet pairs that will be discovered. The survey is therefore requesting the first two (spring/fall) and last two bulge observing seasons over the full time-baseline of the mission. This is also optimal for general astrometry in the bulge fields but may pose problems for understanding long-term variations in the PSF (on timescales of a year or so) prior to the end of the mission.


Appendix B: Typical Astrometry-Related Queries to Object Catalog

These queries were submitted to the archive working group as part of their “20-questions” use case development.

  • AWG-1: Give me positions, IR magnitudes, PMs, and associated uncertainties of all stars in the HLS or EML survey within a color-magnitude box/isochrone cutout.

  • AWG-1a: Also return LSST optical magnitudes, PMs, and associated uncertainties for the selected stars.

  • AWG-2: Return positions, magnitudes, PMs, distances, and associated uncertainties of stars in a specified field that LSST identifies as standard-candle variable stars (e.g., RR Lyrae).

  • AWG-3: Run a group finder (or other analysis software) I provide on the above data.

  • AWG-4: Give me all frames from any observing program and any associated calibration information or data flags that intersect a defined region on the sky (i.e., it would be great to be able to rereduce data from different observations to measure PMs with more frames/longer time baseline).

  • AWG-5: Give me positions, parallaxes, PMs, and associated uncertainties of all Gaia stars within a WFIRST pointing above a mag threshold.

  • AWG-6: Extension of above: Give me predicted positions and uncertainties of Gaia stars within a WFIRST pointing above a mag threshold at some observation date+time.

  • AWG-7: Give me positions, PMs, magnitudes, and associated uncertainties for all XX-type stars within YY pc from the Sun in this ZZ WFIRST pointing.


This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This work presents results from the European Space Agency (ESA) space mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia MultiLateral Agreement (MLA). The Gaia mission website is The Gaia archive website is RES was supported in part by an NSF Astronomy and Astrophysics Postdoctoral Fellowship under Grant No. AST-1400989. The authors thank Charles Shapiro (JPL) for useful discussions, and the anonymous referees for comments that improved the paper.



D. Spergel, N. Gehrels and J. Breckinridge, “WFIRST-2.4: what every astronomer should know,” (2013). Google Scholar


Gaia Collaboration, “The Gaia mission,” Astron. Astrophys., 595 A1 (2016). AAEJAF 0004-6361 Google Scholar


Gaia Collaboration, “Gaia data release 2. Summary of the contents and survey properties,” Astron. Astrophys., 616 A1 (2018). AAEJAF 0004-6361 Google Scholar


A. M. Ghez et al., “Measuring distance and properties of the Milky Way’s central supermassive black hole with Stellar orbits,” Astrophys. J., 689 1044 –1062 (2008). ASJOAB 0004-637X Google Scholar


A. Gould et al., “WFIRST Ultra-precise astrometry II: asteroseismology,” J. Korean Astron. Soc., 48 93 –104 (2015). CHACDE 0253-3065 Google Scholar


Jr. M. W. Hosek et al., “The arches cluster: extended structure and tidal radius,” Astrophys. J., 813 27 (2015). ASJOAB 0004-637X Google Scholar


A. W. McConnachie, “The observed properties of dwarf galaxies in and around the local group,” Astron. J., 144 4 (2012). Google Scholar


N. Kallivayalil et al., “Third-epoch magellanic cloud proper motions. I. Hubble Space Telescope/WFC3 data and orbit implications,” Astrophys. J., 764 161 (2013). ASJOAB 0004-637X Google Scholar


S. T. Sohn et al., “The space motion of Leo I: Hubble Space telescope proper motion and implied orbit,” Astrophys. J., 768 139 (2013). ASJOAB 0004-637X Google Scholar


S. T. Sohn et al., “Space motions of the dwarf spheroidal galaxies Draco and sculptor based on HST proper motions with a ∼10 yr time baseline,” Astrophys. J., 849 93 (2017). ASJOAB 0004-637X Google Scholar


S. T. Sohn et al., “Absolute Hubble Space Telescope proper motion (HSTPROMO) of distant Milky Way globular clusters: Galactocentric space velocities and the Milky Way mass,” Astrophys. J., 862 52 (2018). ASJOAB 0004-637X Google Scholar


Gaia Collaboration, “Gaia data release 2. Kinematics of globular clusters and dwarf galaxies around the Milky Way,” Astron. Astrophys., 616 A12 (2018). AAEJAF 0004-6361 Google Scholar


J. Li et al., “Detecting the Sagittarius stream with LAMOST DR4 M giants and Gaia DR2,” Astrophys. J., 874 138 (2019). ASJOAB 0004-637X Google Scholar


G. M. Eadie and W. E. Harris, “Bayesian mass estimates of the Milky Way: the dark and light sides of parameter assumptions,” Astrophys. J., 829 108 (2016). ASJOAB 0004-637X Google Scholar


The Theia Collaboration, “Theia: faint objects in motion or the new astrometry frontier,” (2017). Google Scholar


S. T. Sohn et al., “Hubble Space Telescope proper motions of individual stars in stellar streams: orphan, Sagittarius, Lethe, and the new ‘parallel’ stream,” (2016). Google Scholar


N. Kallivayalil et al., “A Hubble astrometry initiative: laying the foundation for the next-generation proper-motion survey of the local group,” (2015). Google Scholar


R. P. van der Marel and N. Kallivayalil, “Third-epoch magellanic cloud proper motions. II. The large magellanic cloud rotation field in three dimensions,” Astrophys. J., 781 121 (2014). ASJOAB 0004-637X Google Scholar


R. E. Sanderson et al., “New views of the distant stellar halo,” Mon. Not. R. Astron. Soc., 470 5014 –5031 (2017). MNRAA4 0035-8711 Google Scholar


J. J. Bochanski et al., “The most distant stars in the Milky Way,” Astrophys. J., 790 L5 (2014). ASJOAB 0004-637X Google Scholar


A. J. Deason et al., “The cold veil of the Milky Way stellar halo,” Mon. Not. R. Astron. Soc., 425 2840 –2853 (2012). MNRAA4 0035-8711 Google Scholar


J. J. Bochanski et al., “Hunting the most distant stars in the Milky Way: methods and initial results,” Astron. J., 147 76 (2014). ANJOAA 0004-6256 Google Scholar


B. Sesar et al., “The 100 kpc distant spur of the Sagittarius stream and the outer Virgo over density, as seen in PS1 RR Lyrae stars,” Astrophys. J., 844 L4 (2017). ASJOAB 0004-637X Google Scholar


B. Diemer and A. V. Kravtsov, “Dependence of the outer density profiles of halos on their mass accretion rate,” Astrophys. J., 789 1 (2014). ASJOAB 0004-637X Google Scholar


V. Belokurov et al., “Precession of the Sagittarius stream,” Mon. Not. R. Astron. Soc., 437 116 –131 (2014). MNRAA4 0035-8711 Google Scholar


W. R. Brown, “Hypervelocity stars,” Annu. Rev. Astron. Astrophys., 53 15 –49 (2015). ARAAAJ 0066-4146 Google Scholar


J. S. Bullock and K. V. Johnston, “Tracing galaxy formation with Stellar Halos. I. Methods,” Astrophys. J., 635 931 –949 (2005). ASJOAB 0004-637X Google Scholar


R. A. Ibata, G. F. Lewis and N. F. Martin, “Feeling the pull: a study of natural Galactic accelerometers. I. Photometry of the delicate stellar stream of the Palomar 5 globular cluster,” Astrophys. J., 819 1 (2016). ASJOAB 0004-637X Google Scholar


J. Bovy, D. Erkal and J. L. Sanders, “Linear perturbation theory for tidal streams and the small-scale CDM power spectrum,” Mon. Not. R. Astron. Soc., 466 628 –668 (2017). MNRAA4 0035-8711 Google Scholar


A. G. Riess et al., “Parallax beyond a kiloparsec from spatially scanning the wide field camera 3 on the Hubble Space Telescope,” Astrophys. J., 785 161 (2014). ASJOAB 0004-637X Google Scholar


S. Casertano et al., “Parallax of galactic cepheids from spatially scanning the wide field camera 3 on the Hubble Space Telescope: the case of SS Canis Majoris,” Astrophys. J., 825 11 (2016). ASJOAB 0004-637X Google Scholar


P. Melchior, D. Spergel and A. Lanz, “In the crosshair: astrometric exoplanet detection with WFIRST’s diffraction spikes,” (2018). Google Scholar


D. Pourbaix and A. Jorissen, “Re-processing the Hipparcos transit data and intermediate astrometric data of spectroscopic binaries. I. Ba, CH and Tc-poor S stars,” Astron. Astrophys. Suppl. Ser., 145 161 –183 (2000). AAESB9 0365-0138 Google Scholar


A. R. Wetzel et al., “Reconciling Dwarf galaxies with ΛCDM cosmology: simulating a realistic population of satellites around a Milky Way-mass galaxy,” Astrophys. J., 827 L23 (2016). ASJOAB 0004-637X Google Scholar


B. Paczyński and K. Z. Stanek, “Galactocentric distance with the optical gravitational Lensing experiment and HIPPARCOS red clump stars,” Astrophys. J., 494 L219 –L222 (1998). ASJOAB 0004-637X Google Scholar


P. Marigo et al., “Evolution of asymptotic giant branch stars. II. Optical to far-infrared isochrones with improved TP-AGB models,” Astron. Astrophys., 482 883 –905 (2008). AAEJAF 0004-6361 Google Scholar


R. E. Sanderson et al., “Synthetic Gaia surveys from the FIRE cosmological simulations of Milky-Way-mass galaxies,” (2018). Google Scholar


D. J. Schlegel, D. P. Finkbeiner and M. Davis, “Maps of dust infrared emission for use in estimation of reddening and cosmic microwave background radiation foregrounds,” Astrophys. J., 500 525 –553 (1998). ASJOAB 0004-637X Google Scholar


A. G. A. Brown, “Pygaia: python toolkit for basic gaia data simulation, manipulation, and analysis,” Google Scholar


M. Romero-Gómez et al., “The analysis of realistic stellar Gaia mock catalogues - I. Red clump stars as tracers of the central bar,” Mon. Not. R. Astron. Soc., 447 218 –233 (2015). MNRAA4 0035-8711 Google Scholar


W. Clarkson et al., “Stellar proper motions in the galactic bulge from deep Hubble Space Telescope ACS WFC photometry,” Astrophys. J., 684 1110 –1142 (2008). ASJOAB 0004-637X Google Scholar


D. Minniti, M. Zoccali, “The galactic bulge: a review,” Formation and Evolution of Galaxy Bulges, 245 323 –332 Cambridge University Press, Cambridge (2008). Google Scholar


S. Vásquez et al., “3D kinematics through the X-shaped Milky Way bulge,” Astron. Astrophys., 555 A91 (2013). AAEJAF 0004-6361 Google Scholar


J. Hosek et al., “The unusual initial mass function of the arches cluster,” Astrophys. J., 870 44 (2019). ASJOAB 0004-637X Google Scholar


A. Stolte et al., “The proper motion of the arches cluster with keck laser-guide star adaptive optics,” Astrophys. J., 675 1278 –1292 (2008). ASJOAB 0004-637X Google Scholar


W. I. Clarkson et al., “Proper motions of the arches cluster with Keck laser guide star adaptive optics: the first kinematic mass measurement of the arches,” Astrophys. J., 751 132 (2012). ASJOAB 0004-637X Google Scholar


E. Agol and M. Kamionkowski, “X-rays from isolated black holes in the Milky Way,” MNRAS, 334 553 –562 (2002). Google Scholar


J. R. Lu et al., “A search for stellar-mass black holes via astrometric microlensing,” Astrophys. J., 830 41 (2016). ASJOAB 0004-637X Google Scholar


E. Carretta et al., “Na-O anticorrelation and HB. VII. The chemical composition of first and second-generation stars in 15 globular clusters from GIRAFFE spectra,” Astron. Astrophys., 505 117 –138 (2009). AAEJAF 0004-6361 Google Scholar


A. P. Milone et al., “The infrared eye of the wide-field camera 3 on the Hubble Space telescope reveals multiple main sequences of very low mass stars in NGC 2808,” Astrophys. J., 754 L34 (2012). ASJOAB 0004-637X Google Scholar


R. G. Gratton, E. Carretta and A. Bragaglia, “Multiple populations in globular clusters. Lessons learned from the Milky Way globular clusters,” Astron. Astrophys. Rev., 20 50 (2012). AASREB 0935-4956 Google Scholar


A. Bellini et al., “The state-of-the-art HST astro-photometric analysis of the core of ω Centauri. I. The catalog,” Astrophys. J., 842 6 (2017). ASJOAB 0004-637X Google Scholar


A. P. Milone et al., “The HST large programme on ω Centauri - I. Multiple stellar populations at the bottom of the main sequence probed in NIR-Optical,” Mon. Not. R. Astron. Soc., 469 800 –812 (2017). MNRAA4 0035-8711 Google Scholar


G. Piotto et al., “The Hubble Space Telescope UV legacy survey of galactic globular clusters. I. Overview of the project and detection of multiple stellar populations,” Astron. J., 149 91 (2015). ANJOAA 0004-6256 Google Scholar


A. Bellini et al., “The state-of-the-art HST astro-photometric analysis of the core of ω Centauri. III. The main sequence’s multiple populations galore,” Astrophys. J., 844 164 (2017). ASJOAB 0004-637X Google Scholar


G. Piotto et al., “Metallicities on the double main sequence of ω Centauri imply large helium enhancement,” Astrophys. J., 621 777 –784 (2005). ASJOAB 0004-637X Google Scholar


G. Piotto et al., “A triple main sequence in the globular cluster NGC 2808,” Astrophys. J., 661 L53 –L56 (2007). ASJOAB 0004-637X Google Scholar


J. E. Norris and G. S. Da Costa, “The giant branch of Omega Centauri. IV. Abundance Patterns based on echelle spectra of 40 red giants,” Astrophys. J., 447 680 (1995). ASJOAB 0004-637X Google Scholar


C. I. Johnson and C. A. Pilachowski, “Chemical abundances for 855 giants in the globular cluster Omega Centauri (NGC 5139),” Astrophys. J., 722 1373 –1410 (2010). ASJOAB 0004-637X Google Scholar


E. Carretta et al., “Detailed abundances of a large sample of giant stars in M 54 and in the Sagittarius nucleus,” Astron. Astrophys., 520 A95 (2010). AAEJAF 0004-6361 Google Scholar


A. P. Milone et al., “The Hubble Space Telescope UV legacy survey of galactic globular clusters - II. The seven stellar populations of NGC 7089 (M2)*,” Mon. Not. R. Astron. Soc., 447 927 –938 (2015). MNRAA4 0035-8711 Google Scholar


A. Bellini et al., “Hubble Space Telescope proper motion (HSTPROMO) catalogs of galactic globular clusters. I. Sample selection, data reduction, and NGC 7078 results,” Astrophys. J., 797 115 (2014). ASJOAB 0004-637X Google Scholar


H. B. Richer et al., “A dynamical signature of multiple stellar populations in 47 Tucanae,” Astrophys. J., 771 L15 (2013). ASJOAB 0004-637X Google Scholar


J. Anderson and R. P. van der Marel, “New limits on an intermediate-mass black hole in Omega Centauri. I. Hubble Space Telescope photometry and proper motions,” Astrophys. J., 710 1032 –1062 (2010). ASJOAB 0004-637X Google Scholar


A. Bellini et al., “The Hubble Space Telescope UV legacy survey of galactic globular clusters: the internal kinematics of the multiple stellar populations in NGC 2808,” Astrophys. J., 810 L13 (2015). ASJOAB 0004-637X Google Scholar


A. D’Ercole et al., “Formation and dynamical evolution of multiple stellar generations in globular clusters,” Mon. Not. R. Astron. Soc., 391 825 –843 (2008). MNRAA4 0035-8711 Google Scholar


A. Bellini et al., “The HST large programme on ω Centauri. II. Internal kinematics,” Astrophys. J., 853 86 (2018). ASJOAB 0004-637X Google Scholar


M. Trenti and R. van der Marel, “No energy equipartition in globular clusters,” Mon. Not. R. Astron. Soc., 435 3272 –3282 (2013). MNRAA4 0035-8711 Google Scholar


P. Bianchini et al., “A novel look at energy equipartition in globular clusters,” Mon. Not. R. Astron. Soc., 458 3644 –3654 (2016). MNRAA4 0035-8711 Google Scholar


M. Libralato et al., “Hubble Space Telescope proper motion (HSTPROMO) catalogs of galactic globular cluster. VI. Improved data reduction and internal-kinematic analysis of NGC 362,” Astrophys. J., 861 99 (2018). ASJOAB 0004-637X Google Scholar


L. R. Bedin et al., “Color-magnitude diagram and luminosity function of M4 near the hydrogen-burning limit,” Astrophys. J., 560 L75 –L78 (2001). ASJOAB 0004-637X Google Scholar


A. Dieball et al., “Deep near-IR observations of the globular cluster M4: hunting for brown Dwarfs,” Astrophys. J., 817 48 (2016). ASJOAB 0004-637X Google Scholar


A. Bellini et al., “2MASS/Gaia as an absolute astrometric reference frame,” Baltimore (2017). Google Scholar


E. Nelan et al., “WFIRST-AFTA guide stars and the 2MASS point source catalog,” Baltimore (2015). Google Scholar


L. Colangeli, “Extended life for ESA’s science missions,” Google Scholar


J. H. J. de Bruijne, “Science performance of Gaia, ESA’s space-astrometry mission,” Astrophys. Space Sci., 341 31 –41 (2012). APSSBE 0004-640X Google Scholar


F. Arenou et al., “Gaia data release 2. Catalogue validation,” Astron. Astrophys., 616 A17 (2018). AAEJAF 0004-6361 Google Scholar


L. Lindegren et al., “Gaia data release 2. The astrometric solution,” Astron. Astrophys., 616 A2 (2018). AAEJAF 0004-6361 Google Scholar


J. Anderson and I. R. King, “An improved distortion solution for the Hubble Space Telescope’s WFPC2,” Publ. Astron. Soc. Pac., 115 113 –131 (2003). PASPAU 0004-6280 Google Scholar


J. Anderson and I. R. King, “Multi-filter PSFs and distortion corrections for the HRC,” Baltimore (2004). Google Scholar


J. Anderson et al., “Ground-based CCD astrometry with wide field imagers. I. Observations just a few years apart allow decontamination of field objects from members in two globular clusters,” Astron. Astrophys., 454 1029 –1045 (2006). AAEJAF 0004-6361 Google Scholar


J. Anderson and I. R. King, “PSFs, photometry, and astronomy for the ACS/WFC,” Baltimore (2006). Google Scholar


A. Bellini and L. R. Bedin, “Astrometry and photometry with HST WFC3. I. Geometric distortion corrections of F225W, F275W, F336W bands of the UVIS channel,” Publ. Astron. Soc. Pac., 121 1419 –1428 (2009). PASPAU 0004-6280 Google Scholar


A. Bellini and L. R. Bedin, “Ground-based CCD astrometry with wide field imagers. IV. An improved geometric-distortion correction for the blue prime-focus camera at the LBT,” Astron. Astrophys., 517 A34 (2010). AAEJAF 0004-6361 Google Scholar


A. Bellini, J. Anderson and L. R. Bedin, “Astrometry and photometry with HST WFC3. II. improved geometric-distortion corrections for 10 filters of the UVIS channel,” Publ. Astron. Soc. Pac., 123 622 –637 (2011). PASPAU 0004-6280 Google Scholar


M. Libralato et al., “Ground-based astrometry with wide field imagers. V. Application to near-infrared detectors: HAWK-I@VLT/ESO,” Astron. Astrophys., 563 A80 (2014). AAEJAF 0004-6361 Google Scholar


M. Libralato et al., “High-precision astrometry with VVV - I. An independent reduction pipeline for VIRCAM@VISTA,” Mon. Not. R. Astron. Soc., 450 1664 –1673 (2015). MNRAA4 0035-8711 Google Scholar


D. Massari et al., “Astrometry with MCAO: HST-GeMS proper motions in the globular cluster NGC 6681,” Astron. Astrophys., 595 L2 (2016). AAEJAF 0004-6361 Google Scholar


E. Dalessandro et al., “GeMS/GSAOI photometric and astrometric performance in dense stellar fields,” Astrophys. J., 833 111 (2016). ASJOAB 0004-637X Google Scholar


A. Bellini, “Will Gaia be precise enough to solve for the geometric distortion of the WFI?,” Baltimore (2018). Google Scholar


W. E. Harris, “A catalog of parameters for globular clusters in the Milky Way,” Astron. J., 112 1487 (1996). ANJOAA 0004-6256 Google Scholar


J. Anderson and I. R. King, “Toward high-precision astrometry with WFPC2. I. Deriving an accurate point-spread function,” Publ. Astron. Soc. Pac., 112 1360 –1382 (2000). PASPAU 0004-6280 Google Scholar


T. Hardy, C. Willot and J. Pazder, “Intra-pixel response of the new jwst infrared detector arrays,” Proc. SPIE, 9154 91542D (2014). PSISDG 0277-786X Google Scholar


L. Dressel, “Wide field camera 3 instrument handbook, Version 11.0,” Baltimore (2019). Google Scholar


D. Gruen et al., “Characterization and correction of charge-induced pixel shifts in DECam,” J. Instrum., 10 C05032 (2015). Google Scholar


N. Rowlands, C. Midwinter and G. Warner, “The impact of the brighter-fatter effect on the performance of the JWST fine guidance sensor,” Proc. SPIE, 10709 107091K (2018). PSISDG 0277-786X Google Scholar


A. A. Plazas et al., “Laboratory measurement of the brighter-fatter effect in an H2RG infrared detector,” Publ. Astron. Soc. Pac., 130 065004 (2018). PASPAU 0004-6280 Google Scholar


C. Shapiro et al., “Precision projector laboratory: detector characterization with an astronomical emulation testbed,” (2018). Google Scholar


D. Spergel et al., “Wide-field infrarred survey telescope-astrophysics focused telescope assets WFIRST-AFTA 2015 report,” (2015). Google Scholar

Biographies of the authors are not available.

© The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Received: 19 December 2018; Accepted: 24 September 2019; Published: 29 October 2019


CHORUS Article. This article was made freely available starting 28 October 2020

Back to Top