3.1. Decoding the physics of galaxy evolution by mapping the cosmic web

Galaxy properties should correlate with the underlying dark matter: the masses of the halos, their spins, and their positions within the cosmic web. Today, we have only approximate estimates of the stellar-mass/halo-mass ratio, with rather large uncertainties, from abundance matching, clustering, and WL (Behroozi et al. Reference Behroozi2013; 368 Kravtsov et al. Reference Kravtsov2013; Hudson et al. Reference Hudson2015). Direct clustering measurements allow us to estimate the stellar-mass/halo-mass ratio as a function of galaxy properties, which is impossible to do via abundance matching. It is also clear that the stellar-mass/halo-mass ratio does not capture the entire influence of dark matter on galaxy properties, as the poorly understood ‘conformity’ of galaxy properties is seen on scales much larger than the radius of even the most massive collapsed halos (Hearin, Behroozi, & van den Bosch Reference Hearin, Behroozi and van den Bosch2016).

Spectroscopic redshifts over large areas are required to connect galaxy properties [e.g. stellar masses and star-formation rates (SFRs)] to the underlying dark matter halo masses and environments that are key to understanding galaxy formation physics. This has been demonstrated by the remarkable success of the Sloan Digital Sky Survey (most recently, Abolfathi et al. Reference Abolfathi2017) at z ∼ 0. The ATLAS surveys will extend the redshift baseline of all SDSS-like galaxy science to z ∼ 3 and in favourable cases to z ∼ 7 and beyond (Figures 6 and 7). Typical halo mass limits in Figure 6 are derived starting from empirical average star-formation histories (SFHs) as a function of halo mass and redshift from Behroozi et al. (Reference Behroozi2013). These SFHs are post-processed using FSPS (Conroy, Gunn, & White Reference Conroy, Gunn and White2009) to compute observer-frame F160W (H-band) apparent magnitudes, assuming Charlot and Fall (Reference Charlot and Fall2000) dust and a Chabrier (Reference Chabrier2003) initial mass function.

Figure 6. ATLAS will resolve key galaxy formation physics over unprecedented halo mass and redshift ranges. This figure shows limiting halo masses for each of the ATLAS Surveys as a function of redshift. Typical H-band (F160W) apparent magnitudes calculated using FSPS (Conroy et al. Reference Conroy, Gunn and White2009) from average star-formation histories in Behroozi et al. (Reference Behroozi2013), assuming Charlot and Fall (Reference Charlot and Fall2000) dust.

Figure 7. ATLAS Surveys will allow direct measurements of the stellar mass—halo mass relationship (top-left), spectroscopic group catalogues (top-right), environmental densities (bottom-left), and dark matter accretion rates (bottom-right). See text for details on how halo mass and redshift ranges were estimated.

We provide example halo mass and redshift ranges for four galaxy formation science cases in Figure 7. The stellar mass—halo mass relationship (e.g. Moster et al. Reference Moster, Naab and White2013; Behroozi et al. Reference Behroozi2013) measures the efficiency with which galaxies turn incoming gas into stars, and is a key probe of the strength of feedback from stars and supermassive black holes. The spectroscopic clustering measurements provided by ATLAS will directly measure this relationship for massive halos (bias b ≳ 1, Figure 7). Additional constraints will be provided by group catalogues, wherein halo masses are measured for individual galaxies based on spectroscopically identified satellite galaxy counts. These catalogues will also provide information about whether a given galaxy is a satellite or not; this is important for interpreting whether galaxy properties arise mainly from internal processes or instead from interaction with a larger neighbor. ATLAS will provide direct stellar mass–halo mass constraints out to z ≳ 7 and group catalogues to z ≳ 3, assuming a minimum of three satellites per group (Figure 7).

Besides halo mass, many open questions in galaxy formation concern how halo assembly affects galaxy properties and supermassive black hole activity. Halo assembly is not directly measurable, but many assembly properties (e.g. concentration, spin, mass accretion rate) correlate strongly with environmental density (e.g. Lee et al. Reference Lee2017). ATLAS will be able to measure environmental densities to z ≳ 3 (Figure 7), with galaxies in median-density environments having at least five neighbors within a 4 Mpc by 2 000 km/s cylinder even at the highest redshifts. Finally, ATLAS will be able to measure average dark matter accretion rates for galaxies via the detection of the splashback radius (a.k.a., turnaround radius) of their satellites as in More et al. (Reference More2016) out to z ∼ 5 (Figure 7). Here, we conservatively assume that a stack of 20 000 galaxies is sufficient to detect the splashback feature; More et al. (Reference More2016) achieved a > 6σ detection using only 8 000 galaxies.

At redshifts z > 1, most (but not all) galaxies were forming stars rapidly, and their spectra feature strong emission lines. With its 1–4 μm spectral range, ATLAS can observe the strongest optical rest-frame emission lines at high redshift (typically Hα at 0.5 < z < 5, and [O iii] at 1 < z < 7) without the spectral gaps that result from OH airglow contamination and water vapour absorption and low SNR from thermal backgrounds that are unavoidable from the ground. This continuous spectral coverage translates to continuous redshift sampling and very simple selection functions that cannot be achieved with ground-based observations. As the Universe ages, a growing number of galaxies cease star formation and evolve quiescently. This ‘quenching’ is related to both galaxy mass and environment, with quiescent galaxies dominating the denser regions of galaxy clusters and groups today. ATLAS will measure redshifts for these galaxies too, detecting stellar absorption lines and spectral breaks thanks to the minimal sky background for slit spectroscopy in space.

Today and in the next decade, ground-based spectroscopic surveys are measuring cosmic structure with outstanding statistics at 0 < z < 1, with more limited forays to higher redshifts, almost exclusively for strong emission-line galaxies. The emission-line flux limit of the ATLAS Wide Survey reaches 15 times fainter than the characteristic luminosity L* (Hα) at z = 2.23 (Sobral et al. Reference Sobral2013; Pozzetti et al. Reference Pozzetti2016) and L* ([O iii]) at z = 3.24 (Khostovan et al. Reference Khostovan2015), ensuring excellent sampling of large-scale structure. The continuum limits of the Wide Survey reach 2–3.5 times fainter than optical rest-frame continuum L* at z = 2–3.5 (e.g. Marchesini et al. Reference Marchesini2007), while the Medium and Deep Surveys reach 3.3–8.3 times fainter still, ensuring that ATLAS can measure absorption line and continuum break redshifts even for non-star-forming galaxies that might inhabit the densest cosmic structures. The continuum sensitivities of the Medium and Deep Surveys are tuned to match L* in the UV rest-frame for Lyman break galaxies at z ≃ 4 and 6, respectively (e.g. Finkelstein Reference Finkelstein2016). The Medium and Deep Surveys will uniquely map 3D clustering at 5 < z < 7 via the [O iii] emission line, which is observed to become stronger at higher redshifts, where galaxy metallicities are lower, and with absorption line redshifts for L* galaxies. The ‘wedding cake’ of ATLAS spectroscopic surveys will achieve outstanding spatial sampling of structure, from small to cosmological scales, and from z = 7 to 1. This is demonstrated in Appendix A, which gives the measurement errors on the galaxy 2-point correlation function ξ(r) as a function of the comoving separation r for the ATLAS Medium and Wide galaxy surveys.

ATLAS Probe galaxy surveys will obtain galaxy evolution data at high z that is equivalent to that from GAMA (Driver et al. Reference Driver2009) and SDSSFootnote ^h at low z, revolutionising our understanding of how galaxies have formed and evolved in the cosmic web of dark matter.

3.2. Emergence of galaxy properties

ATLAS will measure key diagnostics of galaxy properties through emission and absorption line spectroscopy. Its spectral range covers the most important optical rest-frame nebular emission lines at 1 < z < 5, including Hα, Hβ4861, [O iii], [N ii]6549,6583, and [S ii]6716,6731, with sufficient spectral resolution to resolve Hα+[N ii] and the [S ii] doublet. These lines trace star formation, gas excitation, metal abundance, electron density, and dust extinction, offering a wealth of information about the interstellar medium in distant galaxies. The sensitivity of the ATLAS Wide Survey ensures that both stronger and weaker lines can be measured in most typical, L* galaxies at 2 < z < 3.5. The ATLAS Medium and Deep Surveys will measure line diagnostics for galaxies at higher redshifts and/or lower luminosities. Moreover, in all surveys, spectra can be co-added for vast numbers of galaxies, grouped by other properties such as mass, SFR, excitation, and local environment density, in order to measure still weaker features, such as the auroral lines (e.g. [O iii] 4363, [O ii] 7320, 7330) needed for direct electron temperature and metallicity determinations. The ATLAS Medium and Deep Surveys will measure stellar continuum and absorption lines for faint galaxies (individually and in stacked bins), providing diagnostics of stellar population ages, SFHs, and chemical abundances. Observations of these features from terrestrial telescopes are limited by the Earth’s atmosphere, leading to highly incomplete and non-continuous spectral coverage and redshift range. JWST will obtain outstanding NIR spectra for faint galaxies, but only ATLAS, with its wide field and high multiplex, and its dedicated use for large surveys, can obtain spectra for hundreds of millions of galaxies over large cosmic volumes, needed to tie the evolving properties of galaxies to the context of their environments.

3.3. Black holes, AGN, and feedback

The distribution function of galaxy stellar masses is strikingly different from that of dark matter halos. Its shape implies a characteristic mass at which galaxies most effectively convert gas into stars. At lower and higher masses, ‘feedback’ is invoked to prevent gas from cooling onto galaxies, or to expel gas, thus suppressing star-formation efficiency. The most massive galaxies today are predominantly quiescent, with little active star formation, and with dominant bulges that universally contain supermassive black holes. These, when fuelled, can power active nuclei which may dump tremendous energy into their environments. AGN feedback, perhaps coupled with environmental effects, may be the dominant process regulating star formation and growth for massive galaxies. Like cosmic star formation, the luminous output of AGN and (by inference) the black hole growth rate, peaked at high redshift. ATLAS spectroscopy at 1–4 μm will be uniquely suited to identifying vast samples of AGN over a wide range of luminosity and redshift, using standard nebular excitation diagnostics (e.g. Baldwin, Phillips, & Terlevich Reference Baldwin, Phillips and Terlevich1981) and through detection of high-ionisation emission lines (e.g. [Ne v], He ii, C iv, N v). ATLAS will connect AGN activity to local and large-scale environment with exquisite statistical accuracy that is only possible today in the local universe, and [O iii] luminosities will provide a critical measure of the distribution of accretion luminosities and black hole growth rates. ATLAS spectroscopy will be highly complementary to surveys with Athena, the next major X-ray observatory, and uniquely suitable for redshift measurement of the most distant AGN candidates identified by Euclid, WFIRST, or Athena.

ATLAS will also enable spectroscopic redshifts to be measured for the majority of radio sources detected in wide-area surveys with the ASKAP such as EMU (Norris et al. Reference Norris2011). EMU will detect tens of millions of star-forming galaxies out to z ∼ 3 and several million quasars out to the earliest cosmic times. Although photometric redshifts of many of these galaxies will be possible through a combination of LSST and Euclid data, characterisation of these objects, including metallicities and black hole masses and their location within the cosmic web, will require spectroscopic data which ATLAS will be able to measure.

3.4. Reionisation and cosmic structure

Current observations indicate that the intergalactic medium (IGM) completed its transition from neutral to ionised by z = 6.5. This process is poorly understood: it is usually presumed that star-forming galaxies were responsible, but there is little evidence that sufficient ionising radiation escapes from early galaxies to accomplish this. Reionisation may have been highly inhomogeneous as well, with expanding bubbles driven by strongly clustered young galaxies that are highly biased tracers of dark matter structure. In coming decades, new radio facilities (LOFAR, HERA, MeerKAT, SKA) will map (at least statistically) the distribution of neutral hydrogen in the epoch of reionisation; ATLAS will provide essential complementary information about the spatial distribution of the (potentially) ionising galaxies themselves over the same sky areas and redshift ranges. This requires accurate spectroscopy deep enough to detect z = 7 galaxies over very wide sky areas—exactly what ATLAS will deliver. At 5 < z < 7, ATLAS will characterise the clustering of early galaxies with hard ionising spectra that produce [O iii] emission, a signature of the low-metallicity population that may be the main driver of IGM reionisation. While the redshift range falls mostly below the end of reionisation, an accurate measurement of 3D clustering will strongly constrain theoretical models that can then be extrapolated to higher redshifts. At z > 7.2, ATLAS can observe Lyman α (Lyα). There is already evidence that Lyα may be inhomogeneously suppressed by the neutral IGM at those redshifts (e.g. Tilvi et al. Reference Tilvi2014), and that its escape may correlate with galaxy overdensities that can more effectively ionise large IGM volumes (Castellano et al. Reference Castellano2016). The ATLAS Medium and Deep Surveys can target galaxies at z ≃ 7–8 selected from deep WFIRST Guest Observer science programmes over many square degrees, and measure (or set severe limits on) Lyα emission for statistical correlation with 21 cm surveys of the neutral IGM in the same volumes. JWST cannot survey areas wide enough to correlate large-scale structure with HI surveys, while the WFIRST slitless spectroscopy does not have suitable sensitivity for this. The ATLAS Wide Survey may also discover hundreds of examples of exceedingly rare, exceedingly luminous Lyα galaxies like ‘CR7’ (Sobral et al. Reference Sobral2015), which may also exhibit HeII 1640 emission, a signature expected from ionisation by primordial Population III stars (but see Bowler et al. Reference Bowler2017; Sobral et al. Reference Sobral2019).

3.5. The circumgalactic medium (CGM)

The CGM is a key, but poorly understood component of the cosmic web. Containing the bulk of baryons, it is the gas reservoir that cools to form stars and which is heated and enriched by feedback. The nature of the CGM (density, temperature, metallicity) as a function of galaxy properties is crucial to understanding galaxy formation, especially in the phase of rapid growth at z > 2. The high source density and wide spectral band of ATLAS offer a powerful probe of the CGM at these redshifts. By stacking spectra of many background galaxies, we can measure average CGM absorption around foreground galaxies with known redshifts as a function of impact parameter, mass, SFR, metallicity, local environment, and other properties. In the ATLAS Wide Survey, we can measure equivalent widths of 0.005 Å (3σ) for Ca H&K absorptions (using background galaxies at z > 2) and 0.010 Å for Mg ii (at z > 2.6). In addition, we can reach a 5σ sensitivity of 8 × 10⁻²² erg s⁻¹cm⁻² arcsec⁻² for detecting Hα emission from gaseous halos at z ∼ 2, probing the CGM to radii >100 kpc.

3.6. Protoclusters

ATLAS Medium Survey will have the survey sensitivity and volume to discover forming clusters and protoclusters, all the way from z = 2 to 6 and thus to unveil the crucial phases of cluster formation and galaxy formation in cluster environments. Based on cosmological N-body simulations, cluster at z > 2 are extended enough to utilise the wide-field high spectroscopic multiplex capability of ATLAS (see e.g. Orsi et al. Reference Orsi, Fanidakis, Lacey and Baugh2016). In addition, the ATLAS spectral resolution of R = 1 000 is sufficient for obtaining useful measurements of cluster galaxy velocities to probe cluster dynamics. ATLAS will provide a unique probe of clusters and their cosmic web environments.

The first collapsed structures hosted within a single dark matter halo with masses of few times of 10¹³ M_⊙ are expected to be in place in reasonable numbers at z ∼ 2–3, with space densities ranging from 10⁻⁵ to 10⁻⁶Mpc⁻³, depending on the exact mass cut, corresponding to several per square degree over 2 < z < 3 (so on average one in each ATLAS FoV). ATLAS will be able to obtain in-depth spectroscopic identifications of a dozen members down to low levels in the mass function in less than 1 h observations (line fluxes down to few 10⁻¹⁸ erg s⁻¹ cm⁻²). Many hundreds will be found, for example, over a systematic survey of 100 deg² by ATLAS Medium. Pointed observations of similar objects found by Euclid/Athena and/or other facilities will be very efficient. Observations of the first structures (z = 1.99, Gobat et al. Reference Gobat2013; z = 2.50, Wang et al. Reference Wang2016) show that these objects host spectacular action: from strong star-formation activity to quenching, from AGN feedback to the morphological transformation of galaxies and the formation of ellipticals. At those redshifts, the AGN activity will become much more prominent, following in parallel the rise of SFR and gas content in the Universe. In these conditions, interactions between the ICM and the cluster’s environment through the inflow of pristine cold gas, the outflows originated from stellar and AGN’s winds, and the deposition of warm plasma will shape the baryon distribution, the metal content and the energy budget of the structures (e.g. Valentino et al. Reference Valentino2015, Reference Valentino2016).

At higher redshifts, the first lower mass group-like halos (masses of 10¹³ M_⊙ or below) as well as protoclusters (larger scale overdensities not yet collapsed) will be within reach of spectroscopic identification by ATLAS. Groups of 10¹³ M_⊙ have an expected density of 1 per tens of ATLAS fields at 5 < z < 7, but several dozens could be found over the 100 square degree ATLAS Medium Survey. A basic characterisation with identification of most massive galaxies (few 10⁹ M_⊙) will require a few hours of integration. Detailed characterisation down the stellar mass function will require many dozens of hours of observations with ATLAS (fluxes of a few 10⁻¹⁹ erg s⁻¹ cm⁻²). Observing such objects to z ∼ 6 will be key to explore the early phases of environmental effects on galaxy formation and evolution as well as the beginning of structure formation (see e.g. Orsi et al. Reference Orsi, Fanidakis, Lacey and Baugh2016; Izquierdo et al. Reference Izquierdo-Villalba, Orsi, Bonoli, Lacey, Baugh and Griffin2018).

3.7. Galaxy kinematics

With the spectral resolution of R = 1 000, ATLAS Probe will allow us to obtain accurate kinematic measurements that will be essential in several galaxy evolution cases. For example:

• Absorption line velocity dispersions, which enable the measurement of scaling relations and, in combination with galaxy sizes estimated from imaging, dynamical masses of spheroidal/quiescent galaxies. Belli et al. (Reference Belli2017) and Van Dokkum, Kriek, & Franx (Reference Van Dokkum, Kriek and Franx2009) have used such measurements to demonstrate the evolution towards higher stellar velocity dispersion at fixed mass at the highest redshifts where such measurements have been performed to date (i.e. z ∼ 2).

• Emission-line virial kinematics, which can be interpreted in concert with WFIRST imaging and structural properties to infer dynamical masses (Price et al. Reference Price2016; Alcorn et al. Reference Alcorn2018). Such dynamical masses can be compared with inferred baryonic masses to probe the amount of dark matter within galaxy effective radii. It will be possible to infer baryonic masses from the sum of gas masses estimated from dust-corrected SFRs, assuming the Schmidt–Kennicutt relation, and stellar masses estimated from WFIRST galaxy photometry. At the same time, such comparisons place constraints on the allowed range of stellar initial mass functions (Price et al. Reference Price2016).

• Velocity shifts of interstellar absorption lines (e.g. Na i 5890,5896) relative to galaxy systemic velocity, to track outflows/inflows. Rest-optical nebular emission lines such as Hα and [O iii] 5007 can be used to establish the galaxy systemic velocity, as they are dominated by the surface-brightness weighted ensemble-averaged emission from H ii regions in galaxies. Over a wide range of redshifts, kinematic evidence for gas flows around galaxies has been assembled based on the small differences in redshift measured for interstellar absorption features relative to the systemic redshift. Most kinematic evidence of this type suggests outflows (e.g. Shapley et al. Reference Shapley2003; Steidel et al. Reference Steidel2010; Weiner et al. Reference Weiner2009; Chen et al. Reference Chen2010), but there are rare instances of kinematic evidence of infalling gas (Rubin et al. Reference Rubin2012; Martin et al. Reference Martin2012).

• Probe the dark matter content of disc galaxies as a function of redshift for 0.5 < z ≲ 7, by obtaining galaxy rotation curves from stacking thousands of spectra of galaxies with similar colours in redshift bins, scaled using stellar disc sizes measured by WFIRST. This approach is independent of rotation curve models, and shown to be unbiased (Tiley et al. Reference Tiley2019). We will use WFIRST imaging to select nearly-edge-on disc galaxies at 0.5 < z ≲ 7, and obtain several spectra along the major axis out to a few arsec for each galaxy. This will allow us to construct a rotation curve for each galaxy. The spectra from many galaxies, after being scaled using the observed disc size, can be stacked to obtain a well sampled rotation curve (Tiley et al. Reference Tiley2019). ATLAS Probe can obtain disc galaxy rotation curves out to z = 7 this way, since it can detect both [O iii] and [O ii] lines (both strong emission lines) out to z = 7. At the highest redshifts (z ≫ 3), the very compact sizes of galaxies will likely make the measurement of galaxy rotation curves challenging, but ATLAS Probe measurements will at least set interesting limits. ATLAS complements the ground-based projects by extending the redshift reach of rotation curve measurements to high redshifts not possible from the ground.

In addition, ATLAS Probe galaxy kinematic measurements can be used to study galaxy pairs and the evolution of merger fraction and merger rate, measure virial masses of clusters and groups based on galaxy velocity dispersion, and better deblend spectral lines, leading to better accuracy in the measurement of equivalent widths, metal abundances (e.g. ISM absorptions), stellar lines, etc.

3.8. Physics enabled by emission lines

A wide range of physical information is encoded in a galaxy’s spectrum. The main strong optical lines are particularly well-studied and start to shift into the observer frame of ATLAS at z > 0.5. Figure 8 shows when various key diagnostic features become accessible to ATLAS with emission lines shown as lines and absorption features as shaded regions.

Figure 8. The wavelength at which different spectroscopic features fall as a function of redshift, for the ATLAS wavelength range of 1–4 μm. The thinner lines denotes emission lines as indicated in the figure. Also indicated with broader bands are the D4000 and B2640 absorption features as well as the region including the Mg ii absorption lines.

The optical emission lines ([N ii]6548,6584, Hα, [S ii] 6717,6731, Hβ, [O iii]4959,5007 and [O ii]3727 in particular) have a long history of being used to estimate SFRs, gas-phase metallicities and ionisation conditions. At 1.7 < z < 5.0, all of these strong, optical rest-frame lines are accessible to ATLAS. Metallicity and ionisation conditions can be studied in large samples of galaxies over a wide range of cosmic epoch, drastically reducing the current systematic uncertainties that are caused by patching together different indicators at different redshifts. The hardness of the ionising spectrum can be traced by weaker optical lines, such as [O i]6300, while in the UV the [C iii]1909, [O iii] 1661,1666, and C iv 1548,1550 doublets and the He ii 1640 line have been shown to be present frequently in low-mass star-forming galaxies with low metallicity/high redshift (e.g. Senchyna et al. Reference Senchyna2017; Stark et al. Reference Stark2015) but they are also prominent tracers of AGN activity (Feltre et al. Reference Feltre, Charlot and Gutkin2016).

Finally, the Lyα line offers a potential cornucopia of information on the properties of galaxies. It is related to the SFR in galaxies, modulated by resonant scattering and dust absorption through the neutral interstellar and circumgalactic medium. The shape of the Lyα profile can be used to constrain the distribution and kinematics of outflowing and inflowing gas around galaxies (Verhamme, Schaerer, & Maselli Reference Verhamme, Schaerer and Maselli2006; Steidel et al. Reference Steidel2010). Recent research on low-redshift galaxies with Lyman continuum escape has shown that the separation of the blue and red peaks in the emission line also correlate with the escape fraction of ionising photons, f _esc (Verhamme et al. Reference Verhamme2017; Izotov et al. Reference Izotov2018). At f _esc > 0.05, the typical separation is <400 km/s which is marginally resolvable by ATLAS with R = 1 000. When combined with a study of the systematic changes in the overall line profiles, this offers ATLAS the possibility to trace changes in the escape fraction with redshift at the very highest redshifts.

4.1. Weighing dark-matter-dominated filaments in the cosmic web

A key prediction of the cold dark matter model is the existence of the ‘cosmic web’, a network of low-density filaments connecting dark matter halos. Spectroscopic redshifts are required to locate the galaxy groups and clusters that are connected by filaments in redshift space. The uncertainty associated with photometric redshifts scatters the true physical pairs and smears the filamentary structure (see Figure 5). Epps & Hudson (Reference Epps and Hudson2017) recently demonstrated that the stacked dark-matter-dominated filament connecting pairs of massive galaxies can be detected when spectroscopic redshifts along with WL data are available for most of the source galaxies lensed by the massive galaxy pairs. ATLAS Probe spectroscopy and WFIRST imaging together provide the ideal data set for such detections. The detection of these filaments by ATLAS Probe over a significant cosmic volume will test the cold dark matter model for structure formation, and provide key insight into large-scale structure in the Universe.

Epps & Hudson (Reference Epps and Hudson2017) used the overlap in the CFHTLenS imaging and BOSS spectroscopic data over 105 deg² (Miyatake et al. Reference Miyatake2015) to select a sample of ∼20 400 luminous red galaxies (LRGs), from which they constructed a catalogue of LRG pairs by selecting pairs that were separated by Δz _spec < 0.002 (∼5 h ⁻¹ Mpc in comoving separation) and by 6 to 10 h ⁻¹ Mpc in the transverse direction. This yielded a sample of ∼23 000 pairs of LRGs, with a mean separation of ∼8.23 h ⁻¹ Mpc at a mean redshift of 0.42, and a mean stellar mass of 10^11.3 M_⊙ (with the expected halo mass of 10^13.04 M_⊙, corresponding to galaxy groups). Epps & Hudson (Reference Epps and Hudson2017) constructed shear and convergence maps by stacking all the LRG pairs in their sample, and obtained a 5σ detection of the stacked filament connecting the LRG pairs.

The Canada–France Imaging SurveyFootnote ⁱ (CFIS; Ibata et al. Reference Ibata2017) is a CFHT Large Program that has been allocated 271 nights from 2017 to 2021. It will provide imaging to complement the spectroscopy from BOSS and DESI, to enable dark matter filament measurements over a wide sky area at z ≲ 0.8. ATLAS Probe complements this by enabling dark matter filament measurements over 2 000 deg² at 0.5 < z < 3.0, with much higher galaxy number densities (see Figure 9).

Figure 10 shows the estimated uncertainties of filament mass measurements for slices in redshift z of width indicated by the horizontal error bars (so averaging over larger redshift range at higher z). The large blue cross is the CFHTLenS result from Epps & Hudson (Reference Epps and Hudson2017). Orange symbols at low z are forecasts for CFIS (imaging) with BOSS (spectroscopy), and green represents conservative forecasts for ATLAS Probe (spectroscopy) with WFIRST (imaging). Note that the signal drops at high z due to structure growth and the noise increases at high z due to lack of source galaxies behind the higher lens redshifts.

Figure 10. Conservative estimate of the measurement errors of filament mass from ATLAS Wide, for a 6–10 Mpc h ⁻¹ cosmic web filament between 10¹³ M_⊙ halos. See text for details.

We have assumed that filaments with mass ∼1.5 × 10²³ M_⊙ connect halos that are ∼10¹³ M_⊙ at z ∼ 0.4, and have a number density of ∼3 × 10⁻⁴(h Mpc⁻¹)³ (BOSS-LRG-like in density but the halos do not have to be passive, just massive). To derive galaxy bias b(z), we assumed constant clustering, i.e. b(z)G(z) is constant, where G(z) is the linear growth factor. Since the galaxies have b = 2 at z ∼ 0.43, we find b(z) = b(0.43)G(0.43)/G(z). The filament mass is assumed to scale with redshift like the three point correlation function.

The mass estimate uncertainties in Figure 10 are conservative because it should also be possible to detect filaments between lower mass halos. In future work, we will carry out calculations based on realistic simulations to increase the fidelity of our estimates, and study the measurement of filament mass at higher redshifts.

The massive galaxies that mark the galaxy groups can be identified through colour selection in the WFIRST HLS imaging data. We expect that ∼70% of the source galaxy redshifts will be measured by ATLAS Wide. The lens galaxies are expected to be massive galaxies, some of which may be passive ellipticals with very low emission-line flux. However, the LRGs are very bright. Zhai et al. (Reference Zhai2017) showed that the LRGs from eBOSS (0.6 < z < 0.9) all have i ≲ 22, with the majority at i < 21, with very mild redshift evolution. Since ATLAS will reach a continuum depth of AB = 23 at 3σ in 5 000 s (the exposure time per target for ATLAS Wide), we expect to obtain redshifts for nearly all of the LRGs used for the dark matter filament mass measurements.

The dark matter filament mass measurements by ATLAS Probe are synergistic with galaxy evolution studies using the three-tiered ATLAS Probe galaxy surveys. Using the GAMA spectroscopic survey in the redshift range 0.03 ≤ z ≤ 0.25, Kraljic et al. (Reference Kraljic2018) have shown that the cosmic web plays an important role in shaping galaxy properties. Key galaxy properties (such as stellar mass, dust corrected colour, and specific SFR) are found to be correlated with galaxy distances to the 3D cosmic web features such as nodes, filaments, and walls. ATLAS Probe is unique in providing samples of different galaxy types, i.e. with different bias factors. These enable the ATLAS Probe galaxies to trace the cosmic web in unprecedented details and precision. ATLAS Probe measurements of dark matter filaments can be used to constrain models of the cosmic web, and enable the calibration of the reconstruction of the underlying dark matter distribution. This will enable transformative progress in our understanding of galaxy evolution and the nature of dark matter.

4.2. Differentiating dark energy and modification of general relativity

The observed cosmic acceleration is one of the biggest unsolved mysteries in cosmology today (Riess et al. Reference Riess1998; Perlmutter et al. Reference Perlmutter1999). We do not even know if it is due to an unknown energy component in the Universe (i.e. dark energy), or a modification of General Relativity (i.e. modified gravity). The measurement of both the cosmic expansion history and the growth history of cosmic large-scale structure is required to solve this mystery (see e.g. Guzzo et al. Reference Guzzo2008; Wang Reference Wang2008). WL and galaxy clustering both probe cosmic expansion history, and probe different ways in which gravity can be modified. The ATLAS Wide Survey will provide the ultimate data set to remove systematic effects in our quest to discover the cause for cosmic acceleration.

(i) Multi-tracer and ultra-high-density BAO/RSD. Galaxy clustering data from 3D distributions of galaxies is the most robust probe of cosmic acceleration. The baryon acoustic oscillation (BAO) measurements provide a direct measurement of cosmic expansion history H(z) and angular diameter distance D_A(z) (Blake & Glazebrook Reference Blake and Glazebrook2003; Seo & Eisenstein Reference Seo and Eisenstein2003), and the redshift-space distortions (RSD) enable measurement of the growth rate of cosmic large-scale structure f_g(z) (Guzzo et al. Reference Guzzo2008; Wang Reference Wang2008). The ATLAS Wide Survey will obtain spectroscopic redshifts of 183M galaxies over 2 000 deg², ∼12 times higher in galaxy number density compared to WFIRST. These provide multiple galaxy tracers of BAO/RSD (red galaxies, different emission-line-selected galaxies, and WL shear-selected galaxies) over 0.5 < z < 4, with each at high number densities. These enable robust modelling of BAO/RSD (e.g. the removal of the nonlinear effects via the reconstruction of the linear density field), and significantly tightens constraints on dark energy and modified gravity by evading the cosmic variance when used as multi-tracers (McDonald & Seljak Reference McDonald and Seljak2009). ATLAS measures H(z), D_A(z), and f_g(z) over the wide redshift range of 0.5 < z < 4 (see Figure 11), with high precision over 0.5 < z < 3.5, improving the constraints from WFIRST by a factor of three or more at 2 < z < 3 (beyond the reach of Euclid), and extending constraints to 3 < z < 4, beyond the reach of both WFIRST and Euclid. The 14 000 deg² ELG survey from the ground-based project DESI (Aghamousa et al. Reference Aghamousa2016) is complementary to the space-based surveys, as shown in Figure 11.

Figure 11. ATLAS vs. WFIRST and Euclid comparison of the measurement errors from BAO/RSD on the cosmic expansion history H(z), angular diameter distance D_A(z), and the growth rate of cosmic large-scale structure f_g(z) = β(z)b(z) [where β(z) is the linear redshift-space distortion parameter, and b(z) is the bias factor]. The BOSS survey and the future 14 000 deg² ELG survey from the ground-based project DESI (Aghamousa et al. Reference Aghamousa2016) are also shown for comparison. H(z) and D_A(z) are scaled with the BAO scale s (the comoving sound horizon at the drag epoch), which is measured precisely by Planck.

We have followed the methodology of Wang, Chuang, & Hirata (Reference Wang, Chuang and Hirata2013) in deriving the constraints shown in Figure 11, with 100% nonlinear effects, a cut-off of k _max = 0.25h/Mpc, and peculiar velocity dispersion of 290 km/s. We assumed a bias function b(z) = 0.84G(0)/G(z), following the fitting of the Mostek et al. (Reference Mostek2013) results from Aghamousa et al. (Reference Aghamousa2016), where G(z) is the growth factor. This corresponds to constant clustering. It is appropriate for the majority of the WFIRST WL source galaxies, most of which are emission-line objects at some level (though the lines are fainter than the ones found by slitless spectroscopy on Euclid and WFIRST). For the slitless surveys, we assumed an efficiency of 75% for WFIRST, and 40% for Euclid, since the signal-to-noise thresholds for WFIRST and Euclid are 7σ and 3.5σ, respectively. The galaxy number density for ATLAS is shown in the left panel of Figure 9 (normalised for a total of 183M galaxies). For WFIRST and Euclid, we use the forecasts shown in Merson et al. (Reference Merson2018). The galaxy number density for DESI is from Aghamousa et al. (Reference Aghamousa2016).

(ii) 3D WL with spectroscopic redshifts. The ATLAS Wide Survey replaces photometric redshifts with spectroscopic redshifts for ∼70% of the lensed galaxies in the WFIRST HLS WL sample. This eliminates the photo-z calibration ladder as a major source of systematic uncertainty in WL. Furthermore, one can identify every pair of galaxies that are near each other in 3D space for 70% of the sample, dramatically suppressing the systematic error from intrinsic galaxy alignments. This improves the measurement of the growth rate index (used to parametrise the gravity model) by ∼50%. In future work, we will investigate the cosmological constraints from the ATLAS spectroscopic subsample of the WFIRST WL sample.

(iii) Joint analysis of weak lensing and galaxy clustering. Less than 5% of the galaxies in the WFIRST WL sample will have WFIRST spectroscopy; ∼70% will have ATLAS spectroscopy from the Wide Survey at z > 0.5. This super data set of both WL shear and 3D galaxy clustering data for the same 183M galaxies over 0.5 < z < 4 enables a straightforward joint analysis of the data. This facilitates the precise modelling of the bias between galaxies and matter b(z), and provides the ultimate measurement of H(z) and f _g(z). We expect this to result in definitive measurements of dark energy and modified gravity for 0.5 < z < 4, with the removal of photo-z errors as the primary systematic for WL, and detailed modelling of bias systematic for BAO/RSD.

(iv) Type Ia Supernovae (SNe Ia). ATLAS Probe will obtain the host galaxy redshifts of nearly all 30 000 SNe Ia that will have light curves measured by the LSST DDF and WFIRST SN surveys, over the redshift range of 0.2 < z < 2.0. This will provide a powerful measurement of H(z) that is independent of cosmic large-scale structure. The WFIRST and LSST SN surveys will cover ∼40 deg² each. There are currently no plans in place to obtain host galaxy redshifts for LSST SNe with z > 0.5. For WFIRST, there may be an Integrated Field Channel on board to obtain spectra of SNe and hosts; however, for z > 1.0, the number of spectra acquired will be ∼10% of the full discovery sample (Hounsell et al. Reference Hounsell2018). Other plans, like a joint effort with Prime Focus Spectrograph, will obtain host galaxy redshifts to z ∼ 1.0, but a small and largely biased sample for z > 1.0. Multiple studies (e.g. Sullivan et al. Reference Sullivan2010) have shown that there is a relation between SN luminosity and host galaxy properties, so a large effort should be made to create an unbiased sample of host galaxies to as high a redshift as possible.

With host galaxy redshifts up to z = 2.0 from ATLAS Probe, the statistical distance precision of WFIRST and LSST SNe should be sub-1% out to z = 2.0 (Hounsell et al. Reference Hounsell2018). This is comparable to the constraints shown in Figure 11 from the ATLAS measurement errors from BAO/RSD. Without the host galaxy redshifts, SN analyses will be forced to rely on photometric redshifts, which for the full SN sample degrades the cosmological precision by introducing a series of systematic uncertainties and reducing the statistical precision as well.

(v) Clusters. The abundance of mass-calibrated galaxy clusters provides a complementary measurement of cosmic expansion history and growth rate of large-scale structure. In the ATLAS Wide Survey, we will carry out spectroscopy for the 40 000 clusters with M > 10¹⁴ M_⊙ expected to be found by WFIRST HLS imaging (Spergel et al. Reference Spergel2015). The ATLAS/WFIRST cluster sample will have mass accurately measured by the deep 3D WL data, and cluster membership precisely determined by spectroscopic redshifts. This will provide a powerful crosscheck to constraints from BAO/RSD and 3D WL, and a unique data set of spectroscopic clusters for studying cluster astrophysics.

(vi) Higher-order statistics. The very high number density galaxy samples from ATLAS Probe provide the ideal data set for studying higher-order statistics of galaxy clustering, which will boost the power to probe dark energy and gravity by a significant factor. While the use of galaxy clustering 2-point statistics is now standard in cosmology, the use of the 3-point statistics is still limited due to a number of technical challenges. Since the galaxy 3-point statistics provides additional information to that from the galaxy 2-point statistics, the combination of these is needed to optimally extract the cosmological information from galaxy clustering data (see e.g. Gagrani & Samushia Reference Gagrani and Samushia2017).

(vii) Calibration of photometric redshifts. In addition to providing redshifts for ∼70% of the WFIRST HLS lensing sample at z > 0.5, ATLAS should greatly improve photometric redshifts for those WFIRST galaxies which it does not obtain a robust redshift determination for. This will be true for two reasons. First, the redshifts ATLAS obtains will all be helpful for improving our knowledge of the relationship between galaxy colours/magnitudes and redshift (or, equivalently, for constraining models of galaxy SED evolution); the Medium Survey should be particularly well-suited for this, especially if the survey area is distributed amongst multiple fields to mitigate the effects of sample/cosmic variance (in contrast, variance in the 1 deg² deep field will be sufficient to limit its use for this purpose). If the survey is appropriately designed, ATLAS should greatly exceed the requirements for a Stage IV photometric redshift training sample established in Newman et al. (Reference Newman2015).

Second, the large sample of galaxy redshifts over a large volume provided by the ATLAS Wide Survey would provide a high-precision calibration of photometric redshifts, even in the presence of systematic incompleteness, via cross-correlation between ATLAS galaxies with secure redshifts and photometric-only samples (Newman Reference Newman2008). Gains would be particularly great at z ≳ 1.5, where the DESI and 4MOST samples will be limited to QSOs (and hence relatively sparse). The resulting calibration would ensure that accurate dark energy inference can be done even if spectroscopic redshift samples are systematically incomplete. We note that not only WFIRST but other Stage IV (or future) photometric dark energy experiments (e.g. LSST and the WL imaging survey on Euclid) would benefit from the improved photometric redshift training and calibration provided by ATLAS.

4.3. Particle cosmology: inflationary physics, neutrino mass, and particle dark matter

Probe inflationary physics. Current observational data are consistent with inflation having occurred when the Universe was 10⁻³⁴ s old. However, we have very little insight on the physics of inflation. While the detection of primordial gravity waves by future CMB observations would provide the definitive proof that inflation occurred, the combination of CMB data from Planck and the ultra high-precision matter power spectrum P(k, z) (with 0.5 < z < 4) measured from the ATLAS Wide Survey will characterise the primordial matter power spectrum P(k)_in with definitive precision and accuracy. The measured P(k)_in allows us to probe inflationary physics, e.g. whether inflation occurred in one stage, or in multiple stages, and whether inflation occurred smoothly, or in spikes (Wang, Spergel, & Strauss Reference Wang, Spergel and Strauss1999). These help determine the particle physics required for inflation, e.g. properties of the inflaton, including its coupling to other particles, which will help in its identification. We expect ATLAS to provide an order of magnitude improvement in measuring P(k)_in over other planned galaxy redshift surveys due to its extremely high galaxy number density over a broad redshift coverage, opening up in particular the high-redshift window where the nonlinearities are less important and the modelling more robust.

Measure the total mass of neutrinos. Solar and atmospheric neutrino data indicate that neutrinos have a small but nonzero mass (see e.g. Lesgourgues & Pastor Reference Lesgourgues and Pastor2006). This has profound implications for particle physics and cosmology. Measuring the tiny mass of neutrinos is a major enterprise in experimental particle physics, but cosmological observations will likely provide the most stringent constraints on the total mass of neutrino species, and possibly be sensitive enough to the mass hierarchy (Carbone et al. Reference Carbone, Verde, Wang and Cimatti2011). By covering λ > 2m and obtaining spectra for faint galaxies, the ATLAS galaxy surveys probe structures at early times and is not limited to the most massive/biased objects. ATLAS thus has the potential to be the least affected by systematic errors in modelling nonlinear structure evolution, allowing us to measure P(k, z) (0.5 < z < 4) down to very small scales. This will allow us to measure the total mass of neutrinos to high precision, as the effect of massive neutrinos increases at small scales. The ATLAS Medium Survey will provide a unique new angle to measure neutrino mass by providing precise measurement of galaxy clustering on the scales of 1–50 h ⁻¹ Mpc for 0.5 < z < 7. This will help tighten the ATLAS neutrino mass measurement by probing the effect of massive neutrinos to even earlier times. If the neutrino mass m_v ≥ 0.06eV, ATLAS Probe will be able to detect it with a significance exceeding 3σ (Ballardini et al., in preparation). Because ATLAS Probe will have the tightest possible control of systematic uncertainties in both galaxy clustering and WL, we expect it to lead to the most accurate measurement of m_v within the next few decades. In particular, ATLAS Probe data will enable the modelling of galaxy bias on small scales via the joint analysis of 2-point galaxy statistics with 3-point galaxy statistics and WL data, removing the main systematic uncertainty in the neutrino mass measurement.

Constrain axions as dark matter. Axion is the leading alternative to WIMPs (weakly interacting massive particles) as a particle candidate for dark matter. Axions are motivated by known particle physics, as it solves the strong CP problem. If dark matter were made of axions, or a mixture of axions and WIMPS, they would give rise to specific features in the cosmic large-scale structure (Hlozek et al. Reference Hlozek, Grin, Marsh and Ferreira2015; Cedeno, Gonzalez-Morales, & Urena-Lopez Reference Cedeno, Gonzalez-Morales and Urena-Lopez2017), which can be measured to exquisite statistical precisions and accuracy using ATLAS Wide data. We will explore these constraints in future work.

5.1. The 3D structure of the hidden Milky Way

Currently, we know more about the structure and the spatially resolved SFHs of galaxies in the Local Group (and even beyond) than we do about our own Galaxy. Results from Gaia will revolutionise our understanding of the 3D Milky Way, but Gaia is limited, especially in the inner Galaxy, due to it being an optical survey. Figure 12 illustrates the extinction challenge. For the inner Galactic plane (|l| < 65° and |b| < 1°), 98.5% of the Galactic plane has a minimum extinction of 2.5 mag in the J band. This drops to only 10.4% of the Galactic plane in the K band. The infrared range of ATLAS opens up large regions of the Galaxy for spectral investigation.

Figure 12. A Spitzer/IRAC colour image of the inner Galactic plane/bulge from 3.6 μm (blue) to 8.0 μm (red) showing the dramatic reduction of extinction from the J band (1.2 μm; blue contours) to the K band (2.2 μm; red contours) to the Spitzer 3.6 and 4.5 μm bands (background image) using the extinction map of Gonzalez et al. (Reference Gonzalez2012). Extinction in the infrared bands compared to the visual are A _band/A _V=0.27, 0.17, 0.10, 0.06, and 0.03 for the J, H, K, L (3.4 μm), and M (4.8 μm) bands, respectively (Draine Reference Draine and Draine2011). Stars with 30 magnitudes of visual extinction in the midplane will only suffer one magnitude of extinction for the longest wavelengths covered by ATLAS. The respective footprints of WFIRST and ATLAS are shown as are all of the Galactic plane fields that have been imaged by HST.

ATLAS will acquire the spectra in a strip of ±1° centred on the Galactic Equator, corresponding to the approximate area of the family of Spitzer/IRAC 3.6 and 4.5 μm surveys: GLIMPSE-I (|l| = 10°–65°), GLIMPSE-II (|l| < 10°); GLIMPSE-360 (|l| > 65°) (Churchwell et al. Reference Churchwell2009; Hora et al. Reference Hora2007; Carey et al. Reference Carey2008; Majewski et al. Reference Majewski2007; Whitney et al. Reference Whitney2008).Footnote ^j These cover both the inner and outer Galaxy, and thus very different regions in terms of spiral structure, star formation, evolved stars and interstellar dust. As an example, Figure 13 illustrates some of the 720 deg² footprint of the Spitzer/GLIMPSE surveys. A future Galactic Plane imaging survey with WFIRST will improve our knowledge of the stellar populations, especially in the inner Galaxy where GLIMPSE is several magnitudes shallower than in the outer Galaxy, and would provide the ideal source catalogue for ATLAS.

Figure 13. Left: GLIMPSE image. Right: Magnitude distribution at 3.6 μm of the several hundred million GLIMPSE catalogue sources. Reprinted from Meade et al. (Reference Meade2014), ‘GLIMPSE360 Data Description GLIMPSE360: Completing the Spitzer Galactic Plane Survey’, https://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/doc/glimpse360_dataprod_v1.5.pdf.

The ATLAS Galactic Plane Survey will represent a formidable follow-up to the GLIMPSE surveys, and complement shallower existing Galactic plane surveys: the NIR UKIDSS (JHK) and optical IPHAS (r, i, and Hα), as well as the even larger PanSTARRS optical, WISE 3–23 μm mid-infrared, and SCUBA-2 sub-mm survey, the FCRAO surveys of CO in the outer Galaxy, the Herschel Hi-GAL survey and Planck, plus all future major optical/infrared surveys like Gaia, ZTF, LSST, and WFIRST. ATLAS spectra will be crucial to interpreting WFIRST photometry in particular, in which photospheric temperatures and extinction effects will be almost completely degenerate due to the filter choices. ATLAS spectra will also complement the planned SDSS-V spectroscopic survey programmes covering brighter Milky Way objects (H < 11.5 mag).

The ATLAS Galactic Plane Survey will probe galactic structure, study star-forming regions, and measure interstellar extinction.

(a) Galactic structure. By taking advantage of the low extinction in the 1–4 μm region and the relatively high angular resolution and sensitivity of ATLAS compared to Spitzer, a spectroscopic survey of the inner Galaxy can produce unique information on the bar(s) of the Galaxy, the nuclear region, the stellar disc, and spiral arms. Existing and planned photometric surveys can provide useful results using star counts techniques and statistical analysis of colour–colour and colour magnitude diagrams. However, only spectroscopic information allows firm derivation of the extinction and effective temperature of the sources, and therefore their position in the HR diagram. Gaia will determine 1 (10)% distances for 10⁶ (10⁸) stars out to 2.5 (25) kpc, as faint as V ∼ 19 mag or K ∼ 11–19 mag, depending on intrinsic colour. Not only will ATLAS be able to provide spectral information on all of the (nearby) stars observed by Gaia in the optical, it will be able to see many additional stars that are not detectable in the optical. The calibration of ‘spectroscopic parallax’ using Gaia for nearby stars will allow us to use improved spectroscopic parallax for the much more distant stars that will be detected by ATLAS. Hogg, Eilers, & Rix (Reference Hogg, Eilers and Rix2018) show an example of this, using Gaia observations of nearby red giants to improve distance estimates of APOGEE selected red giants across the Galactic disc. Reconstructing the 3D structure of the Galaxy will allow progress on a number of fundamental questions, e.g. regarding the scale-length of the Galactic disc, whether the stellar warp and the gas warp coincide, and existence of stellar streams across the Galactic plane. ATLAS will address these issues by reaching disc regions well beyond the extinction-limited range sampled by SDSS, PanSTARRS, Gaia, and LSST.

(b) Star formation.Spitzer/GLIMPSE surveys produced an unprecedented picture of star formation in our Galaxy by unveiling hundreds of new star-forming regions. ATLAS will quantify the SFR of the Galaxy, its variation with Galactocentric radius, and its association with various dynamical features in the Galaxy, testing theories of star formation both on a global scale and at the molecular cloud level. Moreover, as the metallicity declines in the Outer Galaxy (reaching values intermediate between that of the LMC and in the inner Galactic disc), ATLAS will allow measurement of variations in the SFR representative of typical conditions at z = 2–3, the peak of the cosmic star-formation efficiency. ATLAS will characterise Young Stellar Objects (YSO) and their surrounding environments, analysing the energy budget and spatial correlation between ionisation (through recombination lines), photodissociation (through PAH emission), and high-velocity outflows (from shocked H ₂). The specific early type (<B2) stars that excite Galactic Plane H ii regions out to 10 kpc will be identified via their distinguishing 1–4 μm spectroscopic features, revealing important details regarding the spatially resolved SFH of the Galaxy over the past <50–100 Myr. Among lower mass Young Stellar Objects (left panels of Figure 14), accretion and outflow rates can be measured across individual star-forming regions, probing accretion histories and quantifying burst vs. steady accretion scenarios.

Figure 14. Example published spectra of: young stars (left panels out to 2.5 μm only; Connelley & Greene Reference Connelley and Greene2010), low-mass stars (middle panels, covering 1–4 μm; Cushing, Rayner, & Vacca Reference Cushing, Rayner and Vacca2005), and evolved objects (right panels out to 2.5 μm only; Cooper et al. Reference Cooper2013). Figures reprinted with permission from the authors.

(c) Interstellar extinction. By measuring the extinction law over tens of millions of lines of sight, ATLAS will provide a unique dataset to study in detail the variation of the extinction curve out to the edge of the Galactic disc. Variations with Galactic longitude have been identified and attributed to small variations in ISM density, mean grain size, or disc metallicity gradient. Variations of chemical composition of dust grains reflect the abundance/depletion of metals in the ISM, and hence the cooling mechanisms that control the efficiency of star formation. Characterising the variations in infrared extinction is important for enabling extinction corrections of Galactic plane sources from photometric surveys and understanding the nature of dust grains. A particular puzzle from Spitzer photometric studies was the fact that extinction in the 3–8 μm range was unchanging (flat) in contradiction to expectation from models by Draine & Li (Reference Draine and Li2001). It has been hypothesised that this indicates a population of micron-sized grains which accounts for 15% of the total dust mass (Wang et al. Reference Wang, Li and Jiang2015). The spectra from ATLAS of stars with known intrinsic spectra, e.g. red clump giants, will allow one to characterise how the infrared extinction varies with line of sight and search for deviations from the flat extinction curve indicated by current photometric data.

5.2. Low-mass stars, brown dwarfs, and exoplanets

ATLAS Wide Survey (the spectroscopic follow-up to the WFIRST HLS imaging), with its high Galactic latitude, complements the ATLAS Galactic Plane Survey. It will probe the low-mass stellar content of the Milky Way.

Brown dwarfs cool through the L-T-Y spectral sequence as they age (Burrows et al. Reference Burrows1997; Kirkpatrick Reference Kirkpatrick2005; Cushing et al. Reference Cushing2011). These objects provide important information on the shape and low-mass cut-off of the field mass function, and they also stand as proxies for exoplanets, with their atmospheres having similar effective temperatures but without the complicating effect of an irradiating host star. Based on the theoretical predictions of Baraffe et al. (Reference Baraffe, Chabrier, Barman, Allard and Hauschildt2003), at the 3σ limit of J ≈ 23.5 mag AB ATLAS will be able to detect 5-Gyr-old field brown dwarfs as low in mass as 10M _Jup (T _eff ∼ 300K, an early-Y dwarf) at 10 pc and 35M _Jup (T _eff ∼ 700 K, late-T) at 100 pc. For a robust, well-characterised sample in which ATLAS spectra of SNR = 100 are required, corresponding to J ≈ 18.5 mag AB, these limits are 35M _Jup (T _eff ∼ 700 K, late-T) at 10 pc and 75M _Jup (T _eff ∼ 1800 K, mid-L) at 100 pc. In terms of spectral type, ATLAS can detect and characterise mid-L dwarfs to 250 pc and mid-T dwarfs to 100 pc at J = 21.0 mag AB (SNR ∼ 10).

For low-mass stars, ATLAS can detect and characterise sources at the distance of the Galactic Bulge. These are important for a variety of Milky Way research topics (see section above) as well as providing valuable data for the interpretation of microlensing light curves when these objects are lensed by an intervening system. Spectra with SNR > 100 are obtainable at a distance of 8.0 kpc for types earlier than ∼F8 V and with SNR > 10 for types earlier than ∼K7 V. Detections at the SNR = 3 limit are possible down to a type of ∼M3 V.

5.2.1. The mass function at low metallicities

Low-mass stars and brown dwarfs are brightest in the NIR (see Figure 14, middle), but the nearest examples to the Sun almost all have metallicities near solar. Is the mass function at lowest masses dependent upon the composition of the gas from which the objects formed or is it invariant to metallicity, as early theoretical investigations suggest (Bate Reference Bate2014)? Gizis & Reid (Reference Gizis and Reid1999) determined that the local M star population is comprised of 99.8% dwarfs with [M/H] ∼0.0 and only 0.2% subdwarfs with [M/H] <−0.5. Therefore, to find large collections of these metal-poor objects with which to study population statistics, deeper surveys in the NIR are needed.

The WFIRST HLIS provides a seed list with which to answer this question. Subdwarf candidates will be selected photometrically using the WFIRST data themselves, with ATLAS providing spectroscopic verification or refutation. The ATLAS data will be further used to establish both the spectral type and metallicity class of those that are verified. These, combined with WFIRST HLIS photometric measurements, will establish space densities of these objects with which we can determine mass functions of the stellar subdwarf populations as a function of metallicity.

Only a few dozen late-M and early-L subdwarfs have been uncovered to date in the immediate Solar Neighbourhood, thanks primarily to 2MASS, which probes to J ∼ 17 mag AB. For sources with J < 21.0 mag AB (spectra with SNR > 10), ATLAS follow-up of candidates selected from the 2 000-deg² HLIS can reach a volume ∼15 times larger than the all-sky survey of 2MASS. The spectroscopic hallmark of these objects falls within the ATLAS observing window; specifically, strong collision-induced absorption by H₂ is the dominant shaper of the continuum, as shown in Figure 15(left).

Figure 15. Left: The 1–2.5 μm spectra of two normal late-M dwarfs (top), two metal-poor late-M/early-L subdwarfs (middle), and two even more metal-poor late-M ultrasubdwarfs (bottom) from Kirkpatrick et al. (Reference Kirkpatrick2010, Reference Kirkpatrick2016). For metal-poor objects, the broad trough caused by collision-induced absorption by H₂ flattens out the entire near-infrared spectrum, with the strongest effect at K-band. This flattening along with the concomitant bluer J−K colours are the hallmarks of low metallicity. In this plot, regions of telluric absorption are shown by the dark grey bands if the atmospheric transmission is less than 20% and by light grey bands if the atmospheric transmission is between 20% and 80% (Rayner et al. Reference Rayner, Cushing and Vacca2009).. Spectra are normalised at 1.28 μm and offsets added to separate the spectra vertically. Right: A Hammer projection in Galactic coordinates showing the sky distribution of known L, T, and Y dwarfs (black points), taken from Dwarf Archives (see http:/dwarfarchives.org) and the List of All UltraCool Dwarfs (see https://jgagneastro.wordpress.com/list-of-ultracool-dwarfs/). Note the lack of brown dwarfs in a wide band centred along the Galactic plane.

Rarer still are late-L and T subdwarfs that populate the brown dwarf side of the low-metallicity mass function. Less than a dozen of these objects are known near the Sun (Kirkpatrick et al. Reference Kirkpatrick2014, Reference Kirkpatrick2016; Pinfield et al. Reference Pinfield2014) so identifying more of these objects in the ATLAS Wide Survey would enable a mass function analysis below the substellar break. As more of these objects are uncovered, the width of the temperature gap between the stellar and substellar sequences at a particular metallicity can be measured, offering an independent measurement of the age of these populations.

5.3. Evolved stars and other objects

Evolved stars. Because of their absolute luminosity, evolved stars can be traced across a large fraction of the Galactic disc. At 3 μm, ATLAS Galactic Plane Survey will detect Red Giant Clump stars standard candles with L _mag = −1.75 (0.71 in AB mag) at 10kpc with A_v = 30 mag, reaching the heavily obscured regions of the Galactic centre or the outer edges of the Milky Way over at least the two outer quadrants. Clusters of Red Supergiants, ∼100 times brighter, can be studied via spectroscopy across the entire Galactic disc. Late-type stars tend to produce significant amount of dust and may have colours similar to YSOs and galaxies; IR spectroscopy is needed to disentangle their nature (Figure 14, right). The more massive AGB stars (initial mass 4–8 M_⊙) are sufficiently short lived that they can be used to trace the spiral arms of the Galaxy. Their number and distribution provides a fossil record of the recent history of star formation in the Milky Way.

Other objects. ATLAS will provide an exceptional database for studying a wide variety of sources: bubbles and young SN remnants, infrared dark clouds, heavily embedded high-mass clusters, planetary nebulae, old SN remnants, variable stars including outbursting YSOs, extremely metal-poor halo stars, and even galaxies and clusters hiding beyond regions of high extinction, identified as such by their spectra.

6.1. KBO model

To estimate the 1–4 μm brightness of a non-volatile rich KBO, we avoid the specific compositional modelling approach presented in Parker et al. (Reference Parker2015), but rather, turn to other well studied planetesimals as examples of what we might find in the Kuiper Belt. Our best analogy for the IR behaviour of a KBO is the spectrum of Phoebe (see Figure 16), which we believe has originated from the same primordial population that also populated the Kuiper Belt (Nesvorný, Vokrouhlický, & Morbidelli Reference Nesvorný, Vokrouhlický and Morbidelli2007; Levison et al. 1491 Reference Levison, Morbidelli, Van Laerhoven, Gomes and Tsiganis2008). Phoebe has likely avoided temperatures that would cause water–ice to become volatile, and therefore likely presents a spectral behaviour similar to that of KBOs in the water dominated 3–4 μm region. While the optical colour is significantly less red than many KBOs, possibly due to increased volatile loss of the organic compounds that possibly contribute to the development of a red refractory or organic coating on KBOs (Sheppard Reference Sheppard2010), Phoebe still exhibits water–ice absorption depths and 1–4 μm colours compatible with small KBOs (Fraser & Brown Reference Fraser and Brown2018). As far as available observations suggest, Phoebe appears as an excellent KBO analogue. For our modelling, we adopt the spectrum of moderate water absorption depths, but we point out that from the range of spectra exhibited on Phoebe, the global variations in 1–4 μm reflectance exhibited by KBOs are possibly as large as exhibited across Pheobe.

Figure 16. (a) Spectra of the different spectra units found on Phoebe’s surface. These spectra are from Fraser & Brown (Reference Fraser and Brown2018).. ‘Poor’, ‘Moderate’, and ‘Rich’ represent different levels in the water absorption depth. In our observable KBO density estimates, we adopt the moderate spectrum which exhibits water–ice absorption typical of most KBOs. (b) NIR colours of KBOs gathered with the Wide Field Camera 3 on the Hubble Space Telescope (black points; Fraser, Brown, & Glass 2015)). Simulated colours of Phoebe of the spectra presented in panel (a) are shown by the coloured points. The red point is that drawn from the moderate spectrum shown in panel (a).

KBOs are known to exhibit two compositional classes (e.g. Fraser & Brown Reference Fraser and Brown2012; Peixinho et al. Reference Peixinho, Delsanti, Guilbert-Lepoutre, Gafeira and Lacerda2012), with distinctly different optical colours. We model the 1–4 μm detectability of each class using two separate optical spectral slopesFootnote ^k: 10% for the neutral class, and 30% for the red class (Peixinho, Delsanti, & Doressoundiram Reference Peixinho, Delsanti and Doressoundiram2015; Pike et al. Reference Pike2017; Fraser et al. Reference Fraser2017). We model the 1–4 μm reflectance as approximately flat, as KBO spectra are known to be roughly neutral for λ ≳ 1 μm. In estimating the KBO 1–4 μm flux then, we simply scale the relative reflectance of the Phoebe spectrum at 1 μm according to the optical spectral slopes of spectral class. We also adopt an intrinsic neutral:red population ratio of 7:1 (Schwamb et al. in preparation).

6.2. ATLAS Solar System Survey

ATLAS Solar System survey is a pointed survey that targets individual known and tracked objects, in order to maximise the overall number of KBOs observed. We consider the AB = 22.5 (3σ) depth achieved in 2 500 s exposure time, since a spectrum of sufficient quality is gathered that not only enables the detection of water–ice absorption, but also allows the differentiation between water–ice absorption and hydroxyl absorption. Phoebe’s 3 μm hydroxyl absorption feature is a ∼50% absorption feature, while the ice feature is much shallower, with depth of only ∼10%. Thus, the latter sets our SNR requirement at roughly SNR > 10. From resolution considerations, the instrumental spectra can be binned to a lower resolution by a factor of ∼15 while still preserving sufficient spectral resolution to detect each of these broad-band features. With the above colour considerations, translating Phoebe’s 2.5 micron reflectance to r-band, the single exposure limiting magnitudes are r′ = 22.7 and r′ = 23.2 for the neutral and red classes, respectively.

At the AB = 22.5 (3σ) continuum depth limit, the number of detectable KBOs over the entire sky is ≳3 000. The vast majority of these will be observed and tracked by the LSST, allowing direct pointed observations by ATLAS for these targets. To observe each target with 2 500 s exposure would require roughly 2 083 h of on-target observations. Assuming 95% observing efficiency, a dedicated KBO ATLAS survey that is brightness complete to AB = 22.5 would require roughly 0.25 observing years.

6.3. KBOs from ATLAS Wide Survey

Now we estimate the number of KBOs that will be observed in the ATLAS Wide Survey. For this, we consider the (stacked) 3σ depth of 23.4 AB mags of ATLAS Wide. Each field in ATLAS Wide will be observed for 10 000 s (10×1 000 s exposures). Since the KBO targets are known and small in number, we will keep these on the target lists for ATLAS Wide, and stack the spectra gathered over 10 exposures of each KBO target in the wide survey. The stacked depths for the neutral and red classes become r′ = 23.6 and r′ = 23.1, respectively. These depths represent the practical limit achievable by the ATLAS Probe. Conveniently, all KBOs with r′ < 23.6 will have already been detected and tracked by the LSST, providing fully sufficient ephemerides for direct acquisition of spectra with ATLAS.

To determine the number of KBOs that will be observed in the ATLAS Wide Survey, we execute a survey simulation. We consider the nominal WFIRST wide survey region as presented in Holler.Footnote ^l As a Kuiper Belt model, we make use of the Canada France Ecliptic Plane Survey (CFEPS) model v7 (Petit et al. Reference Petit2011), which is the most up to date outer Solar System model. We note that this model does not include certain highly populated Kuiper Belt dynamical classes, such as the scattered disc and objects in the 5:2 mean-motion resonance with Neptune (Gladman, Marsden, & Vanlaerhoven Reference Gladman, Marsden, Vanlaerhoven, Barucci, Boehnhardt, Cruikshank, Morbidelli and Dotson2008). As such, detectability estimates drawn from this model likely represent lower bounds on the true detectable number of targets.

From the CFEPS model, at the stacked depths of ATLAS Wide, the number of KBOs detected will be ≃250 and ≃20 for the neutral and red classes, respectively.

8.1. Assumptions and requirements

We assume that the 2 048 × 1 080 DMD projects onto a 4 096 × 4 096 detector with 2× magnification, and that spectra are dispersed along detector columns (see Section 10.3 and Figure 24). To avoid overlap of adjacent spectra, targets must be separated by at least one column of ‘closed’ DMDs. Thus, a maximum of 1024 DMD columns (per quadrant) may be used to observe targets. At low spectral resolution, the projected length of spectra on the detectors is short enough that more than one target per DMD column can be observed without spectral overlap.

Spectra are assumed to span one octave in wavelength, e.g. from 1 to 2 microns, or 2 to 4 microns. We consider a range of instrument spectroscopic resolutions, R ≡ λ/Δλ, from R = 400–1 600, defined at the mid-point of the spectral range, with the instrument producing spectra with constant dispersion in nm/pixel. Note that the baseline ATLAS design assumes R = 1 000. The dispersion scale (nm/pixel) and total spectral length on the detector (pixels) is set assuming that the resolution element is sampled by 2.5 pixels in the dispersion direction. The total spectral length in detector pixels is n _spec = 5R/3 e.g. n _spec = 1 667 pixels for R = 1 000. The maximum spectral length that ensures that the full spectral range falls onto the detector for any DMD mirror is 1 936 pixels. Therefore, for spectral resolutions R > 1 162, spectra for targets located in the top rows of DMD mirrors will not fall entirely on the detector. In the present simulations, target positions are limited to the portion of the DMD where full spectral coverage is possible.

For spectral resolution R < 1 162, targets may be placed anywhere within a DMD column, whereas at higher resolution the top rows must be avoided in order for the entire dispersed spectrum to fall on the detector. The maximum column multiplex, the maximum number of galaxies that can be targeted per DMD mirror column without spectral overlap in the dispersion direction, is as follows:

(1) \begin{equation} M^C_{\rm max} = 1 + \mathrm{int}[j_{\rm max} \times 6 / (5R)]\end{equation}

where j _max is the maximum DMD row allowed for target placement, i.e. j _max = 1 080 for R < 1 162, and j _max = 2 048−(5R/6) for R > 1 162. Table 4 reports the maximum column multiplex ( $M^C_{\rm max}$ ) and total DMD multiplex ( $M^T_{\rm max} = 1\,024 M^C_{\rm max}$ ) for various ranges in spectral resolution R. In practice, the achieved multiplex is generally lower and depends on the underlying density of potential targets. When that density is high, there is a greater chance that potential targets will be found in the right places to fit several spectra into a given DMD column.

Table 4. Maximum target multiplex as a function of spectral resolution

We further assume that targets should be fairly well centred on a DMD mirror in order to be observed spectroscopically. If not, a significant fraction of their light may be lost outside the mirror and not directed to the spectrograph. Placing a target near the mirror edge may also have deleterious consequences for the effective PSF of the spectrum recorded on the detector. In the present simulations, we require that a target falls within the central 25% surface area of a DMD mirror; i.e. if the fractional x, y position of a source within the mirror ranges from 0 to 1 in each axis, a target can only be selected if it falls in the range 0.25 < (x, y) < 0.75.

The mirror centring criterion and the use of alternate DMD columns to avoid spectral overlap lead naturally to an observing strategy that cycles through a dithered pattern of four telescope pointings. Defining the angular coordinates of an arbitrary pointing as (0,0), the four pointings would have coordinates of (0,0), (0, 0.5θ _DMD), (0.5 θ _DMD, 0) and (0.5 θ _DMD, 0.5 θ _DMD) where θ _DMD is the angular size of a DMD mirror projected on the sky (0.754 arcsec for f/2.5). Additionally, two exposures will be taken at each pointing with a different set of micro-mirrors turned on to allow observation of galaxies in DMD columns that were previously excluded in order to avoid overlap of adjacent spectra. Half-step dithering in x and y combined with the two exposures per pointing ensure that all galaxies will fall within the central 25% surface area of one of the DMD mirrors. A minimum of 8 exposures will then be taken, more if there are still targets to observe, most probably in that case in multiples of 8 exposures.

8.2. Target densities

We assume that a specified number (N _targ) of potential targets are randomly distributed (i.e. with no clustering) over the DMD FOV. The present simulations consider target densities in the range log N _targ = 3.4 to 4.8 per DMD FOV. For a 1.5 m telescope with f/2.5, the DMD FOV is 0.0969 deg², and the simulations correspond to galaxy surface density from 10^4.4 to 10^5.8 deg⁻². From NIR (1.6 micron, F160W) source counts in the GOODS-S CANDELS-Deep region (Figure 18), these cumulative densities would approximately correspond to magnitude limits from m _F160W = 21.5 to 27. Actual ATLAS targeting strategies could preselect ranges of photometric redshifts, etc., so for the present purposes we represent the simulation results as a function of the number of potential targets per DMD FOV, rather than by limiting magnitude.

Figure 18. 1.6 micron (F160W) source counts from the CANDELS ultra deep region of GOODS-S. Upper panel: differential source counts, N(m), with no correction for incompleteness. The raw counts become significantly incomplete at m _F160W > 27 mag. Lower panel: cumulative source counts, N(<m).

8.3. Simulating target selection

In the present simulations, targets are selected for spectroscopy following a very simple procedure that is not optimised in any formal sense. Starting at DMD column i = 0 (zero indexed), all galaxies that fall within that column and which meet the mirror centring criterion are identified. A viable target falling in the lowest row (index j) is selected for observation. Then, the j index is incremented by the spectral length in DMD mirror units (n _spec/2), in order to avoid selecting any other galaxies whose spectra would overlap that of the previous object. If there are viable targets that fall within the remaining rows of that DMD column, the object in the lowest-j row is also selected as a target, and the process is repeated until the maximum DMD row (j = 1 079) is reached, or, for higher spectral resolutions, the maximum row for which the full spectral range would be projected onto the detector (see above). Next, if one or more targets are selected for a given DMD column, then the column index is incremented by two, i.e. an empty column of mirrors is left as a buffer, as described above. However, if no targets are selected in a given DMD column, then the column index is incremented by one. This process is repeated over all DMD columns, until targets have been selected over the full DMD FOV. In practice, the target selection process will be optimised to remove selection biases and to measure the sky through a small number of slits distributed across the field.

10.1. Preliminary optical design

Our preliminary optical design is based on the design originally developed by Content (Reference Content2008) and Content et al. (Reference Content2008) for the SPACE mission proposal (Cimatti et al. Reference Cimatti2009) and the following pre-Phase-A and Phase-A studies of the ESA Euclid mission, which resulted from the merger of the SPACE and DUNE proposals. Euclid did not adopt the concept of slit mask MOS preferring the scientifically less ambitious slitless spectroscopy approach that can be more simply added to its NIR camera.

The ATLAS telescope is a 1.5-m f/11.2 modified Ritchey–Chrétien design similar to the SPACE telescope design. The telescope focal plane falls on a four-faceted pyramid mirror, a solution adopted, e.g. by the WFPC2 camera on the HST. A full view of the ATLAS Probe instrument is shown in Figure 21. The 4 faces of the pyramid mirror act as pickoff mirrors (Figure 22, left) sending the light to a fore-optics system. The 4-petals shape of the field allows for almost perfect tessellation of the sky, except for 2 small squares totalling only 0.07% of the area, due to the aspect ratio of the DMD which is different from the ideal 2:1 ratio (full field coverage will be achieved through dithering). Our approach also permits the split of the instrument into 4 identical and separable blocks, which helps simplify the assembly, integration and testing.

Figure 21. A full view of the preliminary optical design for the ATLAS Probe instrument. The large grey circle is the back of the primary.

The 3 mirrors fore-optics (Figure 22, right) receive the f/11.2 telescope beam and release a ≃f/2.5 × f/2.3 beam on the DMD, accommodating for the 12° micro-mirror tilt along the diagonal. The DMD acts as an addressable multi-slit mask, reflecting to the spectrograph the light of the science targets (+12°, ‘on’ beam) and sending the rest of the field (−12°, ‘off’ beam) to a light dump.

Figure 22. Preliminary optical design for the ATLAS Probe instrument. Left: Pickoff mirrors in focal plane; positions and shapes are approximate. Right: Shade layout of the fore-optics.

The spectrograph (Figure 23) includes a four mirror collimator followed by a beam splitter reflecting (1–2.1 μm) and transmitting (2.1–4 μm) light from the beam to a short and long-wavelength channel, respectively. Similar beam splitters (with cut-off at 2.4 micron) have been adopted, e.g. by NIRCam on JWST. We take advantage of the relatively modest resolving power R = 1 000 to make use of prisms instead of gratings as dispersive element. Prisms provide a series of advantages over gratings:

Figure 23. Preliminary optical design for the ATLAS Probe instrument. Shade layout of the spectrograph; the short wavelength camera is on the opposite side.

1. (1) They have a much better transmission especially at the extreme wavelengths of the spectra.

2. (2) They can have a wavelength range larger than an octave between the extreme wavelengths. In this preliminary design the short wavelength camera has a wavelength range from 1 to 2.1 μm, with the possibility of further increasing the total wavelength range up to 0.9–5 μm in a future trade-off study.

3. (3) The resolution is more uniform as a function of wavelength.

4. (4) There is no need for a sorting order filter since there are no orders.

5. (5) There are no same-bandwidth parasite orders, especially 0 and 2, cross-contaminating spectra aligned with each other. These would be unusually strong with our bandwidths where the longest wavelength is the double of the shortest.

Each prism directs the beam toward its camera which is made of only one mirror, an unusually low number. This is obtained by making the beam on the prism non-telecentric and correcting the resulting aberrations with the collimator and camera.

For target acquisition it is necessary to directly image the field; as originally proposed for SPACE/Euclid, direct imaging can be achieved through an achromatic grism where grating and prism dispersions cancel each other, with excellent image quality. ATLAS has two cameras, but only one is needed for target acquisition. Since the ATLAS short wavelength camera has a wavelength range only slightly different from Euclid, we envision a slightly modified design of the Euclid achromatic grism.

Our preliminary design already has excellent image quality in the spectrograph. The Gaussian equivalent full width at half maximum (GEFWHM) is about 1/2 of a DMD mirror image. The image quality of the fore-optics on the DMD is 1/2 of a DMD mirror, satisfactory at this point of the preliminary design.

To attain the performances described above we make use of aspheric surfaces, some of them quite complex but still well inside the limits of modern manufacturing.

Table 6 lists the main parameters of our system. Figure 24 shows the focal plane layout, with the DMD reimaged onto the lower part of the detector and the detector’s upper part left ‘dark’ to capture the full spectrum of each source, regardless of its position in the reimaged field.

Table 6. ATLAS main optical parameters

Figure 24. Layout of the focal planes for ATLAS Probe. Left: Focal plane projected on the DMD. Right: Focal plane project on the IR detector. The DMD is reimaged onto the lower part of the detector, and the detector’s upper part is left ‘dark’ to capture the full spectrum of each source, regardless of its position in the reimaged field. The moderate lengths of the R ≃ 1 000 spectra (1 333 pixels for 2 pixels, or 1 667 pixels for 2.5 pixels per spectral resolution element sampling) allows for stacking multiple spectra on the same column, increasing source multiplexing. In order for the full spectrum to be captured by the IR detector for any target position in the field, the maximum spectrum length is 1 936 pixels, corresponding to about 2.9 pixels per resolution element.

10.2 Detectors and ASICs

Our baseline detector is the Teledyne H4RG-10, the same type currently under development for WFIRST. The long-wavelength cut-off of our spectroscopic channels (2.1 and 4 μm) match the standard ∼2.35 μm and ∼5.37 μm cut-off of WFIRST and JWST devices, respectively. Fine-tuning the Hg vs. Cd stoichiometric ratio can further reduce dark current, allowing warmer operating temperatures. H4RG arrays tested by WFIRST have typical QE > 90% in the 0.8–2.35 μm range, readout noise ∼15e in Double Correlated Mode and mean dark current <0.01 e s⁻¹ pixel⁻¹. JWST devices have similar performance. The architecture of the multiplexer allows multiple non-destructive reads of each pixel during a single exposure (‘sampling up the ramp’); this enables mitigation against cosmic rays and reduction in read noise: typical readout noise with 16 samples drops to about 5e. Teledyne has developed SIDECAR ASICs (Application Specific Integrated Circuit) to manage all aspects of FPA operation and output digitisation in cold environment. By keeping analogue signal paths as short as possible to reduce noise and output capacitance loading, ASICS improve power consumption, speed, weight and performance. A Teledyne ASIC device is currently driving HST/ACS and several Teledyne ASICs will soon fly on three out of four JWST instruments; we will adopt such devices for ATLAS.

10.3. DMDs and controller

The DMD is the core of our system. We base our design on the 2k CINEMA model of 2 048 × 1 080 micro-mirrors, 13.7 μm on a side. DMDs have been successfully used on ground-based spectrographs like RITMOS (Meyer et al. Reference Meyer2004) and IRMOS (MacKenty et al. Reference MacKenty2006). A new DMD-based spectrograph, SAMOS, is under construction of the SOAR telescope in Chile (Robberto et al. Reference Robberto2016). NASA has funded a Strategic Astrophysics Technology (SAT) programme (PI Ninkov, ATLAS Probe core team member) to raise the TRL level of DMDs to TRL5-6 before the 2020 decadal survey. DMDs have successfully passed proton and heavy-ion irradiation testing (Fourspring et al. Reference Fourspring2013; Travinsky et al. Reference Travinsky2016). Following NASA General Environmental Verification Standard (GEVS), the team performed random vibration, sine burst, and mechanical shock testing of manufacturer-sealed DMDs (Vorobiev et al. Reference Vorobiev2016) as well as of devices re-windowed for better UV and IR capabilities (Quijada et al. Reference Quijada2016). These GEVS tests suggest that DMDs are robust and insensitive to the potential vibroacoustic environments experienced during launch. Low-temperature testing of DMDs were also performed, the main concern being micro-mirror stiction. Tests at RIT have shown that temperatures as low as 130 K do not affect the performance of DMDs. More recent data obtained at JHU at 80 K confirm these results. A general overview of the results obtained by these test campaigns has been recently presented (Travinsky et al. Reference Travinsky2017).

These findings place DMDs between TRL levels 5 and 6 (note that level 6 requires testing against the specific environment of a mission, thus is only achievable after a mission has been identified, e.g. ATLAS at L2). DMDs are normally controlled by commercially available boards based on the DLP Discovery 4100 chipset from Texas Instruments. The chipset is not designed to operate in a cryogenic environment or in space. ATLAS core team members at JHU have recently developed custom electronics to operate Cinema DMDs at cryogenic temperatures. For ATLAS, we intend to produce a version of this system based on rad-hard components, suitable for operations at L2.

Micro Shutter Arrays (MSA) may represent an alternative to DMDs. However, their limited size implies that one would need 142 devices similar to those on JWST/NIRSpec to obtain the same number of selectable slits achieved with four DMDs. Also, the filling factor ≲70% of their clear apertures greatly complicates operations, since a large fraction of the visible targets ends up partially occulted. If future developments will lead to the construction of MSAs with ∼30 times more apertures, clear-aperture ratio larger than 90%, and cosmetic quality comparable to those of DMDs, we will evaluate their use in the ATLAS Probe design.