The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We investigate the nature of the relation among stellar mass, star formation rate, and gas-phase metallicity (the M∗-SFR-Z relation) at high redshifts using a sample of 260 star-forming galaxies at z∼2.3 from the MOSDEF survey. We present an analysis of the high-redshift M∗-SFR-Z relation based on several emission-line ratios for the first time. We show that a M∗-SFR-Z relation clearly exists at z∼2.3. The strength of this relation is similar to predictions from cosmological hydrodynamical simulations. By performing a direct comparison of stacks of z∼0 and z∼2.3 galaxies, we find that z∼2.3 galaxies have ∼0.1 dex lower metallicity at fixed M∗ and SFR. In the context of chemical evolution models, this evolution of the M∗-SFR-Z relation suggests an increase with redshift of the mass-loading factor at fixed M∗, as well as a decrease in the metallicity of infalling gas that is likely due to a lower importance of gas recycling relative to accretion from the intergalactic medium at high redshifts. Performing this analysis simultaneously with multiple metallicity-sensitive line ratios allows us to rule out the evolution in physical conditions (e.g., N/O ratio, ionization parameter, and hardness of the ionizing spectrum) at fixed metallicity as the source of the observed trends with redshift and with SFR at fixed M∗ at z∼2.3. While this study highlights the promise of performing high-order tests of chemical evolution models at high redshifts, detailed quantitative comparisons ultimately await a full understanding of the evolution of metallicity calibrations with redshift.

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R = 3000-3650) rest-frame optical spectra (∼3700-7000 ) for ∼1500 galaxies at in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals:, down to fixed (F160W) magnitudes of 24.0, 24.5, and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [O ii], Hβ, [O iii], H, [N ii], and [S ii]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-ii H and K, Mgb, 4000 break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ∼80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass () and star formation rate. The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

Astronomical software is now a fact of daily life for all hands-on members of our community. Purpose-built software for data reduction and modeling tasks becomes ever more critical as we handle larger amounts of data and simulations. However, the writing of astronomical software is unglamorous, the rewards are not always clear, and there are structural disincentives to releasing software publicly and to embedding it in the scientific literature, which can lead to significant duplication of effort and an incomplete scientific record. We identify some of these structural disincentives and suggest a variety of approaches to address them, with the goals of raising the quality of astronomical software, improving the lot of scientist-authors, and providing benefits to the entire community, analogous to the benefits provided by open access to large survey and simulation datasets. Our aim is to open a conversation on how to move forward. We advocate that: (1) the astronomical community consider software as an integral and fundable part of facility construction and science programs; (2) that software release be considered as integral to the open and reproducible scientific process as are publication and data release; (3) that we adopt technologies and repositories for releasing and collaboration on software that have worked for open-source software; (4) that we seek structural incentives to make the release of software and related publications easier for scientist-authors; (5) that we consider new ways of funding the development of grass-roots software; (6) and that we rethink our values to acknowledge that astronomical software development is not just a technical endeavor, but a fundamental part of our scientific practice.

We present the stellar mass (M∗)gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts (0.5 ≤ z ≤ 0.7) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given M∗ of this magnitude-limited (R ≳ 24 AB) sample are representative of normal starforming galaxies. For masses below 109 Mo, our sample of 237 galaxies is ∼10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at 108 Mo < M∗ < 1011 Mo: 12 + log (O/H) = (5.83 ± 0.19) +(0.30 ± 0.02) log (M∗/Mo). At 109 Mo < M∗ < 1010.5 Mo, our MZR shows agreement with others measured at similar redshifts in the literature. Our power-law slope is, however, shallower than the extrapolation of the MZRs of others to masses below 109 Mo. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the fundamental metallicity relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as Mo decreases. Our result implies that either the scatter of the baryonic accretion rate (δṀ? ) or the scatter of the M∗-Mhalo relation (δSHMR) increases as M∗ decreases. Moreover, our measure of scatter at z = 0.7 appears consistent with that found for local galaxies. This lack of redshift evolution constrains models of galaxy evolution to have both δṀ and δSHMR remain unchanged from z = 0.7 to z = 0.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

We report the discovery of a candidate galaxy with a photo-z of z ~ 12 in the first epoch of the James Webb Space Telescope (JWST) Cosmic Evolution Early Release Science Survey. Following conservative selection criteria, we identify a source with a robust zphot = 11.8-0.2+0.3 (1σ uncertainty) with mF200W = 27.3 and ≥7σ detections in five filters. The source is not detected at λ < 1.4 μm in deep imaging from both Hubble Space Telescope (HST) and JWST and has faint ~3σ detections in JWST F150W and HST F160W, which signal a Lyα break near the red edge of both filters, implying z ~ 12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2σ lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z probability distribution function favoring z > 11. All data-quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (rh = 340 ± 14 pc). Maisie's Galaxy has log M∗/M⊙ ~ 8.5 and is highly star-forming (log sSFR~-8.2 yr-1), with a blue rest- UV color (β~-2.5) indicating little dust, though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions that smoothly decline with increasing redshift. Should follow-up spectroscopy validate this redshift, our universe was already aglow with galaxies less than 400 Myr after the Big Bang.