The origin of the correlations between mass, morphology, quenched fraction, and formation history in galaxies is difficult to define, primarily due to the uncertainties in galaxy star-formation histories. Star-formation histories are better constrained for higher redshift galaxies, observed closer to their formation and quenching epochs. Here we use "non-parametric" star-formation histories and a nested sampling method to derive constraints on the formation and quenching timescales of quiescent galaxies at $0.7<z<2.5$. We model deep HST grism spectroscopy and photometry from the CLEAR (CANDELS Lyman$-\alpha$ Emission at Reionization) survey. The galaxy formation redshifts, $z_{50}$ (defined as the point where they had formed 50\% of their stellar mass) range from $z_{50}\sim 2$ (shortly prior to the observed epoch) up to $z_{50} \simeq 5-8$. \editone{We find that early formation redshifts are correlated with high stellar-mass surface densities, $\log \Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) >$10.25, where $\Sigma_1$ is the stellar mass within 1~pkpc (proper kpc). Quiescent galaxies with the highest stellar-mass surface density, $\log\Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) > 10.25$, } show a \textit{minimum} formation redshift: all such objects in our sample have $z_{50} > 2.9$. Quiescent galaxies with lower surface density, $\log \Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) = 9.5 - 10.25$, show a range of formation epochs ($z_{50} \simeq 1.5 - 8$), implying these galaxies experienced a range of formation and assembly histories. We argue that the surface density threshold $\log\Sigma_1/(M_\odot\ \mathrm{kpc}^{-2})>10.25$ uniquely identifies galaxies that formed in the first few Gyr after the Big Bang, and we discuss the implications this has for galaxy formation models.

Redshift $z=9--10$ object selection is the effective limit of Hubble Space Telescope imaging capability, even when confirmed with Spitzer. If only a few photometry data points are available, it becomes attractive to add criteria based on their morphology in these J- and H-band images. One could do so through visual inspection, a size criterion, or alternate morphometrics. We explore a vetted sample of BoRG $z\sim9$ and $z\sim10$ candidate galaxies and the object rejected by Morishita+ (2018) to explore the utility of a size criterion in z=9-10 candidate selection. A stringent, PSF-corrected effective radius criterion ($r_e<0\farcs3$) would result in the rejection of 65-70\% of the interlopers visually rejected by Morishita+. It may also remove up to $\sim20$\% of bona-fide brightest ($L>>L^*$) z=9 or 10 candidates from a BoRG selected sample based on the Mason+ (2015) luminosity functions, assuming the Holwerda+ (2015) $z\sim9$ size-luminosity relation. We argue that including a size constraint in lieu of a visual inspection may serve in wide-field searches for these objects in e.g. EUCLID or HST archival imaging with the understanding that some brightest ($L>>L^*$) candidates may be missed. The sizes of the candidates found by Morishita+ (2018) follow the expected size distribution of $z\sim9$ for bright galaxies, consistent with the lognormal in Shibuya+ (2015) and single objects. Two candidates show high star-formation surface density ($\Sigma_{SFR} > 25 M_\odot/kpc^2$) and all merit further investigation and follow-up observations.

We study the global dust and (molecular) gas content in the Lyman Alpha Reference Sample (LARS), i.e. 14 local star-forming galaxies. We characterize their interstellar medium and relate newly derived properties to quantities relevant for Ly$\alpha$ escape. We observed LARS galaxies with Herschel/PACS, SOFIA/FIFI-LS, the IRAM 30m telescope and APEX, targeting far-infrared (FIR) continuum and emission lines of [C II]158$\mu$m, [O I]63$\mu$m, [O III]88$\mu$m and low-J CO lines. Using Bayesian methods we derive dust model parameters and estimate total gas masses for all LARS galaxies, taking into account a metallicity-dependent gas-to-dust ratio. Star formation rates were estimated from FIR, [C II]158$\mu$m and [O I]63$\mu$m luminosities. LARS covers a wide dynamic range in the derived properties, with FIR-based star formation rates from $\sim$0.5-100 $M_{\odot}\ yr^{-1}$, gas fractions between $\sim$15-80% and gas depletion times ranging from a few hundred Myr up to more than 10 Gyr. The distribution of LARS galaxies in the $\Sigma_{gas}$ vs. $\Sigma_{SFR}$ (Kennicutt-Schmidt plane) is thus quite heterogeneous. However, we find that LARS galaxies with the longest gas depletion times, i.e. relatively high gas surface densities ($\Sigma_{gas}$) and low star formation rate densities ($\Sigma_{SFR}$), have by far the highest Ly$\alpha$ escape fraction. A strong $\sim$linear relation is found between Ly$\alpha$ escape fraction and the total gas (HI+H$_2$) depletion time. We argue that the Ly$\alpha$ escape in those galaxies is driven by turbulence in the star-forming gas that shifts the Ly$\alpha$ photons out of resonance close to the places where they originate. We further report on an extreme [C II]158$\mu$m excess in LARS 5, corresponding to $\sim$14$\pm$3% of the FIR luminosity, i.e. the most extreme [C II]-to-FIR ratio observed in a non-AGN galaxy to date.

Dust lanes bisect the plane of a typical edge-on spiral galaxy as a dark optical absorption feature. Their appearance is linked to the gravitational stability of spiral disks; the fraction of edge-on galaxies that displays a dust lane is a direct indicator of the typical vertical balance between gravity and turbulence; a balance struck between the energy input from star-formation and the gravitational pull into the plane of the disk. Based on morphological classifications by the Galaxy~Zoo project on the Kilo-Degree Survey (KiDS) imaging data in the Galaxy and Mass Assembly (GAMA) fields, we explore the relation of dust lanes to the galaxy characteristics, most of which were determined using the {\sc magphys} spectral energy distribution fitting tool: stellar mass, total and specific star-formation rates, and several parameters describing the cold dust component. We find that the fraction of dust lanes does depend on the stellar mass of the galaxy; they start to appear at $M^* \sim 10^9 M_\odot$. A dust lane also implies strongly a dust mass of at least $10^5 M_\odot$, but otherwise does not correlate with cold dust mass parameters of the {\sc magphys} spectral energy distribution analysis, nor is there a link with star-formation rate, specific or total. Dust lane identification does not depend on disk ellipticity (disk thickness) or Sersic profile but correlates with bulge morphology; a round bulge favors dust lane votes. The central component along the line of sight that produces the dust lane is not associated with either one of the components fit by {\sc magphys}, the cold diffuse component or the localized, heated component in HII regions, but a mix of these two.

We assemble an unbiased sample of 29 galaxies with [O II] $\lambda 3727$ and/or [O III] $\lambda 5007$ detections at $z < 0.15$ from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Pilot Survey (HPS). HPS finds galaxies without pre-selection based on their detected emission lines via integral field spectroscopy. Sixteen of these objects were followed up with the second-generation, low resolution spectrograph (LRS2) on the upgraded Hobby-Eberly Telescope. Oxygen abundances were then derived via strong emission lines using a Bayesian approach. We find most of the galaxies fall along the mass-metallicity relation derived from photometrically selected star forming galaxies in the Sloan Digital Sky Survey (SDSS). However, two of these galaxies have low metallicity (similar to the very rare green pea galaxies in mass-metallicity space). The star formation rates of this sample fall in an intermediate space between the SDSS star forming main sequence and the extreme green pea galaxies. We conclude that spectroscopic selection fills part of the mass-metallicity-SFR phase space that is missed in photometric surveys with pre-selection like SDSS, i.e., we find galaxies that are actively forming stars but are faint in continuum. We use the results of this pilot investigation to make predictions for the upcoming unbiased, large spectroscopic sample of local line emitters from HETDEX. With the larger HETDEX survey we will determine if galaxies selected spectroscopically without continuum brightness pre-selection have metallicities that fall on a continuum that bridges typical star forming and rarer, more extreme systems like green peas.

Support for early career astronomers who are just beginning to explore astronomy research is imperative to increase retention of diverse practitioners in the field. Since 2010, Astrobites has played an instrumental role in engaging members of the community -- particularly undergraduate and graduate students -- in research. In this white paper, the Astrobites collaboration outlines our multi-faceted online education platform that both eases the transition into astronomy research and promotes inclusive professional development opportunities. We additionally offer recommendations for how the astronomy community can reduce barriers to entry to astronomy research in the coming decade.

The James Webb Space Telescope (JWST) will observe several stars for long cumulative durations while pursuing exoplanets as primary science targets for both Guaranteed Time Observations (GTO) and very likely General Observer (GO) programs. Here we argue in favor of an automatic default parallels program to observe e.g., using the F200W/F277W filters or grism of NIRCAM/NIRISS in order to find high redshift (z >> 10) galaxies, cool red/brown dwarf-sub-stellar objects, Solar System objects, and observations of serendipitous planetary transits. We argue here the need for automated exploratory astrophysical observations with unused JWST instruments during these long duration exoplanet observations. Randomized fields that are observed in parallel mode reduce errors due to cosmic variance more effectively than single continuous fields of a typical wedding cake observing strategy (Trenti & Stiavelli 2008). Hence, we argue that the proposed automated survey will explore a unique and rich discovery space in high redshift Universe, Galactic structure, and Solar System. We show that the GTO and highly-probable GO target list of exoplanets covers the Galactic disk/halo and high redshift Universe, mostly well out of the plane of the disk of the Milky Way. Exposure times are of the order of the CEERS GTO medium deep survey in a single filter, comparable to CANDELS in HST's surveys and deep fields. The area covered by NIRISS and NIRCam combined could accumulate to a half square degree surveyed.

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

We present the Super Eight galaxies - a set of very luminous, high-redshift ($7.1<z<8.0$) galaxy candidates found in Brightest of Reionizing Galaxies (BoRG) Survey fields. The original sample includes eight galaxies that are $Y$-band dropout objects with $H$-band magnitudes of $m_H<25.5$. Four of these objects were originally reported in Calvi et al. 2016. Combining new Hubble Space Telescope (HST) WFC3/F814W imaging and $Spitzer$ IRAC data with archival imaging from BoRG and other surveys, we explore the properties of these galaxies. Photometric redshift fitting places six of these galaxies in the redshift range of $7.1<z<8.0$, resulting in three new high-redshift galaxies and confirming three of the four high-redshift galaxy candidates from Calvi et al. 2016. We calculate the half-light radii of the Super Eight galaxies using the HST F160W filter and find that the Super Eight sizes are in line with typical evolution of size with redshift. The Super Eights have a mean mass of log(M$_*$/M$_\odot$) $\sim10$, which is typical for sources in this luminosity range. Finally, we place our sample on the UV $z\sim8$ luminosity function and find that the Super Eight number density is consistent with other surveys in this magnitude and redshift range.

We study young star-forming clumps on physical scales of 10-500 pc in the Lyman-Alpha Reference Sample (LARS), a collection of low-redshift (z = 0.03-0.2) UV-selected star-forming galaxies. In each of the 14 galaxies of the sample, we detect clumps for which we derive sizes and magnitudes in 5 UV-optical filters. The final sample includes $\sim$1400 clumps, of which $\sim$600 have magnitude uncertainties below 0.3 in all filters. The UV luminosity function for the total sample of clumps is described by a power-law with slope $\alpha = -2.03^{+0.11}_{-0.13}$. Clumps in the LARS galaxies have on average $\Sigma_{SFR}$ values higher than what observed in HII regions of local galaxies and comparable to typical SFR densities of clumps in z = 1-3 galaxies. We derive the clumpiness as the relative contribution from clumps to the UV emission of each galaxy, and study it as a function of galactic-scale properties, i.e. $\Sigma_{SFR}$ and the ratio between rotational and dispersion velocities of the gas ($v_s/\sigma_0$). We find that in galaxies with higher $\Sigma_{SFR}$ or lower $v_s/\sigma_0$, clumps dominate the UV emission of their host systems. All LARS galaxies with Ly$\alpha$ escape fractions larger than 10% have more than 50% of the UV luminosity from clumps. We tested the robustness of these results against the effect of different physical resolutions. At low resolution, the measured clumpiness appears more elevated than if we could resolve clumps down to single clusters. This effect is small in the redshift range covered by LARS, thus our results are not driven by the physical resolution.

The galaxy catalogs generated from low-resolution emission line surveys often contain both foreground and background interlopers due to line misidentification, which can bias the cosmological parameter estimation. In this paper, we present a method for correcting the interloper bias by using the joint-analysis of auto- and cross-power spectra of the main and the interloper samples. In particular, we can measure the interloper fractions from the cross-correlation between the interlopers and survey galaxies, because the true cross-correlation must be negligibly small. The estimated interloper fractions, in turn, remove the interloper bias in the cosmological parameter estimation. For example, in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) low-redshift ($z<0.5$) [O II] $\lambda3727${\AA} emitters contaminate high-redshift ($1.9<z<3.5$) Lyman-$\alpha$ line emitters. We demonstrate that the joint-analysis method yields a high signal-to-noise ratio measurement of the interloper fractions while only marginally increasing the uncertainties in the cosmological parameters relative to the case without interlopers. We also show the same is true for the high-latitude spectroscopic survey of Wide-Field Infrared Survey Telescope (WFIRST) mission where contamination occurs between the Balmer-$\alpha$ line emitters at lower redshifts ($1.1<z<1.9$) and Oxygen ([O III] $\lambda5007${\AA}) line emitters at higher redshifts ($1.7<z<2.8$).

U-band observations with the LSST have yet to be fully optimized in cadence. The straw man survey design is a simple coverage of the medium-deep-fast survey. Here we argue that deep coverage of the four deep drilling fields (XMM-LSS, ECDFS, ELAIS-S1 and COSMOS) has a much higher scientific return, given that these are also the target of the Southern Hemisphere's Square Kilometer Array Pathfinder, the MeerKAT specifically, deep radio observations.

We present a full data analysis of the pure-parallel Hubble Space Telescope (HST) imaging observations in the Brightest of Reionizing Galaxies Survey (BoRG[z9]) in Cycle 22. The medium-deep exposures with five HST/WFC3IR+UVIS filter bands from 79 independent sightlines (~370 arcmin^2) provide the least biased determination of number density for z>9 bright galaxies against cosmic variance. After a strict two-step selection for candidate galaxies, including dropout color and photometric redshift analyses, and revision of previous BoRG candidates, we identify one source at z~10 and two sources at z~9. The z~10 candidate shows evidence of line-of-sight lens magnification (mu~1.5), yet it appears surprisingly luminous (MUV ~ -22.6\pm0.3 mag), making it one of the brightest candidates at z > 8 known (~ 0.3 mag brighter than the z = 8.68 galaxy EGSY8p7, spectroscopically confirmed by Zitrin and collaborators). For z ~ 9 candidates, we include previous data points at fainter magnitudes and find that the data are well fitted by a Schechter luminosity function with alpha ~ -2.1, MUV ~ -21.5 mag, and log phi ~ -4.5 Mpc^-3mag^-1, for the first time without fixing any parameters. The inferred cosmic star formation rate density is consistent with unaccelerated evolution from lower redshift.

We use deep \textit{Hubble Space Telescope} spectroscopy to constrain the metallicities and (\editone{light-weighted}) ages of massive ($\log M_\ast/M_\odot\gtrsim10$) galaxies selected to have quiescent stellar populations at $1.0<z<1.8$. The data include 12--orbit depth coverage with the WFC3/G102 grism covering $\sim$ $8,000<\lambda<11,500$~\AA\, at a spectral resolution of $R\sim 210$ taken as part of the CANDELS Lyman-$\alpha$ Emission at Reionization (CLEAR) survey. At $1.0<z<1.8$, the spectra cover important stellar population features in the rest-frame optical. We simulate a suite of stellar population models at the grism resolution, fit these to the data for each galaxy, and derive posterior likelihood distributions for metallicity and age. We stack the posteriors for subgroups of galaxies in different redshift ranges that include different combinations of stellar absorption features. Our results give \editone{light-weighted ages of $t_{z \sim 1.1}= 3.2\pm 0.7$~Gyr, $t_{z \sim 1.2}= 2.2\pm 0.6$~Gyr, $t_{z\sim1.3}= 3.1\pm 0.6$~Gyr, and $t_{z\sim1.6}= 2.0 \pm 0.6$~Gyr, \editone{for galaxies at $z\sim 1.1$, 1.2, 1.3, and 1.6. This} implies that most of the massive quiescent galaxies at $1<z<1.8$ had formed $>68$\% of their stellar mass by a redshift of $z>2$}. The posteriors give metallicities of \editone{$Z_{z\sim1.1}=1.16 \pm 0.29$~$Z_\odot$, $Z_{z\sim1.2}=1.05 \pm 0.34$~$Z_\odot$, $Z_{z\sim1.3}=1.00 \pm 0.31$~$Z_\odot$, and $Z_{z\sim1.6}=0.95 \pm 0.39$~$Z_\odot$}. This is evidence that massive galaxies had enriched rapidly to approximately Solar metallicities as early as $z\sim3$.

AIMS: We aim to evaluate the near-infrared colors of brown dwarfs as observed with four major infrared imaging space observatories: the Hubble Space Telescope (HST), the James Webb Space Telescope (JWST), the Euclid mission, and the WFIRST telescope. METHODS: We used the SPLAT SPEX/ISPEX spectroscopic library to map out the colors of the M-, L-, and T-type dwarfs. We have identified which color-color combination is optimal for identifying broad type and which single color is optimal to then identify the subtype (e.g., T0-9). We evaluated each observatory separately as well as the narrow-field (HST and JWST) and wide-field (Euclid and WFIRST) combinations. RESULTS: The Euclid filters perform poorly typing brown dwarfs and WFIRST performs only marginally better, despite a wider selection of filters. WFIRST's W146 and F062 combined with Euclid's Y-band discriminates somewhat better between broad brown dwarf categories. However, subtyping with any combination of Euclid and WFIRST observations remains uncertain due to the lack of medium or narrow-band filters. We argue that a medium band added to the WFIRST filter selection would greatly improve its ability to preselect brown dwarfs its imaging surveys. CONCLUSIONS: The HST filters used in high-redshift searches are close to optimal to identify broad stellar type. However, the addition of F127M to the commonly used broad filter sets would allow for unambiguous subtyping. An improvement over HST is one of two broad and medium filter combinations on JWST: pairing F140M with either F150W or F162M discriminates very well between subtypes.

The Lyman Alpha Reference Sample (LARS) of 14 star-forming galaxies offers a wealth of insight into the workings of these local analogs to high-redshift star-forming galaxies. The sample has been well-studied in terms of LyA and other emission line properties, such as HI mass, gas kinematics, and morphology. We analyze deep surface photometry of the LARS sample in UBIK broadband imaging obtained at the Nordic Optical Telescope and the Canada-France-Hawaii Telescope, and juxtaposition their derived properties with a sample of local high-redshift galaxy analogs, namely, with blue compact galaxies (BCGs). We construct radial surface brightness and color profiles with both elliptical and isophotal integration, as well as RGB images, deep contours, color maps, a burst fraction estimate, and a radial mass-to-light ratio profile for each LARS galaxy. Standard morphological parameters like asymmetry, clumpiness, the Gini and M20 coefficients are [...] analyzed, as well as isophotal asymmetry profiles for each galaxy. [...] We compare the LARS to the properties of the BCG sample and highlight the differences. Several diagnostics indicate that the LARS galaxies have highly disturbed morphologies even at the level of the faintest isophotes [...]. The ground-based photometry [...] reveals previously unexplored isophotes [...]. The burst fraction estimate suggests a spatially more extended burst region in LARS than in the BCGs. [...] The galaxies in the LARS sample appear to be in earlier stages of a merger event compared to the BCGs. Standard morphological diagnostics like asymmetry, clumpiness, Gini and M20 coefficients cannot separate the two samples, although an isophotal asymmetry profile successfully captures the average difference in morphology. These morphological diagnostics do not show any correlation with the equivalent width or the escape fraction of Lyman Alpha. [abridged]

We examine the dust geometry and Ly{\alpha} scattering in the galaxies of the Lyman Alpha Reference Sample (LARS), a set of 14 nearby (0.02 < $z$ < 0.2) Ly{\alpha} emitting and starbursting systems with Hubble Space Telescope Ly{\alpha}, H{\alpha}, and H{\beta} imaging. We find that the global dust properties determined by line ratios are consistent with other studies, with some of the LARS galaxies exhibiting clumpy dust media while others of them show significantly lower Ly{\alpha} emission compared to their Balmer decrement. With the LARS imaging, we present Ly{\alpha}/H{\alpha} and H{\alpha}/H{\beta} maps with spatial resolutions as low as $\sim$ 40 pc, and use these data to show that in most galaxies, the dust geometry is best modeled by three distinct regions: a central core where dust acts as a screen, an annulus where dust is distributed in clumps, and an outer envelope where Ly{\alpha} photons only scatter. We show that the dust that affects the escape of Ly{\alpha} is more restricted to the galaxies' central regions, while the larger Ly{\alpha} halos are generated by scattering at large radii. We present an empirical modeling technique to quantify how much Ly{\alpha} scatters in the halo, and find that this "characteristic" scattering distance correlates with the measured size of the Ly{\alpha} halo. We note that there exists a slight anti-correlation between the scattering distance of Ly{\alpha} and global dust properties.

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.

We compare the rest-frame ultraviolet and rest-frame optical morphologies of 2 < z < 3 star-forming galaxies in the GOODS-S field using Hubble Space Telescope WFC3 and ACS images from the CANDELS, GOODS, and ERS programs. We show that the distribution of sizes and concentrations for 1.90 < z < 2.35 galaxies selected via their rest-frame optical emission-lines are statistically indistinguishable from those of Lyman-alpha emitting systems found at z ~ 2.1 and z ~ 3.1. We also show that the z > 2 star-forming systems of all sizes and masses become smaller and more compact as one shifts the observing window from the UV to the optical. We argue that this offset is due to inside-out galaxy formation over the first ~ 2 Gyr of cosmic time.

Differentiating between active galactic nuclei (AGN) activity and star formation in z ~ 2 galaxies is difficult because traditional methods, such as line ratio diagnostics, change with redshift while multi-wavelength methods (X-ray, radio, IR) are sensitive to only the brightest AGN. We have developed a new method for spatially resolving emission lines in HST/WFC3 G141 grism spectra and quantifying AGN activity through the spatial gradient of the [O III]/H$\beta$ line ratio. Through detailed simulations, we show that our novel line-ratio gradient approach identifies ~ sim 40% more low-mass and obscured AGN than obtained by classical methods. Based on our simulations, we developed a relationship that maps stellar mass, star formation rate, and measured [O III]/H$\beta$ gradient to AGN Eddington ratio. We apply our technique to previously studied stacked samples of galaxies at z ~2 and find that our results are consistent with these studies. Using this gradient method will also be able to inform other galaxy evolution science, such as inside-out quenching and metallicity gradients, and will be widely applicable to future spatially resolved JWST data.

We present a Bayesian approach to the redshift classification of emission-line galaxies when only a single emission line is detected spectroscopically. We consider the case of surveys for high-redshift Lyman-alpha-emitting galaxies (LAEs), which have traditionally been classified via an inferred rest-frame equivalent width (EW) greater than 20 angstrom. Our Bayesian method relies on known prior probabilities in measured emission-line luminosity functions and equivalent width distributions for the galaxy populations, and returns the probability that an object in question is an LAE given the characteristics observed. This approach will be directly relevant for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which seeks to classify ~10^6 emission-line galaxies into LAEs and low-redshift [O II] emitters. For a simulated HETDEX catalog with realistic measurement noise, our Bayesian method recovers 86% of LAEs missed by the traditional EW > 20 angstrom cutoff over 2 < z < 3, outperforming the EW cut in both contamination and incompleteness. This is due to the method's ability to trade off between the two types of binary classification error by adjusting the stringency of the probability requirement for classifying an observed object as an LAE. In our simulations of HETDEX, this method reduces the uncertainty in cosmological distance measurements by 14% with respect to the EW cut, equivalent to recovering 29% more cosmological information. Rather than using binary object labels, this method enables the use of classification probabilities in large-scale structure analyses. It can be applied to narrowband emission-line surveys as well as upcoming large spectroscopic surveys including Euclid and WFIRST.

We compare the physical and morphological properties of z ~ 2 Lyman-alpha emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ~ 2 optical emission line galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5 < logM < 10.5), size (0.5 < R < 3.0 kpc), and star-formation rate (~1 < SFR < 100). Remarkably, a comparison of the most commonly used physical and morphological parameters -- stellar mass, half-light radius, UV slope, star formation rate, ellipticity, nearest neighbor distance, star formation surface density, specific star formation rate, [O III] luminosity, and [O III] equivalent width -- reveals no statistically significant differences between the populations. This suggests that the processes and conditions which regulate the escape of Ly-alpha from a z ~ 2 star-forming galaxy do not depend on these quantities. In particular, the lack of dependence on the UV slope suggests that Ly-alpha emission is not being significantly modulated by diffuse dust in the interstellar medium. We develop a simple model of Ly-alpha emission that connects LAEs to all high-redshift star forming galaxies where the escape of Ly-alpha depends on the sightline through the galaxy. Using this model, we find that mean solid angle for Ly-alpha escape is 2.4+/-0.8 steradians; this value is consistent with those calculated from other studies.

We explore the evolution of the Stellar Mass-Star Formation Rate-Metallicity Relation using a set of 256 COSMOS and GOODS galaxies in the redshift range 1.90 < z < 2.35. We present the galaxies' rest-frame optical emission-line fluxes derived from IR-grism spectroscopy with the Hubble Space Telescope and combine these data with star formation rates and stellar masses obtained from deep, multi-wavelength (rest-frame UV to IR) photometry. We then compare these measurements to those for a local sample of galaxies carefully matched in stellar mass (7.5 < log(M*/Msol) < 10.5) and star formation rate (-0.5 < log(SFR) < 2.5 in Msol yr^-1). We find that the distribution of z ~ 2.1 galaxies in stellar mass-SFR-metallicity space is clearly different from that derived for our sample of similarly bright (L_H\b{eta} > 3 . 10^40 ergs s^-1) local galaxies, and this offset cannot be explained by simple systematic offsets in the derived quantities. At stellar masses above ~10^9 Msol and star formation rates above ~10 Msol yr^-1, the z ~ 2.1 galaxies have higher oxygen abundances than their local counterparts, while the opposite is true for lower-mass, lower-SFR systems.

We derive the mean wavelength dependence of stellar attenuation in a sample of 239 high redshift (1.90 < z < 2.35) galaxies selected via Hubble Space Telescope (HST) WFC3 IR grism observations of their rest-frame optical emission lines. Our analysis indicates that the average reddening law follows a form similar to that derived by Calzetti et al. for local starburst galaxies. However, over the mass range 7.2 < log M/Msolar < 10.2, the slope of the attenuation law in the UV is shallower than that seen locally, and the UV slope steepens as the mass increases. These trends are in qualitative agreement with Kriek & Conroy, who found that the wavelength dependence of attenuation varies with galaxy spectral type. However, we find no evidence of an extinction "bump" at 2175 A in any of the three stellar mass bins, or in the sample as a whole. We quantify the relation between the attenuation curve and stellar mass and discuss its implications.

We use 317,000 emission-line galaxies from the Sloan Digital Sky Survey to investigate line-ratio selection of active galactic nuclei (AGNs). In particular, we demonstrate that "star formation dilution" by HII regions causes a significant bias against AGN selection in low-mass, blue, star-forming, disk-dominated galaxies. This bias is responsible for the observed preference of AGNs among high-mass, green, moderately star-forming, bulge-dominated hosts. We account for the bias and simulate the intrinsic population of emission-line AGNs using a physically-motivated Eddington ratio distribution, intrinsic AGN narrow line region line ratios, a luminosity-dependent Lbol/L[OIII] bolometric correction, and the observed Mbh-sigma relation. These simulations indicate that, in massive (log(M*/Msun) > 10) galaxies, AGN accretion is correlated with specific star formation rate but is otherwise uniform with stellar mass. There is some hint of lower black hole occupation in low-mass (log(M*/Msun) < 10) hosts, although our modeling is limited by uncertainties in measuring and interpreting the velocity dispersions of low-mass galaxies. The presence of star formation dilution means that AGNs contribute little to the observed strong optical emission lines (e.g., [OIII] and Ha) in low-mass and star-forming hosts. However the AGN population recovered by our modeling indicates that feedback by typical (low- to moderate-accretion) low-redshift AGNs has nearly uniform efficiency at all stellar masses, star formation rates, and morphologies. Taken together, our characterization of the observational bias and resultant AGN occupation function suggest that AGNs are unlikely to be the dominant source of star formation quenching in galaxies, but instead are fueled by the same gas which drives star formation activity.

The Hobby-Eberly Dark Energy Experiment pilot survey identified 284 [O II] 3727 emitting galaxies in a 169 square-arcminute field of sky in the redshift range 0 < z < 0.57. This line flux limited sample provides a bridge between studies in the local universe and higher-redshift [O II] surveys. We present an analysis of the star formation rates (SFRs) of these galaxies as a function of stellar mass as determined via spectral energy distribution fitting. The [O II] emitters fall on the "main sequence" of star-forming galaxies with SFR decreasing at lower masses and redshifts. However, the slope of our relation is flatter than that found for most other samples, a result of the metallicity dependence of the [O II] star formation rate indicator. The mass specific SFR is higher for lower mass objects, supporting the idea that massive galaxies formed more quickly and efficiently than their lower mass counterparts. This is confirmed by the fact that the equivalent widths of the [O II] emission lines trend smaller with larger stellar mass. Examination of the morphologies of the [O II] emitters reveals that their star formation is not a result of mergers, and the galaxies' half-light radii do not indicate evolution of physical sizes.

We use near-IR grism spectroscopy from the Hubble Space Telescope to examine the strength of [Ne~III] 3869 relative to H-beta, [O~II] 3727 and [O~III] 5007 in 236 low mass (7.5 < log (Mstar/Msolar) < 10.5) star-forming galaxies in the redshift range 1.90 < z < 2.35. By stacking the data by stellar mass, we show that the [Ne~III]/[O~II] ratios of the z ~ 2 universe are marginally higher than those seen in a comparable set of local SDSS galaxies, and that [Ne~III]/[O~III] is enhanced by ~0.2 dex. We consider the possible explanations for this ~4-sigma result, including higher oxygen depletion out of the gas-phase, denser H~II regions, higher production of Ne22 via Wolf-Rayet stars, and the existence of a larger population of X-ray obscured AGN at z ~ 2 compared to z ~ 0. None of these simple scenarios, alone, are favored to explain the observed line ratios. We conclude by suggesting several avenues of future observations to further explore the mystery of enhanced [Ne~III] emission.

We present the results from a VLT/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyman-alpha emitters (LAEs) at $z$ = 2.1 - 2.5 in the COSMOS and GOODS-N fields discovered from the HETDEX Pilot Survey. We detect rest-frame optical nebular lines (H$\alpha$ and/or [OIII]$\lambda$5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas-mass fraction, and Ly$\alpha$ velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Ly$\alpha$ relative to the systemic redshift ranging between +85 and +296 km s$^{-1}$ with a mean of +180 km s$^{-1}$. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Ly$\alpha$ velocity offsets show a moderate correlation with the measured star formation rate (2.5$\sigma$), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass ($\lesssim$ 1.5$\sigma$). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Ly$\alpha$ photons, we find no signature of selective quenching of resonantly scattered Ly$\alpha$ photons. However, we also find no evidence that a clumpy ISM is enhancing the Ly$\alpha$ equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density as well as less dust, easing the escape of Ly$\alpha$ photons over that in continuum-selected star-forming galaxies.

We use Hubble Space Telescope near-IR grism spectroscopy to examine the H-beta line strengths of 260 star-forming galaxies in the redshift range 1.90 < z < 2.35. We show that at these epochs, the H-beta star formation rate (SFR) is a factor of ~1.8 higher than what would be expected from the systems' rest-frame UV flux density, suggesting a shift in the standard conversion between these quantities and star formation rate. We demonstrate that at least part of this shift can be attributed to metallicity, as H-beta is more greatly enhanced in systems with lower oxygen abundance. This offset must be considered when measuring the star formation rate history of the universe. We also show that the relation between stellar and nebular extinction in our z ~ 2 sample is consistent with that observed in the local universe.

We use broadband photometry extending from the rest-frame UV to the near-IR to fit the individual spectral energy distributions (SEDs) of 63 bright (L(Ly-alpha) > 10^43 ergs/s) Ly-alpha emitting galaxies (LAEs) in the redshift range 1.9 < z < 3.6. We find that these LAEs are quite heterogeneous, with stellar masses that span over three orders of magnitude, from 7.5 < log M < 10.5. Moreover, although most LAEs have small amounts of extinction, some high-mass objects have stellar reddenings as large as E(B-V) ~0.4. Interestingly, in dusty objects the optical depths for Ly-alpha and the UV continuum are always similar, indicating that Ly-alpha photons are not undergoing many scatters before escaping their galaxy. In contrast, the ratio of optical depths in low-reddening systems can vary widely, illustrating the diverse nature of the systems. Finally, we show that in the star formation rate (SFR)-log mass diagram, our LAEs fall above the "main-sequence" defined by z ~ 3 continuum selected star-forming galaxies. In this respect, they are similar to sub-mm-selected galaxies, although most LAEs have much lower mass.