Recent studies show that a universal relation between black-hole (BH) growth and stellar mass ($M_\bigstar$) or star formation rate (SFR) is an oversimplification of BH-galaxy co-evolution, and that morphological and structural properties of host galaxies must also be considered. Particularly, a possible connection between BH growth and host-galaxy compactness was identified among star-forming (SF) galaxies. Utilizing $\approx 6300$ massive galaxies with $I_{\rm 814W}~<~24$ at $z$ $<$ 1.2 in the COSMOS field, we perform systematic partial-correlation analyses to investigate how sample-averaged BH accretion rate ($\rm \overline{BHAR}$) depends on host-galaxy compactness among SF galaxies, when controlling for morphology and $M_\bigstar$ (or SFR). The projected central surface-mass density within 1 kpc, $\Sigma_{1}$, is utilized to represent host-galaxy compactness in our study. We find that the $\rm \overline{BHAR}$-$\Sigma_{1}$ relation is stronger than either the $\rm \overline{BHAR}$-$M_\bigstar$ or $\rm \overline{BHAR}$-SFR relation among SF galaxies, and this $\rm \overline{BHAR}$-$\Sigma_{1}$ relation applies to both bulge-dominated galaxies and galaxies that are not dominated by bulges. This $\rm \overline{BHAR}$-$\Sigma_{1}$ relation among SF galaxies suggests a link between BH growth and the central gas density of host galaxies on the kpc scale, which may further imply a common origin of the gas in the vicinity of the BH and in the central $\sim$ kpc of the galaxy. This $\rm \overline{BHAR}$-$\Sigma_{1}$ relation can also be interpreted as the relation between BH growth and the central velocity dispersion of host galaxies at a given gas content, indicating the role of the host-galaxy potential well in feeding BHs.

We present the discovery of the optical afterglow and host galaxy of the {\it Swift} short-duration gamma-ray burst, GRB\,181123B. Observations with Gemini-North starting at $\approx 9.1$~hr after the burst reveal a faint optical afterglow with $i\approx25.1$~mag, at an angular offset of 0.59 $\pm$ 0.16$''$ from its host galaxy. Using $grizYJHK$ observations, we measure a photometric redshift of the host galaxy of $z = 1.77^{+0.30}_{-0.17}$. From a combination of Gemini and Keck spectroscopy of the host galaxy spanning 4500-18000~\AA , we detect a single emission line at 13390~\AA, inferred as H$\beta$ at $z = 1.754 \pm 0.001$ and corroborating the photometric redshift. The host galaxy properties of GRB\,181123B are typical to those of other SGRB hosts, with an inferred stellar mass of $\approx 1.7 \times 10^{10}\,M_{\odot}$, mass-weighted age of $\approx 0.9$~Gyr and optical luminosity of $\approx 0.9L^{*}$. At $z=1.754$, GRB\,181123B is the most distant secure SGRB with an optical afterglow detection, and one of only three at $z>1.5$. Motivated by a growing number of high-$z$ SGRBs, we explore the effects of a missing $z>1.5$ SGRB population among the current {\it Swift} sample on delay time distribution models. We find that log-normal models with mean delay times of $\approx 4-6$~Gyr are consistent with the observed distribution, but can be ruled out to $95\%$ confidence with an additional $\approx1-5$~{\it Swift} SGRBs recovered at $z>1.5$. In contrast, power-law models with $\propto$ $t^{-1}$ are consistent with the redshift distribution and can accommodate up to $\approx30$ SGRBs at these redshifts. Under this model, we predict that $\approx 1/3$ of the current {\it Swift} population of SGRBs is at $z>1$. The future discovery or recovery of existing high-$z$ SGRBs will provide significant discriminating power on their delay time distributions, and thus their formation channels.

The primary method for inferring the stellar mass ($M_*$) of a galaxy is through spectral energy distribution (SED) modeling. However, the technique rests on assumptions such as the galaxy star formation history and dust attenuation law that can severely impact the accuracy of derived physical properties from SED modeling. Here, we examine the effect that the assumed star formation history (SFH) has on the stellar properties inferred from SED fitting by ground truthing them against mock observations of high-resolution cosmological hydrodynamic galaxy formation simulations. Classically, SFHs are modeled with simplified parameterized functional forms, but these forms are unlikely to capture the true diversity of galaxy SFHs and may impose systematic biases with under-reported uncertainties on results. We demonstrate that flexible nonparametric star formation histories outperform traditional parametric forms in capturing variations in galaxy star formation histories, and as a result, lead to significantly improved stellar masses in SED fitting. We find a decrease in the average bias of 0.4 dex with a delayed-$\tau$ model to a bias of just under 0.05 dex for the nonparametric model. Similarly, using nonparametric star formation histories in SED fitting result in increased accuracy in recovered galaxy star formation rates (SFRs) and stellar ages.

We present SPECULATOR - a fast, accurate, and flexible framework for emulating stellar population synthesis (SPS) models for predicting galaxy spectra and photometry. For emulating spectra, we use principal component analysis to construct a set of basis functions, and neural networks to learn the basis coefficients as a function of the SPS model parameters. For photometry, we parameterize the magnitudes (for the filters of interest) as a function of SPS parameters by a neural network. The resulting emulators are able to predict spectra and photometry under both simple and complicated SPS model parameterizations to percent-level accuracy, giving a factor of $10^3$-$10^4$ speed up over direct SPS computation. They have readily-computable derivatives, making them amenable to gradient-based inference and optimization methods. The emulators are also straightforward to call from a GPU, giving an additional order-of-magnitude speed-up. Rapid SPS computations delivered by emulation offers a massive reduction in the computational resources required to infer the physical properties of galaxies from observed spectra or photometry and simulate galaxy populations under SPS models, whilst maintaining the accuracy required for a range of applications.

There has been a long-standing factor-of-two tension between the observed star formation rate density and the observed stellar mass buildup after $z\sim2$. Recently we have proposed that sophisticated panchromatic SED models can resolve this tension, as these methods infer systematically higher masses and lower star formation rates than standard approaches. In a series of papers we now extend this analysis and present a complete, self-consistent census of galaxy formation over $0.2 < z < 3$ inferred with the \texttt{Prospector} galaxy SED-fitting code. In this work, Paper I, we present the evolution of the galaxy stellar mass function using new mass measurements of $\sim$10$^5$ galaxies in the 3D-HST and COSMOS-2015 surveys. We employ a new methodology to infer the mass function from the observed stellar masses: instead of fitting independent mass functions in a series of fixed redshift intervals, we construct a continuity model that directly fits for the redshift evolution of the mass function. This approach ensures a smooth picture of galaxy assembly and makes use of the full, non-Gaussian uncertainty contours in our stellar mass inferences. The resulting mass function has higher number densities at a fixed stellar mass than almost any other measurement in the literature, largely owing to the older stellar ages inferred by \texttt{Prospector}. The stellar mass density is $\sim$50% higher than previous measurements, with the offset peaking at $z\sim1$. The next two papers in this series will present the new measurements of star-forming main sequence and the cosmic star formation rate density, respectively.

Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST.

We present BVRI and unfiltered light curves of 93 Type Ia supernovae (SNe Ia) from the Lick Observatory Supernova Search (LOSS) follow-up program conducted between 2005 and 2018. Our sample consists of 78 spectroscopically normal SNe Ia, with the remainder divided between distinct subclasses (three SN 1991bg-like, three SN 1991T-like, four SNe Iax, two peculiar, and three super-Chandrasekhar events), and has a median redshift of 0.0192. The SNe in our sample have a median coverage of 16 photometric epochs at a cadence of 5.4 days, and the median first observed epoch is ~4.6 days before maximum B-band light. We describe how the SNe in our sample are discovered, observed, and processed, and we compare the results from our newly developed automated photometry pipeline to those from the previous processing pipeline used by LOSS. After investigating potential biases, we derive a final systematic uncertainty of 0.03 mag in BVRI for our dataset. We perform an analysis of our light curves with particular focus on using template fitting to measure the parameters that are useful in standardising SNe Ia as distance indicators. All of the data are available to the community, and we encourage future studies to incorporate our light curves in their analyses.

We report the serendipitous detection of two 3 mm continuum sources found in deep ALMA Band 3 observations to study intermediate redshift galaxies in the COSMOS field. One is near a foreground galaxy at 1.3", but is a previously unknown dust-obscured star-forming galaxy (DSFG) at probable $z_{CO}=3.329$, illustrating the risk of misidentifying shorter wavelength counterparts. The optical-to-mm spectral energy distribution (SED) favors a grey $\lambda^{-0.4}$ attenuation curve and results in significantly larger stellar mass and SFR compared to a Calzetti starburst law, suggesting caution when relating progenitors and descendants based on these quantities. The other source is missing from all previous optical/near-infrared/sub-mm/radio catalogs ("ALMA-only"), and remains undetected even in stacked ultradeep optical ($>29.6$ AB) and near-infrared ($>27.9$ AB) images. Using the ALMA position as a prior reveals faint $SNR\sim3$ measurements in stacked IRAC 3.6+4.5, ultradeep SCUBA2 850$\mu$m, and VLA 3GHz, indicating the source is real. The SED is robustly reproduced by a massive $M^*=10^{10.8}$M$_\odot$ and $M_{gas}=10^{11}$M$_\odot$, highly obscured $A_V\sim4$, star forming $SFR\sim300$ M$_{\odot}$yr$^{-1}$ galaxy at redshift $z=5.5\pm$1.1. The ultrasmall 8 arcmin$^{2}$ survey area implies a large yet uncertain contribution to the cosmic star formation rate density CSFRD(z=5) $\sim0.9\times10^{-2}$ M$_{\odot}$ yr$^{-1}$ Mpc$^{-3}$, comparable to all ultraviolet-selected galaxies combined. These results indicate the existence of a prominent population of DSFGs at $z>4$, below the typical detection limit of bright galaxies found in single-dish sub-mm surveys, but with larger space densities $\sim3 \times 10^{-5}$ Mpc$^{-3}$, higher duty cycles $50-100\%$, contributing more to the CSFRD, and potentially dominating the high-mass galaxy stellar mass function.

This manuscript describes the public release of the Hubble Legacy Fields (HLF) project photometric catalog for the extended GOODS-South region from the Hubble Space Telescope (HST) archival program AR-13252. The analysis is based on the version 2.0 HLF data release that now includes all ultraviolet (UV) imaging, combining three major UV surveys. The HLF data combines over a decade worth of 7475 exposures taken in 2635 orbits totaling 6.3 Msec with the HST Advanced Camera for Surveys Wide Field Channel (ACS/WFC) and the Wide Field Camera 3 (WFC3) UVIS/IR Channels in the greater GOODS-S extragalactic field, covering all major observational efforts (e.g., GOODS, GEMS, CANDELS, ERS, UVUDF and many other programs; see Illingworth et al 2019, in prep). The HLF GOODS-S catalogs include photometry in 13 bandpasses from the UV (WFC3/UVIS F225W, F275W and F336W filters), optical (ACS/WFC F435W, F606W, F775W, F814W and F850LP filters), to near-infrared (WFC3/IR F098M, F105W, F125W, F140W and F160W filters). Such a data set makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range from high resolution mosaics that are largely contiguous. Here, we describe a photometric analysis of 186,474 objects in the HST imaging at wavelengths 0.2--1.6$\mu$m. We detect objects from an ultra-deep image combining the PSF-homogenized and noise-equalized F850LP, F125W, F140W and F160W images, including Gaia astrometric corrections. SEDs were determined by carefully taking the effects of the point-spread function in each observation into account. All of the data presented herein are available through the HLF website (https://archive.stsci.edu/prepds/hlf/).

We present and analyse a new tidal disruption event (TDE), AT2017eqx at redshift z=0.1089, discovered by Pan-STARRS and ATLAS. The position of the transient is consistent with the nucleus of its host galaxy; it peaks at a luminosity of $L \approx 10^{44}$ erg s$^{-1}$; and the spectrum shows a persistent blackbody temperature $T \gtrsim 20,000$ K with broad H I and He II emission. The lines are initially centered at zero velocity, but by 100 days the H I lines disappear while the He II develops a blueshift of $\gtrsim 5,000$ km s$^{-1}$. Both the early- and late-time morphologies have been seen in other TDEs, but the complete transition between them is unprecedented. The evolution can be explained by combining an extended atmosphere, undergoing slow contraction, with a wind in the polar direction becoming visible at late times. Our observations confirm that a lack of hydrogen a TDE spectrum does not indicate a stripped star, while the proposed model implies that much of the diversity in TDEs may be due to the observer viewing angle. Modelling the light curve suggests AT2017eqx resulted from the complete disruption of a solar-mass star by a black hole of $\sim 10^{6.3} M_\odot$. The host is another quiescent, Balmer-strong galaxy, though fainter and less centrally concentrated than most TDE hosts. Radio limits rule out a relativistic jet, while X-ray limits at 500 days are among the deepest for a TDE at this phase.

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$.

The UVJ color-color diagram is a popular and efficient method to distinguish between quiescent and star-forming galaxies through their rest-frame $U-V$ vs. $V-J$ colors. Here we explore the information content of this color-color space using the Bayesian inference machine Prospector. We fit the same physical model to two datasets: (i) UVJ fluxes alone, and (ii) full UV-mid IR (MIR) broadband SEDs from the 3D-HST survey. Notably this model uses both nonparametric SFHs and a flexible dust attenuation curve, both of which have the potential to `break' the typical correlations observed in UVJ color-color space. Instead, these fits confirm observed trends between UVJ colors and observed galaxy properties, including specific star formation rate (sSFR), dust attenuation, stellar age, and stellar metallicity. They also demonstrate that UVJ colors do not, on their own, constrain stellar age or metallicity; the observed trends in the UVJ diagram are instead driven by galaxy scaling relationships and thus will evolve with cosmological time. We also show that UVJ colors 'saturate' below $\log(\mathrm{sSFR/yr}^{-1})\lesssim -10.5$, i.e. changing sSFR no longer produces substantial changes in UVJ colors. We show that far-UV and/or MIR fluxes continue to correlate with sSFR down to low sSFRs and can be used in color-color diagrams to efficiently target galaxies with much lower levels of ongoing star formation. We provide selection criteria in these new color-color spaces as a function of desired sample sSFR.

Galaxy observations are influenced by many physical parameters: stellar masses, star formation rates (SFRs), star formation histories (SFHs), metallicities, dust, black hole activity, and more. As a result, inferring accurate physical parameters requires high-dimensional models which capture or marginalize over this complexity. Here we re-assess inferences of galaxy stellar masses and SFRs using the 14-parameter physical model Prospector-$\alpha$ built in the Prospector Bayesian inference framework. We fit the photometry of 58,461 galaxies from the 3D-HST catalogs at $0.5 < z < 2.5$. The resulting stellar masses are $\sim0.1-0.3$ dex larger than the fiducial masses while remaining consistent with dynamical constraints. This change is primarily due to the systematically older SFHs inferred with Prospector. The SFRs are $\sim0.1-1+$ dex lower than UV+IR SFRs, with the largest offsets caused by emission from "old" ($t>100$ Myr) stars. These new inferences lower the observed cosmic star formation rate density by $\sim0.2$ dex and increase the observed stellar mass growth by $\sim 0.1$ dex, finally bringing these two quantities into agreement and implying an older, more quiescent Universe than found by previous studies at these redshifts. We corroborate these results by showing that the Prospector-$\alpha$ SFHs are both more physically realistic and are much better predictors of the evolution of the stellar mass function. Finally, we highlight examples of observational data which can break degeneracies in the current model; these observations can be incorporated into priors in future models to produce new & more accurate physical parameters.

In Paper I (Safarzadeh & Berger 2019) we studied the determination of the delay time distribution (DTD) of binary neutron stars (BNS) through scaling relations between halo/stellar mass and the star formation history (SFH) of galaxies hosting gravitational wave (GW) events in the local universe. Here we explore how a detailed reconstruction of the individual SFHs of BNS merger host galaxies can improve on the use of the scaling relations. We use galaxies from the Galaxy and Mass Assembly (GAMA) survey, which is mass complete at $M_*>10^9$ M$_\odot$ in the redshift range $0.05<z<0.08$. We use the reconstructed SFHs derived from the Prospector code, for two distinct sets of priors (favoring continuous and bursty SFHs), and convolve those with power law DTDs characterized by an index $\Gamma$ and a minimum delay time $t_{\rm min}$. We find that with this approach $\mathcal{O}(100)-\mathcal{O}(300)$ host galaxies are required to constrain the DTD parameters, with the number depending on the choice of SFH prior and on the parameters of the true DTD. We further show that using only the host galaxies of BNS mergers, as opposed to the full population of potential host galaxies in the relevant cosmic volume, leads to a minor bias in the recovered DTD parameters. The required host galaxy sample size is nearly an order of magnitude smaller relative to the approach of using scaling relations, and we expect such a host galaxy sample to be collected within a decade or two, prior to the advent of third-generation GW detectors.

Nonparametric star formation histories (SFHs) have long promised to be the `gold standard' for galaxy spectral energy distribution (SED) modeling as they are flexible enough to describe the full diversity of SFH shapes, whereas parametric models rule out a significant fraction of these shapes {\it a priori}. However, this flexibility is not fully constrained even with high-quality observations, making it critical to choose a well-motivated prior. Here, we use the SED-fitting code \texttt{Prospector} to explore the effect of different nonparametric priors by fitting SFHs to mock UV-IR photometry generated from a diverse set of input SFHs. First, we confirm that nonparametric SFHs recover input SFHs with less bias and return more accurate errors than do parametric SFHs. We further find that, while nonparametric SFHs robustly recover the overall shape of the input SFH, the primary determinant of the size and shape of the posterior star formation rate (SFR) as a function of time is the choice of prior, rather than the photometric noise. As a practical demonstration, we fit the UV-IR photometry of $\sim$6000 galaxies from the GAMA survey and measure inter-prior scatters in mass (0.1 dex), SFR$_{100\; \mathrm{Myr}}$ (0.8 dex), and mass-weighted ages (0.2 dex), with the bluest star-forming galaxies showing the most sensitivity. An important distinguishing characteristic for nonparametric models is the characteristic timescale for changes in SFR(t). This difference controls whether galaxies are assembled in bursts or in steady-state star formation, corresponding respectively to (feedback-dominated/accretion-dominated) models of galaxy formation and to (larger/smaller) confidence intervals derived from SED-fitting. High-quality spectroscopy has the potential to further distinguish between these proposed models of SFR(t).

We present optical photometry and spectroscopy of SN 2016iet, an unprecedented Type I supernova (SN) at $z=0.0676$ with no obvious analog in the existing literature. The peculiar light curve has two roughly equal brightness peaks ($\approx -19$ mag) separated by 100 days, and a subsequent slow decline by 5 mag in 650 rest-frame days. The spectra are dominated by emission lines of calcium and oxygen, with a width of only $3400$ km s$^{-1}$, superposed on a strong blue continuum in the first year, and with a large ratio of $L_{\rm [Ca\,II]}/L_{\rm [O\,I]}\approx 4$ at late times. There is no clear evidence for hydrogen or helium associated with the SN at any phase. We model the light curves with several potential energy sources: radioactive decay, central engine, and circumstellar medium (CSM) interaction. Regardless of the model, the inferred progenitor mass near the end of its life (i.e., CO core mass) is $\gtrsim 55$ M$_\odot$ and up to $120$ M$_\odot$, placing the event in the regime of pulsational pair instability supernovae (PPISNe) or pair instability supernovae (PISNe). The models of CSM interaction provide the most consistent explanation for the light curves and spectra, and require a CSM mass of $\approx 35$ M$_\odot$ ejected in the final decade before explosion. We further find that SN 2016iet is located at an unusually large offset ($16.5$ kpc) from its low metallicity dwarf host galaxy ($Z\approx 0.1$ Z$_\odot$, $M\approx 10^{8.5}$ M$_\odot$), supporting the PPISN/PISN interpretation. In the final spectrum, we detect narrow H$\alpha$ emission at the SN location, likely due to a dim underlying galaxy host or an H II region. Despite the overall consistency of the SN and its unusual environment with PPISNe and PISNe, we find that the inferred properties of SN\,2016iet challenge existing models of such events.

Parametric models for galaxy star-formation histories (SFHs) are widely used, though they are known to impose strong priors on physical parameters. This has consequences for measurements of the galaxy stellar-mass function (GSMF), star-formation-rate density (SFRD) and star-forming main sequence (SFMS). We investigate the effects of the exponentially declining, delayed exponentially declining, lognormal and double power law SFH models using BAGPIPES. We demonstrate that each of these models imposes strong priors on specific star-formation rates (sSFRs), potentially biasing the SFMS, and also imposes a strong prior preference for young stellar populations. We show that stellar mass, SFR and mass-weighted age inferences from high-quality mock photometry vary with the choice of SFH model by at least 0.1, 0.3 and 0.2 dex respectively. However the biases with respect to the true values depend more on the true SFH shape than the choice of model. We also demonstrate that photometric data cannot discriminate between SFH models, meaning it is important to perform independent tests to find well-motivated priors. We finally fit a low-redshift, volume-complete sample of galaxies from the Galaxy and Mass Assembly (GAMA) Survey with each model. We demonstrate that our stellar masses and SFRs at redshift, $z\sim0.05$ are consistent with other analyses. However, our inferred cosmic SFRDs peak at $z\sim0.4$, approximately 6 Gyr later than direct observations suggest, meaning our mass-weighted ages are significantly underestimated. This makes the use of parametric SFH models for understanding mass assembly in galaxies challenging. In a companion paper we consider non-parametric SFH models.

Using the \prospector\ spectral energy distribution (SED) fitting code, we analyze the properties of 19 Extreme Emission Line Galaxies (EELGs) identified in the bluest composite SED in the \zfourge\ survey at $2.5 \leq z \leq 4$. \prospector\ includes a physical model for nebular emission and returns probability distributions for stellar mass, stellar metallicity, dust attenuation, and nonparametric star formation history (SFH). The EELGs show evidence for a starburst in the most recent 50 Myr, with the median EELG having a specific star formation rate (sSFR) of 4.6 Gyr$^{-1}$ and forming 15\% of its mass in this short time. For a sample of more typical star-forming galaxies (SFGs) at the same redshifts, the median SFG has a sSFR of 1.1 Gyr$^{-1}$ and forms only $4\%$ of its mass in the last 50 Myr. We find that virtually all of our EELGs have rising SFHs, while most of our SFGs do not. From our analysis, we hypothesize that many, if not most, star-forming galaxies at $z \geq 2.5$ undergo an extreme H$\beta$+$[\hbox{{\rm O}\kern 0.1em{\sc iii}}]$ emission line phase early in their lifetimes. In a companion paper, we obtain spectroscopic confirmation of the EELGs as part of our {\sc MOSEL} survey. In the future, explorations of uncertainties in modeling the UV slope for galaxies at $z>2$ are needed to better constrain their properties, e.g. stellar metallicities.

We present COSMOS-Drift And SHift (DASH), a Hubble Space Telescope WFC3 imaging survey of the COSMOS field in the H_160 filter. The survey comprises 456 individual WFC3 pointings corresponding to an area of 0.49 deg^2 (0.66 deg^2 when including archival data) and reaches a 5 point-source limit of H_160 =25.1 (0".3 aperture). COSMOS-DASH is the widest HST/WFC3 imaging survey in H_160 filter, tripling the extragalactic survey area in the near-infrared at HST resolution. We make the reduced H_160 mosaic available to the community. We use this dataset to measure the sizes of 162 galaxies with log(M_star/M_sun) > 11.3 at 1.5 < z < 3.0, and augment this sample with 748 galaxies at 0.1 < z < 1.5 using archival ACS imaging. We find that the median size of galaxies in this mass range changes with redshift as r_eff = (10.4+/-0.4)(1 +z)^(0.65+/-0.05) kpc. Separating the galaxies into star forming and quiescent galaxies using their restframe U-V and V-J colors, we find no statistical difference between the median sizes of the most massive star-forming and quiescent galaxies at z = 2.5: they are 4.9+/-0.9 kpc and 4.3 +/-0.3 kpc respectively. However, we do find a significant difference in the S`ersic index between the two samples, such that massive quiescent galaxies have higher central densities than star forming galaxies. We extend the size-mass analysis to lower masses by combining it with the 3D-HST/CANDELS sample of van der Wel et al. (2014), and derive empirical relations between size, mass, and redshift. Fitting a relation of the form r_eff = A m_star^a, m_star = M_star/5x10^10 M_sun and r_eff in kpc, we find log A = -0.25 log (1 + z) + 0.79 and a = -0.13 log(1 + z) + 0.27. We also provide relations for the subsamples of star forming and quiescent galaxies. Our results confirm previous studies that were based on smaller samples or ground-based imaging.

Forward modeling of the full galaxy SED is a powerful technique, providing self-consistent constraints on stellar ages, dust properties, and metallicities. However, the accuracy of these results is contingent on the accuracy of the model. One significant source of uncertainty is the contribution of obscured AGN, as they are relatively common and can produce substantial mid-IR (MIR) emission. Here we include emission from dusty AGN torii in the Prospector SED-fitting framework, and fit the UV-IR broadband photometry of 129 nearby galaxies. We find that 10% of the fitted galaxies host an AGN contributing $>$10% of the observed galaxy MIR luminosity. We demonstrate the necessity of this AGN component in the following ways. First, we compare observed spectral features to spectral features predicted from our model fit to the photometry. We find that the AGN component greatly improves predictions for observed H$\alpha$ and H$\beta$ luminosities, as well as mid-infrared Akari and Spitzer/IRS spectra. Second, we show that inclusion of the AGN component changes stellar ages and SFRs by up to a factor of 10, and dust attenuations by up to a factor of 2.5. Finally, we show that the strength of our model AGN component correlates with independent AGN indicators, suggesting that these galaxies truly host AGN. Notably, only 46% of the SED-detected AGN would be detected with a simple MIR color selection. Based on these results, we conclude that SED models which fit MIR data without AGN components are vulnerable to substantial bias in their derived parameters.

We present the properties of NGC 4993, the host galaxy of GW170817, the first gravitational wave (GW) event from the merger of a binary neutron star (BNS) system and the first with an electromagnetic (EM) counterpart. We use both archival photometry and new optical/near-IR imaging and spectroscopy, together with stellar population synthesis models to infer the global properties of the host galaxy. We infer a star formation history peaked at $\gtrsim 10$ Gyr ago, with subsequent exponential decline leading to a low current star formation rate of 0.01 M$_{\odot}$ yr$^{-1}$, which we convert into a binary merger timescale probability distribution. We find a median merger timescale of $11.2^{+0.7}_{-1.4}$ Gyr, with a 90% confidence range of $6.8-13.6$ Gyr. This in turn indicates an initial binary separation of $\approx 4.5$ R$_{\odot}$, comparable to the inferred values for Galactic BNS systems. We also use new and archival $Hubble$ $Space$ $Telescope$ images to measure a projected offset of the optical counterpart of $2.1$ kpc (0.64$r_{e}$) from the center of NGC 4993 and to place a limit of $M_{r} \gtrsim -7.2$ mag on any pre-existing emission, which rules out the brighter half of the globular cluster luminosity function. Finally, the age and offset of the system indicates it experienced a modest natal kick with an upper limit of $\sim 200$ km s$^{-1}$. Future GW$-$EM observations of BNS mergers will enable measurement of their population delay time distribution, which will directly inform their viability as the dominant source of $r$-process enrichment in the Universe.

At redshift z=0.03, the recently-discovered SN 2017egm is the nearest Type I superluminous supernova (SLSN) to date, and first near the center of a massive spiral galaxy (NGC 3191). Using SDSS spectra of NGC 3191, we find a metallicity ~2 Z$_\odot$ at the nucleus and ~1.3 Z$_\odot$ for a star forming region at a radial offset similar to SN 2017egm. Archival radio-to-UV photometry reveals a star-formation rate ~15 M$_\odot$ yr$^{-1}$ (with ~70% dust-obscured), which can account for a Swift X-ray detection, and stellar mass ~$10^{10.7}$ M$_\odot$. We model the early UV-optical light curves with a magnetar central-engine model, using the Bayesian light curve fitting tool MOSFiT. The fits indicate ejecta mass 2-4 M$_\odot$, spin period 4-6 ms, magnetic field (0.7-1.7)$\times 10^{14}$G, and kinetic energy 1-2 $\times10^{51}$ erg. These parameters are consistent with the overall distributions for SLSNe, modeled by Nicholl et al (2017), although the derived mass and spin are towards the low end, possibly indicating enhanced loss of mass and angular momentum before explosion. This has two implications: (i) SLSNe can occur at solar metallicity, although with a low fraction ~10%; and (ii) metallicity has at most a modest effect on their properties. Both conclusions are in line with results for long gamma-ray bursts. Assuming a monotonic rise gives an explosion date MJD $57889\pm1$. However, a short-lived excess in the data relative to the best-fitting models may indicate an early-time `bump'. If confirmed, SN 2017egm would be the first SLSN with a spectrum during the bump-phase; this shows the same O II lines seen at maximum light, which may be an important clue for explaining these bumps.

Broadband photometry of galaxies measures an unresolved mix of complex stellar populations, gas, and dust. Interpreting these data is a challenge for models: many studies have shown that properties derived from modeling galaxy photometry are uncertain by a factor of two or more, and yet answering key questions in the field now requires higher accuracy than this. Here, we present a new model framework specifically designed for these complexities. Our model, Prospector-$\alpha$, includes dust attenuation and re-radiation, a flexible attenuation curve, nebular emission, stellar metallicity, and a 6-component nonparametric star formation history. The flexibility and range of the parameter space, coupled with MCMC sampling within the Prospector inference framework, is designed to provide unbiased parameters and realistic error bars. We assess the accuracy of the model with aperture-matched optical spectroscopy, which was excluded from the fits. We compare spectral features predicted solely from fits to the broadband photometry to the observed spectral features. Our model predicts H$\alpha$ luminosities with a scatter of $\sim$0.18 dex and an offset of $\sim$0.1 dex across a wide range of morphological types and stellar masses. This agreement is remarkable, as the H$\alpha$ luminosity is dependent on accurate star formation rates, dust attenuation, and stellar metallicities. The model also accurately predicts dust-sensitive Balmer decrements, spectroscopic stellar metallicities, PAH mass fractions, and the age- and metallicity-sensitive features D$_{\mathrm{n}}$4000 and H$\delta$. Although the model passes all these tests, we caution that we have not yet assessed its performance at higher redshift or the accuracy of recovered stellar masses.

[Abridged] We present observations of PS16dtm, a luminous transient that occurred at the nucleus of a known Narrow-line Seyfert 1 galaxy hosting a 10$^6$ M$_\odot$ black hole. The transient was previously claimed to be a Type IIn SLSN due to its luminosity and hydrogen emission lines. The light curve shows that PS16dtm brightened by about two magnitudes in ~50 days relative to the archival host brightness and then exhibited a plateau phase for about 100 days followed by the onset of fading in the UV. During the plateau PS16dtm showed no color evolution, maintained a blackbody temperature of 1.7 x 10$^4$ K, and radiated at approximately $L_{Edd}$ of the SMBH. The spectra exhibit multi-component hydrogen emission lines and strong FeII emission, show little evolution with time, and closely resemble the spectra of NLS1s while being distinct from those of Type IIn SNe. Moreover, PS16dtm is undetected in the X-rays to a limit an order of magnitude below an archival X-ray detection of its host galaxy. These observations strongly link PS16dtm to activity associated with the SMBH and are difficult to reconcile with a SN origin or any known form of AGN variability, and therefore we argue that it is a TDE in which the accretion of the stellar debris powers the rise in the continuum and excitation of the pre-existing broad line region, while providing material that obscures the X-ray emitting region of the pre-existing AGN accretion disk. A detailed TDE model fit to the light curve indicates that PS16dtm will remain bright for several years; we further predict that the X-ray emission will reappear on a similar timescale as the accretion rate declines. Finally, we place PS16dtm in the context of other TDEs and find that TDEs in AGN galaxies are an order of magnitude more efficient and reach Eddington luminosities, likely due to interaction of the stellar debris with the pre-existing accretion disk.

Recent surveys have identified a seemingly ubiquitous population of galaxies with elevated [OIII]/H$\beta$ emission line ratios at $z > 1$, though the nature of this phenomenon continues to be debated. The [OIII]/H$\beta$ line ratio is of interest because it is a main component of the standard diagnostic tools used to differentiate between active galactic nuclei (AGN) and star-forming galaxies, as well as the gas-phase metallicity indicators $O_{23}$ and $R_{23}$. Here, we investigate the primary driver of increased [OIII]/H$\beta$ ratios by median-stacking rest-frame optical spectra for a sample of star-forming galaxies in the 3D-HST survey in the redshift range $z\sim1.4-2.2$. Using $N = 4220$ star-forming galaxies, we stack the data in bins of mass and specific star formation rates (sSFR) respectively. After accounting for stellar Balmer absorption, we measure [OIII]$\lambda5007$\AA/H$\beta$ down to $\mathrm{M} \sim 10^{9.2} \ \mathrm{M_\odot}$ and sSFR $\sim 10^{-9.6} \ \mathrm{yr}^{-1}$, more than an order of magnitude lower than previous work at similar redshifts. We find an offset of $0.59\pm0.05$ dex between the median ratios at $z\sim2$ and $z\sim0$ at fixed stellar mass, in agreement with existing studies. However, with respect to sSFR, the $z \sim 2$ stacks all lie within 1$\sigma$ of the median SDSS ratios, with an average offset of only $-0.06\pm 0.05$. We find that the excitation properties of galaxies are tightly correlated with their sSFR at both $z\sim2$ and $z\sim0$, with a relation that appears to be roughly constant over the last 10 Gyr of cosmic time.

We present a new technique for wide and shallow observations using the near-infrared channel of Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). Wide-field near-IR surveys with HST are generally inefficient, as guide star acquisitions make it impractical to observe more than one pointing per orbit. This limitation can be circumvented by guiding with gyros alone, which is possible as long as the telescope has three functional gyros. The method presented here allows us to observe mosaics of eight independent WFC3-IR pointings in a single orbit by utilizing the fact that HST drifts by only a very small amount in the 25 seconds between non-destructive reads of unguided exposures. By shifting the reads and treating them as independent exposures the full resolution of WFC3 can be restored. We use this "drift and shift" (DASH) method in the Cycle 23 COSMOS-DASH program, which will obtain 456 WFC3 $H_{160}$ pointings in 57 orbits, covering an area of 0.6 degree$^2$ in the COSMOS field down to $H_{160} = 25$. When completed, the program will more than triple the area of extra-galactic survey fields covered by near-IR imaging at HST resolution. We demonstrate the viability of the method with the first four orbits (32 pointings) of this program. We show that the resolution of the WFC3 camera is preserved, and that structural parameters of galaxies are consistent with those measured in guided observations.

We explore the presence of non-stellar rest-frame near-IR ($2-5 \ \mu \mathrm{m}$) emission in galaxies at $z \sim 1$. Previous studies identified this excess in relatively small samples and suggested that such non-stellar emission, which could be linked to the $3.3 \ \mu \mathrm{m}$ polycyclic aromatic hydrocarbons feature or hot dust emission, is associated with an increased star formation rate (SFR). In this Letter, we confirm and quantify the presence of an IR excess in a significant fraction of galaxies in the 3D-HST GOODS catalogs. By constructing a matched sample of galaxies with and without strong non-stellar near-IR emission, we find that galaxies with such emission are predominantly star-forming galaxies. Moreover, star-forming galaxies with an excess show increased mid- and far-IR and H$\alpha$ emission compared to other star-forming galaxies without. While galaxies with a near-IR excess show a larger fraction of individually detected X-ray active galactic nuclei (AGNs), an X-ray stacking analysis, together with the IR-colors and H$\alpha$ profiles, shows that AGNs are unlikely to be the dominant source of the excess in the majority of galaxies. Our results suggest that non-stellar near-IR emission is linked to increased SFRs and is ubiquitous among star-forming galaxies. As such, the near-IR emission might be a powerful tool to measure SFRs in the era of the James Webb Space Telescope.

We present reduced data and data products from the 3D-HST survey, a 248-orbit HST Treasury program. The survey obtained WFC3 G141 grism spectroscopy in four of the five CANDELS fields: AEGIS, COSMOS, GOODS-S, and UDS, along with WFC3 $H_{140}$ imaging, parallel ACS G800L spectroscopy, and parallel $I_{814}$ imaging. In a previous paper (Skelton et al. 2014) we presented photometric catalogs in these four fields and in GOODS-N, the fifth CANDELS field. Here we describe and present the WFC3 G141 spectroscopic data, again augmented with data from GO-1600 in GOODS-N. The data analysis is complicated by the fact that no slits are used: all objects in the WFC3 field are dispersed, and many spectra overlap. We developed software to automatically and optimally extract interlaced 2D and 1D spectra for all objects in the Skelton et al. (2014) photometric catalogs. The 2D spectra and the multi-band photometry were fit simultaneously to determine redshifts and emission line strengths, taking the morphology of the galaxies explicitly into account. The resulting catalog has 98,663 measured redshifts and line strengths down to $JH_{IR}\leq 26$ and 22,548 with $JH_{IR}\leq 24$, where we comfortably detect continuum emission. Of this sample 5,459 galaxies are at $z>1.5$ and 9,621 are at $0.7<z<1.5$, where H$\alpha$ falls in the G141 wavelength coverage. Based on comparisons with ground-based spectroscopic redshifts, and on analyses of paired galaxies and repeat observations, the typical redshift error for $JH_{IR}\leq 24$ galaxies in our catalog is $\sigma_z \approx 0.003 \times (1+z)$, i.e., one native WFC3 pixel. The $3\sigma$ limit for emission line fluxes of point sources is $1.5\times10^{-17}$ ergs s$^{-1}$ cm$^{-2}$. We show various representations of the full dataset, as well as individual examples that highlight the range of spectra that we find in the survey.

We present a study of photometric redshift accuracy in the 3D-HST photometric catalogs, using 3D-HST grism redshifts to quantify and dissect trends in redshift accuracy for galaxies brighter than $H_{F140W}<24$ with an unprecedented and representative high-redshift galaxy sample. We find an average scatter of $0.0197\pm0.0003(1+z)$ in the Skelton et al. (2014) photometric redshifts. Photometric redshift accuracy decreases with magnitude and redshift, but does not vary monotonically with color or stellar mass. The 1-$\sigma$ scatter lies between $0.01-0.03$(1+z) for galaxies of all masses and colors below $z<2.5$ (for $H_{F140W}{<}24$), with the exception of a population of very red ($U-V > 2$), dusty star-forming galaxies for which the scatter increases to $\sim0.1(1+z)$. Although the overall photometric redshift accuracy for quiescent galaxies is better than for star-forming galaxies, scatter depends more strongly on magnitude and redshift than on galaxy type. We verify these trends using the redshift distributions of close pairs and extend the analysis to fainter objects, where photometric redshift errors further increase to $\sim0.046(1+z)$ at $H_{F160W}=26$. We demonstrate that photometric redshift accuracy is strongly filter-dependent and quantify the contribution of multiple filter combinations. We evaluate the widths of redshift probability distribution functions and find that error estimates are underestimated by a factor of $\sim1.1-1.6$, but that uniformly broadening the distribution does not adequately account for fitting outliers. Finally, we suggest possible applications of these data in planning for current and future surveys and simulate photometric redshift performance in the LSST, DES, and combined DES and VHS surveys.

In this paper we study a key phase in the formation of massive galaxies: the transition of star forming galaxies into massive (M_stars~10^11 Msun), compact (r_e~1 kpc) quiescent galaxies, which takes place from z~3 to z~1.5. We use HST grism redshifts and extensive photometry in all five 3D-HST/CANDELS fields, more than doubling the area used previously for such studies, and combine these data with Keck MOSFIRE and NIRSPEC spectroscopy. We first confirm that a population of massive, compact, star forming galaxies exists at z~2, using K-band spectroscopy of 25 of these objects at 2.0<z<2.5. They have a median NII/Halpha ratio of 0.6, are highly obscured with SFR(tot)/SFR(Halpha)~10, and have a large range of observed line widths. We infer from the kinematics and spatial distribution of Halpha that the galaxies have rotating disks of ionized gas that are a factor of ~2 more extended than the stellar distribution. By combining measurements of individual galaxies, we find that the kinematics are consistent with a nearly Keplerian fall-off from V_rot~500 km/s at 1 kpc to V_rot~250 km/s at 7 kpc, and that the total mass out to this radius is dominated by the dense stellar component. Next, we study the size and mass evolution of the progenitors of compact massive galaxies. Even though individual galaxies may have had complex histories with periods of compaction and mergers, we show that the population of progenitors likely followed a simple inside-out growth track in the size-mass plane of d(log r_e) ~ 0.3 d(log M_stars). This mode of growth gradually increases the stellar mass within a fixed physical radius, and galaxies quench when they reach a stellar density or velocity dispersion threshold. As shown in other studies, the mode of growth changes after quenching, as dry mergers take the galaxies on a relatively steep track in the size-mass plane.

It is well established that (1) star-forming galaxies follow a relation between their star formation rate (SFR) and stellar mass (M$_{\star}$), the "star-formation sequence", and (2) the SFRs of galaxies correlate with their structure, where star-forming galaxies are less concentrated than quiescent galaxies at fixed mass. Here, we consider whether the scatter and slope of the star-formation sequence is correlated with systematic variations in the Sersic indices, $n$, of galaxies across the SFR-M$_{\star}$ plane. We use a mass-complete sample of 23,848 galaxies at $0.5<z<2.5$ selected from the 3D-HST photometric catalogs. Galaxy light profiles parameterized by $n$ are based on Hubble Space Telescope CANDELS near-infrared imaging. We use a single SFR indicator empirically-calibrated from stacks of Spitzer/MIPS 24$\mu$m imaging, adding the unobscured and obscured star formation. We find that the scatter of the star-formation sequence is related in part to galaxy structure; the scatter due to variations in $n$ at fixed mass for star-forming galaxies ranges from 0.14$\pm$0.02 dex at $z\sim2$ to 0.30$\pm$0.04 dex at $z<1$. While the slope of the log(SFR)-log(M$_{\star}$) relation is of order unity for disk-like galaxies, galaxies with $n>2$ (implying more dominant bulges) have significantly lower SFR/M$_{\star}$ than the main ridgeline of the star-formation sequence. These results suggest that bulges in massive $z\sim2$ galaxies are actively building up, where the stars in the central concentration are relatively young. At $z<1$, the presence of older bulges within star-forming galaxies lowers global SFR/M$_{\star}$, decreasing the slope and contributing significantly to the scatter of the star-formation sequence.

We present Ha maps at 1kpc spatial resolution for star-forming galaxies at z~1, made possible by the WFC3 grism on HST. Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 2676 galaxies. By creating deep stacked Halpha (Ha) images, we reach surface brightness limits of 1x10^-18\erg\s\cm^2\arcsec^2, allowing us to map the distribution of ionized gas out to >10kpc for typical L* galaxies at this epoch. We find that the spatial extent of the Ha distribution increases with stellar mass as r(Ha)[kpc]=1.5(Mstars/10^10Msun)^0.23. Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect, however, is mass dependent with r(Ha)/r(stars)=1.1(M/10^10Msun)^0.054, such that at low masses r(Ha)~r(stars). We map the Ha distribution as a function of SFR(IR+UV) and find evidence for `coherent star formation' across the SFR-M plane: above the main sequence, Ha is enhanced at all radii; below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at fixed mass. At high masses (10^10.5<M/Msun<10^11), above the main sequence, Ha is particularly enhanced in the center, plausibly building bulges and/or supermassive black holes. Below the main sequence, the star forming disks are more compact and a strong central dip in the EW(Ha), and the inferred specific star formation rate, appears. Importantly though, across the entirety of the SFR-M plane, the absolute star formation rate as traced by Ha is always centrally peaked, even in galaxies below the main sequence.

Interpreting observations of distant galaxies in terms of constraints on physical parameters - such as stellar mass, star-formation rate (SFR) and dust optical depth - requires spectral synthesis modelling. We analyse the reliability of these physical parameters as determined under commonly adopted `classical' assumptions: star-formation histories assumed to be exponentially declining functions of time, a simple dust law and no emission-line contribution. Improved modelling techniques and data quality now allow us to use a more sophisticated approach, including realistic star-formation histories, combined with modern prescriptions for dust attenuation and nebular emission (Pacifici et al. 2012). We present a Bayesian analysis of the spectra and multi-wavelength photometry of 1048 galaxies from the 3D-HST survey in the redshift range 0.7<z<2.8 and in the stellar mass range 9<log(M/Mo)<12. We find that, using the classical spectral library, stellar masses are systematically overestimated (~0.1 dex) and SFRs are systematically underestimated (~0.6 dex) relative to our more sophisticated approach. We also find that the simultaneous fit of photometric fluxes and emission-line equivalent widths helps break a degeneracy between SFR and optical depth of the dust, reducing the uncertainties on these parameters. Finally, we show how the biases of classical approaches can affect the correlation between stellar mass and SFR for star-forming galaxies (the `Star-Formation Main Sequence'). We conclude that the normalization, slope and scatter of this relation strongly depend on the adopted approach and demonstrate that the classical, oversimplified approach cannot recover the true distribution of stellar mass and SFR.

We examine the connection between the observed star-forming sequence (SFR $\propto$ $M^{\alpha}$) and the observed evolution of the stellar mass function between $0.2 < z < 2.5$. We find the star-forming sequence cannot have a slope $\alpha$ $\lesssim$ 0.9 at all masses and redshifts, as this would result in a much higher number density at $10 < \log(\mathrm{M/M_{\odot}}) < 11$ by $z=1$ than is observed. We show that a transition in the slope of the star-forming sequence, such that $\alpha=1$ at $\log(\mathrm{M/M_{\odot}})<10.5$ and $\alpha=0.7-0.13z$ ({Whitaker} {et~al.} 2012) at $\log(\mathrm{M/M_{\odot}})>10.5$, greatly improves agreement with the evolution of the stellar mass function. We then derive a star-forming sequence which reproduces the evolution of the mass function by design. This star-forming sequence is also well-described by a broken-power law, with a shallow slope at high masses and a steep slope at low masses. At $z=2$, it is offset by $\sim$0.3 dex from the observed star-forming sequence, consistent with the mild disagreement between the cosmic SFR and recent observations of the growth of the stellar mass density. It is unclear whether this problem stems from errors in stellar mass estimates, errors in SFRs, or other effects. We show that a mass-dependent slope is also seen in other self-consistent models of galaxy evolution, including semi-analytical, hydrodynamical, and abundance-matching models. As part of the analysis, we demonstrate that neither mergers nor hidden low-mass quiescent galaxies are likely to reconcile the evolution of the mass function and the star-forming sequence. These results are supported by observations from {Whitaker} {et~al.} (2014).

We constrain the slope of the star formation rate ($\log\Psi$) to stellar mass ($\log\mathrm{M_{\star}}$) relation down to $\log(\mathrm{M_{\star}/M_{\odot}})=8.4$ ($\log(\mathrm{M_{\star}/M_{\odot}})=9.2$) at $z=0.5$ ($z=2.5$) with a mass-complete sample of 39,106 star-forming galaxies selected from the 3D-HST photometric catalogs, using deep photometry in the CANDELS fields. For the first time, we find that the slope is dependent on stellar mass, such that it is steeper at low masses ($\log\mathrm{\Psi}\propto\log\mathrm{M_{\star}}$) than at high masses ($\log\mathrm{\Psi}\propto(0.3-0.6)\log\mathrm{M_{\star}}$). These steeper low mass slopes are found for three different star formation indicators: the combination of the ultraviolet (UV) and infrared (IR), calibrated from a stacking analysis of Spitzer/MIPS 24$\mu$m imaging; $\beta$-corrected UV SFRs; and H$\alpha$ SFRs. The normalization of the sequence evolves differently in distinct mass regimes as well: for galaxies less massive than $\log(\mathrm{M_{\star}/M_{\odot}})<10$ the specific SFR ($\Psi/\mathrm{M_{\star}}$) is observed to be roughly self-similar with $\Psi/\mathrm{M_{\star}}\propto(1+z)^{1.9}$, whereas more massive galaxies show a stronger evolution with $\Psi/\mathrm{M_{\star}}\propto(1+z)^{2.2-3.5}$ for $\log(\mathrm{M_{\star}/M_{\odot}})=10.2-11.2$. The fact that we find a steep slope of the star formation sequence for the lower mass galaxies will help reconcile theoretical galaxy formation models with the observations. The results of this study support the analytical conclusions of Leja et al. (2014).

The dense interiors of massive galaxies are among the most intriguing environments in the Universe. In this paper we ask when these dense cores were formed and determine how galaxies gradually assembled around them. We select galaxies that have a stellar mass >3x10^10 Msun inside r=1 kpc out to z=2.5, using the 3D-HST survey and data at low redshift. Remarkably, the number density of galaxies with dense cores appears to have decreased from z=2.5 to the present. This decrease is probably mostly due to stellar mass loss and the resulting adiabatic expansion, with some contribution from merging. We infer that dense cores were mostly formed at z>2.5, consistent with their largely quiescent stellar populations. While the cores appear to form early, the galaxies in which they reside show strong evolution: their total masses increase by a factor of 2-3 from z=2.5 to z=0 and their effective radii increase by a factor of 5-6. As a result, the contribution of dense cores to the total mass of the galaxies in which they reside decreases from ~50% at z=2.5 to ~15% at z=0. Because of their early formation, the contribution of dense cores to the total stellar mass budget of the Universe is a strong function of redshift. The stars in cores with M_1kpc>3x10^10 Msun make up ~0.1% of the stellar mass density of the Universe today but 10%-20% at z~2, depending on their IMF. The formation of these cores required the conversion of ~10^11 Msun of gas into stars within ~1 kpc, while preventing significant star formation at larger radii.

The 3D-HST and CANDELS programs have provided WFC3 and ACS spectroscopy and photometry over ~900 square arcminutes in five fields: AEGIS, COSMOS, GOODS-North, GOODS-South, and the UKIDSS UDS field. All these fields have a wealth of publicly available imaging datasets in addition to the HST data, which makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range. In this paper we describe a photometric analysis of the CANDELS and 3D-HST HST imaging and the ancillary imaging data at wavelengths 0.3um to 8um. Objects were selected in the WFC3 near-IR bands, and their SEDs were determined by carefully taking the effects of the point spread function in each observation into account. A total of 147 distinct imaging datasets were used in the analysis. The photometry is made available in the form of six catalogs: one for each field, as well as a master catalog containing all objects in the entire survey. We also provide derived data products: photometric redshifts, determined with the EAZY code, and stellar population parameters determined with the FAST code. We make all the imaging data that were used in the analysis available, including our reductions of the WFC3 imaging in all five fields. 3D-HST is a spectroscopic survey with the WFC3 and ACS grisms, and the photometric catalogs presented here constitute a necessary first step in the analysis of these grism data. All the data presented in this paper are available through the 3D-HST website.

Spectroscopic + photometric redshifts, stellar mass estimates, and rest-frame colors from the 3D-HST survey are combined with structural parameter measurements from CANDELS imaging to determine the galaxy size-mass distribution over the redshift range 0<z<3. Separating early- and late-type galaxies on the basis of star-formation activity, we confirm that early-type galaxies are on average smaller than late-type galaxies at all redshifts, and find a significantly different rate of average size evolution at fixed galaxy mass, with fast evolution for the early-type population, R_eff ~ (1+z)^-1.48, and moderate evolution for the late-type population, R_eff ~ (1+z)^-0.75. The large sample size and dynamic range in both galaxy mass and redshift, in combination with the high fidelity of our measurements due to the extensive use of spectroscopic data, not only fortify previous results, but also enable us to probe beyond simple average galaxy size measurements. At all redshifts the slope of the size-mass relation is shallow, R_eff ~ M_star^0.22, for late-type galaxies with stellar mass >3x10^9 M_sol, and steep, R_eff M_star^0.75, for early-type galaxies with stellar mass >2x10^10 M_sol. The intrinsic scatter is <~0.2 dex for all galaxy types and redshifts. For late-type galaxies, the logarithmic size distribution is not symmetric, but skewed toward small sizes: at all redshifts and masses a tail of small late-type galaxies exists that overlaps in size with the early-type galaxy population. The number density of massive (~10^11 M_sol), compact (R_eff < 2 kpc) early-type galaxies increases from z=3 to z=1.5-2 and then strongly decreases at later cosmic times.

Most massive galaxies are thought to have formed their dense stellar cores at early cosmic epochs. However, cores in their formation phase have not yet been observed. Previous studies have found galaxies with high gas velocity dispersions or small apparent sizes but so far no objects have been identified with both the stellar structure and the gas dynamics of a forming core. Here we present a candidate core in formation 11 billion years ago, at z=2.3. GOODS-N-774 has a stellar mass of 1.0x10^11 Msun, a half-light radius of 1.0 kpc, and a star formation rate of 90[+45-20]Msun/yr. The star forming gas has a velocity dispersion 317+-30 km/s, amongst the highest ever measured. It is similar to the stellar velocity dispersions of the putative descendants of GOODS-N-774, compact quiescent galaxies at z~2 and giant elliptical galaxies in the nearby Universe. Galaxies such as GOODS-N-774 appear to be rare; however, from the star formation rate and size of the galaxy we infer that many star forming cores may be heavily obscured, and could be missed in optical and near-infrared surveys.

We present direct observational evidence for star formation quenching in galaxy groups in the redshift range 0<z<2.5. We utilize a large sample of nearly 6000 groups, selected by fixed cumulative number density from three photometric catalogs, to follow the evolving quiescent fractions of central and satellite galaxies over roughly 11 Gyr. At z~0, central galaxies in our sample range in stellar mass from Milky Way/M31 analogs (M=6.5x10^10 M\solar) to nearby massive ellipticals (M=1.5x10^11 M\solar). Satellite galaxies in the same groups reach masses as low as twice that of the Large Magellanic Cloud (M=6.5x10^9 M\solar). Using statistical background subtraction, we measure the average rest-frame colors of galaxies in our groups and calculate the evolving quiescent fractions of centrals and satellites over seven redshift bins. Our analysis shows clear evidence for star formation quenching in group halos, with a different quenching onset for centrals and their satellite galaxies. Using halo mass estimates for our central galaxies, we find that star formation shuts off in centrals when typical halo masses reach between 10^12 and 10^13 M\solar, consistent with predictions from the halo quenching model. In contrast, satellite galaxies in the same groups most likely undergo quenching by environmental processes, whose onset is delayed with respect to their central galaxy. Although star formation is suppressed in all galaxies over time, the processes that govern quenching are different for centrals and satellites. While mass plays an important role in determining the star formation activity of central galaxies, quenching in satellite galaxies is dominated by the environment in which they reside.

We present Keck/MOSFIRE K-band spectroscopy of the first mass-selected sample of galaxies at $z\sim2.3$. Targets are selected from the 3D-HST Treasury survey. The six detected galaxies have a mean [NII]$\lambda$6584/H$\alpha$ ratio of $0.27\pm0.01$, with a small standard deviation of 0.05. This mean value is similar to that of UV-selected galaxies of the same mass. The mean gas-phase oxygen abundance inferred from the [NII]/H$\alpha$ ratios depends on the calibration method, and ranges from 12+log(O/H)$_{gas}=8.57$ for the {Pettini} & {Pagel} (2004) calibration to 12+log(O/H)$_{gas}= 8.87$ for the {Maiolino} {et~al.} (2008) calibration. Measurements of the stellar oxygen abundance in nearby quiescent galaxies with the same number density indicate 12+log(O/H)$_{stars}= 8.95$, similar to the gas-phase abundances of the $z\sim2.3$ galaxies if the {Maiolino} {et~al.} (2008) calibration is used. This suggests that these high-redshift star forming galaxies may be progenitors of today's massive early-type galaxies. The main uncertainties are the absolute calibration of the gas-phase oxygen abundance and the incompleteness of the $z\sim2.3$ sample: the galaxies with detected H$\alpha$ tend to be larger and have higher star formation rates than the galaxies without detected H$\alpha$, and we may still be missing the most dust-obscured progenitors.

We follow the structural evolution of star forming galaxies (SFGs) like the Milky Way by selecting progenitors to z~1.3 based on the stellar mass growth inferred from the evolution of the star forming sequence. We select our sample from the 3D-HST survey, which utilizes spectroscopy from the HST WFC3 G141 near-IR grism and enables precise redshift measurements for our sample of SFGs. Structural properties are obtained from Sersic profile fits to CANDELS WFC3 imaging. The progenitors of z=0 SFGs with stellar mass M=10^{10.5} Msun are typically half as massive at z~1. This late-time stellar mass assembly is consistent with recent studies that employ abundance matching techniques. The descendant SFGs at z~0 have grown in half-light radius by a factor of ~1.4 since z~1. The half-light radius grows with stellar mass as r_e M^{0.29}. While most of the stellar mass is clearly assembling at large radii, the mass surface density profiles reveal ongoing mass growth also in the central regions where bulges and pseudobulges are common features in present day late-type galaxies. Some portion of this growth in the central regions is due to star formation as recent observations of H-alpha maps for SFGs at z~1 are found to be extended but centrally peaked. Connecting our lookback study with galactic archeology, we find the stellar mass surface density at R=8 kpc to have increased by a factor of ~2 since z~1, in good agreement with measurements derived for the solar neighborhood of the Milky Way.

The Fundamental Plane (FP) is an empirical relation between the size, surface brightness, and velocity dispersion of early-type galaxies. This relation has been studied extensively for early-type galaxies in the local universe to constrain galaxy formation mechanisms. The evolution of the zeropoint of this plane has been extended to high redshifts to study the luminosity evolution of massive galaxies, under the assumption of structural homology. In this work, we assess this assumption by replacing surface brightness with stellar mass density and present the evolution of the "mass FP" for massive, quiescent galaxies since z~2. By accounting for stellar populations, we thereby isolate and trace structural and dynamical evolution. Despite the observed dramatic evolution in the sizes and morphologies of massive galaxies since z~3, we find that quiescent galaxies lie on the mass FP out to z~2. In contrast with ~1.4 dex evolution in the luminosity FP, average residuals from the z~0 mass FP are less than ~0.15 dex since z~2. Assuming the Hyde & Bernardi (2009) mass FP slope, we find that this minimal offset scales as (1+z)^{-0.095+/-0.043}. This result lends credence to previous studies that derived luminosity evolution from the FP. Therefore, despite their compact sizes and suggestions that massive galaxies are more disk-like at z~2, the relationship between their dynamics and structural properties are consistent with local early-type galaxies. Finally, we find no strong evidence for a tilt of the mass FP relative to the Virial plane, but emphasize the need for full models including selection biases to fully investigate this issue.

Galaxies with the mass of the Milky Way dominate the stellar mass density of the Universe but it is uncertain how and when they were assembled. Here we study progenitors of these galaxies out to z=2.5, using data from the 3D-HST and CANDELS Treasury surveys. We find that galaxies with present-day stellar masses of log(M)~10.7 built ~90% of their stellar mass since z=2.5, with most of the star formation occurring before z=1. In marked contrast to the assembly history of massive elliptical galaxies, mass growth is not limited to large radii: the mass in the central 2 kpc of the galaxies increased by a factor of 3.2+-0.8 between z=2.5 and z=1. We therefore rule out simple models in which bulges were fully assembled at high redshift and disks gradually formed around them. Instead, bulges (and black holes) likely formed in lockstep with disks, through bar instabilities, migration, or other processes. We find that after z=1 the growth in the central regions gradually stopped and the disk continued to build, consistent with recent studies of the gas distributions in z~1 galaxies and the properties of many spiral galaxies today.

We present a statistical study of the environments of massive galaxies in four redshift bins between z=0.04 and z=1.6, using data from the Sloan Digital Sky Survey (SDSS) and the NEWFIRM Medium Band Survey (NMBS). We measure the projected radial distribution of galaxies in cylinders around a constant number density selected sample of massive galaxies and utilize a statistical subtraction of contaminating sources. Our analysis shows that massive primary galaxies typically live in group halos and are surrounded by 2 to 3 satellites with masses more than one-tenth of the primary galaxy mass. The cumulative stellar mass in these satellites roughly equals the mass of the primary galaxy itself. We further find that the radial number density profile of galaxies around massive primaries has not evolved significantly in either slope or overall normalization in the past 9.5 Gyr. A simplistic interpretation of this result can be taken as evidence for a lack of mergers in the studied groups and as support for a static evolution model of halos containing massive primaries. Alternatively, there exists a tight balance between mergers and accretion of new satellites such that the overall distribution of galaxies in and around the halo is preserved. The latter interpretation is supported by a comparison to a semi-analytic model, which shows a similar constant average satellite distribution over the same redshift range.

A central challenge in observational studies of galaxy formation is how to associate progenitor galaxies with their descendants at lower redshifts. One promising approach is to link galaxies at fixed number density, rather than fixed luminosity or mass. This method is effective if stellar mass rank order is broadly conserved through cosmic time. In this paper, we use the Guo et al. (2011) semi-analytical model to analyze under what circumstances this assumption is valid in the context of a cosmological simulation. Specifically, we select progenitor galaxies at a constant number density and compare the stellar mass evolution of their descendants to the evolution at a constant number density. The median stellar mass of the descendants increases by a factor of four (0.6 dex) from z = 3 to z = 0. Constant number density selection reproduces this to within 40% (0.15 dex) over a wide range of number densities. We show that the discrepancy primarily results from scatter in the stellar mass growth rates and merging. After applying simple, observationally-based corrections for these processes, the discrepancy is reduced to 12% (0.05 dex). We conclude that number density selection can be used to predict the median descendant mass of high-redshift progenitor galaxies. The main uncertainty in this study is that semi-analytical models do not reproduce the observed mass evolution of galaxies, which makes the quantitative aggregate effects of star formation, merging, and quenching on the rank order of galaxies somewhat uncertain.

The assembly of galaxies can be described by the distribution of their star formation as a function of cosmic time. Thanks to the WFC3 grism on HST it is now possible to measure this beyond the local Universe. Here we present the spatial distribution of Halpha emission for a sample of 54 strongly star-forming galaxies at z~1 in the 3D-HST Treasury survey. By stacking the Halpha emission we find that star formation occurred in approximately exponential distributions at z~1, with median Sersic index of n=1.0+-0.2. The stacks are elongated with median axis ratios of b/a=0.58+-0.09 in Halpha, consistent with (possibly thick) disks at random orientation angles. Keck spectra obtained for a subset of eight of the galaxies show clear evidence for rotation, with inclination-corrected velocities of 90 to 330 km/s. The most straightforward interpretation of our results is that star formation in strongly star-forming galaxies at z~1 generally occurred in disks. The disks appear to be "scaled-up" versions of nearby spiral galaxies: they have EW(Halpha)~100 Angstroms out to the solar orbit and they have star formation surface densities above the threshold for driving galactic scale winds.

We present 3D-HST, a near-infrared spectroscopic Treasury program with the Hubble Space Telescope for studying the processes that shape galaxies in the distant Universe. 3D-HST provides rest-frame optical spectra for a sample of ~7000 galaxies at 1<z<3.5, the epoch when 60% of all star formation took place, the number density of quasars peaked, the first galaxies stopped forming stars, and the structural regularity that we see in galaxies today must have emerged. 3D-HST will cover 3/4 (625 sq.arcmin) of the CANDELS survey area with two orbits of primary WFC3/G141 grism coverage and two to four parallel orbits with the ACS/G800L grism. In the IR these exposure times yield a continuum signal-to-noise of ~5 per resolution element at H~23.1 and a 5sigma emission line sensitivity of 5x10-17 erg/s/cm2 for typical objects, improving by a factor of ~2 for compact sources in images with low sky background levels. The WFC3/G141 spectra provide continuous wavelength coverage from 1.1-1.6 um at a spatial resolution of ~0."13, which, combined with their depth, makes them a unique resource for studying galaxy evolution. We present the preliminary reduction and analysis of the grism observations, including emission line and redshift measurements from combined fits to the extracted grism spectra and photometry from ancillary multi-wavelength catalogs. The present analysis yields redshift estimates with a precision of sigma(z)=0.0034(1+z), or sigma(v)~1000 km/s. We illustrate how the generalized nature of the survey yields near-infrared spectra of remarkable quality for many different types of objects, including a quasar at z=4.7, quiescent galaxies at z~2, and the most distant T-type brown dwarf star known. The CANDELS and 3D-HST surveys combined will provide the definitive imaging and spectroscopic dataset for studies of the 1<z<3.5 Universe until the launch of the James Webb Space Telescope.