In the $\Lambda$CDM paradigm the Galactic stellar halo is predicted to harbor the accreted debris of smaller systems. To identify these systems, the H3 Spectroscopic Survey, combined with $Gaia$, is gathering 6D phase-space and chemical information in the distant Galaxy. Here we present a comprehensive inventory of structure within 50 kpc from the Galactic center using a sample of 5684 giants at $|b|>40^{\circ}$ and $|Z|>2$ kpc. We identify known structures including the high-$\alpha$ disk, the in-situ halo (disk stars heated to eccentric orbits), Sagittarius (Sgr), $Gaia$-Sausage-Enceladus (GSE), the Helmi Streams, Sequoia, and Thamnos. Additionally, we identify the following new structures: (i) Aleph ([Fe/H]$=-0.5$), a low eccentricity structure that rises a surprising 10 kpc off the plane, (ii, iii) Arjuna ([Fe/H]$=-1.2$) and I'itoi ([Fe/H]$<-2$), which comprise the high-energy retrograde halo along with Sequoia, and (iv) Wukong ([Fe/H]$=-1.6$), a prograde phase-space overdensity chemically distinct from GSE. For each structure we provide [Fe/H], [$\alpha$/Fe], and orbital parameters. Stars born within the Galaxy are a major component at $|Z|\sim$2 kpc ($\approx$60$\%$), but their relative fraction declines sharply to $\lesssim$5$\%$ past 15 kpc. Beyond 15 kpc, $>$80$\%$ of the halo is built by two massive ($M_{\star}\sim10^{8}-10^{9}M_{\odot}$) accreted dwarfs: GSE ([Fe/H]$=-1.2$) within 25 kpc, and Sgr ([Fe/H]$=-1.0$) beyond 25 kpc. This explains the relatively high overall metallicity of the halo ([Fe/H]$\approx-1.2$). We attribute $\gtrsim$95$\%$ of the sample to one of the listed structures, pointing to a halo built entirely from accreted dwarfs and heating of the disk.

The archaeological record of stars in the Milky Way opens a uniquely detailed window into the early formation and assembly of galaxies. Here we use 11,000 main-sequence turn-off stars with well-measured ages, [Fe/H], [$\alpha$/Fe], and orbits from the H3 Survey and Gaia to time the major events in the early Galaxy. Located beyond the Galactic plane, $1\lesssim |Z|/\rm kpc \lesssim4$, this sample contains three chemically distinct groups: a low metallicity population, and low-$\alpha$ and high-$\alpha$ groups at higher metallicity. The age and orbit distributions of these populations show that: 1) the high-$\alpha$ group, which includes both disk stars and the in-situ halo, has a star-formation history independent of eccentricity that abruptly truncated $8.3\pm0.1$ Gyr ago ($z\simeq1$); 2) the low metallicity population, which we identify as the accreted stellar halo, is on eccentric orbits and its star formation truncated $10.2.^{+0.2}_{-0.1}$ Gyr ago ($z\simeq2$); 3) the low-$\alpha$ population is primarily on low eccentricity orbits and the bulk of its stars formed less than 8 Gyr ago. These results suggest a scenario in which the Milky Way accreted a satellite galaxy at $z\approx2$ that merged with the early disk by $z\approx1$. This merger truncated star formation in the early high-$\alpha$ disk and perturbed a fraction of that disk onto halo-like orbits. The merger enabled the formation of a chemically distinct, low-$\alpha$ disk at $z\lesssim1$. The lack of any stars on halo-like orbits at younger ages indicates that this event was the last significant disturbance to the Milky Way disk.

The protagonists of cosmic reionization remain elusive. Faint star-forming galaxies are leading candidates because they are numerous and may have significant ionizing photon escape fractions ($f_{esc}$). Here we update this picture via an empirical model that successfully predicts latest observations (e.g., the drop in star-formation density at z>8). We generate an ionizing spectrum for each galaxy in our model and constrain $f_{esc}$ using latest measurements of the reionization timeline (e.g., Ly$\alpha$ damping of quasars and galaxies at z>7). Assuming a constant $f_{esc}$, we find $M_{UV}$<-13.5 galaxies need $f_{esc}=0.21^{+0.06}_{-0.04}$ to complete reionization. The inferred IGM neutral fraction is [0.9, 0.5, 0.1] at z=[8.2, 6.8, 6.2]$\pm$0.2, i.e., the bulk of reionization transpires in 300 Myrs. Inspired by the emergent sample of Lyman Continuum (LyC) leakers that overwhelmingly displays higher-than-average star-formation surface density ($\Sigma$), we propose a model relating $f_{esc}$ to $\Sigma$ and find $f_{esc}\propto\Sigma^{0.4\pm0.1}$. Since $\Sigma$ falls by ~2.5 dex between z=8 and z=0, our model explains the humble upper limits on $f_{esc}$ at lower redshifts and its required evolution to ~0.2 at z>6. Within this model, strikingly, <5% of galaxies with $M_{UV}$<-18 (the `oligarchs') account for >80% of the reionization budget. In fact, faint sources ($M_{UV}$>-16) must be relegated to a limited role to ensure high neutral fractions at z=7-8. Shallow faint-end slopes of the UV luminosity function ($\alpha$>-2) and/or $f_{esc}$ distributions skewed toward bright galaxies produce the required late and rapid reionization. We predict LyC leakers like COLA1 (z=6.6, $f_{esc}$~30%, $M_{UV}$=-21.5) become increasingly common towards z~6 and that the drivers of reionization do not lie hidden across the faint-end of the luminosity function, but are already known to us. (abridged)

The $100^\circ$-long thin stellar stream in the Milky Way halo, GD-1, has an ensemble of features that may be due to dynamical interactions. Using high-resolution MMT/Hectochelle spectroscopy we show that a spur of GD-1-like stars outside of the main stream are kinematically and chemically consistent with the main stream. In the spur, as in the main stream, GD-1 has a low intrinsic radial velocity dispersion, $\sigma_{V_r}\lesssim1\,\rm km\,s^{-1}$, is metal-poor, $\rm [Fe/H]\approx-2.3$, with little $\rm [Fe/H]$ spread and some variation in $\rm [\alpha/Fe]$ abundances, which point to a common globular cluster progenitor. At a fixed location along the stream, the median radial velocity offset between the spur and the main stream is smaller than $0.5\,\rm km\,s^{-1}$, comparable to the measurement uncertainty. A flyby of a massive, compact object can change orbits of stars in a stellar stream and produce features like the spur observed in GD-1. In this scenario, the radial velocity of the GD-1 spur relative to the stream constrains the orbit of the perturber and its current on-sky position to $\approx5,000\,\rm deg^2$. The family of acceptable perturber orbits overlaps the stellar and dark-matter debris of the Sagittarius dwarf galaxy in present-day position and velocity. This suggests that GD-1 may have been perturbed by a globular cluster or an extremely compact dark-matter subhalo formerly associated with Sagittarius.

The Galactic stellar halo is predicted to have formed at least partially from the tidal disruption of accreted dwarf galaxies. This assembly history should be detectable in the orbital and chemical properties of stars. The H3 Survey is obtaining spectra for 200,000 stars, and, when combined with Gaia data, is providing detailed orbital and chemical properties of Galactic halo stars. Unlike previous surveys of the halo, the H3 target selection is based solely on magnitude and Gaia parallax; the survey therefore provides a nearly unbiased view of the entire stellar halo at high latitudes. In this paper we present the distribution of stellar metallicities as a function of Galactocentric distance and orbital properties for a sample of 4232 kinematically-selected halo giants to 100 kpc. The stellar halo is relatively metal-rich, [Fe/H]=-1.2, and there is no discernable metallicity gradient over the range $6<R_{\rm gal}<100$ kpc. However, the halo metallicity distribution is highly structured including distinct metal-rich and metal-poor components at $R_{\rm gal}<10$ kpc and $R_{\rm gal}>30$ kpc, respectively. Metal-poor stars with [Fe/H]$<-2$ are a small population of the halo at all distances and orbital categories. We associate the "in-situ" stellar halo with stars displaying thick-disk chemistry on halo-like orbits; such stars are confined to $|z|<10$ kpc. The majority of the stellar halo is resolved into discrete features in orbital-chemical space, suggesting that the bulk of the stellar halo formed from the accretion and tidal disruption of dwarf galaxies. (ABRIDGED)

The timing argument provides a lower limit on the mass of the Milky Way. We find, using a sample of 32 stars at $R > 60$ kpc drawn from the H3 Spectroscopic Survey and mock catalogs created from published numerical simulations, that M$_{200} > 0.91\times 10^{12}$ M$_\odot$ with 90% confidence. We recommend using this limit to refine the allowed prior mass range in more complex and sophisticated statistical treatments of Milky Way dynamics. The use of such a prior would have significantly reduced many previously published uncertainty ranges. Our analysis suggests that the most likely value of M$_{200}$ is $\sim 1.4 \times 10^{12}$ M$_\odot$, but establishing this as the Milky Way mass requires a larger sample of outer halo stars and a more complete analysis of the inner halo stars in H3. The imminent growth in the sample of outer halo stars due to ongoing and planned surveys will make this possible.

Modelling reionization often requires significant assumptions about the properties of ionizing sources. Here, we infer the total output of hydrogen-ionizing photons (the ionizing emissivity, $\dot{N}_\textrm{ion}$) at $z=4-14$ from current reionization constraints, being maximally agnostic to the properties of ionizing sources. We use a Bayesian analysis to fit for a non-parametric form of $\dot{N}_\textrm{ion}$, allowing us to flexibly explore the entire prior volume. We infer a declining $\dot{N}_\textrm{ion}$ with redshift at $z>6$, which can be used as a benchmark for reionization models. Model-independent reionization constraints from the CMB optical depth and Ly$\alpha$ and Ly$\beta$ forest dark pixel fraction produce $\dot{N}_\textrm{ion}$ evolution ($d\log_{10}\dot{N}_\textrm{ion}/dz|_{z=6\rightarrow8} = -0.31\pm0.35$ dex) consistent with the declining UV luminosity density of galaxies, assuming constant ionizing photon escape fraction and efficiency. Including measurements from Ly$\alpha$ damping of galaxies and quasars produces a more rapid decline: $d\log_{10}\dot{N}_\textrm{ion}/dz|_{z=6\rightarrow8} =-0.44\pm0.22$ dex, steeper than the declining galaxy luminosity density (if extrapolated beyond $M_\mathrm{UV} \lesssim -13$), and constrains the mid-point of reionization to $z = 6.93\pm0.14$.

Modern theories of galaxy formation predict that the Galactic stellar halo was hierarchically assembled from the accretion and disruption of smaller systems. This hierarchical assembly is expected to produce a high degree of structure in the combined phase and chemistry space; this structure should provide a relatively direct probe of the accretion history of our Galaxy. Revealing this structure requires precise 3D positions (including distances), 3D velocities, and chemistry for large samples of stars. The Gaia satellite is delivering proper motions and parallaxes for >1 billion stars to G~20. However, radial velocities and metallicities will only be available to G~15, which is insufficient to probe the outer stellar halo (>10 kpc). Moreover, parallaxes will not be precise enough to deliver high-quality distances for stars beyond ~10 kpc. Identifying accreted systems throughout the stellar halo therefore requires a large ground-based spectroscopic survey to complement Gaia. Here we provide an overview of the H3 Stellar Spectroscopic Survey, which will deliver precise stellar parameters and spectrophotometric distances for 200,000 stars to r=18. Spectra are obtained with the Hectochelle instrument at the MMT, which is configured for the H3 Survey to deliver resolution R~23,000 spectra covering the wavelength range 5150A-5300A. The survey is optimized for stellar halo science and therefore focuses on high Galactic latitude fields (|b|>30 deg.), sparsely sampling 15,000 sq. degrees. Targets are selected on the basis of Gaia parallaxes, enabling very efficient selection of bone fide halo stars. The survey began in the Fall of 2017 and has collected 88,000 spectra to-date. All of the data, including the derived stellar parameters, will eventually be made publicly available via the survey website: h3survey.rc.fas.harvard.edu.

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy.

We present a color-based method for identifying red giants based on Pan-STARRS grz and WISE W1 and W2 photometry. We utilize a subsample of bright stars with precise parallaxes from Gaia DR2 to verify that the color-based selection reliably separates dwarfs from giants. The selection is conservative in the sense that contamination is small (~30%) but not all giants are included (the selection primarily identifies K giants). The color-based selection can be applied to stars brighter than $W1\approx16$, more than two magnitudes fainter than techniques relying on shallower 2MASS photometry. Many streams and clouds are visible in the resulting sky maps, especially when binned by Gaia DR2 proper motions, including the Sagittarius stream, the Hercules-Aquila Cloud, the Eastern Banded Structure, Monoceros, and the Virgo Overdensity. In addition to the characterization of new and known stellar streams, we expect that this method for selecting red giants will enable detailed analysis of the diffuse stellar halo to distances exceeding 100 kpc.

We use deep F275W imaging from the $Hubble$ Deep UV Legacy Survey (HDUV) and G280 grism spectroscopy from $HST$/WFC3, along with new and archival optical spectra from Keck/DEIMOS, to search for candidate ionizing sources in the GOODS-N field at z ~ 2.5 - 3. Spectroscopic identification of our UV-selected sources are 99% complete to F275W = 25.5 in the region of the UV imaging, and we identify 6 potential ionizing galaxies or AGNs at z ~ 3. By far the brightest of these is a z = 2.583 AGN that totally dominates the ionizing flux in the region, with a specific ionizing volume emissivity at $912~ \AA$ of $\epsilon_{912}$ = $8.3^{27}_{1.4} \times 10^{24}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$. Based on our spectroscopic data, we find four candidates are contaminated by foreground galaxies at z ~ 0.5 - 0.7. At $\epsilon_{912}$ = $2.2^{7.2}_{0.4} \times 10^{23}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$, the remaining candidate galaxy's contribution to the ionizing background lies well below the flux required to ionize the intergalactic medium at z ~ 2.5 - 3, consistent with previous observations that show AGNs provide the bulk of the ionizing background at these redshifts.

We present the Hubble Deep UV Legacy Survey (HDUV), a 132 orbit imaging program with the WFC3/UVIS camera onboard the Hubble Space Telescope (HST). The HDUV extends and builds on the few previous UV imaging surveys in the two GOODS/CANDELS-Deep fields to provide deep images over a total area of ~100 arcmin2 in the two filters F275W and F336W. Our release also includes all the F275W imaging data taken by the CANDELS survey, which were aligned using a novel approach and combined with the HDUV survey data. By reaching depths of 27.5-28.0 mag (5sigma, in 0.4" apertures), these are the deepest high-resolution UV data over such a large area taken to date. Such unique UV imaging enables a wide range of science by the community. Among the main goals of the HDUV survey are: (1) provide a complete sample of faint star-forming galaxies at z~1-3, (2) constrain the ionizing photon escape fraction from galaxies at z~2-3, and (3) track the build-up of bulges and the disappearance of clumpy disk galaxies through reliable internal stellar population properties at sub-kpc resolution out to z~3. The addition of the HDUV data further enhances the legacy value of the two GOODS/CANDELS-Deep fields, which now include deep 11-band HST imaging as well as very deep ancillary data from X-ray to radio, enabling unique multi-wavelength studies. Here, we provide an overview of the survey design, describe the data reduction, and highlight a few basic analyses on the images which are released to the community as high level science products via the Mikulski Archive for Space Telescopes (MAST).

Recent work has suggested extreme [OIII] emitting star-forming galaxies are important to reionization. Relatedly, [OIII]/[OII] has been put forward as an indirect estimator of the Lyman Continuum (LyC) escape fraction ($f_{esc}$) at $z\gtrsim4.5$ when the opaque IGM renders LyC photons unobservable. Using deep archival U-band (VLT/VIMOS) imaging of a recently confirmed overdensity at $z\sim3.5$ we calculate tight constraints on $f_{esc}$ for a sample (N=73) dominated by extreme [OIII] emitters. We find no Lyman Continuum signal ($f_{esc}^{rel} < 6.3^{+0.7}_{-0.7} \%$ at $1\sigma$) in a deep U-band stack of our sample (31.98 mag at 1$\sigma$). This constraint is in agreement with recent studies of star-forming galaxies spanning $z\sim1-4$ that have found very low average $f_{esc}$. Despite the galaxies in our study having an estimated average rest-frame EW([OIII]$\lambda5007$)$\sim400\AA$ and [OIII]/[OII]$\sim 4$ from composite SED-fitting, we find no LyC detection, which brings into question the potential of [OIII]/[OII] as an effective probe of the LyC--a majority of LyC emitters have [OIII]/[OII]$>3$, but we establish here that [OIII]/[OII]$>3$ does not guarantee significant LyC leakage for a population. Since even extreme star-forming galaxies are unable to produce the $f_{esc}\sim10-15\%$ required by most theoretical calculations for star-forming galaxies to drive reionization, there must either be a rapid evolution of $f_{esc}$ between $z\sim3.5$ and the Epoch of Reionization, or hitherto observationally unstudied sources (e.g. ultra-faint low-mass galaxies with $\log(M/M_\odot)\sim7-8.5$) must make an outsized contribution to reionization.

We present six galaxies at z~2 that show evidence of Lyman continuum (LyC) emission based on the newly acquired UV imaging of the Hubble Deep UV legacy survey (HDUV) conducted with the WFC3/UVIS camera on the Hubble Space Telescope (HST). At the redshift of these sources, the HDUV F275W images partially probe the ionizing continuum. By exploiting the HST multi-wavelength data available in the HDUV/GOODS fields, models of the UV spectral energy distributions, and detailed Monte Carlo simulations of the intergalactic medium absorption, we estimate the absolute ionizing photon escape fractions of these galaxies to be very high -- typically >60% (>13% for all sources at 90% likelihood). Our findings are in broad agreement with previous studies that found only a small fraction of galaxies to show high escape fraction. These six galaxies comprise the largest sample yet of LyC leaking candidates at z~2 whose inferred LyC flux has been cleanly observed at HST resolution. While three of our six candidates show evidence of hosting an active galactic nucleus (AGN), two of these are heavily obscured and their LyC emission appears to originate from star-forming regions rather than the central nucleus. This suggests an AGN-aided pathway for LyC escape from these sources. Extensive multi-wavelength data in the GOODS fields, especially the near-IR grism spectra from the 3D-HST survey, enable us to study the candidates in detail and tentatively test some recently proposed indirect methods to probe LyC leakage -- namely, the [OIII]/[OII] line ratio and the H$\beta-$UV slope diagram. High-resolution spectroscopic followup of our candidates will help constrain such indirect methods which are our only hope of studying $f_{esc}$ at z~5-9 in the fast-approaching era of the James Webb Space Telescope.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.