Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7µm and 3.0-5.0µm, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed highredshift galaxies. We observed 46 WFC3/IR pointings (∼200 arcmin 2 ) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ∼ 6 − 10 candidates, including the brightest known at z ∼ 6, and the most spatially-resolved distant lensed arc known at z ∼ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ∼ 10 candidates from z ∼ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO Cycle 1 call for proposals.

PypeIt is a Python package for semi-automated reduction of astronomical, spectroscopic data. Its algorithms build on decades-long development of previous data reduction pipelines by the developers (Bernstein, Burles, & Prochaska, 2015;Bochanski et al., 2009). The reduction procedure -including a complete list of the input parameters and available functionality -is provided as online documentation hosted by Read the Docs, which is regularly updated. In what follows, we provide a brief description of the algorithms, but refer the interested reader to the online documentation for the most up-to-date information.Release v1.0.3 serves the following spectrographs: Gemini/GNIRS, Gemini/GMOS, Gemini/FLAMINGOS 2, Lick/Kast, Magellan/MagE, Magellan/Fire, MDM/OSMOS, Keck/DEIMOS (600ZD, 830G, 1200G), Keck/LRIS, Keck/MOSFIRE (Y, J, K gratings tested), Keck/NIRES, Keck/NIRSPEC (low-dispersion), LBT/Luci-I, Luci-II, LBT/MODS (beta), NOT/ALFOSC (grism4), VLT/X-Shooter (VIS, NIR), VLT/FORS2 (300I, 300V).

We investigate the impact of local environment on the galaxy stellar mass function (SMF) spanning a wide range of galaxy densities from the field up to dense cores of massive galaxy clusters. Data are drawn from a sample of eight fields from the Observations of Redshift Evolution in Large-Scale Environments (ORELSE) survey. Deep photometry allow us to select mass-complete samples of galaxies down to 10 9 M . Taking advantage of >4000 secure spectroscopic redshifts from ORELSE and precise photometric redshifts, we construct 3-dimensional density maps between 0.55 < z < 1.3 using a Voronoi tessellation approach. We find that the shape of the SMF depends strongly on local environment exhibited by a smooth, continual increase in the relative numbers of high-to low-mass galaxies towards denser environments. A straightforward implication is that local environment proportionally increases the efficiency of (a) destroying lower-mass galaxies and/or (b) growth of higher-mass galaxies. We also find a presence of this environmental dependence in the SMFs of star-forming and quiescent galaxies, although not quite as strongly for the quiescent subsample. To characterize the connection between the SMF of field galaxies and that of denser environments we devise a simple semi-empirical model. The model begins with a sample of ≈10 6 galaxies at z start =5 with stellar masses distributed according to the field. Simulated galaxies then evolve down to z final =0.8 following empirical prescriptions for star-formation, quenching, and galaxy-galaxy merging. We run the simulation multiple times, testing a variety of scenarios with differing overall amounts of merging. Our model suggests that a large number of mergers are required to reproduce the SMF in dense environments. Additionally, a large majority of these mergers would have to occur in intermediate density environments (e.g. galaxy groups).

Using ∼5000 spectroscopically-confirmed galaxies drawn from the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey we investigate the relationship between color and galaxy density for galaxy populations of various stellar masses in the redshift range 0.55 ≤ z ≤ 1.4. The fraction of galaxies with colors consistent with no ongoing star formation ( f q ) is broadly observed to increase with increasing stellar mass, increasing galaxy density, and decreasing redshift, with clear differences observed in f q between field and group/cluster galaxies at the highest redshifts studied. We use a semi-empirical model to generate a suite of mock group/cluster galaxies unaffected by environmentally-specific processes and compare these galaxies at fixed stellar mass and redshift to observed populations to constrain the efficiency of environmentally-driven quenching (Ψ convert ). High-density environments from 0.55 ≤ z ≤ 1.4 appear capable of efficiently quenching galaxies with log(M * /M ⊙ ) > 10.45. Lower stellar mass galaxies also appear efficiently quenched at the lowest redshifts studied here, but this quenching efficiency is seen to drop precipitously with increasing redshift. Quenching efficiencies, combined with simulated group/cluster accretion histories and results on the star formation rate-density relation from a companion ORELSE study, are used to constrain the average time from group/cluster accretion to quiescence and the elapsed time between accretion and the inception of the quenching event. These timescales were constrained to be t convert = 2.4 ± 0.3 and t delay = 1.3 ± 0.4 Gyr, respectively, for galaxies with log(M * /M ⊙ ) > 10.45 and t convert = 3.3 ± 0.3 and t delay = 2.2 ± 0.4 Gyr for lower stellar mass galaxies. These quenching efficiencies and associated timescales are used to rule out certain environmental mechanisms as being the primary processes responsible for transforming the star-formation properties of galaxies over this 4 Gyr window in cosmic time.