We present a new constraint on the Lyman Continuum (LyC) escape fraction at z~1.3. We obtain deep, high sensitivity far-UV imaging with the Advanced Camera for Surveys (ACS) Solar Blind Channel (SBC) on the Hubble Space Telescope (HST), targeting 11 star-forming galaxies at 1.2<z<1.4. The galaxies are selected from the 3D-HST survey to have high H$\alpha$ equivalent width (EW) with EW > 190 \AA, low stellar mass (M* < 10^10 M_sun) and U-band magnitude of U<24.2. These criteria identify young, low metallicity star bursting populations similar to the primordial star-forming galaxies believed to have reionized the universe. We do not detect any LyC signal (with S/N >3) in the individual galaxies or in the stack in the far-UV images. We place $3\sigma$ limits on the relative escape fraction of individual galaxies to be f_{esc,rel}<[0.10-0.22] and a stacked $3\sigma$ limit of f_{esc,rel}<0.07. Comparing to the confirmed LyC emitters from the literature, the galaxies in our sample span similar ranges of various galaxy properties including stellar mass, dust attenuation, and star formation rate (SFR). In particular, we compare the distribution of H$\alpha$ and [OIII] EWs of confirmed LyC emitters and non-detections including the galaxies in this study. Finally, we discuss if a dichotomy seen in the distribution of H$\alpha$ EWs can perhaps distinguish the LyC emitters from the non-detections.

Future surveys focusing on understanding the nature of dark energy (e.g., Euclid and WFIRST) will cover large fractions of the extragalactic sky in near-IR slitless spectroscopy. These surveys will detect a large number of galaxies that will have only one emission line in the covered spectral range. In order to maximize the scientific return of these missions, it is imperative that single emission lines are correctly identified. Using a supervised machine-learning approach, we classified a sample of single emission lines extracted from the WFC3 IR Spectroscopic Parallel survey (WISP), one of the closest existing analogs to future slitless surveys. Our automatic software integrates a SED fitting strategy with additional independent sources of information. We calibrated it and tested it on a "gold" sample of securely identified objects with multiple lines detected. The algorithm correctly classifies real emission lines with an accuracy of 82.6%, whereas the accuracy of the SED fitting technique alone is low (~50%) due to the limited amount of photometric data available (<=6 bands). While not specifically designed for the Euclid and WFIRST surveys, the algorithm represents an important precursor of similar algorithms to be used in these future missions.

The mutually complementary Euclid and Roman galaxy redshift surveys will use Halpha- and [OIII]-selected emission line galaxies as tracers of the large scale structure at $0.9 \lesssim z \lesssim 1.9$ (Halpha) and $1.5 \lesssim z \lesssim 2.7$ ([OIII]). It is essential to have a reliable and sufficiently precise knowledge of the expected numbers of Halpha-emitting galaxies in the survey volume in order to optimize these redshift surveys for the study of dark energy. Additionally, these future samples of emission-line galaxies will, like all slitless spectroscopy surveys, be affected by a complex selection function that depends on galaxy size and luminosity, line equivalent width, and redshift errors arising from the misidentification of single emission-line galaxies. Focusing on the specifics of the Euclid survey, we combine two slitless spectroscopic WFC3-IR datasets -- 3D-HST+AGHAST and the WISP survey -- to construct a Euclid-like sample that covers an area of 0.56 deg$^2$ and includes 1277 emission line galaxies. We detect 1091 ($\sim$3270 deg$^{-2}$) Halpha+[NII]-emitting galaxies in the range $0.9\leq z \leq 1.6$ and 162 ($\sim$440 deg$^{-2}$) [OIII]$\lambda$5007-emitters over $1.5\leq z \leq 2.3$ with line fluxes $\geq 2 \times 10^{-16}$ erg s$^{-1}$ cm$^{-2}$. The median of the Halpha+[NII] equivalent width distribution is $\sim$250\r{A}, and the effective radii of the continuum and Halpha+[NII] emission are correlated with a median of $\sim$0.38" and significant scatter ($\sigma \sim $0.2"$-$0.35"). Finally, we explore the prevalence of redshift misidentification in future Euclid samples, finding potential contamination rates of $\sim$14-20% and $\sim$6% down to $2\times 10^{-16}$ and $6 \times 10^{-17}$ erg s$^{-1}$ cm$^{-2}$, respectively, though with increased wavelength coverage these percentages drop to nearly zero.

Looking for bright galaxies born in the early universe is fundamental to investigating the Epoch of Reionization, the era when the first stars and galaxies ionized the intergalactic medium. We utilize Hubble Space Telescope pure parallel imaging to select galaxy candidates at a time 500 to 650 million years after the Big Bang, which corresponds to redshifts z ~ 8-10. These data come from the Brightest of Reionizing Galaxies Survey (BoRG) Cycle 22 dataset, which consists of pure-parallel imaging in ~ 90 different lines of sight that sum up to an area of ~ 420 arcmin^2. This survey uses five filters and has the advantage (compared to the Cycle 21 BoRG program) of including imaging in the JH140 band, covering continuous wavelengths from the visible to near-infrared (lambda = 0.35um - 1.7um). This allows us to perform reliable selection of galaxies at z>8 using the photometric redshift technique. We use these galaxy candidates to constrain the bright end of the rest-frame ultraviolet luminosity function in this epoch. These candidates are excellent targets for follow-up observations, particularly with the James Webb Space Telescope.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

We present the results of Galaxy Nurseries project, which was designed to enable crowdsourced analysis of slitless spectroscopic data by volunteer citizen scientists using the Zooniverse online interface. The dataset was obtained by the WFC3 Infrared Spectroscopic Parallel (WISP) Survey collaboration and comprises NIR grism (G102 and G141) and direct imaging. Volunteers were instructed to evaluate indicated spectral features and decide whether it was a genuine emission line or more likely an artifact. Galaxy Nurseries was completed in only 40 days, gathering 414,360 classifications from 3003 volunteers for 27,333 putative emission lines. The results of Galaxy Nurseries demonstrate the feasibility of identifying genuine emission lines in slitless spectra by citizen scientists. Volunteer responses for each subject were aggregated to compute $f_{\mathrm{Real}}$, the fraction of volunteers who classified the corresponding emission line as "Real". To evaluate the accuracy of volunteer classifications, their aggregated responses were compared with independent assessments provided by members of the WISP Survey Science Team (WSST). Overall, there is a broad agreement between the WSST and volunteers' classifications, although we recognize that robust scientific analyses typically require samples with higher purity and completeness than raw volunteer classifications provide. Nonetheless, choosing optimal threshold values for $f_{\mathrm{Real}}$ allows a large fraction of spurious lines to be vetoed, substantially reducing the timescale for subsequent professional analysis of the remaining potential lines.

We present the results of a systematic search for Lyman-alpha emitters (LAEs) at $6 \lesssim z \lesssim 7.6$ using the HST WFC3 Infrared Spectroscopic Parallel (WISP) Survey. Our total volume over this redshift range is $\sim 8 \times10^5$ Mpc$^3$, comparable to many of the narrowband surveys despite their larger area coverage. We find two LAEs at $z=6.38$ and $6.44$ with line luminosities of L$_{\mathrm{Ly}\alpha} \sim 4.7 \times 10^{43}$ erg s$^{-1}$, putting them among the brightest LAEs discovered at these redshifts. Taking advantage of the broad spectral coverage of WISP, we are able to rule out almost all lower-redshift contaminants. The WISP LAEs have a high number density of $7.7\times10^{-6}$ Mpc$^{-3}$. We argue that the LAEs reside in Mpc-scale ionized bubbles that allow the Lyman-alpha photons to redshift out of resonance before encountering the neutral IGM. We discuss possible ionizing sources and conclude that the observed LAEs alone are not sufficient to ionize the bubbles.

We present the first results from our survey of the star-forming complex W3, combining VRI photometry with multiobject spectroscopy to identify and characterize the high-mass stellar population across the region. With 79 new spectral classifications, we bring the total number of spectroscopically-confirmed O- and B-type stars in W3 to 105. We find that the high-mass slope of the mass function in W3 is consistent with a Salpeter IMF, and that the extinction toward the region is best characterized by an Rv of approximately 3.6. B-type stars are found to be more widely dispersed across the W3 giant molecular cloud (GMC) than previously realized: they are not confined to the high-density layer (HDL) created by the expansion of the neighboring W4 HII region into the GMC. This broader B-type population suggests that star formation in W3 began spontaneously up to 8--10 Myr ago, although at a lower level than the more recent star formation episodes in the HDL. In addition, we describe a method of optimizing sky subtraction for fiber spectra in regions of strong and spatially-variable nebular emission.

Upcoming space-based surveys such as Euclid and WFIRST-AFTA plan to measure Baryonic Acoustic Oscillations (BAOs) in order to study dark energy. These surveys will use IR slitless grism spectroscopy to measure redshifts of a large number of galaxies over a significant redshift range. In this paper, we use the WFC3 Infrared Spectroscopic Parallel Survey (WISP) to estimate the expected number of Halpha (Ha) emitters observable by these future surveys. WISP is an ongoing HST slitless spectroscopic survey, covering the 0.8-1.65micron wavelength range and allowing the detection of Ha emitters up to z~1.5 and [OIII] emitters to z~2.3. We derive the Ha-[OIII] bivariate line luminosity function for WISP galaxies at z~1 using a maximum likelihood estimator that properly accounts for uncertainties in line luminosity measurement, and demonstrate how it can be used to derive the Ha luminosity function from exclusively fitting [OIII] data. Using the z~2 [OIII] line luminosity function, and assuming that the relation between Ha and [OIII] luminosity does not change significantly over the redshift range, we predict the Ha number counts at z~2 - the upper end of the redshift range of interest for the future surveys. For the redshift range 0.7<z<2, we expect ~3000 galaxies/deg^2 for a flux limit of 3x10^{-16} ergs/s/cm^2 (the proposed depth of Euclid galaxy redshift survey) and ~20,000 galaxies/deg^2 for a flux limit of ~10^{-16} ergs/s/cm^2 (the baseline depth of WFIRST galaxy redshift survey).

We present near infrared polarimetry data of background stars shining through a selection of starless cores taken in the $K$ band, probing visual extinctions up to $A_{V} \sim 48$. We find that $P_K/{\tau _K}$ continues to decline with increasing $A_{V}$ with a power law slope of roughly -0.5. Examination of published submillimeter (submm) polarimetry of starless cores suggests that by $A_{V} \gtrsim 20$ the slope for $P$ vs. $\tau$ becomes $\sim -1$, indicating no grain alignment at greater optical depths. Combining these two data sets, we find good evidence that, in the absence of a central illuminating source, the dust grains in dense molecular cloud cores with no internal radiation source cease to become aligned with the local magnetic field at optical depths greater than $A_V \sim 20$. A simple model relating the alignment efficiency to the optical depth into the cloud reproduces the observations well.

An N-body hybrid simulation, integrating both massive and tracer particles, of a Galactic disk is used to study the stellar phase space distribution or velocity distributions in different local neighborhoods. Pattern speeds identified in Fourier spectrograms suggest that two-armed and three-armed spiral density waves, a bar and a lopsided motion are coupled in this simulation, with resonances of one pattern lying near resonances of other patterns. We construct radial and tangential (uv) velocity distributions from particles in different local neighborhoods. More than one clump is common in these local velocity distributions regardless of the position in the disk. Features in the velocity distribution observed at one galactic radius are also seen in nearby neighborhoods (at larger and smaller radii) but with shifted mean v values. This is expected if the v velocity component of a clump sets the mean orbital galactic radius of its stars. We find that gaps in the velocity distribution are associated with the radii of kinks or discontinuities in the spiral arms. These gaps also seem to be associated with Lindblad resonances with spiral density waves and so denote boundaries between different dominant patterns in the disk. We discuss implications for interpretations of the Milky Way disk based on local velocity distributions. Velocity distributions created from regions just outside the bar's Outer Lindblad resonance and with the bar oriented at 45 degrees from the Sun-Galactic center line more closely resemble that seen in the solar neighborhood (triangular in shape at lower uv and with a Hercules like stream) when there is a strong nearby spiral arm, consistent with the observed Centaurus Arm tangent, just interior to the solar neighborhood.

We present the first observations from a large-scale survey of nearby (z < 1) long-duration gamma-ray burst (LGRB) host galaxies, which consist of eight rest-frame optical spectra obtained at Keck and Magellan. Along with two host galaxy observations from the literature, we use optical emission line diagnostics to determine metallicities, ionization parameters, young stellar population ages, and star formation rates. We compare the LGRB host environments to a variety of local and intermediate-redshift galaxy populations, as well as the newest grid of stellar population synthesis and photoionization models generated with the Starburst99/Mappings codes. With these comparisons we investigate whether the GRB host galaxies are consistent with the properties of the general galaxy population, and therefore whether they may be used as reliable tracers of star formation. We find that LGRB host galaxies generally have low-metallicity ISM environments out to z ~ 1. The ISM properties of our GRB hosts, including metallicity, ionization parameter, and young stellar population age, are significantly different from the general galaxy population, host galaxies of nearby broad-lined Type Ic supernovae, and nearby metal-poor galaxies.

The morphology of the outer rings of early-type spiral galaxies is compared to integrations of massless collisionless particles initially in nearly circular orbits. Particles are perturbed by a quadrupolar gravitational potential corresponding to a growing and secularly evolving bar. We find that outer rings with R1R2 morphology and pseudorings are exhibited by the simulations even though they lack gaseous dissipation. Simulations with stronger bars form pseudorings earlier and more quickly than those with weaker bars. We find that the R1 ring, perpendicular to the bar, is fragile and dissolves after a few bar rotation periods if the bar pattern speed increases by more than ~ 8%, bar strength increases (by >~ 140%) after bar growth, or the bar is too strong (Q_T>0.3). If the bar slows down after formation, pseudoring morphology persists and the R2 ring perpendicular to the bar is populated due to resonance capture. The R2 ring remains misaligned with the bar and increases in ellipticity as the bar slows down. The R2 ring becomes scalloped and does not resemble any ringed galaxies if the bar slows down more than 3.5% suggesting that bars decrease in strength before they slow down this much. We compare the morphology of our simulations to B-band images of 9 ringed galaxies from the Ohio State University Bright Spiral Galaxy Survey, and we find a reasonable match in morphologies to R1R2' pseudorings seen within a few bar rotation periods of bar formation. Some of the features previously interpreted in terms of dissipative models may be due to transient structure associated with recent bar growth and evolution.