Many of the fundamental physical constants in Physics, as a discipline, are measured to exquisite levels of precision. The fundamental constants that define Cosmology, however, are largely determined via a handful of independent techniques that are applied to even fewer datasets. The history of the measurement of the Hubble Constant (H0), which serves to anchor the expansion history of the Universe to its current value, is an exemplar of the difficulties of cosmological measurement; indeed, as we approach the centennial of its first measurement, the quest for H0 still consumes a great number of resources. In this white paper, we demonstrate how the approaching era of Extremely Large Telescopes (ELTs) will transform the astrophysical measure of H0 from the limited and few into a fundamentally new regime where (i) multiple, independent techniques are employed with modest use of large aperture facilities and (ii) 1% or better precision is readily attainable. This quantum leap in how we approach H0 is due to the unparalleled sensitivity and spatial resolution of ELT's and the ability to use integral field observations for simultaneous spectroscopy and photometry, which together permit both familiar and new techniques to effectively by-pass the conventional 'ladder' framework to minimize total uncertainty. Three independent techniques are discussed -- (i) standard candles via a two-step distance ladder applied to metal, poor stellar populations, (ii) standard clocks via gravitational lens cosmography, and (iii) standard sirens via gravitational wave sources -- each of which can reach 1% with relatively modest investment from 30-m class facilities.

We present in this study a first analysis of the astrometric error budget of absolute astrometry relative to back- ground galaxies using adaptive optics. We use for this analysis multi-conjugated adaptive optics (MCAO) images obtained with GeMS/GSAOI at Gemini South. We find that it is possible to obtain 0.3 mas reference precision in a random field with 1 hour on source using faint background galaxies. Systematic errors are correctable below that level, such that the overall error is approximately 0.4 mas. Because the reference sources are extended, we find it necessary to correct for the dependency of the PSF centroid on the used aperture size, which would oth- erwise cause an important bias. This effect needs also to be considered for Extremely Large Telescopes (ELTs). When this effect is corrected, ELTs have the potential to measure proper motions of dwarfs galaxies around M31 with 10 km/s accuracy over a baseline of 5 years.

The standard picture of galaxy formation motivates the decomposition of the Milky Way into 3--4 stellar populations with distinct kinematic and elemental abundance distributions: the thin disk, thick disk, bulge, and stellar halo. To test this idea, we construct a Gaussian mixture model (GMM) for both simulated and observed stars in the Solar neighborhood, using measured velocities and iron abundances (i.e., an augmented Toomre diagram) as the distributions to be decomposed. We compare results for the Gaia-APOGEE DR16 crossmatch catalog of the Solar neighborhood with those from a suite of synthetic Gaia-APOGEE crossmatches constructed from FIRE-2 cosmological simulations of Milky Way-mass galaxies. We find that in both the synthetic and real data, the best-fit GMM uses five independent components, some of whose properties resemble the standard populations predicted by galaxy formation theory. Two components can be identified unambiguously as the thin disk and another as the halo. However, instead of a single counterpart to the thick disk, there are three intermediate components with different age and alpha abundance distributions (although these data are not used to construct the model). We use decompositions of the synthetic data to show that the classified components indeed correspond to stars with different origins. By analogy with the simulated data, we show that our mixture model of the real Gaia-APOGEE crossmatch distinguishes the following components: (1) a classic thin disk of young stars on circular orbits (46%), (2) thin disk stars heated by interactions with satellites (22%), (3, 4) two components representing the velocity asymmetry of the alpha-enhanced thick disk (27%), and (5) a stellar halo consistent with early, massive accretion (4%).

We present the first release of the MaNGA Stellar Library (MaStar), which is a large, well-calibrated, high-quality empirical library covering the wavelength range 3622-10354 Å at a resolving power of R ∼ 1800. The spectra were obtained using the same instrument as used by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, by piggybacking on the Sloan Digital Sky Survey (SDSS-IV)/Apache Point Observatory Galaxy Evolution Experiment 2-N (APOGEE-2N) observations. Compared to previous empirical libraries, the MaStar library will have a higher number of stars and a more comprehensive stellar-parameter coverage, especially of cool dwarfs, low-metallicity stars, and stars with different [α/Fe], achieved by a sophisticated target-selection strategy that takes advantage of stellar-parameter catalogs from the literature. This empirical library will provide a new basis for stellar-population synthesis and is particularly well suited for stellar-population analysis of MaNGA galaxies. The first version of the library contains 8646 high-quality per-visit spectra for 3321 unique stars. Compared to photometry, the relative flux calibration of the library is accurate to 3.9% in g-r, 2.7% in r-i, and 2.2% in i-z. The data are released as part of SDSS Data Release 15. We expect the final release of the library to contain more than 10,000 stars.

Type Ia supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of them explode. Here, we present infant-phase detections of SN 2018aoz from a very low brightness of −10.5 AB absolute magnitude, revealing a hitherto unseen plateau in the B band that results in a rapid redward colour evolution between 1.0 and 12.4 hours after the estimated epoch of first light. The missing B-band flux is best explained by line-blanket absorption from Fe-peak elements in the outer 1% of the ejected mass. The observed B − V colour evolution of the supernova also matches the prediction from an over-density of Fe-peak elements in the same outer 1% of the ejected mass, whereas bluer colours are expected from a purely monotonic distribution of Fe-peak elements. The presence of excess nucleosynthetic material in the extreme outer layers of the ejecta points to enhanced surface nuclear burning or extended subsonic mixing processes in some normal type Ia SN explosions.

SN 2018aoz is a Type Ia SN with a B-band plateau and excess emission in infant-phase light curves ≲1 day after the first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN 2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show that the SN is intermediate between two subtypes of normal Type Ia: core normal and broad line. The excess emission may be attributable to the radioactive decay of surface iron-peak elements as well as the interaction of ejecta with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on Hα and He i favor a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O i] and He i disfavors a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN 2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe ii] and [Ni ii]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1D models are incompatible with the infant-phase excess emission, B max - V max color, and weak strength of nebular-phase [Ca ii]. Although the explosion processes of SN 2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear to be normal after ∼1 day.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is a dual-hemisphere, near-infrared (NIR), spectroscopic survey with the goal of producing a chemodynamical mapping of the Milky Way. The targeting for APOGEE-2 is complex and has evolved with time. In this paper, we present the updates and additions to the initial targeting strategy for APOGEE-2N presented in Zasowski et al. (2017). These modifications come in two implementation modes: (i) "Ancillary Science Programs"competitively awarded to Sloan Digital Sky Survey IV PIs through proposal calls in 2015 and 2017 for the pursuit of new scientific avenues outside the main survey, and (ii) an effective 1.5 yr expansion of the survey, known as the Bright Time Extension (BTX), made possible through accrued efficiency gains over the first years of the APOGEE-2N project. For the 23 distinct ancillary programs, we provide descriptions of the scientific aims, target selection, and how to identify these targets within the APOGEE-2 sample. The BTX permitted changes to the main survey strategy, the inclusion of new programs in response to scientific discoveries or to exploit major new data sets not available at the outset of the survey design, and expansions of existing programs to enhance their scientific success and reach. After describing the motivations, implementation, and assessment of these programs, we also leave a summary of lessons learned from nearly a decade of APOGEE-1 and APOGEE-2 survey operations. A companion paper, F. Santana et al. (submitted; AAS29036), provides a complementary presentation of targeting modifications relevant to APOGEE-2 operations in the Southern Hemisphere.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (2014 July-2017 July). This is the third data release for SDSS-IV, and the 15th from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 data cubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g., stellar and gas kinematics, emission-line and other maps) from the MaNGA Data Analysis Pipeline, and a new data visualization and access tool we call "Marvin." The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper, we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials, and examples of data use. Although SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

This paper documents the 16th data release (DR16) from the Sloan Digital Sky Surveys (SDSS), the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the Southern Hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey and new data from the SPectroscopic IDentification of ERosita Survey programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).