The Galaxy Evolution Probe (GEP) is a concept for a mid- and far-infrared space observatory to measure key properties of large samples of galaxies with large and unbiased surveys. GEP will attempt to achieve zodiacal light and Galactic dust emission photon background-limited observations by utilizing a 6-K, 2.0-m primary mirror and sensitive arrays of kinetic inductance detectors (KIDs). It will have two instrument modules: a 10 to 400 μm hyperspectral imager with spectral resolution R = λ / Δλ ≥ 8 (GEP-I) and a 24 to 193 μm, R = 200 grating spectrometer (GEP-S). GEP-I surveys will identify star-forming galaxies via their thermal dust emission and simultaneously measure redshifts using polycyclic aromatic hydrocarbon emission lines. Galaxy luminosities derived from star formation and nuclear supermassive black hole accretion will be measured for each source, enabling the cosmic star formation history to be measured to much greater precision than previously possible. Using optically thin far-infrared fine-structure lines, surveys with GEP-S will measure the growth of metallicity in the hearts of galaxies over cosmic time and extraplanar gas will be mapped in spiral galaxies in the local universe to investigate feedback processes. The science case and mission architecture designed to meet the science requirements is described, and the KID and readout electronics state of the art and needed developments are described. This paper supersedes the GEP concept study report cited in it by providing new content, including: a summary of recent mid-infrared KID development, a discussion of microlens array fabrication for mid-infrared KIDs, and additional context for galaxy surveys. The reader interested in more technical details may want to consult the concept study report.
Galaxy Evolution Probe in Context
Tracing the mass assembly history of galaxies is an essential component of understanding the origins of the Hubble sequence and of using galaxies as cosmological probes of dark matter and dark energy. Observational and theoretical work over the past three decades have established a baseline framework for galaxy assembly in a cosmological context. These studies have shown that stellar and black hole mass assembly varies substantially with redshift, increasing by more than an order of magnitude between the local universe and ,1–3 though the behavior at higher redshifts is uncertain.4 Environment and large-scale structure profoundly affect galaxy assembly, with factors such as local galaxy density and dark matter halo properties known to play key roles in shaping the galaxy mass function.5–8 Finally, there is a complex and subtle relationship between star formation and active galactic nucleus (AGN) activity,9 including both positive and negative feedback effects.10–13 Infrared surveys from facilities including IRAS,14 ISO,15 AKARI,16 Herschel,17 WISE,18 and Spitzer19 have played fundamental roles in these studies, since a substantial fraction of all galaxies over the history of the universe were permeated with dust through their most active evolutionary stages.
Building upon decades of observational and theoretical work, collectively a data-driven, self-consistent model for galaxy evolution that starts from cosmology, incorporates stellar dynamics and evolution, and includes interstellar processes, star formation, and supermassive black hole growth is within reach. Achieving this goal requires addressing the questions that have arisen from previous and current generations of multiwavelength surveys. These include: what was the role of feedback from black holes and stars themselves in regulating star-formation? How did galaxies’ external environments and internal contents influence their evolutionary trajectories? When were the Universe’s heavy elements formed by stars in galaxies and how did they escape into the circumgalactic and intergalactic medium?
Answering these questions will require large panchromatic surveys measuring bulk properties of hundreds of thousands to millions of galaxies over most of cosmic history paired with detailed high angular and spectral resolution studies of gas and star formation in individual galaxies. Such large samples are needed to precisely disentangle the effects of redshift and environment in driving galaxy assembly. High-resolution observations of representative and outlier galaxies identified in surveys trace detailed physical processes on a galaxy-by-galaxy basis, enabling astrophysical processes to be mapped onto cosmological processes.
New mid- and far-infrared observations enabled by the rapid advances in infrared detectors and technology20 are a vital component of these next-generation surveys. They require sensitivity to detect Milky Way-type galaxies at (), prior to when most stellar mass had been assembled. Star formation rates and supermassive black hole accretion rates should be measured over a full range of cosmic environments, from isolated field galaxies, to galaxies in groups and in massive clusters. Measuring the heating of dust and gas by star formation and AGN is needed to understand how stellar and supermassive black hole growth were linked over cosmic time. Spectral mapping capability is needed for unbiased spatial–spectral surveys, line luminosities and intensity mapping, and line mapping of nearby galaxies to measure gas column densities, ionization parameters, and metallicities with tracers unaffected by dust obscuration. Crucially, galaxies detected by their dust continuum emission must have measured redshifts so their epochs and luminosity distances are known and their multiwavelength counterparts can be identified.
The Galaxy Evolution Probe (GEP) is a NASA Astrophysics Probe concept that capitalizes on new detector capability to address these questions with a powerful mid- and far-infrared toolset. GEP will measure star-formation rates and detect AGN even under conditions of heavy dust extinction. It will measure supermassive black hole accretion rates to address the connection between the masses of stellar populations and supermassive black holes. The same observations will measure metallicities with extinction-free tracers to observe growth of metals over the last two-third of the Universe’s age. In nearby galaxies, GEP will observe feedback between star-formation, AGN, and the interstellar medium to understand the processes that regulate star-formation. Mapping nearby galaxies and the Galactic interstellar medium will reveal the energy balance by measuring the total interstellar material mass, ionization state, and the local radiation field using fine-structure transitions of ions, polycyclic aromatic hydrocarbon (PAH) molecules, and the mid-infrared dust continuum.
GEP (Fig. 1) is designed with a 2.0-m primary mirror that will be cooled to 6 K to enable sensitivity limited by photon shot noise from foreground astrophysical sources: zodiacal dust emission and Galactic dust emission. There is one scientific instrument with two modules: an imager, GEP-I, and a dispersive spectrometer, GEP-S. GEP-I has 23 photometric bands distributed on the focal plane: 18 resolution bands from 10 to designed to measure redshifts with PAHs and mid-infrared spectral energy distributions, and five resolution bands from 95 to to measure dust spectral energy distributions encompassing the peak to beyond . GEP-S utilizes four long-slit grating spectrometers with spectral resolution from 24 to . Both modules baseline arrays of kinetic inductance detectors (KIDs) cooled to 100 mK by a multistage adiabatic demagnetization refrigerator (ADR) backed by a hybrid Joule–Thomson and Stirling cryocooler, which will also cool the telescope and coupling optics. GEP launch is targeted for January 1, 2029, with a planned mission duration of 4 years at Earth–Sun L2 (Table 1).
Basic GEP parameters.
GEP is optimized for large, multitiered surveys for galaxies detected by their mid- and far-infrared emission from dust, PAHs, and atomic fine-structure lines. The GEP reference mission includes two types of surveys: photometric hyperspectral surveys with GEP-I and spectroscopic surveys with GEP-S. The GEP-I survey areas will be 3, 30, and 300 square degrees and an all-sky survey (Fig. 2). The spectral surveys with -level redshift precision will cover a range of low- and high-ionization atomic fine-structure lines. Spectral surveys will consist of “blind” surveys utilizing a long-slit configuration, intensity mapping, and follow-up, deep pointed observations of galaxies identified in the GEP-I surveys, and regions of the Milky Way and nearby galaxies. GEP will require 4 years to do these surveys; however, it utilizes no expendable cryogens, and a Guest Observer mission can be envisioned after the surveys are completed.
GEP shares common elements with two other cryogenic far-infrared space missions under study: SPICA, an ESA-JAXA collaboration in a similar cost class to GEP, and Origins Space Telescope (OST), a NASA concept study for a future flagship infrared mission. SPICA’s design features a 2.5-m, 6- to 8-K telescope and a spectrometer using sensitive bolometer arrays, although it does not have the 23-band, moderate spectral resolution hyperspectral imaging that is central to GEP’s architecture. OST is designed with a 5.9-m, 4.5-K telescope; its earliest possible realization is well into the 2030 decade.
GEP was described in the concept study report submitted to NASA,21 and its science case was briefly summarized in a short white paper submitted to the 2020 Astronomy Decadal Survey.22 The initial design concept for the mission and the optical design were described by Glenn et al.,23 and the cryogenic design was described by Moore et al.24 This paper summarizes the design and science and provides new content: new engineering design, a summary of recent mid-infrared KID development, a discussion of microlens array fabrication for mid-infrared KIDs, and additional context for galaxy surveys.
GEP has two overarching science goals. The first is to map the history of galaxy growth through star formation and accretion by supermassive black holes and to characterize the relationship between these processes. The second is to measure the growth of metals in galaxies and the changing of star-formation interstellar medium environments in galaxies over cosmic time. These science goals are broken down into specific objectives, which are translated into requirements that drive the GEP design in the concept study science traceability matrix.21
Broadly, the star formation began in the first billion years of the Universe, rose to a peak or broad plateau around , and then declined sharply.2 Despite substantial success over the past decades in understanding the average star formation and supermassive black hole accretion rates with redshift, large uncertainties and significant questions about galaxy evolution still exist, including: star formation rates derived from infrared galaxy surveys of limited size have not probed low luminosities25 (essentially, below ) and are limited by sample variance, redshifts of large samples of far-infrared continuum-detected galaxies are uncertain, and extinction correction in rest-frame ultraviolet observations is substantial and sometimes uncertain.2 GEP is designed to address these concerns with large, deep infrared surveys for star-forming galaxies to meet its first objective: to measure the coevolution of star formation and supermassive black hole growth in galaxies. GEP-I will yield star formation rates from far-infrared luminosities26,27 and separate AGN and star formation contributions to total infrared luminosities via mid-infrared spectral shapes,28 even in cases of heavy extinction because the dust emission is almost always optically thin. The 23-band GEP-I photometry of the mid-infrared PAH features will yield redshifts for the galaxies. PAHs have been detected spectroscopically with Spitzer at high redshifts, e.g., and 2.9629 and ,30 in galaxies with and without prominent AGN. From mid- and far-infrared fine-structure lines, GEP-S will measure star formation rates31,32 and supermassive black hole accretion rates.32
GEP will access the early epochs of galaxy growth () by utilizing the brightening from gravitational lensing. It is clear that much of the star formation and black hole growth in massive galaxies since reionization occurred in dusty regions, and there is considerable uncertainty about how much star formation rate density may be missed in the ultraviolet census of that early phase of galaxy evolution.2 Wide-area GEP-I surveys taking advantage of the brightening provided by gravitational lensing will address this problem with lensed high- galaxies.
Understanding the role of accreting supermassive black holes in galaxy evolution requires infrared observations because these processes occur on small scales at the centers of galaxies, in most cases obscured by the very gas and dust that is accreting. A principal objective of GEP is to identify obscured AGN in galaxies and relate their accretion luminosities to their star formation rates. GEP will identify and quantify luminosities of dust-obscured AGN via mid-infrared spectral signatures, including the continuum shape with GEP-I and high-ionization fine-structure atomic transitions with GEP-S, such as [Ne V] at rest-frame wavelengths of 14.3 and and [O IV] at . GEP complements x-ray detection of AGN because x-ray observations can miss Compton-thick AGN or underestimate accretion rates.
The masses of supermassive black holes in the centers of modern-day galaxies are well correlated with galaxies’ bulge masses.33,34 This has led to the now commonly adopted hypothesis that a feedback loop exists in which AGN activity governs the rate of star-formation in massive galaxies or at least in galactic bulges.35 Theoretical models10,36–38 invoke AGN feedback as a primary mechanism to explain the observed distribution of galaxy masses today. Without AGN feedback, models are unable to explain the low ratio of galactic stellar mass to halo mass for high-mass galaxies.39 However, the efficiency of AGN feedback for regulating star formation has not been established empirically and remains controversial.12,40 GEP’s goal is to determine how feedback from buried accreting black holes was related to the decline of star-formation. GEP will assess the role of feedback in galaxies by searching for faint wings or asymmetries in the velocity profiles of bright fine-structure lines in stacked spectra of AGN and starburst samples. Such line-wing signatures have been measured interferometrically in CO (e.g., Ref. 13) as well as in the optical (e.g., Ref. 41). They can have velocities in excess of , particularly for the AGN sources. Given the modest spectral resolving power of GEP-S ( or ), GEP’s far-IR measurement will be accomplished with a combination of (a) high signal-to-noise ratio ( or greater) in the main line core for either the high-redshift spectral stacks or individual nearby galaxy spectra and (b) knowledge of the spectral response function of the instrument, with fidelity to 1 part in 500 or better. The excellent raw sensitivity of GEP-S, combined with the fact that we are considering the brightest lines in the spectrum, will satisfy the first condition. For the second, a combination of laboratory measurements preflight and measurements of narrow Milky Way line-emitting standards should provide ample opportunity for accurate spectral calibration. In nearby galaxies, GEP will directly obtain a spatially resolved view of feedback and its effects with a spectroscopic study of galactic outflows and fountains in local galaxies in various atomic fine-structure emission lines (e.g., [C II], [N II], and [O III]).
GEP is designed not only to measure the star formation rates of galaxies over cosmic time but also the evolving conditions of star formation regions. Specifically, the measurements of fine-structure lines, including [C II], [N II], and [O III], will be used to infer the mass and density of the interstellar medium gas and the hardness of the ambient radiation, which has implications for the stellar initial mass function, and the density of H II regions from which the interstellar pressure can be inferred. The resulting relations between metallicity, star-formation rate, and other galaxy properties will inform models of galactic winds by placing constraints on the presence of gas-phase metals at their source.
Metallicity represents the integrated effects of star-formation, inflow, and outflow from galaxies. A key objective of GEP is to track the buildup of heavy elements in galaxies over the peak epoch of star-formation utilizing spectroscopic surveys of metals in atomic gas (and secondarily lower-resolution spectra of PAHs and dust). Metallicities of galaxies have not been measured beyond the local universe with the extinction-free probes afforded by far-infrared atomic fine-structure lines. GEP will measure the metallicity in galaxies down to galaxies at using the nitrogen-to-oxygen ratio measured with [OIII] and [NIII] fine-structure lines. This ratio measures metallicity through its connection to stellar processing since nitrogen is a secondary nucleosynthesis product and comes on later in the nucleosynthesis process. The nitrogen/oxygen approach has been used with optical lines but the far-IR lines offer the dual benefit of dust immunity and temperature insensitivity.42 Alternatively, for the highest redshifts, metallicities will be measured using a diagnostic formed from the mid-IR tracers of neon and sulfur. This has been shown to track metallicity in photoionization models and has been validated in a careful study of local galaxies that have Spitzer mid-IR spectroscopy in hand.43
Mission Architecture: Thermal Design and Observing Modes
GEP is designed as a Class B mission. A NASA Class B mission is a mission whose loss would be highly impactful to national science objectives and therefore that only low risks to mission success will be tolerated. Less risk is tolerated for Class A missions and more for Classes C and D. GEP is based on the Ball Aerospace BCP2000 reference bus, of Kepler heritage, customized to meet the science and mission requirements of GEP. The 4-year survey program is divided into a GEP-I campaign and a GEP-S campaign, from which 350 TB of observational data will be downlinked. Mission details can be found in the concept study report.21
GEP’s thermal system employs multiple passive and active stages to meet the temperature intercept requirements of the instruments and optical assembly (Fig. 3). A continuous multistage ADR provides cooling for the detectors at 100 mK, with a 1-K thermal intercept to reduce thermal noise and parasitic loads. A hybrid Joule–Thomson/Stirling cryocooler intercepts heat at 4 K from the ADR and from the cryogenic amplifiers and parasitic loads. The hybrid cooler has an 18 K intercept to cool a second stage of amplifiers for the detector signals and an active shield.
The sunshield assembly consists of three passively cooled reflective shields with a total area of and a active shield located above the sunshield under the focal plane boxes. The substrate of all four shields is an internally self-supporting 1-mm-thick M55J laminate with an areal density of . The reflective coating on the three warmer, passively cooled sunshields consists of a thin (0.005 in.) layer of aluminized Kapton adhered to the surface. The bottom surface of the actively cooled sunshield is also aluminized Kapton; however, the top surface is a 1-mm thick, high purity aluminum thermal spreader layer. Fiberglass composite struts provide intershield supports at the perimeter to accept launch loads. The passively cooled shields intercept conducted loads from the bipods and harnesses and radiative loads from the Sun. With these nondeployable sunshields, the spacecraft can tilt , maintaining all cryogenic components in the shadow cone and thus preserve thermal system operability. The cooler and ADR electronics dissipate heat at ambient temperature (), along with all bus components, which is radiated by of radiators mounted on the bus behind the solar panels.
The GEP-I and GEP-S instrument modules will observe one at a time and use the same readout electronics. They share a scan survey observing mode, where mapping is performed as the spacecraft slews at . The scan survey mode will be used for all the GEP-I surveys, for GEP-S’s 1.5 and 100 square degree surveys, and for mapping of nearby galaxies. GEP-S will also have a pointed observation mode in which a chopping mirror, with a total throw of , modulates the signal for noise mitigation. One 5-h DSN pass per day will be required to downlink 0.24 TB, during which the all-sky survey will be conducted by spinning the spacecraft while the high-gain antenna is pointed at Earth. The intent is to do the entire all-sky survey in this mode while spinning at the same rate as the data taken in the nominal scan mode. If it is determined that the scan rate must be different for the all-sky mode than for nominal scanning, a separate calibration will be required.
The GEP payload includes the optical telescope assembly, the GEP instrument comprised of a hyperspectral imager module (GEP-I) and four spectrometer modules (GEP-S), and a payload thermal subsystem (Fig. 4).
The fundamental optical design requirement is to collect light with a 2.0-m diameter primary mirror and to form an focus. An unobscured three-mirror astigmat (TMA) is an oft-used configuration that is well suited to the first-order optical and mechanical requirements of this system. The powered mirrors are all conic shapes with parent surfaces that have mutual tilts and decentrations to reduce wavefront error. The baseline for the GEP primary, secondary, and tertiary and chopping mirrors is unactuated silicon carbide (SiC), as was flown on the Herschel Space Observatory and Gaia. GEP’s optical design and simulated performance are described in detail in Ref. 23.
The five instrument interface planes are nearly coplanar at the common focal surface of the TMA, with minor differences in final focus to minimize wavefront error in each channel. Based on the sizes of the focal-plane array and the spectrometer slits and enclosures, the optimized field of view is , in excess of what is required for the instrument modules. The centers of the GEP-S slits are in the plane of symmetry of the TMA so that they can be untilted with respect to the central ray to each spectrometer. Stray light suppression is achieved with a pupil stop and a field stop.
To meet the extragalactic confusion requirement (Sec. 5.3), the primary mirror is required to be diffraction limited at . However, the root-mean-square (RMS) wavefront error across the field of view in the optical design is diffraction limited at the minimum wavelength of . The 100% encircled energy diameter is , much less than the 3.43″ detector pixel size (see Sec. 4.2).
GEP-I Hyperspectral Imager
GEP-I is designed to obtain repeated measurements in each of the 23 wavebands to build survey depth by continuously scanning areas of sky covered by each survey. Each waveband occupies the same amount of focal plane area, 0.002 square degrees, half of which is occupied by detectors and half of which is allotted for bandpass filter mounting in the current design configuration (Figs. 5). Bands 1 to 18 have spectral resolution , whereas bands 19 to 23 are . It is likely that the mid-infrared (e.g., 10 to ) bandpass filters could be replaced with linear-variable filters, resulting a spectral resolution of at least , which will be the subject of a future investigation.
GEP-I will have 25,735 KIDs. GEP I’s optical design performance is diffraction limited at ; however, the primary mirror is specified to be diffraction limited at , corresponding to (FWHM) beam size. The primary reason for this is that the smallest KID pixel size we expect to be able to fabricate without exceeding the readout bandwidth is (the number of detectors is inversely proportional to their physical area), which corresponds to 3.43″. With 3″ beams and 3.43″ pixels, the full angular resolution of the telescope will not be taken advantage of for GEP-I bands 1 to 9 (up to ), although all flux density will be recovered. For bands 9 to 23, the pixel size will be less than the beam FWHM. Should greater bandwidth become feasible through improvements in data acquisition and computing speed (a very likely development), the pixel sizes can be reduced to less than the beam FWHM.
GEP-S Long-Slit Spectrometer
GEP’s spectrometer was designed to meet the science requirements calling for observing mid- and far-infrared atomic fine-structure lines from galaxies over a range of redshifts. Specifically, the [Ne V] line starting at for AGN identification and the [O I] line at for probing star-forming galaxies at or beyond the cosmic peak. The entire bandwidth should be available to identify spectral lines for galaxies of unknown redshift. Sufficient spectral resolution is required to achieve good sensitivity through dispersion of the astrophysical background photons. Spectral resolution meets these requirements.
GEP-S’s design is implemented with four diffraction gratings, each with fractional bandwidth of 1.6, identified as bands 1 to 4 (Table 1), with a total of 24,640 KIDs. The ray trace diagram and CAD model of band 4 is shown in Fig. 6. Each band, or diffraction grating module, is comprised of an enclosure, a slit, a collimator, a grating operated in first order, a focusing mirror, and an array of KIDs. The long slit lengths enable spectral mapping. As with GEP-I, the shortest-wavelength bands 1 and 2 are placed nearest to the center of the field of view where the optical performance is the best. A chopping mirror [Fig. 4(a)] is included for detector modulation for staring GEP-S observations.
Kinetic Inductance Detectors and Readout Electronics
KIDs have been baselined for the GEP-I and GEP-S focal plane arrays. KIDs are high-quality factor superconducting resonators whose resonances vary with absorbed light through modulation of kinetic inductance by pair breaking.44 They are in wide use in recent and imminent ground-based and suborbital optical, far-infrared, submillimeter, and millimeter-wave instruments (Fig. 7).
Two other detector array technologies were considered for GEP: transition-edge sensors (TESs) and Si:As IBCs (impurity band conductors). KIDs are not as technologically mature as TESs but their simple fabrication and focal plane electronics make them advantageous, partially because TESs require complex hybridization with SQUID readouts. Like KIDs, TESs have not been demonstrated in the 10 to wavelength regime. Si:As IBCs are technologically mature but they are unusable beyond , so KIDs (or TESs) would still be needed for 28 to . However, pixel pitches and operating temperature requirements for IBCs (7 K) and KIDs (100 mK) are very different, and the divergent optical and cryogenic requirements would effectively necessitate two instruments each for GEP-I and GEP-S, which would drive up complexity, cost, and risk. Thus, KIDs have been adopted as the baseline detector technology for all of the wavebands of GEP-I and GEP-S.
GEP-I and GEP-S will utilize backilluminated, lumped-element, microlens-coupled aluminum KIDs. Each KID is a superconducting thin-film microresonator, comprised of an absorbing meandered inductor and an interdigitated capacitor deposited on a silicon substrate. Light is concentrated onto the inductors by a microlens array. The geometry of the inductor controls the absorption characteristics as a function of wavelength and polarization. Absorption of photons changes the kinetic inductance, producing a small shift in the resonant frequency and resonance depth, which are detected in the transmitted probe signal. To enable frequency-multiplexed readout with a minimal bandwidth requirement, the interdigitated capacitors are unique for each pixel, yielding different resonant frequencies separated by several resonance line widths. Submillimeter-wavelength KIDs with sensitivity sufficient for GEP-I have been demonstrated with a different architecture.45 NEP improvement, not just different operational wavelengths, is needed for GEP-S.
GEP KIDs are based on the MAKO design49 and adapted to work at shorter wavelengths. The MAKO design will work down to a wavelength of , but below that modification will be required: the wire grid absorber linewidth used in the MAKO detectors is . Scaling this design with wavelength would require prohibitively narrow lines for , even if using an ultraviolet stepper for film deposition. Our approach is to interrupt the absorber line with a meander that increases the resistance per unit length (and therefore the absorption efficiency), while simultaneously introducing a distributed capacitance to compensate for the accompanying inductance to tune the impedance.63 An impedance match to silicon may then be maintained with much wider and easier-to-fabricate lines.
This absorber design is shown in Fig. 8. For radiation, the linewidth is 200 nm and the unit cell is square, repeated continuously in the circular envelope of the absorber with electrical continuity meandered vertically. Simulations indicate that this geometry achieves a peak single-polarization absorptivity of just over 70% near for to 0.8 Ω, bracketing the range expected for a 40-nm Al film. Perido et al.63 presented dark measurements of the KID array shown in Fig. 8. The frequency noise was dominated by two-level system noise, with typical values at 10 Hz of at to 200 mK, consistent with the expected scaling. Perido et al.63 also presented a resonant dual-polarization design for which simulations yield absorption of 70% in each polarization.
The pixel-to-pixel spacing is , dominated in area by the resonators’ interdigitated capacitors. Efficiently coupling the inductive absorbers to free-space radiation will be achieved with microlens arrays, and large feedhorn arrays are currently impractical for short mid-infrared wavelengths because manufacturing tolerances are too small. For minimal absorption by the substrates, the arrays should be fabricated with silicon for wavelengths and germanium for 10 to because of the silicon absorption features in the 14 to range. Microlens arrays could either be fabricated separately and bonded to back-illuminated KID arrays64 or etched into the wafer frontside before or after KID fabrication. There are examples of silicon microlens fabrication in the literature that could achieve the micron-level tolerances required for operation down to ,65–71 although those involving micromachining and laser etching would likely be too slow for arrays of detectors. A promising approach is the deposition of Fresnel “zone plate” microlens arrays on the wafer frontsides, which have already been demonstrated at for antenna coupling.72 This approach has the virtue of fabrication simplicity but it rejects 50% of radiation. If noise-equivalent powers (NEPs) must be improved by a factor of , three-dimensional (3-D) hemispherical lenses or Fresnel lenses could be etched on the wafer frontsides instead of depositing Fresnel zones.
For all GEP focal planes, KIDs are organized into groups of detectors that are spread across a 0.6 to 1.6 GHz band and read out using electronics illustrated in Fig. 9. The choice of a 1.1-GHz center frequency resulted from a trade study in which smaller pixels were favored by the optical design but larger pixels reduced the readout frequency and bandwidth, and thus power dissipation. The readout electronics generate analog waveforms using RF-DACs that are transmitted to the cold focal plane, exciting all 1500 resonators. The 1-GHz bandwidth return signal from all 1500 KIDs is digitized with RF-ADCs and digitally channelized with sufficient resolution to separate the individual KID frequencies. The KID readout scheme, initially demonstrated in 2006,73 has now been implemented in various forms for ground-based and balloon-borne instruments.50,74–82
The GEP-I and GEP-S modules will share the same set of readout electronics (Fig. 9): they have separate, nonsimultaneous observing modes. Microwave switches on the readout lines enable selection of GEP-S or GEP-I for readout. There are 24 parallel readout channels available from three RF readout boards with eight channels each. The eight-channel RF readout cards share a payload electronics chassis with a chopping mirror driver card and two clock, processor, memory, and power cards for dual-string redundancy. The total estimated power consumption for the readout electronics is 484 W.
KID Development and Readout Outlook
For a launch date as early as January 1, 2029, technologies required for GEP must be at or above TRL 6 in 2025. Mid-infrared KIDs are currently at TRL 3, with the proof-of-concepts having been demonstrated. To reach TRL 6, they must be validated at a component level in a relevant environment and in a subsystem operational environment. In this case, the latter indicates integration with optics and readout electronics. Here, we review the current state-of-the-art and itemize the technological improvements that must be made at a component level. This will have to followed by demonstrations of sensitivity, yield, and sufficiently low cosmic ray susceptibility for arrays while read out with flight-like electronics. Of the existing demonstrations shown in Fig. 7, the SPACEKIDs results45 are the closest to meeting the GEP requirements, which are listed in Table 2. By comparing the capability gap between GEP KID focal plane requirements and SPACEKIDs performance, specific advances required to achieve TRL 6 can be identified.
Detector requirements for GEP versus achieved for SPACEKIDs.45 Notes: Minimum tile sizes for GEP-I/GEP-S shown; actual arrays could be multiples thereof. Tiles with 12×120=1440 pixel format are envisioned for GEP-I bands 1 to 18. GEP-S bands 1 and 2 assume arrays with 112×70 format, which could consist of tiles with 28×35=980 pixels. The dynamic range is specified for a 1-Jy calibration source, e.g., an asteroid.45,83 Techniques to mitigate electrical and optical crosstalk, and cosmic ray susceptibility, have been demonstrated.45,84,85
To date, all KID instrument readout systems use field-programmable gate arrays (FPGAs)88 for the digital channelization. This channelization is conceptually equivalent to an FFT and sometimes is implemented as such.74,80 First-generation KID readouts50,75–77 typically used Virtex-589 hardware to process a 500-MHz bandwidth using about 50 W, across which up to 4000 KIDs could be multiplexed.50 Similarly, the CORE mission study90 concluded that a 1-GHz bandwidth could be read with 50 W using existing space-grade, TRL-9 Virtex-5 FPGAs.91 Newer-generation FPGAs could have lower power consumption,92 (Fig. 11), thus the laboratory figure of 50 W per 1 GHz channel is conservative, and there is strong interest in the use of advanced FPGAs in space.98–103 Radiation test results on late-generation FPGAs have been positive.91,104–106 As a result, Xilinx’s 20-nm KU060 FPGA (used in the SMURF readout electronics developed at SLAC)82 will be available as a space-rated product by late 2020.91 This bodes well for even more advanced options, such as the new 16-nm FinFET Xilinx RF-system on chip (RFSoC), which integrates eight ADCs, eight DACs, and considerable FPGA logic into a single chip; the entire GEP readout could potentially be reduced to one or two such chips. A board with this chip has been released;107 initial power dissipation estimates are well below the GEP conservative assumption of 25 W per 1 GHz channel. Alternatively, mixed-signal application-specific integrated circuits that integrate ADCs and signal processing have been developed for similar applications.108
GEP Surveys, Simulations, and Science
GEP is designed to do large, unbiased surveys. Its overarching goal is to provide a legacy dataset with broad utility for studying the evolution of star-formation, interstellar medium, and black hole growth in galaxies. Parameters of GEP’s surveys in the context of other surveys are shown in Figs. 12 and 13 with complementary surveys and the astrophysics that they probe.
The GEP-I survey program is optimized to sample a comprehensive range of redshifts and galaxy luminosities (Table 3). A combination of four depths and areas will sample low redshift and rare, luminous galaxies, faint, high-redshift galaxies, and intermediate redshift and luminosity galaxies. Large numbers of galaxies will be detected in each tier, preventing sample variance and small number statistics from hampering conclusions on galaxy evolution with redshift (Fig. 14). It is anticipated that the GEP-I surveys (except the all-sky survey) will be centered on and divided between the north and south ecliptic poles to minimize the photon backgrounds from primarily zodiacal dust and secondarily Galactic dust. This will provide sky coverage overlap with Euclid and Nancy Grace Roman Space Telescope (NGRST) surveys, which will provide near-infrared photometry, morphologies, and stellar masses for galaxies detected by GEP.
GEP galaxy surveys and yields obtainable with background-limited sensitivities.
Three types of spectroscopic measurements with GEP-S will complement GEP-I photometric surveys: (1) individual observations of specific galaxies identified in the GEP-I surveys to provide precise redshifts and to validate PAH redshifts and to obtain measurements of the far-infrared emission lines. (2) Unbiased spectroscopic surveys obtained by rastering GEP-S on the sky. (3) Spectral maps of nearby galaxies. GEP-S will perform a deep spectroscopic survey over 1.5 square degrees and a wide spectroscopic survey over 100 square degrees (Table 3). The spectral survey datasets will detect galaxies the far-infrared fine-structure transitions (and the continuum, when binned). The wide survey will be used to stack GEP spectra on the NGRST and/or Euclid grism sources to provide high signal-to-noise ratio average galaxy spectra in bins (Table 4). The GEP-S spectral mapping speed is shown in Fig. 15, where it is seen that GEP would be a substantial improvement of the state of the art with times to survey sky regions in the far-infrared six orders of magnitude faster than previous observatories.
GEP stacking of IR galaxy datasets obtainable with background-limited sensitivities.
Theoretical Framework and Simulations
To quantify the science yield of the GEP-I survey program, a mock survey was constructed using a combination of the Millennium N-body simulation127 to provide the distribution of large-scale structure and the Galacticus semianalytic model128 to populate that simulation with galaxies based on physical models. For each dark matter halo in the simulation volume, the star-formation rate and black hole accretion rate of the galaxy were computed. Bolometric infrared luminosities from star-formation were estimated as () ergs ,129 and those due to AGN activity as (where is the radiative efficiency computed by Galacticus from the black hole spin and accretion rate, and is the black hole accretion rate). Infrared spectra with the corresponding AGN fraction were then assigned to galaxies using the models of Dale et al.28 and normalized to give the computed total infrared bolometric luminosity. In the Dale et al.28 models, the distribution of dust heating intensities (and, therefore, temperatures) is parameterized in terms of , the exponent of the power-law distribution of heating intensities. Values of for model galaxies were drawn from a Gaussian distribution with mean, , and dispersion, , with no redshift dependence. A range of different possible values of and was explored, with the above values selected as those which result in the closest match between our simulated and observed numbers counts of galaxies over the 24- and wavelength range.130,131 These spectra were then used to compute broadband luminosities for model galaxies in each GEP-I band. Finally, a light cone from this synthetic catalog was extracted, corresponding to a 4-square degree area, from to 3, and observed fluxes were determined in all GEP-I bands for each galaxy.
The galaxy number counts estimated in this way are lower limits because these (and other) galaxy evolution models generically underpredict the number of very high luminosity galaxies132 (). Although we attempted to tune the Galacticus model number counts to match observations from Spitzer and Herschel by judiciously choosing from the Dale et al.28 spectral library, at high flux densities (i.e., at ) the model counts are almost an order of magnitude too low compared with observations.131 The agreement is much better at 24 and , coming close to matching observations.130 Some of the bright observed galaxies that are not accounted for in the models likely result from gravitational lensing, although lensing is unlikely to account for the majority of the discrepancy. We adopt our model predictions for the detection rates with the understanding that they are likely conservative.
Noise was added to the simulated maps corresponding to the expected map depths after they were convolved with GEP beam sizes. A small region of the maps stepping through the bands is animated in Fig. 16.
Extragalactic Source Confusion
Extragalactic source confusion arises when images of galaxies—the point spread function for unresolved galaxies—overlap. It happens at far-infrared wavelengths where diffraction-limited beams can be several arcseconds or larger. Confusion “noise” is the uncertainty in the extraction of a given source’s flux due to the presence of a background of other sources which cannot accurately be subtracted from the signal. Analytically, the confusion noise () is just the standard deviation of the flux from beam to beam, and for extraction of sources from a map, a “confusion limit” arises at this value. This limiting flux below which extractions are unreliable typically corresponds to that for which density of sources is 1 per several () beams, depending on the source count relation. This basic paradigm holds true when referring to extracting sources from a map without prior information, but there are well-developed methods to employ prior higher-resolution datasets to extract source parameters in a larger far-IR beam. An excellent example is use of Spitzer IRAC and MIPS positions to extract (statistically) fluxes in the 70 and maps.133 In this section, we assess the important issue of source confusion with the GEP-I datasets, but first note that spectral surveys with GEP-S should not suffer appreciably from source confusion.
For spectral mapping with GEP-S, source confusion should not limit the extraction of line intensities. The key point is that the bright mid- and far-IR spectral lines are sparse and well-separated spectrally, and they form an unambiguous template for redshift identification. Thus, each line can be conclusively identified and measured, even with multiple line emitting sources in the beam. For a 3-D spatial–spectral survey, the analog of source counts for two-dimensional is line counts, the number of line emitters above a given flux level in a typical spatial–spectral bin of the survey. This issue was studied for the OST, which also features a wideband, moderate-resolving-power spectrograph; Ref. 134 is a good reference. As an example they predict that at (which is approximately GEP-S long-wavelength range limit), the density of lines with fluxes above (corresponding to a 25-h integration with GEP-S) per beam-bin is 1/150 for Origins. Correcting this for the larger GEP beam and slightly smaller resolving power results in an estimate of one source per 12 beam-bins for GEP-S at this wavelength. This is approaching but not exceeding a practical confusion limit. Shorter wavelengths of course are more forgiving with the smaller beam.
We estimate the extragalactic confusion noise to be expected for GEP-I observations by considering confusion noise measurements from previous observations. The confusion noise flux density was measured—or upper limits were placed—with Spitzer 70 and 135,136 and Herschel 70, 100, and .131,137,138 We start by deriving an empirical relationship for confusion noise as a function of telescope aperture diameter to assess the dependence of confusion on telescope size and to check for consistency between Spitzer and Herschel observations, which had different wavebands and aperture sizes (0.85 and 3.5 m, respectively). Using the number counts models of Béthermin et al.,139 the flux density corresponding to 15 beams per source (where galaxies are not be substantially blended) scales approximately as the inverse of aperture diameter squared (, slightly steeper for apertures smaller than 1.7 m, Fig. 17). Confusion noise measurements in the literature cited above scaled by this relationship are consistent, which validates the scaling relation for interpolation to the GEP 2.0 m aperture for estimation of confusion noise. Critically, the scaling also shows that aperture diameters below 2.0 m will be increasingly susceptible to confusion. Since the confusion flux density drops slowly with telescope diameter beyond 2.0 m, 2.0 m represents something of a “sweet spot” in the trade between confusion and cost.
A comparison of the expected GEP-I survey map depths for a 2.0-m aperture with the confusion measurements cited in the literature already mentioned yields the following:
In short, GEP-I will likely have significant confusion noise at and longer, but not at shorter wavelengths. GEP-I will have to integrate deeper than the confusion noise in the two deepest surveys for PAH redshifts of galaxies with the wavebands at and below, which will not be limited by confusion noise. In addition, monochromatic fluctuation “probability of deflection” P(D) analyses show that it is possible to constrain galaxy populations meaningfully with observations deeper than the confusion noise level.140,141 A polychromatic P(D) analysis with GEP-I observations covering the redshifted PAH features would be extremely powerful; it would yield precise galaxy number counts and redshifts statistically for the ensemble, and therefore luminosity functions, and tightly constrain galaxy evolution models. Furthermore, using cross-identification with counterparts at shorter wavelengths, galaxy properties can be measured even when there is source confusion. For example, Labbé et al.142 showed that contamination by confusion can be reduced a factor of six with short-wavelength prior-based photometry, and positional and flux priors from ancillary data could enable improved spectral energy distribution extraction. In addition, because the detector count in wavebands longer than (GEP-I band 19) is small (4.2) compared with shorter wavelengths, the cost of retaining them is merited for measurements of far-infrared luminosities, for nearby galaxy science that do not require observations as deep, and for P(D)-type fluctuation analyses.
For wavelengths , corresponding to GEP-I bands 17 to 23, photometry of partly blended sources can still be extracted reliably. The Next Generation (X)Cross Identification (XID+) code was developed to estimate flux densities accurately from confusion-limited Herschel imaging;143 we apply it here to GEP-I. The performance of XID+ was quantified by measuring the differences between fitted and true galaxy flux densities as a function of the distance to the nearest galaxy neighbor using the simulations described in Sec. 5.2. Positional priors were assumed from the mid-infrared wavebands and the fitted galaxy far-infrared flux densities were compared to their true values. Figure 18 shows that down to galaxy separations of the beam size—below the classical confusion “limit”—galaxy flux densities can be deblended with small fractional errors and little or no bias in most cases. Using the deblended sources, the map flux densities were reproduced with small residuals, although with poorer performance at the map edges because of missing galaxies outside the field of view. This will not be a concern for the majority of galaxies detected in GEP-I surveys, for which the minimum size is 3 square degrees.
GEP-I Redshift Measurements
To measure the cosmic star formation history of the Universe, GEP requires redshifts with precision to , a requirement that our simulations show is achieved with margin. In general, the redshift uncertainty is set by the width of the GEP-I bands relative to the widths of the PAH features.
Expected GEP-I redshift precision was estimated by adding noise to the model spectra according to each of the map depths given by the photon backgrounds. The mock spectra were binned into the GEP-I bands and were calculated by comparison to the spectral model and two other different spectral models: since the spectra of galaxies will not be known a priori, the comparison models were used to ascertain the uncertainties incurred by having a spectrum different from the model. Only the first 11 GEP-I bands (10 to ) were used for the redshift measurements because the steeply rising mid-infrared dust spectra influence the redshifts and the dust temperatures will not be known a priori. A sample GEP-I redshift chi-squared distribution is shown in Fig. 19. The nominal model had strong PAH emission features and the alternate comparison models had: (1) strong PAH emission features but cooler dust (hence a more slowly rising spectrum with wavelength) and (2) hot dust that substantially overwhelmed the PAH features above . The results were as follows:
A comprehensive set of simulations is required to quantify redshift measurements as a function of GEP-I spectral resolving power, AGN fraction and other sources of weak PAH emission (e.g., low metallicity), galaxy-to-galaxy variability, and instrument parameters, such as nonideal bandpass shapes and noise. These simulations likely represent best-case scenarios. These effects will be explored as a subsequent effort.
GEP Redshift Measurements for Galaxies with Strong AGN
GEP has two means for measuring redshifts. One, GEP’s spectrometer GEP-S will provide redshifts for AGN-dominated galaxies via (low- and) high-ionization fine-structure lines. Two, even in the instances of weak PAH-to-continuum emission in strongly AGN-dominated galaxies, the rest-frame silicate absorption feature will still enable redshift measurements. The harsh environments immediately around AGNs can in theory destroy PAHs,144,145 but observationally, even galaxies with AGN often display these features, although they are weaker. Detailed Spitzer IRS spectra of hundreds of galaxies with AGN activity found that some PAH features are seen in most cases,146 showing that PAHs can survive in some of these environments. Li147 summarized the PAH observations in AGN by noting that the 6.2, 7.7, and PAH features, although suppressed with respect to the band (perhaps because smaller PAHs are more susceptible to the hard radiation field of an AGN) nevertheless can be present. Figure 5 in Ref. 147 of the relative PAH line strengths shows the influence of radiation hardness, as well as the limits of current modeling in Seyfert galaxies, many of which show strong neutral PAH emission.
From a theoretical perspective, coauthors of this paper have examined the role of AGN activity in a galaxy mid-IR SEDs in another paper148 using a Gadget-2 simulation of an idealized (noncosmological) major merger of two identical disk galaxies.149 The galaxies had initial halo and baryonic masses of and , and a central black hole mass of . The star formation and feedback were modeled as described in Refs. 150 and 149. The radiative transfer code SUNRISE151 was used to compute SEDs for seven viewing angles every 10 Myr throughout the simulation run, and by varying the efficiency of the AGN it was effectively turned off or raised to 100% (for Eddington-limited Bondi-Hoyle accretion). The ratio of the two resultant SEDs reveals the influence of the pure AGN (see also Fig. 9 in Ref. 152). The results presented in Fig. 1 of Ref. 148 show that the silicate absorption feature is present and will enable redshift measurements even when galaxies’ luminosities are strongly dominated by AGN.
Other physical effects will also influence the strengths of these mid-IR PAHs features, for example, low metallicity in high-, young galaxies, or destructive supernovae activity.147 Although weak PAH features make photometric redshift measurements more challenging, it is precisely the complexity and variability of these features that ultimately enable diagnosis of the physical conditions present.
Galaxy Luminosity Functions, Star Formation History, Separating AGN, and Clustering
GEP’s science objectives require luminosity functions over a range of redshifts to observe the changes in galaxy formation and the build-up of stellar mass over cosmic time. Faint-end (below ) mid- and far-infrared luminosity function slopes have not been measured above , and there is disagreement about the faint-end slopes even at . GEP-I’s surveys will detect hundreds of millions of galaxies and measure redshifts for millions of them for derivation of luminosity functions. With GEP-I surveys, faint-end slopes below for , below ( at , and below at , will be measured. Figure 19 shows current observational determinations of the infrared bolometric luminosity function and compares these to a sample luminosity function derived from our mock GEP catalogs.
One of GEP’s primary goals is to quantify the star formation history of the Universe by utilizing very large statistical samples of galaxies that overcome Poisson statistics and sample variance from large-scale structure. How well will it do? Figure 14 showed that Poisson statistics and sample variance will be overcome and Fig. 20 demonstrated that luminosity functions will be well quantified even below at . How do these translate into measurement of the history of cosmic star formation rate density? Figure 21 shows that the star formation history will be quantified with uncertainties an order of magnitude below the state of the art to at least . Uncertainties will likely improve substantially beyond also; however, that was not included in the simulations (Sec. 5.2). Because it maps large areas, GEP will also be effective at detecting gravitationally lensed galaxies. This will enable it to serendipitously detect sub- galaxies at with a tail out to higher redshifts, perhaps lensed galaxies in the all-sky survey and 50 in the deeper 300-square degree survey.
Another GEP goal is to identify embedded AGN and quantify their contributions to the infrared luminosities of galaxies. The mid-infrared part of the spectrum provides excellent leverage for this because the dust is hotter—thus the spectra are “bluer”—and the PAHs are more subsumed by the bright dust continuum. Figure 22 shows that a simple principal component analysis with only three components is sufficient to identify AGN and separate their contributions to infrared luminosities from star formation with GEP-I data. The high-ionization lines observed with GEP-S will indicate the presence of AGN and enable the larger catalog of GEP-I galaxies to be calibrated against the fine-structure line luminosities. These observations can be used to assess AGN influence on galactic interstellar media and the coevolution of supermassive black holes and star formation.
Within a cosmology, the clustering of galaxies reveals the masses of dark matter haloes that they occupy. Galaxy mass correlates strongly with halo mass; however, it is possible that mid- and far-infrared luminosities of galaxies correlate more weakly with halo mass because even low-mass galaxies can be temporarily infrared-bright from star formation bursts. Clustering of GEP-I catalogs will strongly test galaxy formation and evolution theories that predict halo occupations as a function of redshift (Fig. 23).
Milky Way and Nearby Galaxies
GEP will map the Milky Way and the Magellanic Clouds with GEP-I during the all-sky survey. The sensitivity required of the telescope and instruments for the deep extragalactic observations will enable high signal-to-noise ratio mid- and far-infrared spectral energy distributions of star-formation regions. The 3.43″ mid-infrared resolution will resolve gradients in dust temperatures and PAH excitation (Fig. 24), thereby probing the radiation fields and chemistry, whereas the far-infrared spectral energy distributions will yield luminosities of embedded stars, dust temperatures, and gas masses.
Interstellar Physical Conditions
Like the multiband imager, GEP-S is designed for rapid surveys but targeting the gas-phase spectral features in the rest-frame mid- and far-infrared. The spectra from GEP-S will chart the cosmic history of the interaction of galaxies’ energy sources (stars and supermassive black holes) with their gas reservoirs. Transitions from low-ionization states, such as [C II], [N II], [O III], and [NeII], can be used to infer the masses of interstellar neutral and ionized gas components, the hardness of the ultraviolet radiation fields and its implications for the stellar initial mass function, and the density of HII regions from which the pressure in the ISM can be inferred. Intermediate ionization transitions of nitrogen and oxygen, [NIII] and [OIII], indicated the degree of stellar processing, a proxy for metallicity154 in galaxies irrespective of the dust and gas temperature uncertainties. Finally, spectral tracers of highly ionized gas, such as [S IV] at , [Ne V] at 14.3 and , and [O IV] at , with ionization potentials of 35 to 97 eV, when combined with bright lines with much lower ionization potential (e.g., [Ne II] at ), indicate the relative amount of heating from young stars and AGN.
To utilize this basic spectroscopic toolkit, there are three types of observations envisioned for GEP-S: (1) “blind” field-filling spectroscopic surveys obtained by rastering GEP-S on the sky, described further below; (2) individual pointed observations of specific galaxies identified in the photometric surveys or other facilities, for both precise redshifts to validate the photometric technique and fluxes in the full suite of far-infrared spectral features; and (3) spectral maps of nearby galaxies to study feedback and interstellar energy balance in detail.
Survey spectroscopy with GEP-S is a particular strength of the GEP long-slit spectrometer architecture, and the study envisions two spectral surveys: a deep survey over 1.5 square degrees and a wide spectroscopic survey over 100 square degrees. The resulting datasets will be used in at least three ways. First, they will detect galaxies directly in the far-infrared fine-structure transitions (and the continuum when binned). Second, the wide survey will be used to stack on the NGRST and/or Euclid grism sources to provide high signal-to-noise ratio average galaxy spectra in bins. Both surveys will be excellent datasets to measure the total cosmic luminosity density in the various far-infrared lines and ratios among integrated line intensities.
The synergy with the coming near-infrared grism surveys will be an especially powerful tool. Comparing the anticipated Euclid and NGRST grism depths155 with GEP depth in the 100-square-degree survey suggests that stacks will access populations down to , respectively, in a total sample of 60 (22) million galaxies, respectively, for Euclid and NGRST. The sensitivity is sufficient to form hundreds of independent bins in quantities, such as star formation rate, stellar mass, and metallicity. The resulting spectra will measure not only the bright fine-structure lines but also the “second-tier” transitions, such as [Ne V] (at 14 and rest frame), which provide dust-immune AGN indicators, and the quadrupole pure rotational lines () which indicate warm shocked gas.156 In considering the Euclid and Roman datasets, those sources that are very highly obscured may not yield clear redshifts with the GRISM—they of course cannot be included in the stacking analysis. These sources will, however, make excellent candidates for pointed spectroscopy with GEP-S, which will access the bright and dust immune far-IR fine-structure transitions of [OI] and [OIII]. At a redshift of , the 1-h GEP-S sensitivity of is a good match to a source whose less-obscured counterpart is included in the stacking analysis. While a sample of such sources will not have the full statistical power of the stacking analyses, it will nevertheless provide useful and otherwise inaccessible constraints on the properties of these obscured galaxies.
These spectral surveys will also be used for line intensity mapping, a powerful emerging technique in which the clustering of line-emitting galaxies is detected as fluctuations in a 3-D spatial–spectral dataset with redshift encoded as wavelength. The key virtue is that technique is sensitive to all forms of emission, not just the individually detected galaxies, so it offers the potential for absolute aggregate measurements irrespective of a single-source detection limit. GEP-S intensity mapping of mid- and far-infrared fine structure lines—and possibly PAHs—will yield signals in autocorrelation, but it will also be powerful when cross-correlated with rest-frame optical and ultraviolet surveys and ground-based millimeter-wave surveys that separately probe stars and the atomic and molecular interstellar medium.
Probing Feedback from Extraplanar Gas
A hallmark feature of the cryogenic telescope and sensitive GEP-S spectrometer is excellent surface brightness sensitivity. This will provide a breakthrough capability for assessing feedback effects in local galaxies. Feedback in the form of stellar and AGN winds and supernova explosions are believed to eject gas and dust from star-forming galaxy disks, creating a reservoir of low-column-density gas in the outskirts of galaxies. This material exists below the threshold for star formation (hydrogen column density below 1 to ) and has been difficult to measure. GEP-S will be capable of detecting [CII] and [NII] in this material, even when the local density is well below the critical density for excitation. Sensitive maps of the diffuse gas can be obtained with GEP-S around low-redshift star-forming galaxies and AGN. These maps will provide a powerful test of the theories and models of stellar feedback which have evolved to have great detail but need testing by observations.
Future Science Exploration
GEP’s science reach will be far broader than outlined here. For example, monochromatic probability-of-density—P(D)—fluctuation analyses have shown to be effective for constraining galaxy number counts into the confusion noise.140,141 Multicolor P(D) analyses could go much further: distinguishing between galaxy number count models and thereby constraining luminosity functions as a function of redshift to substantially greater depth than for analyses restricted to galaxies whose brightness are above individual detection thresholds. In nearby galaxies, mid- and far-infrared spectral energy distributions are needed for panchromatic spectral energy distribution fitting to constrain stellar populations, nebular conditions, star formation rates, and identify embedded AGN on a galaxy-by-galaxy basis. In the Milky Way, mapping of PAH emission in various environments, such as the vicinities of hot, young stars, and photodissociation regions, will probe the radiation fields and PAH production, excitation, and destruction. Lastly, with the possibility of linear-variable filters with , the optimum resolution for PAH redshifts should be revisited. Then, GEP-I redshift precision should be quantified as a function of galaxy luminosity, redshift, AGN fraction (with strong AGN mid-infrared continuum leading to weak PAH lines), and instrumental parameters, such as spectral resolving power .
Finally, it must be noted that the GEP design reference mission is designed for dedicated surveys. The surveys should be designed with community input and consideration of the rich landscape of imminent multiwavelength galaxy surveys. And, while the design reference mission does not support an open-time phase to manage costs, with no expendable cryogens the expected lifetime of the mission exceeds the planned survey durations. An extended mission with guest-observer opportunities would yield a large volume of science, as the Spitzer extended mission has shown.
The GEP concept was developed to make precise measurements of the star-formation rates, nuclear accretion rates, and interstellar conditions of galaxies over cosmic time and over the full range of cosmic environments. It was designed as a survey mission to obtain large, well-defined samples of galaxies that will not be limited by counting statistics nor by cosmic sample variance. Its science goals will be realized with 2.0 m, cryogenic telescope () and two instrument modules, each with arrays of KIDs. It will utilize established cryogenic, telescope, and bus technologies. GEP will improve the measurement uncertainties of the cosmic star formation rate density by more than an order of magnitude out to a redshift of at least , the current extent of our predictive simulations. Obscured AGN in galactic centers will be identified on a galaxy-by-galaxy basis, yielding new luminosity functions. Critical measurements of metallicities and interstellar conditions will be obtained for galaxies, complementary to probes at other wavelengths. GEP’s ability to provide these powerful new observational capabilities derives from emerging mid- and far-infrared detector technologies that require continued development and nurturing to realize not only GEP but also a future flagship opportunity, such as OST.
This work was supported, in part, by a NASA Astrophysics Probe Concept Study grant to Jason Glenn (NASA Award Number NNX17AJ89G). This research was funded in part by the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the NASA. Jason Glenn thanks JPL and Ball Aerospace for their strong support of the study effort. The GEP team would like to thank many contributors whose important work helped the successful completion of this study but whose roles traditionally do not qualify for inclusion in authorship lists of scientific publications: engineering aides, documentarians, administrators, administrative aides, and reviewers. The authors have no relevant financial interests in the manuscript and no other potential conflicts of interest to disclose.
https://doi.org/10.1146/annurev-astro-081811-125615 ARAAAJ 0066-4146 Google Scholar
https://doi.org/10.1146/annurev-astro-082708-101734 ARAAAJ 0066-4146 Google Scholar
https://doi.org/10.1146/annurev-astro-081811-125521 Google Scholar
https://smd-prod.s3.amazonaws.com/science-red/s3fs-public/atoms/files/GEP_Study_Rpt.pdf Google Scholar
http://surveygizmoresponseuploads.s3.amazonaws.com/fileuploads/623127/5043187/183-0ec368750b918c40030f7019cf8bbaea_GlennJason.pdf Google Scholar
https://www.xilinx.com/support/documentation/selection-guides/ultrascale-fpga-product-selection-guide.pdf Google Scholar
https://indico.esa.int/event/232/contributions/2147/attachments/1800/2100/SEFUW_2018_FPGA_and_SOC_Airbus_Defence_Space.pdf Google Scholar
https://doi.org/10.1088/0067-0049/197/2/36 APJSA2 0067-0049 Google Scholar
https://doi.org/10.3847/1538-4357/aacb79 ASJOAB 0004-637X Google Scholar
https://roman.gsfc.nasa.gov/science/Roman_Reference_Information.html Google Scholar
https://irsa.ipac.caltech.edu/data/SPITZER/docs/mips/mipsinstrumenthandbook Google Scholar
Jason Glenn received his BS degree in physics from the University of New Mexico in 1991 and his PhD in astronomy from the University of Arizona in 1997. Since 2020, he has been a research astrophysicist, stellar, galactic, and extragalactic in the Observational Cosmology Laboratory at the NASA Goddard Space Flight Research Center. Prior to that, he was a professor of astrophysics at the University of Colorado Boulder for 20 years. He is the author of more than 200 journal papers. His current research interests include galaxy evolution, the cool interstellar medium in galaxies, and infrared through millimeter-wave instrument and detector development.
Charles M. Bradford received his doctorate from Cornell in 2001. He held a Millikan postdoctoral fellowship at Caltech from 2001 to 2003 and has been on the science staff at Jet Propulsion Laboratory (JPL) since that time. He enjoys developing and fielding new submillimeter- and millimeter-wave instrumentation, and using spectroscopy to study interstellar medium conditions in galaxies. His current projects include a balloon-borne far-IR spectrometer, an on-chip mm-wave spectrometer, and ultra-sensitive far-IR detectors for cryogenic space missions.
Rashied Amini, PhD, is a system engineer at JPL, California Institute of Technology, working in mission formulation and in research of autonomous technologies. Recently, he was the HabEx and Galaxy Evolution Probe Study Lead, submitted to the 2020 Astrophysics Decadal Survey. As a result of his formulation work, he is interested in supporting the maturation of technologies critical to science exploration. He received his PhD in physics from Washington University in St. Louis.
Lee Armus is a senior staff scientist at IPAC with over 25 years of experience in infrared spectroscopy and imaging. He is currently the lead scientist at the Roman Space Telescope Science Support Center, and formerly the lead for the IRS Instrument Support Team at the Spitzer Science Center. His research is focused on several topics central to the science of the Galaxy Evolution Probe, including galaxy mergers, galactic outflows, and luminous infrared galaxies.
Andrew J. Benson is a staff scientist at the Observatories of the Carnegie Institution for Science. His research focuses on understanding the nature of dark matter and the process of galaxy formation, with a particular emphasis on formulating a coherent picture of the many different aspects of these problems. He has developed a model of dark matter and galaxy formation physics, Galacticus, which blends both analytic understanding and numerical techniques.
Jeremy Darling is a professor at the University of Colorado. He studies osmology, galaxy evolution, massive black holes, and astrophysical approaches to fundamental physics.
Jeanette L. Domber is a program manager for Ball Aerospace Civil Space. She received her PhD in aerospace engineering from the University of Colorado at Boulder.
Sarah J. Lipscy is the deputy director of Civil Space Business Development at Ball Aerospace. She holds a PhD in astronomy and physics from UCLA.
Raphael Shirley is a postdoctoral researcher at the University of Southampton. His research is focused on galaxy evolution and the role of active galactic nuclei in star formation over cosmic time. He has worked on Herschel imaging, multi-wavelength astronomy, and supernova observations. He is now working on image processing pipelines for the Vera C. Rubin Observatory and the VISTA telescope. He has previously held positions at Sussex, Cambridge, and the IAC in Tenerife.
Howard A. Smith is a senior astrophysicist at the Center for Astrophysics, Harvard & Smithsonian (CfA) in Cambridge, with over 350 published scientific articles, and is a member of the Harvard Astronomy Department. His research field is the origins of stars in the Milky Way and other galaxies, with specialization in techniques of infrared spectroscopy and instrument development.
Jonas Zmuidzinas (BS Caltech 1981, PhD Berkeley 1987) joined the Caltech physics faculty in 1989 and currently serves as the Merle Kingsley Professor of Physics and Director of the Caltech Optical Observatories. He has also concurrently held positions at JPL as a senior research scientist (2006 to 2016), director of the JPL Microdevices Laboratory (2007 to 2011), and JPL chief technologist (2011 to 2016). His research focuses on superconducting detectors and devices and their application to astronomy.