Source author record

Marc Rafelski

Marc Rafelski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

48works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

48 published item(s)

preprint2026arXiv

BEACON: JWST NIRCam Pure-parallel Imaging Survey. II. Physical Properties of $z=7-14$ Galaxies

We present photometric properties of 161 galaxy candidates at $z=7-14$ selected from the second data release (DR2) of BEACON, a JWST Cycle 2 pure-parallel NIRCam imaging program. Carefully selected from 36 independent pointings (corresponding to $\sim350$\,arcmin$^2$ sky coverage), and hence with reduced cosmic variance, our galaxy candidates provide an unbiased sample for investigating galaxy properties over a wide range of environments. We measure the physical properties, including UV continuum slope ($β_{\rm UV}$), stellar mass ($M_*$), star formation rate (SFR), and sizes. Our highest redshift galaxy candidate at $z=13.71\pm0.15$ has a remarkably bright UV luminosity of $M_{\rm UV}=-21.19\pm0.08$, making it the brightest galaxy at $z>12$ if spectroscopically confirmed. With an extremely blue UV slope, compact morphology, and high star formation rate surface density ($Σ_{\rm SFR}$), this candidate may have extremely low metallicity, high ionizing photon escape fraction, or contributions from an AGN. Among our multiple independent sightlines, we identify three fields of galaxy number overdensity with $>3σ$ significance. The properties of galaxies in various environments do not exhibit significant differences, implying either that accelerated galaxy evolution in overdense regions is not yet widespread at $z>7$, or that the current constraints are limited by sample size. Our simulations indicate that increasing the sample by an order of magnitude would allow such environmental trends to be robustly confirmed or ruled out, underscoring the importance of future pure-parallel observations.

preprint2024arXiv

Lyman Continuum Emission from AGN at 2.3$\lesssim$z$\lesssim$3.7 in the UVCANDELS Fields

We present the results of our search for Lyman continuum (LyC) emitting AGN at redshifts 2.3$\lesssim$z$\lesssim$4.9 from HST WFC3 F275W observations in the UVCANDELS fields. We also include LyC emission from AGN using HST WFC3 F225W, F275W, and F336W found in the ERS and HDUV data. We performed exhaustive queries of the Vizier database to locate AGN with high quality spectroscopic redshifts. In total, we found 51 AGN that met our criteria within the UVCANDELS and ERS footprints. Of these 51, we find 12 AGN had $\geq$4$σ$ detected LyC flux in the WFC3/UVIS images. Using space- and ground-based data from X-ray to radio, we fit the multi-wavelength photometric data of each AGN to a CIGALE SED and correlate various SED parameters to the LyC flux. KS-tests of the SED parameter distributions for the LyC-detected and non-detected AGN showed they are likely not distinct samples. However, we find that X-ray luminosity, star-formation onset age, and disk luminosity show strong correlations relative to their emitted LyC flux. We also find strong correlation of the LyC flux to several dust parameters, i.e., polar and toroidal dust emission, 6 $μm$ luminosity, and anti-correlation with metallicity and $A_{FUV}$. We simulate the LyC escape fraction ($f_{esc}$) using the CIGALE and IGM transmission models for the LyC-detected AGN and find an average $f_{esc}$$\simeq$18%, weighted by uncertainties. We stack the LyC flux of subsamples of AGN according to the wavelength continuum region in which they are detected and find no significant distinctions in their LyC emission, although our $sub-mm\ detected$ F336W sample shows the brightest stacked LyC flux. These findings indicate that LyC-production and -escape in AGN is more complicated than the simple assumption of thermal emission and a 100% escape fraction. Further testing of AGN models with larger samples than presented here is needed.

preprint2022arXiv

A Massive, Dusty HI-Absorption-Selected Galaxy at $z \approx 2.46$ Identified in a CO Emission Survey

We report a NOrthern Extended Millimeter Array (NOEMA) and Atacama Large Millimeter/submillimeter Array (ALMA) search for redshifted CO emission from the galaxies associated with seven high-metallicity ([M/H] $\geq -1.03$) damped Lyman-$α$ absorbers (DLAs) at $z\approx1.64-2.51$. Our observations yielded one new detection of CO(3-2) emission from a galaxy at $z=2.4604$ using NOEMA, associated with the $z=2.4628$ DLA towards QSO B0201+365. Including previous searches, our search results in detection rates of CO emission of $\approx56^{+38}_{-24}$ % and $\approx11^{+26}_{-9}$ %, respectively, in the fields of DLAs with ${\rm [M/H]}>-0.3$ and ${\rm [M/H]}<-0.3$. Further, the HI-selected galaxies associated with five DLAs with [M/H] $>-0.3$ all have high molecular gas masses, $\gtrsim5\times10^{10}\ {\rm M}_\odot$. This indicates that the highest-metallicity DLAs at $z\approx2$ are associated with the most massive galaxies. The newly-identified $z\approx2.4604$ HI-selected galaxy, DLA0201+365g, has an impact parameter of $\approx7$ kpc to the QSO sightline, and an implied molecular gas mass of $(5.04\pm0.78) \times10^{10}\times(α_{\rm CO}/4.36)\times(r_{31}/0.55)\ {\rm M}_\odot$. Archival Hubble Space Telescope Wide Field and Planetary Camera 2 imaging covering the rest-frame near-ultraviolet (NUV) and far-ultraviolet (FUV) emission from this galaxy yield non-detections of rest-frame NUV and FUV emission, and a $5σ$ upper limit of 2.3 M$_\odot$ yr$^{-1}$ on the unobscured star formation rate (SFR). The low NUV-based SFR estimate, despite the very high molecular gas mass, indicates that DLA0201+365g either is a very dusty galaxy, or has a molecular gas depletion time that is around two orders of magnitude larger than that of star-forming galaxies at similar redshifts.

preprint2022arXiv

Constraining the Size of the Circumgalactic Medium Using the Transverse Autocorrelation Function of C IV Absorbers in Paired Quasar Spectra

The circumgalactic medium (CGM) plays a vital role in the formation and evolution of galaxies, acting as a lifeline between galaxies and the surrounding intergalactic medium (IGM). In this study we leverage a unique sample of quasar pairs to investigate the properties of the CGM with absorption line tomography. We present a new sample of medium resolution Keck/ESI, Magellan/MagE, and VLT/XSHOOTER spectra of 29 quasar pairs at redshift $2 < z < 3$. We supplement the sample with additional spectra of 32 pairs from the literature, creating a catalog of 61 quasar pairs with angular separations between 1.7" and 132.9" and projected physical separations ($r_\perp$) between 14 kpc and 887 kpc. We construct a catalog of 906 metal-line absorption doublets of C IV ($λλ1548, 1550$) with equivalent widths ranging from 6 mÅ $\leq W_{r, 1550} \leq 2053$ mÅ. The best fit linear model to the log-space equivalent width frequency distribution ($\log f(W_r) = m\log(W_{r}) + b$) of the sample yields coefficients of $m=-1.44\pm0.16$ and $b=-0.43\pm0.16$. To constrain the projected extent of C IV, we calculate the transverse autocorrelation function. The flattening of the autocorrelation function at low $r_\perp$ provides a lower limit for the coherence length of the metal enriched CGM - on the order of 200 $h^{-1}$ comoving kpc. This physical size constraint allows us to refine our understanding of the metals in the CGM, where the extent of C IV in the CGM depends on gas flows, feedback, timescale of metal injection and mixing, and the mass of the host galaxies.

preprint2022arXiv

Discovery of a damped Ly$α$ absorber originating in a spectacular interacting dwarf galaxy pair at $z = 0.026$

We present the discovery of neutral gas detected in both damped Ly$α$ absorption (DLA) and HI 21-cm emission outside of the stellar body of a galaxy, the first such detection in the literature. A joint analysis between the Cosmic Ultraviolet Baryon Survey and the MeerKAT Absorption Line Survey reveals an HI bridge connecting two interacting dwarf galaxies (log$(M_{\text{star}}/\text{M}_{\odot}) = 8.5 \pm 0.2$) that host a $z = 0.026$ DLA with log[$N$(HI)/cm$^{-2}$]$ = 20.60 \pm 0.05$ toward the QSO J2339-5523 ($z_{\text{QSO}} = 1.35$). At impact parameters of $d = 6$ and $33$ kpc, the dwarf galaxies have no companions more luminous than $\approx 0.05L_{*}$ within at least $Δv = \pm 300$ km s$^{-1}$ and $d \approx 350$ kpc. HI 21-cm emission is spatially coincident with the DLA at the 2-3$σ$ level per spectral channel over several adjacent beams. However, HI 21-cm absorption is not detected against the radio-bright QSO; if the background UV and radio sources are spatially aligned, the gas is either warm or clumpy (with spin temperature to covering factor ratio $T_{s}/f_{c} > 1880$ K). VLT-MUSE observations demonstrate that the $α$-element abundance of the ionized ISM is consistent with the DLA ($\approx 10$% solar), suggesting that the neutral gas envelope is perturbed ISM gas. This study showcases the impact of dwarf-dwarf interactions on the physical and chemical state of neutral gas outside of star-forming regions. In the SKA era, joint UV and HI 21-cm analyses will be critical for connecting the cosmic neutral gas content to galaxy environments.

preprint2022arXiv

Jansky Very Large Array detections of CO(1-0) emission in HI-absorption-selected galaxies at $z \gtrsim 2$

We report a Karl G. Jansky Very Large Array search for redshifted CO(1-0) emission from three HI-absorption-selected galaxies at $z \approx 2$, identified earlier in their CO(3-2) or CO(4-3) emission. We detect CO(1-0) emission from DLA B1228-113 at $z\approx2.1933$ and DLA J0918+1636 at $z\approx2.5848$; these are the first detections of CO(1-0) emission in high-$z$ HI-selected galaxies. We obtain high molecular gas masses, $\rm M_{mol}\approx10^{11}\times(α_{\rm CO}/4.36)\ M_\odot$, for the two objects with CO(1-0) detections, which are a factor of $\approx1.5-2$ lower than earlier estimates. We determine the excitation of the mid$-J$ CO rotational levels relative to the $J=1$ level, r$_{ J1}$, in HI-selected galaxies for the first time, obtaining r$_{\rm 31}=1.00\pm0.20$ and r$_{\rm 41}=1.03\pm0.23$ for DLA J0918+1636, and r$_{\rm 31}=0.86\pm0.21$ for DLA B1228-113. These values are consistent with thermal excitation of the $J=3,4$ levels. The excitation of the $J=3$ level in the HI-selected galaxies is similar to that seen in massive main-sequence and sub-mm galaxies at $z\gtrsim2$, but higher than that in main-sequence galaxies at $z\approx1.5$; the higher excitation of the galaxies at $z\gtrsim2$ is likely to be due to their higher star-formation rate (SFR) surface density. We use Hubble Space Telescope Wide Field Camera 3 imaging to detect the rest-frame near-ultraviolet emission of DLA B1228-113, obtaining an NUV SFR of $4.44\pm0.47$ M$_{\odot}$ yr$^{-1}$, significantly lower than that obtained from the total infrared luminosity, indicating significant dust extinction in the $z\approx2.1933$ galaxy.

preprint2022arXiv

Lyman Continuum Galaxy Candidates in COSMOS

Star-forming galaxies are the sources likely to have reionized the universe. As we cannot observe them directly due to the opacity of the intergalactic medium at $z\gtrsim5$, we study $z\sim3\text{--}5$ galaxies as proxies to place observational constraints on cosmic reionization. Using new deep \textit{Hubble Space Telescope} rest-frame UV F336W and F435W imaging (30-orbit, $\sim40$~arcmin$^2$, $\sim29\text{--}30$~mag depth at 5$σ$), we attempt to identify a sample of Lyman continuum galaxies (LCGs). These are individual sources that emit ionizing flux below the Lyman break ($<912~\textÅ$). This population would allow us to constrain cosmic reionization parameters such as the number density and escape fraction ($f_{\rm esc}$) of ionizing sources. We compile a comprehensive parent sample that does not rely on the Lyman-break technique for redshifts. We present three new spectroscopic candidates at $z\sim3.7\text{--}4.4$, and 32 new photometric candidates. The high-resolution multi-band HST imaging and new Keck/Low Resolution Imaging Spectrometer (LRIS) redshifts make these promising spectroscopic LCG candidates. Using both a traditional and probabilistic approach, we find the most likely $f_{\rm esc}$ values for the three spectroscopic LCG candidates are $>100\%$, and therefore not physical. We are unable to confirm the true nature of these sources with the best available imaging and direct blue Keck/LRIS spectroscopy. More spectra, especially from the new class of 30 m telescopes, will be required to build a statistical sample of LCGs to place firm observational constraints on cosmic reionization.

preprint2022arXiv

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble program. IV. Calibration of Dust Depletions vs Abundance Ratios in the Milky Way and Magellanic Clouds and Application to Damped Lyman-alpha Systems

The evolution of the metal content of the universe can be tracked through rest-frame UV spectroscopy of damped Ly-$α$ systems (DLAs). Gas-phase abundances in DLAs must be corrected for dust depletion effects, which can be accomplished by calibrating the relation between abundance ratios such as [Zn/Fe] and depletions (the fraction of metals in gas, as opposed to dust). Using samples of gas-phase abundances and depletions in the Milky Way (MW), LMC, and SMC, we demonstrate that the relation between [Zn/Fe] and other abundance ratios does not change significantly between these local galaxies and DLAs, indicating that [Zn/Fe] should trace depletions of heavy elements in those systems. The availability of photospheric abundances in young massive stars, a proxy for the total (gas+dust) metallicity of neutral gas, in the MW LMC, and SMC allows us to calibrate the relation between [Zn/Fe] and depletions in these nearby galaxies. We apply the local calibrations of depletions to DLA systems. We find that the fraction of metals in dust, the dust-to-gas-ratio, and total abundances are 2-5 times lower than inferred from previous depletion calibrations based on MW measurements and a different formalism. However, the trend of dust abundance vs. metallicity remains only slightly sub-linear for all existing depletion calibrations, contrary to what is inferred from FIR, 21 cm, and CO emission in nearby galaxies and predicted by chemical evolution models. Observational constraints on the FIR dust opacity and depletions at metallicities lower than 20\% solar will be needed to resolve this tension.

preprint2022arXiv

Quantifying Feedback from Narrow Line Region Outflows in Nearby Active Galaxies. IV. The Effects of Different Density Estimates on the Ionized Gas Masses and Outflow Rates

Active galactic nuclei (AGN) can launch outflows of ionized gas that may influence galaxy evolution, and quantifying their full impact requires spatially resolved measurements of the gas masses, velocities, and radial extents. We previously reported these quantities for the ionized narrow-line region (NLR) outflows in six low-redshift AGN, where the gas velocities and extents were determined from Hubble Space Telescope long-slit spectroscopy. However, calculating the gas masses required multi-component photoionization models to account for radial variations in the gas densities, which span $\sim$6 orders of magnitude. In order to simplify this method for larger samples with less spectral coverage, we compare these gas masses with those calculated from techniques in the literature. First, we use a recombination equation with three different estimates for the radial density profiles. These include constant densities, those derived from [S II], and power-law profiles based on constant values of the ionization parameter ($U$). Second, we use single-component photoionization models with power-law density profiles based on constant $U$, and allow $U$ to vary with radius based on the [O III]/H$β$ ratios. We find that assuming a constant density of $n_\mathrm{H} =$ 10$^2$ cm$^{-3}$ overestimates the gas masses for all six outflows, particularly at small radii where the outflow rates peak. The use of [S II] marginally matches the total gas masses, but also overestimates at small radii. Overall, single-component photoionization models where $U$ varies with radius are able to best match the gas mass and outflow rate profiles when there are insufficient emission lines to construct detailed models.

preprint2022arXiv

Tracing Lyman-alpha and Lyman Continuum Escape in Galaxies with Mg II Emission

Star-forming galaxies are considered the likeliest source of the H I ionizing Lyman Continuum (LyC) photons that reionized the intergalactic medium at high redshifts. However, above z >~ 6, the neutral intergalactic medium prevents direct observations of LyC. Therefore, recent years have seen the development of indirect indicators for LyC that can be calibrated at lower redshifts and applied in the Epoch of Reionization. Emission from Mg II \ly\ly 2796, 2803 doublet has been proposed as a promising LyC proxy. In this paper, we present new Hubble Space Telescope/Cosmic Origins Spectrograph observations for 8 LyC emitter candidates, selected to have strong Mg II emission lines. We securely detect LyC emission in 50% (4/8) galaxies with 2$σ$ significance. This high detection rate suggests that strong Mg II emitters might be more likely to leak LyC than similar galaxies without strong Mg II. Using photoionization models, we constrain the escape fraction of Mg II as ~ 15 -- 60%. We confirm that the escape fraction of Mg II correlates tightly with that of Lyman-alpha (LyA), which we interpret as an indication that the escape fraction of both species is controlled by resonant scattering in the same low column density gas. Furthermore, we show that the combination of the Mg II emission and dust attenuation can be used to estimate the escape fraction of LyC statistically. These findings confirm that Mg II emission can be adopted to estimate the escape fraction of LyA and LyC in local star-forming galaxies and may serve as a useful indirect indicator at the Epoch of Reionization.

preprint2021arXiv

Spectroscopically Identified Emission Line Galaxy Pairs in the WISP survey

We identify a sample of spectroscopically measured emission line galaxy (ELG) pairs up to z=1.6 from the WFC3 Infrared Spectroscopic Parallels (WISP) survey. WISP obtained slitless, near-infrared grism spectroscopy along with direct imaging in the J and H bands by observing in the pure-parallel mode with the Wide Field Camera Three (WFC3) on the Hubble Space Telescope (HST). From our search of 419 WISP fields covering an area of ~0.5 deg$^{2}$, we find 413 ELG pair systems, mostly Halpha emitters. We then derive reliable star formation rates (SFRs) based on the attenuation-corrected Halpha fluxes. Compared to isolated galaxies, we find an average SFR enhancement of 40%-65%, which is stronger for major pairs and pairs with smaller velocity separations (Delta_v < 300 km/s). Based on the stacked spectra from various subsamples, we study the trends of emission line ratios in pairs, and find a general consistency with enhanced lower-ionization lines. We study the pair fraction among ELGs, and find a marginally significant increase with redshift $f \propto (1+z)^α$, where the power-law index α=0.58$\pm$0.17 from $z\sim$0.2 to $z\sim$1.6. The fraction of Active galactic Nuclei (AGNs), is found to be the same in the ELG pairs as compared to isolated ELGs.

preprint2020arXiv

A Cold, Massive, Rotating Disk Galaxy 1.5 Billion Years after the Big Bang

Massive disk galaxies like the Milky Way are expected to form at late times in traditional models of galaxy formation, but recent numerical simulations suggest that such galaxies could form as early as a billion years after the Big Bang through the accretion of cold material and mergers. Observationally, it has been difficult to identify disk galaxies in emission at high redshift, in order to discern between competing models of galaxy formation. Here we report imaging, with a resolution of about 1.3 kiloparsecs, of the 158-micrometre emission line from singly ionized carbon, the far-infrared dust continuum and the near-ultraviolet continuum emission from a galaxy at a redshift of 4.2603, identified by detecting its absorption of quasar light. These observations show that the emission arises from gas inside a cold, dusty, rotating disk with a rotational velocity of 272 kilometres per second. The detection of emission from carbon monoxide in the galaxy yields a molecular mass that is consistent with the estimate from the ionized carbon emission of about 72 billion solar masses. The existence of such a massive, rotationally supported, cold disk galaxy when the Universe was only 1.5 billion years old favours formation through either cold-mode accretion or mergers, although its large rotational velocity and large content of cold gas remain challenging to reproduce with most numerical simulations.

preprint2020arXiv

ALMA [CII] 158um imaging of an HI-Selected Major Merger at z~4

We present high spatial-resolution (~2kpc) Atacama Large Millimeter/submillimeter Array (ALMA) observations of [CII] 158um and dust-continuum emission from a galaxy at z=3.7978 selected by its strong HI absorption (a damped Lya absorber, DLA) against a background QSO. Our ALMA images reveal a pair of star-forming galaxies separated by ~6kpc (projected) undergoing a major merger. Between these galaxies is a third emission component with highly elevated (2x) [CII] 158um emission relative to the dust continuum, which is likely to arise from stripped gas associated with the merger. This merger of two otherwise-normal galaxies is not accompanied by enhanced star-formation, contrary to mergers detected in most luminosity-selected samples. The DLA associated with the merger exhibits extreme kinematics, with a velocity width for the low-ionization metal lines of Dv90~470km/s that spans the velocity spread revealed in the [CII] 158um emission. We propose that DLAs with high Dv90 values are a signpost of major mergers in normal galaxies at high redshifts, and use the distribution of the velocity widths of metal lines in high-z DLAs to provide a rough estimate the fraction of z>3 galaxies that are undergoing a major merger.

preprint2020arXiv

Discovery of a damped Ly$α$ galaxy at z $\sim$ 3 towards the quasar SDSS J011852+040644

We report the detection of the host galaxy of a damped Ly$α$ system (DLA) with log N(HI) $ [\rm cm^{-2}]$ = $21.0 \pm 0.10$ at $z \approx 3.0091$ towards the background quasar SDSS J011852+040644 using the Palomar Cosmic Web Imager (PCWI) at the Hale (P200) telescope. We detect Ly$α$ emission in the dark core of the DLA trough at a 3.3$σ$ confidence level, with Ly$α$ luminosity of $L_{\rm Lyα}$ $\rm = (3.8 \pm 0.8) \times 10^{42}\ erg\ s^{-1}$, corresponding to a star formation rate of $\gtrsim 2\ \rm M_{\odot}\ yr^{-1}$ (considering a lower limit on Ly$α$ escape fraction $f_{esc}^{Lyα} \sim 2\%$) as typical for Lyman break galaxies at these redshifts. The Ly$α$ emission is blueshifted with respect to the systemic redshift derived from metal absorption lines by $281 \pm 43$ km/s. The associated galaxy is at very small impact parameter of $\lesssim 12 \rm\ kpc$ from the background quasar, which is in line with the observed anticorrelation between column density and impact parameter in spectroscopic searches tracing the large-scale environments of DLA host galaxies.

preprint2020arXiv

Into the Lyα jungle: exploring the circumgalactic medium of galaxies at z ~ 4-5 with MUSE

We present a study of the galaxy environment of 9 strong HI+CIV absorption line systems ($16.2<{\rm log}(N({\rm HI}))<21.2$) spanning a wide range in metallicity at $z\sim4-5$, using MUSE integral field and X-Shooter spectroscopic data collected in a $z\approx 5.26$ quasar field. We identify galaxies within a 250 kpc and $\pm1000$ km s$^{-1}$ window for 6 out of the 9 absorption systems, with 2 of the absorption line systems showing multiple associated galaxies within the MUSE field of view. The space density of Ly$α$ emitting galaxies (LAEs) around the HI and CIV systems is $\approx10-20$ times the average sky density of LAEs given the flux limit of our survey, showing a clear correlation between the absorption and galaxy populations. Further, we find that the strongest CIV systems in our sample are those that are most closely aligned with galaxies in velocity space, i.e. within velocities of $\pm500$ km s$^{-1}$. The two most metal poor systems lie in the most dense galaxy environments, implying we are potentially tracing gas that is infalling for the first time into star-forming groups at high redshift. Finally, we detect an extended Ly$α$ nebula around the $z\approx 5.26$ quasar, which extends up to $\approx50$ kpc at the surface brightness limit of $3.8 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. After scaling for surface brightness dimming, we find that this nebula is centrally brighter, having a steeper radial profile than the average for nebulae studied at $z\sim3$ and is consistent with the mild redshift evolution seen from $z\approx 2$.

preprint2016arXiv

Evolution of Intrinsic Scatter in the SFR-Stellar Mass Correlation at 0.5<z<3

We present estimates of intrinsic scatter in the Star Formation Rate (SFR) - Stellar Mass (M*) correlation in the redshift range 0.5 < z < 3.0 and in the mass range 10^7 < M* < 10^11 Msun. We utilize photometry in the Hubble Ultradeep Field (HUDF12), Ultraviolet Ultra Deep Field (UVUDF) campaigns and CANDELS/GOODS-S. We estimate SFR, M* from broadband Spectral Energy Distributions (SEDs) and the best available redshifts. The maximum depth of the HUDF photometry (F160W 29.9 AB, 5 sigma depth) probes the SFR-M* correlation down to M* ~ 10 ^7 Msun, a factor of 10-100X lower in M* than previous studies, and comparable to dwarf galaxies in the local universe. We find the slope of the SFR-M* relationship to be near unity at all redshifts and the normalization to decrease with cosmic time. We find a moderate increase in intrinsic scatter with cosmic time from 0.2 to 0.4 dex across the epoch of peak cosmic star formation. None of our redshift bins show a statistically significant increase in intrinsic scatter at low mass. However, it remains possible that intrinsic scatter increases at low mass on timescales shorter than ~ 100 Myr. Our results are consistent with a picture of gradual and self-similar assembly of galaxies across more than three orders of magnitude in stellar mass from as low as 10^7 Msun.

preprint2016arXiv

Exploring Damped Lyman-$α$ System Host Galaxies using Gamma-ray Bursts

We present a sample of 45 Damped Lyman-$α$ system (DLA, \nhi\ $ \geq 2 \times 10^{20} {\rm cm}^{-2}$) counterparts (33 detections, 12 upper limits) which host gamma-ray bursts (GRB-DLAs) in order to investigate star-formation and metallicity within galaxies hosting DLAs. Our sample spans $z \sim 2-6$ and is nearly three times larger than any previously detected DLA counterparts survey based on quasar line-of-sight searches (QSO-DLAs). We report star formation rates (SFRs) from rest-frame UV photometry and SED modeling. We find that DLA counterpart SFRs are not correlated with either redshift or HI column density. Thanks to the combination of \hst\ and ground-based observations, we also investigate DLA host star-formation efficiency. Our GRB-DLA counterpart sample spans both higher efficiency and low efficiency star formation regions compared to the local Kennicutt-Schmidt relation, local star formation laws, and $z \sim 3$ cosmological simulations. We also compare the depletion times of our DLA hosts sample to other objects in the local Universe, our sample appears to deviate from the star formation efficiencies measured in local spiral and dwarf galaxies. Furthermore, we find similar efficiencies as local inner disks, SMC, and LBG outskirts. Finally, our enrichment time measurements show a spread of systems with under- and over-abundance of metals which may suggest that these systems had episodic star formation and a metal enrichment/depletion as a result of strong stellar feedback and/or metal inflow/outflow.

preprint2016arXiv

Hubble Space Telescope Observations of Accretion-Induced Star Formation in the Tadpole Galaxy Kiso 5639

The tadpole galaxy Kiso 5639 has a slowly rotating disk with a drop in metallicity at its star-forming head, suggesting that star formation was triggered by the accretion of metal-poor gas. We present multi-wavelength HST WFC3 images of UV through I band plus Halpha to search for peripheral emission and determine the properties of various regions. The head has a mass in young stars of ~10^6 Mo and an ionization rate of 6.4x10^51 s^{-1}, equivalent to ~2100 O9-type stars. There are four older star-forming regions in the tail, and an underlying disk with a photometric age of ~1 Gyr. The mass distribution function of 61 star clusters is a power law with a slope of -1.73+-0.51. Fourteen young clusters in the head are more massive than 10^4 Mo, suggesting a clustering fraction of 30%-45%. Wispy filaments of Halpha emission and young stars extend away from the galaxy. Shells and holes in the head HII region could be from winds and supernovae. Gravity from the disk should limit the expansion of the HII region, although hot gas might escape through the holes. The star formation surface density determined from Halpha in the head is compared to that expected from likely pre-existing and accreted gas. Unless the surface density of the accreted gas is a factor of ~3 or more larger than what was in the galaxy before, the star formation rate has to exceed the usual Kennicutt-Schmidt rate by a factor of >5.

preprint2016arXiv

Observational evidence of a slow downfall of star formation efficiency in massive galaxies during the last 10 Gyr

In this paper we study the causes of the reported mass-dependence of the slope of SFR-M* relation, the so-called "Main Sequence" of star-forming galaxies, and discuss its implication on the physical processes that shaped the star formation history of massive galaxies over cosmic time. We use the CANDELS near-IR imaging from the Hubble Space Telescope to perform the bulge-to-disk decomposition of distant galaxies and measure for the first time the slope of the SFR-Mdisk relation at z=1. We find that this relation follows very closely the shape of the SFR-M* correlation, still with a pronounced flattening at the high-mass end. This is clearly excluding, at least at z=1, the secular growth of quiescent bulges in star-forming galaxies as the main driver for the change of slope of the Main Sequence. Then, by stacking the Herschel data available in the CANDELS field, we estimate the total gas mass and the star formation efficiency at different positions on the SFR-M* relation. We find that the relatively low SFRs observed in massive galaxies (M* > 5e10 Msun) are caused by a decreased star formation efficiency, by up to a factor of 3 as compared to lower stellar mass galaxies, and not by a reduced gas content. The trend at the lowest masses is likely linked to the dominance of atomic over molecular gas. We argue that this stellar-mass-dependent SFE can explain the varying slope of the Main Sequence since z=1.5, hence over 70% of the Hubble time. The drop of SFE occurs at lower masses in the local Universe (M* > 2e10 Msun) and is not present at z=2. Altogether this provides evidence for a slow downfall of the star formation efficiency in massive Main Sequence galaxies. The resulting loss of star formation is found to be rising starting from z=2 to reach a level comparable to the mass growth of the quiescent population by z=1. We finally discuss the possible physical origin of this phenomenon.

preprint2016arXiv

Stellar Mass--Gas-phase Metallicity Relation at $0.5\leq z\leq0.7$: A Power Law with Increasing Scatter toward the Low-mass Regime

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies.

preprint2016arXiv

The Bursty Star Formation Histories of Low-mass Galaxies at $0.4<z<1$ Revealed by Star Formation Rates Measured from H$β$ and FUV

We investigate the burstiness of star formation histories (SFHs) of galaxies at $0.4<z<1$ by using the ratio of star formation rates (SFRs) measured from H$β$ and FUV (1500 Å) (H$β$--to--FUV ratio). Our sample contains 164 galaxies down to stellar mass (M*) of $10^{8.5} M_\odot$ in the CANDELS GOODS-N region, where Team Keck Redshift Survey DEIMOS spectroscopy and HST/WFC3 F275W images from CANDELS and Hubble Deep UV Legacy Survey are available. When the {\it ratio} of H$β$- and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H$β$--to--FUV ratio of our sample increases with M* and SFR. The median ratio is $\sim$0.7 at M*$\sim10^{8.5} M_\odot$ (or SFR$\sim 0.5 M_\odot/yr$) and increases to $\sim$1 at M*$\sim10^{10} M_\odot$ (or SFR $\sim 10 M_\odot/yr$). At M*$<10^{9.5} M_\odot$, our median H$β$--to--FUV ratio is lower than that of local galaxies at the same M*, implying a redshift evolution. Bursty SFH on a timescale of a few tens of megayears on galactic scales provides a plausible explanation of our results, and the importance of the burstiness increases as M* decreases. Due to sample selection effects, our H$β$--to--FUV ratio may be an upper limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.

preprint2016arXiv

The Evolution Of The Faint End Of The UV Luminosity Function During The Peak Epoch Of Star Formation (1<z<3)

[Abridged] We present a robust measurement of the rest-frame UV luminosity function (LF) and its evolution during the peak epoch of cosmic star formation at 1<z<3. We use our deep near ultraviolet imaging from WFC3/UVIS on the Hubble Space Telescope (HST) and existing ACS/WFC and WFC3/IR imaging of three lensing galaxy clusters, Abell 2744 and MACSJ0717 from the Hubble Frontier Field survey and Abell 1689. We use photometric redshifts to identify 780 ultra-faint galaxies with $M_{UV}$<-12.5 AB mag at 1<z<3. From these samples, we identified 5 new, faint, multiply imaged systems in A1689. We compute the rest-frame UV LF and find the best-fit faint-end slopes of $α=-1.56\pm0.04$, $α=-1.72\pm0.04$ and $α=-1.94\pm0.06$ at 1.0<z<1.6, 1.6<z<2.2 and 2.2<z<3.0, respectively. Our results demonstrate that the UV LF becomes steeper from z\sim1.3 to z\sim2.6 with no sign of a turnover down to $M_{UV}=-14$ AB mag. We further derive the UV LFs using the Lyman break "dropout" selection and confirm the robustness of our conclusions against different selection methodologies. Because the sample sizes are so large, and extend to such faint luminosities, the statistical uncertainties are quite small, and systematic uncertainties (due to the assumed size distribution, for example), likely dominate. If we restrict our analysis to galaxies and volumes above > 50% completeness in order to minimize these systematics, we still find that the faint-end slope is steep and getting steeper with redshift, though with slightly shallower (less negative) values ($α=-1.55\pm0.06$, $-1.69\pm0.07$ and $-1.79\pm0.08$ for $z\sim1.3$, 1.9 and 2.6, respectively). Finally, we conclude that the faint star-forming galaxies with UV magnitudes of $-18.5<M_{UV}<-12.5$ covered in this study, produce the majority (55%-60%) of the unobscured UV luminosity density at 1<z<3.

preprint2016arXiv

The HI Content of the Universe over the Past 10 Gyrs

We use the Hubble Space Telescope (HST) archive of ultraviolet (UV) quasar spectroscopy to conduct the first blind survey for damped Ly-alpha absorbers (DLAs) at low redshift (z < 1.6). Our statistical sample includes 463 quasars with spectral coverage spanning a total redshift path, dz = 123.3 or an absorption path, dX = 229.7. Within this survey path, we identify 4 DLAs, defined as absorbers with HI column density N(HI) >= 10^20.3cm-2, which implies an incidence per absorption length, l(X)= 0.017(+0.014-0.008) at a median survey path redshift of z=0.623. While our estimate of l(X) is lower than earlier estimates at z ~ 0 from HI 21cm emission studies, the results are consistent within the measurement uncertainties. Our dataset is too small to properly sample the N(HI) frequency distribution function f(N(HI),X), but the observed distribution agrees with previous estimates at z > 2. Adopting the z > 2 shape of f(N(HI),X), we infer an HI mass density at z ~ 0.6 of rho_HI = 0.25(+0.20-0.12) x 10^8 Msol Mpc-3. This is significantly lower than previous estimates from targeted DLA surveys with the HST, but consistent with results from low-z HI 21cm observations, and suggests that the neutral gas density of the universe has been decreasing over the past 10 Gyrs.

preprint2016arXiv

The Star Formation Rate Efficiency of Neutral Atomic-dominated Hydrogen Gas in the Outskirts of Star Forming Galaxies from z~1 to z~3

Current observational evidence suggests that the star formation rate (SFR) efficiency of neutral atomic hydrogen gas measured in Damped Ly-alpha Systems (DLAs) at z~3 is more than 10 times lower than predicted by the Kennicutt-Schmidt (KS) relation. To understand the origin of this deficit, and to investigate possible evolution with redshift and galaxy properties, we measure the SFR efficiency of atomic gas at z~1, z~2, and z~3 around star-forming galaxies. We use new robust photometric redshifts in the Hubble Ultra Deep Field to create galaxy stacks in these three redshift bins, and measure the SFR efficiency by combining DLA absorber statistics with the observed rest-frame UV emission in the galaxies' outskirts. We find that the SFR efficiency of HI gas at z>1 is ~1-3% of that predicted by the KS relation. Contrary to simulations and models that predict a reduced SFR efficiency with decreasing metallicity and thus with increasing redshift, we find no significant evolution in the SFR efficiency with redshift. Our analysis instead suggests that the reduced SFR efficiency is driven by the low molecular content of this atomic-dominated phase, with metallicity playing a secondary effect in regulating the conversion between atomic and molecular gas. This interpretation is supported by the similarity between the observed SFR efficiency and that observed in local atomic-dominated gas, such as in the outskirts of local spiral galaxies and local dwarf galaxies.

preprint2015arXiv

Clumpy Galaxies in CANDELS. I. The Definition of UV Clumps and the Fraction of Clumpy Galaxies at 0.5<z<3

Although giant clumps of stars are crucial to galaxy formation and evolution, the most basic demographics of clumps are still uncertain, mainly because the definition of clumps has not been thoroughly discussed. In this paper, we study the basic demographics of clumps in star-forming galaxies (SFGs) at 0.5<z<3, using our proposed physical definition that UV-bright clumps are discrete star-forming regions that individually contribute more than 8% of the rest-frame UV light of their galaxies. Clumps defined this way are significantly brighter than the HII regions of nearby large spiral galaxies, either individually or blended, when physical spatial resolution and cosmological dimming are considered. Under this definition, we measure the fraction of SFGs that contain at least one off-center clump (Fclumpy) and the contributions of clumps to the rest-frame UV light and star formation rate of SFGs in the CANDELS/GOODS-S and UDS fields, where our mass-complete sample consists of 3239 galaxies with axial ratio q>0.5. The redshift evolution of Fclumpy changes with the stellar mass (M*) of the galaxies. Low-mass (log(M*/Msun)<9.8) galaxies keep an almost constant Fclumpy of about 60% from z~3.0 to z~0.5. Intermediate-mass and massive galaxies drop their Fclumpy from 55% at z~3.0 to 40% and 15%, respectively, at z~0.5. We find that (1) the trend of disk stabilization predicted by violent disk instability matches the Fclumpy trend of massive galaxies; (2) minor mergers are a viable explanation of the Fclumpy trend of intermediate-mass galaxies at z<1.5, given a realistic observability timescale; and (3) major mergers are unlikely responsible for the Fclumpy trend in all masses at z<1.5. The clump contribution to the rest-frame UV light of SFGs shows a broad peak around galaxies with log(M*/Msun)~10.5 at all redshifts, possibly linked to the molecular gas fraction of the galaxies. (Abridged)

preprint2015arXiv

Predicting the redshift 2 Halpha luminosity function using [OIII] emission line galaxies

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

preprint2015arXiv

Star-forming blue ETGs in two newly discovered galaxy overdensities in the HUDF at z=1.84 and 1.9: unveiling the progenitors of passive ETGs in cluster cores

We present the discovery of two galaxy overdensities in the Hubble Space Telescope UDF: a proto-cluster, HUDFJ0332.4-2746.6 at $z = 1.84 \pm 0.01$, and a group, HUDFJ0332.5-2747.3 at $z =1.90 \pm 0.01$. Assuming viralization, the velocity dispersion of HUDFJ0332.4-2746.6 implies a mass of $M_{200}= (2.2 \pm 1.8) \times 10^{14} M_{\odot}$, consistent with the lack of extended X-ray emission. Neither overdensity shows evidence of a red sequence. About $50\%$ of their members show interactions and/or disturbed morphologies, which are signatures of merger remnants or disk instability. Most of their ETGs have blue colors and show recent star-formation. These observations reveal for the first time large fractions of spectroscopically confirmed star-forming blue ETGs in proto-clusters at $z\approx 2$. These star-forming ETGs are most likely among the progenitors of the quiescent population in clusters at more recent epochs. Their mass-size relation is consistent with that of passive ETGs in clusters at $z\sim0.7-1.5$. If these galaxies are the progenitors of cluster ETGs at these lower redshifts, their size would evolve according to a similar mass-size relation. It is noteworthy that quiescent ETGs in clusters at $z=1.8-2$ also do not show any significant size evolution over this redshift range, contrary to field ETGs. The ETG fraction is $\lesssim 50\%$, compared to the typical quiescent ETG fraction of $\approx 80\%$ in cluster cores at $z< 1$. The fraction, masses, and colors of the newly discovered ETGs imply that other cluster ETGs will be formed/accreted at a later time.

preprint2015arXiv

The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

preprint2015arXiv

The Lyman Continuum Escape Fraction of Low-Mass Star-Forming Galaxies at z~1

To date no direct detection of Lyman continuum emission has been measured for intermediate--redshift z~1 star-forming galaxies . We combine HST grism spectroscopy with GALEX UV and ground--based optical imaging to extend the search for escaping Lyman continuum to a large (~600) sample of z~1 low-mass, moderately star-forming galaxies selected initially on H$α$ emission. The characteristic escape fraction of LyC from SFGs that populate this parameter space remains weakly constrained by previous surveys, but these faint SFGs are assumed to play a significant role in the reionization of neutral hydrogen in the intergalactic medium (IGM) at high redshift (z>6). We do not make an unambiguous detection of escaping LyC radiation from this $z\sim1$ sample, individual non--detections to constrain the absolute Lyman continuum escape fraction, $f_{esc}$<2.1% (3$σ$). We measure upper limits of $f_{esc}$<9.6% from a sample of SFGs selected on high H$α$ equivalent width (EW>200Å), which are thought to be close analogs of high redshift sources of reionization. For reference, we also present an emissivity--weighted escape fraction which is useful as a measurement of the general contribution of the SFGs to the z~1 ionizing UV background. In the discussion, we consider the implications of these intermediate redshift constraints for the reionization of hydrogen in the intergalactic medium at high (z > 6) redshift. If we assume our $z\sim1$ SFGs, for which we measure this emissivity-weighted $f_{esc}$, are analogs to the high redshift sources of reionization, we find is difficult reconcile reionization by faint (M$_{UV}<-13$) SFGs with a low escape fraction ($f_{esc}<$3%), with constraints from independent high redshift observations. If $f_{esc}$ evolves with redshift, reionization by SFGs may be consistent with observations from Planck.

preprint2015arXiv

The Neutral Hydrogen Cosmological Mass Density at z=5

We present the largest homogeneous survey of $z>4.4$ damped Lyman-$α$ systems (DLAs) using the spectra of 163 QSOs that comprise the Giant Gemini GMOS (GGG) survey. With this survey we make the most precise high-redshift measurement of the cosmological mass density of neutral hydrogen, $Ω_{\rm HI}$. At such high redshift important systematic uncertainties in the identification of DLAs are produced by strong intergalactic medium absorption and QSO continuum placement. These can cause spurious DLA detections, result in real DLAs being missed, or bias the inferred DLA column density distribution. We correct for these effects using a combination of mock and higher-resolution spectra, and show that for the GGG DLA sample the uncertainties introduced are smaller than the statistical errors on $Ω_{\rm HI}$. We find $Ω_{\rm HI}=0.98^{+0.20}_{-0.18}\times10^{-3}$ at $\langle z\rangle=4.9$, assuming a 20% contribution from lower column density systems below the DLA threshold. By comparing to literature measurements at lower redshifts, we show that $Ω_{\rm HI}$ can be described by the functional form $Ω_{\rm HI}(z)\propto(1+z)^{0.4}$. This gradual decrease from $z=5$ to $0$ is consistent with the bulk of HI gas being a transitory phase fuelling star formation, which is continually replenished by more highly-ionized gas from the intergalactic medium, and from recycled galactic winds.

preprint2015arXiv

UVUDF: Ultraviolet Through Near-infrared Catalog and Photometric Redshifts of Galaxies in the Hubble Ultra Deep Field

We present photometry and derived redshifts from up to eleven bandpasses for 9927 galaxies in the Hubble Ultra Deep field (UDF), covering an observed wavelength range from the near-ultraviolet (NUV) to the near-infrared (NIR) with Hubble Space Telescope observations. Our Wide Field Camera 3 (WFC3)/UV F225W, F275W, and F336W image mosaics from the ultra-violet UDF (UVUDF) imaging campaign are newly calibrated to correct for charge transfer inefficiency, and use new dark calibrations to minimize background gradients and pattern noise. Our NIR WFC3/IR image mosaics combine the imaging from the UDF09 and UDF12 campaigns with CANDELS data to provide NIR coverage for the entire UDF field of view. We use aperture-matched point-spread function corrected photometry to measure photometric redshifts in the UDF, sampling both the Lyman break and Balmer break of galaxies at z~0.8-3.4, and one of the breaks over the rest of the redshift range. Our comparison of these results with a compilation of robust spectroscopic redshifts shows an improvement in the galaxy photometric redshifts by a factor of two in scatter and a factor three in outlier fraction over previous UDF catalogs. The inclusion of the new NUV data is responsible for a factor of two decrease in the outlier fraction compared to redshifts determined from only the optical and NIR data, and improves the scatter at z<0.5 and at z>2. The panchromatic coverage of the UDF from the NUV through the NIR yields robust photometric redshifts of the UDF, with the lowest outlier fraction available.

preprint2014arXiv

Directly imaging damped Ly-alpha galaxies at z>2. III: The star formation rates of neutral gas reservoirs at z~2.7

We present results from a survey designed to probe the star formation properties of 32 damped Ly-alpha systems (DLAs) at z~2.7. By using the "double-DLA" technique that eliminates the glare of the bright background quasars, we directly measure the rest-frame FUV flux from DLAs and their neighbouring galaxies. At the position of the absorbing gas, we place stringent constraints on the unobscured star formation rates (SFRs) of DLAs to 2-sigma limits of <0.09-0.27 M/yr, corresponding to SFR surface densities <10^(-2.6)-10^(-1.5) M/yr/kpc^2. The implications of these limits for the star formation law, metal enrichment, and cooling rates of DLAs are examined. By studying the distribution of impact parameters as a function of SFRs for all the galaxies detected around these DLAs, we place new direct constraints on the bright end of the UV luminosity function of DLA hosts. We find that <13% of the hosts have SFR>2 M/yr at impact parameters b < (SFR/(M/yr))^(0.8)+6 kpc, differently from current samples of confirmed DLA galaxies. Our observations also disfavor a scenario in which the majority of DLAs arise from bright LBGs at distances 20 < b < 100 kpc. These new findings corroborate a picture in which DLAs do not originate from highly star forming systems that are coincident with the absorbers, and instead suggest that DLAs are associated with faint, possibly isolated, star-forming galaxies. Potential shortcomings of this scenario and future strategies for further investigation are discussed.

preprint2014arXiv

Hubble Space Telescope Grism Spectroscopy of Extreme Starbursts Across Cosmic Time: The Role of Dwarf Galaxies in the Star Formation History of the Universe

Near infrared slitless spectroscopy with the Wide Field Camera 3, onboard the Hubble Space Telescope, offers a unique opportunity to study low-mass galaxy populations at high-redshift ($z\sim$1-2). While most high$-z$ surveys are biased towards massive galaxies, we are able to select sources via their emission lines that have very-faint continua. We investigate the star formation rate (SFR)-stellar mass ($M_{\star}$) relation for about 1000 emission-line galaxies identified over a wide redshift range of $0.3 \lesssim z \lesssim 2.3$. We use the H$_α$ emission as an accurate SFR indicator and correct the broadband photometry for the strong nebular contribution to derive accurate stellar masses down to $M_{\star} \sim 10^{7} M_{\odot}$. We focus here on a subsample of galaxies that show extremely strong emission lines (EELGs) with rest-frame equivalent widths ranging from 200 to 1500 Å. This population consists of outliers to the normal SFR-$M_{\star}$ sequence with much higher specific SFRs ($> 10$ Gyr$^{-1}$). While on-sequence galaxies follow a continuous star formation process, EELGs are thought to be caught during an extreme burst of star formation that can double their stellar mass in less than $100$ Myr. The contribution of starbursts to the total star formation density appears to be larger than what has been reported for more massive galaxies in previous studies. In the complete mass range $8.2 <$ log($M_{\star}/M_{\odot}$) $< 10$ and a SFR lower completeness limit of about 2 $M_{\odot}$ yr$^{-1}$ (10 $M_{\odot}$ yr$^{-1}$) at $z\sim1$ ($z \sim 2$), we find that starbursts having EW$_{rest}$(H$_α$)$>$ 300, 200, and 100 A contribute up to $\sim13$, 18, and 34 %, respectively, to the total SFR of emission-line selected sample at $z\sim1-2$. The comparison with samples of massive galaxies shows an increase in the contribution of starbursts towards lower masses.

preprint2014arXiv

Optical-faint, Far-infrared-bright Herschel Sources in the CANDELS Fields: Ultra-Luminous Infrared Galaxies at z>1 and the Effect of Source Blending

The Herschel very wide-field surveys have charted hundreds of square degrees in multiple far-IR (FIR) bands. While the Sloan Digital Sky Survey (SDSS) is currently the best resource for optical counterpart identifications over such wide areas, it does not detect a large number of Herschel FIR sources and leaves their nature undetermined. As a test case, we studied seven "SDSS-invisible", very bright 250um sources (S_{250} > 55 mJy) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) fields where we have a rich multi-wavelength data set. We took a new approach to decompose the FIR sources, using the near-IR or the optical images directly for position priors. This is an improvement over the previous decomposition efforts where the priors are from mid-IR data that still suffer from the source blending problem in the first place. We found that in most cases the single Herschel sources are made of multiple components that are not necessarily at the same redshifts. Our decomposition succeeded in identifying and extracting their major contributors. We show that these are all ULIRGs at z ~ 1--2 whose high L_IR is mainly due to dust-obscured star formation. Most of them would not be selected as sub-mm galaxies. They all have complicated morphologies indicative of merger or violent instability, and their stellar populations are heterogeneous in terms of stellar masses, ages and formation histories. Their current ULIRG phases are of various degrees of importance in their stellar mass assembly. Our practice provides a promising starting point to develop an automatic routine to reliably study bright Herschel sources.

preprint2014arXiv

Physical Properties of Emission-Line Galaxies at z ~ 2 from Near-Infrared Spectroscopy with Magellan FIRE

We present results from near-infrared spectroscopy of 26 emission-line galaxies at z ~ 2 obtained with the FIRE spectrometer on the Magellan Baade telescope. The sample was selected from the WISP survey, which uses the near-infrared grism of the Hubble Space Telescope Wide Field Camera 3 to detect emission-line galaxies over 0.3 < z < 2.3. Our FIRE follow-up spectroscopy (R~5000) over 1.0-2.5 micron permits detailed measurements of physical properties of the z~2 emission-line galaxies. Dust-corrected star formation rates for the sample range from ~5-100 M_sun yr-1. We derive a median metallicity for the sample of ~0.45 Z_sun, and the estimated stellar masses range from ~10^8.5 - 10^9.5 M_sun. The average ionization parameters measured for the sample are typically much higher than what is found for local star-forming galaxies. We derive composite spectra from the FIRE sample, from which we infer typical nebular electron densities of ~100-400 cm^-3. Based on the location of the galaxies and composite spectra on BPT diagrams, we do not find evidence for significant AGN activity in the sample. Most of the galaxies as well as the composites are offset in the BPT diagram toward higher [O III]/H-beta at a given [N II]/H-alpha, in agreement with other observations of z > 1 star-forming galaxies, but composite spectra derived from the sample do not show an appreciable offset from the local star-forming sequence on the [O III]/H-beta versus [S II]/H-alpha diagram. We infer a high nitrogen-to-oxygen abundance ratio from the composite spectrum, which may contribute to the offset of the high-redshift galaxies from the local star-forming sequence in the [O III]/H-beta versus [N II]/H-alpha diagram. We speculate that the elevated nitrogen abundance could result from substantial numbers of Wolf-Rayet stars in starbursting galaxies at z~2. (Abridged)

preprint2014arXiv

The Rapid Decline in Metallicity of Damped Ly-alpha Systems at z~5

We present evidence that the cosmological mean metallicity of neutral atomic hydrogen gas shows a sudden decrease at z>4.7 down to <Z>=-2.03(+0.09-0.11), which is 6 sigma deviant from that predicted by a linear fit to the data at lower redshifts. This measurement is made possible by the chemical abundance measurements of 8 new damped Ly-alpha (DLA) systems at z>4.7 observed with the Echellette Spectrograph and Imager on the Keck II telescope, doubling the number of measurements at z>4.7 to 16. Possible explanations for this sudden decrease in metallicity include a change in the physical processes that enrich the neutral gas within disks, or an increase of the covering factor of neutral gas outside disks due to a lower ultra-violet radiation field and higher density at high redshift. The later possibility would result in a new population of presumably lower metallicity DLAs, with an increased contribution to the DLA population at higher redshifts resulting in a reduced mean metallicity. Furthermore, we provide evidence of a possible decrease at z>4.7 in the comoving metal mass density of DLAs, which is flat out to z~4.3. Such a decrease is expected, as otherwise most of the metals from star-forming galaxies would reside in DLAs by z~6. While the metallicity is decreasing at high redshift, the contribution of DLAs to the total metal budget of the universe increases with redshift, with DLAs at z~4.3 accounting for ~20% as many metals as produced by Lyman break galaxies.

preprint2014arXiv

The UV Continuum of z > 1 Star-forming Galaxies in the Hubble Ultraviolet UltraDeep Field

We estimate the UV continuum slope, beta, for 923 galaxies in the range 1 < z < 8 in the Hubble Ultradeep Field (HUDF). These data include 460 galaxies at 1 < z < 2 down to an absolute magnitude M_{UV} = -14 (~0.006 L*_{z=1}; 0.02 L*_{z=0}), comparable to dwarf galaxies in the local universe. We combine deep HST/UVIS photometry in F225W, F275W, F336W wavebands (UVUDF) with recent data from HST/WFC3/IR (HUDF12). Galaxies in the range 1 < z < 2 are significantly bluer than local dwarf galaxies. We find their mean (median) values <beta> = -1.382 (-1.830) +/- 0.002 (random) +/- 0.1 (systematic). We find comparable scatter in beta (standard deviation = 0.43) to local dwarf galaxies and 30% larger scatter than z > 2 galaxies. We study the trends of beta with redshift and absolute magnitude for binned sub-samples and find a modest color-magnitude relation, dbeta/dM = -0.11 +/- 0.01 and no evolution in dbeta/dM with redshift. A modest increase in dust reddening with redshift and luminosity, Delta E(B-V) ~ 0.1, and a comparable increase in the dispersion of dust reddening at z < 2, appears likely to explain the observed trends. At z > 2, we find trends that are consistent with previous works; combining our data with the literature in the range 1 < z < 8, we find a color evolution with redshift, dbeta/dz = -0.09 +/-0.01 for low luminosity (0.05 L*_{z=3}), and dbeta/dz = -0.06 +/-0.01 for medium luminosity (0.25 L*_{z=3}) galaxies.

preprint2013arXiv

Low Masses and High Redshifts: The Evolution of the Mass-Metallicity Relation

We present the first robust measurement of the high redshift mass-metallicity (MZ) relation at 10^{8}< M/M_{\sun} < 10^{10}, obtained by stacking spectra of 83 emission-line galaxies with secure redshifts between 1.3 < z < 2.3. For these redshifts, infrared grism spectroscopy with the Hubble Space Telescope Wide Field Camera 3 is sensitive to the R23 metallicity diagnostic: ([OII]3726,3729 + [OIII] 4959,5007)/Hβ. Using spectra stacked in four mass quartiles, we find a MZ relation that declines significantly with decreasing mass, extending from 12+log(O/H) = 8.8 at M=10^{9.8} M_{\sun} to 12+log(O/H)= 8.2 at M=10^{8.2} M_{\sun}. After correcting for systematic offsets between metallicity indicators, we compare our MZ relation to measurements from the stacked spectra of galaxies with M>10^{9.5} M_{\sun} and z~2.3. Within the statistical uncertainties, our MZ relation agrees with the z~2.3 result, particularly since our somewhat higher metallicities (by around 0.1 dex) are qualitatively consistent with the lower mean redshift z=1.76 of our sample. For the masses probed by our data, the MZ relation shows a steep slope which is suggestive of feedback from energy-driven winds, and a cosmological downsizing evolution where high mass galaxies reach the local MZ relation at earlier times. In addition, we show that our sample falls on an extrapolation of the star-forming main sequence (the SFR-M_{*} relation) at this redshift. This result indicates that grism emission-line selected samples do not have preferentially high SFRs. Finally, we report no evidence for evolution of the mass-metallicity-SFR plane; our stack-averaged measurements show excellent agreement with the local relation.

preprint2013arXiv

Predicting Future Space Near-IR Grism Surveys using the WFC3 Infrared Spectroscopic Parallels Survey

We present near-infrared emission line counts and luminosity functions from the HST WFC3 Infrared Spectroscopic Parallels (WISP) program for 29 fields (0.037 deg^2) observed using both the G102 and G141 grisms. Altogether we identify 1048 emission line galaxies with observed equivalent widths greater than 40 Angstroms, 467 of which have multiple detected emission lines. The WISP survey is sensitive to fainter flux levels (3-5x10^-17 ergs/s/cm^2) than the future space near-infrared grism missions aimed at baryonic acoustic oscillation cosmology (1-4x10^-16 ergs/s/cm^2), allowing us to probe the fainter emission line galaxies that the shallower future surveys may miss. Cumulative number counts of 0.7<z<1.5 galaxies reach 10,000 deg^-2 above an H-alpha flux of 2x10^-16 ergs/s/cm^2. H-alpha-emitting galaxies with comparable [OIII] flux are roughly 5 times less common than galaxies with just H-alpha emission at those flux levels. Galaxies with low H-alpha/[OIII] ratios are very rare at the brighter fluxes that future near-infrared grism surveys will probe; our survey finds no galaxies with H-alpha/[OIII] < 0.95 that have H-alpha flux greater than 3x10^-16 ergs/s/cm^2. Our H-alpha luminosity function contains a comparable number density of faint line emitters to that found by the NICMOS near-infrared grism surveys, but significantly fewer (factors of 3-4 less) high luminosity emitters. We also find that our high redshift (z=0.9-1.5) counts are in agreement with the high redshift (z=1.47) narrow band H-alpha survey of HiZELS (Sobral et al. 2013), while our lower redshift luminosity function (z=0.3-0.9) falls slightly below their z=0.84 result. The evolution in both the H-alpha luminosity function from z=0.3--1.5 and the [OIII] luminosity function from z=0.7-2.3 is almost entirely in the L* parameter, which steadily increases with redshift over those ranges.

preprint2013arXiv

The Fundamental Plane of Damped Lyman Alpha Systems

Using a sample of 100 H I - selected damped Lyman alpha (DLA) systems, observed with the High Resolution Echelle Spectrometer on the Keck I telescope, we present evidence that the scatter in the well-studied correlation between the redshift and metallicity of a DLA is largely due to the existence of a mass-metallicity relationship at each redshift. To describe the fundamental relations that exist between redshift, metallicity and mass, we use a fundamental plane description, which is described by the following equation: [M/H]=(-1.9+-0.5)+(0.74+-0.21)logdv_90-(0.32+-0.06)z. Here, we assert that the velocity width, dv_90, which is defined as the velocity interval containing 90% of the integrated optical depth, traces the mass of the underlying dark matter halo. This description provides two significant improvements over the individual descriptions of the mass-metallicity correlation and metallicity-redshift correlation. Firstly, the fundamental equation reduces the scatter around both relationships by about 20%, providing a more stringent constraint on numerical simulations modeling DLAs. Secondly, it confirms that the dark matter halos that host DLAs satisfy a mass-metallicity relationship at each redshift between redshifts 2 through 5.

preprint2013arXiv

Ultra-faint Ultraviolet Galaxies at z~2 Behind the Lensing Cluster Abell 1689: the Luminosity Function, Dust Extinction and Star Formation Rate Density

We identified the z~2 Lyman break galaxies using deep HST ultraviolet (F275W/F336W) imaging of Abell 1689. Because of the imaging depth and the large magnification provided by the cluster, we detect galaxies 100x fainter (-19.5< M_1500 <-13) than previous surveys at this redshift. We are able to calculate the intrinsic sensitivity of the observations as a function of source plane position, allowing determinations of effective volume as a function of luminosity. We fit the faint-end slope of the luminosity function to be alpha = -1.74 +/-0.08, consistent with the values obtained for 2.5 < z < 6. There is no turnover in the luminosity function down to MUV = -13. The trend of increasingly redder UV spectral slopes with luminosity at higher redshifts is observed in our sample, but with redder slopes at all luminosities and average reddening of < E(B - V) >= 0.15. We assume the stars in these galaxies are metal poor (0.2Z_sun) compared to their brighter counterparts (Z_sun), resulting in bluer assumed intrinsic UV slopes and larger derived dust extinction. The total UV luminosity density at z ~ 2 is 4.31x10^26 erg/s/Hz/Mpc^3, more than 70% of which is emitted by galaxies in the luminosity range of our sample. We determine the star formation rate density at z ~ 2 (assuming constant dust extinction correction of 4.2 over all luminosities and a Kroupa IMF) of 0.148 M/yr/Mpc^3, significantly higher than previous determinations because of the additional population of fainter galaxies and the larger dust correction factors.[abridged]

preprint2013arXiv

UVUDF: Ultraviolet Imaging of the Hubble Ultradeep Field with Wide-field Camera 3

We present an overview of a 90-orbit Hubble Space Telescope treasury program to obtain near ultraviolet imaging of the Hubble Ultra Deep Field using the Wide Field Camera 3 UVIS detector with the F225W, F275W, and F336W filters. This survey is designed to: (i) Investigate the episode of peak star formation activity in galaxies at 1<z<2.5; (ii) Probe the evolution of massive galaxies by resolving sub-galactic units (clumps); (iii) Examine the escape fraction of ionizing radiation from galaxies at z~2-3; (iv) Greatly improve the reliability of photometric redshift estimates; and (v) Measure the star formation rate efficiency of neutral atomic-dominated hydrogen gas at z~1-3. In this overview paper, we describe the survey details and data reduction challenges, including both the necessity of specialized calibrations and the effects of charge transfer inefficiency. We provide a stark demonstration of the effects of charge transfer inefficiency on resultant data products, which when uncorrected, result in uncertain photometry, elongation of morphology in the readout direction, and loss of faint sources far from the readout. We agree with the STScI recommendation that future UVIS observations that require very sensitive measurements use the instrument's capability to add background light through a "post-flash". Preliminary results on number counts of UV-selected galaxies and morphology of galaxies at z~1 are presented. We find that the number density of UV dropouts at redshifts 1.7, 2.1, and 2.7 is largely consistent with the number predicted by published luminosity functions. We also confirm that the image mosaics have sufficient sensitivity and resolution to support the analysis of the evolution of star-forming clumps, reaching 28-29th magnitude depth at 5 sigma in a 0.2 arcsecond radius aperture depending on filter and observing epoch.

preprint2012arXiv

Discovery of Three Distant, Cold Brown Dwarfs in the WFC3 Infrared Spectroscopic Parallels Survey

We present the discovery of three late type (>T4) brown dwarfs, including a probable Y dwarf, in the WFC3 Infrared Spectroscopic Parallels (WISP) Survey. We use the G141 grism spectra to determine the spectral types of the dwarfs and derive distance estimates based on a comparison with nearby T dwarfs with known parallaxes. These are the most distant spectroscopically confirmed T/Y dwarfs, with the farthest at an estimated distance of ~400 pc. We compare the number of cold dwarfs found in the WISP survey with simulations of the brown dwarf mass function. The number found is generally consistent with an initial stellar mass function dN/dM \propto M^{-α} with α= 0.0--0.5, although the identification of a Y dwarf is somewhat surprising and may be indicative of either a flatter absolute magnitude/spectral type relation than previously reported or an upturn in the number of very late type brown dwarfs in the observed volume.

preprint2012arXiv

Metallicity Evolution of Damped Lyman-alpha Systems out to z~5

We present chemical abundance measurements for 47 damped Lyman-alpha systems (DLAs), 30 at z>4, observed with the Echellette Spectrograph and Imager and the High Resolution Echelle Spectrometer on the Keck telescopes. HI column densities of the DLAs are measured with Voigt profile fits to the Lyman-alpha profiles, and we find an increased number of false DLA identifications with SDSS at z>4 due to the increased density of the Lyman-alpha forest. Ionic column densities are determined using the apparent optical depth method, and we combine our new metallicity measurements with 195 from previous surveys to determine the evolution of the cosmic metallicity of neutral gas. We find the metallicity of DLAs decreases with increasing redshift, improving the significance of the trend and extending it to higher redshifts, with a linear fit of -0.22+-0.03 dex per unit redshift from z=0.09-5.06. The metallicity 'floor' of ~1/600 solar continues out to z~5, despite our sensitivity for finding DLAs with much lower metallicities. However, this floor is not statistically different from a steep tail to the distribution. We also find that the intrinsic scatter of metallicity among DLAs of ~0.5 dex continues out to z~5. In addition, the metallicity distribution and the alpha/Fe ratios of z>2 DLAs are consistent with being drawn from the same parent population with those of halo stars. It is therefore possible that the halo stars in the Milky Way formed out of gas that commonly exhibits DLA absorption at z>2.

preprint2011arXiv

Star Formation from DLA Gas in the Outskirts of Lyman Break Galaxies at z~3

We present evidence for spatially extended low surface brightness emission around Lyman break galaxies (LBGs) in the V band image of the Hubble Ultra Deep Field, corresponding to the z~3 rest-frame FUV light, which is a sensitive measure of star formation rates (SFRs). We find that the covering fraction of molecular gas at z~3 is not adequate to explain the emission in the outskirts of LBGs, while the covering fraction of neutral atomic-dominated hydrogen gas at high redshift is sufficient. We develop a theoretical framework to connect this emission around LBGs to the expected emission from neutral H I gas i.e., damped Lyman alpha systems (DLAs), using the Kennicutt-Schmidt (KS) relation. Working under the hypothesis that the observed FUV emission in the outskirts of LBGs is from in situ star formation in atomic-dominated hydrogen gas, the results suggest that the SFR efficiency in such gas at z~3 is between factors of 10 and 50 lower than predictions based on the local KS relation. The total SFR density in atomic-dominated gas at z~3 is constrained to be ~10% of that observed from the inner regions of LBGs. In addition, the metals produced by in situ star formation in the outskirts of LBGs yield metallicities comparable to those of DLAs, which is a possible solution to the "Missing Metals" problem for DLAs. Finally, the atomic-dominated gas in the outskirts of galaxies at both high and low redshifts has similar reduced SFR efficiencies and is consistent with the same power law.

preprint2010arXiv

Teaching Astronomy with an Inquiry Activity on Stellar Populations

We describe a new inquiry design aimed at teaching advanced high-school to senior college students the basics of stellar populations. The inquiry is designed to have students come up with their own version of the Hertzsprung-Russell diagram as a tool to understand how stars evolve based on their color, mass, and luminosity. The inquiry makes use of pictures and spectra of stars, which the students analyze and interpret to answer the questions they come up with at the beginning. The students undergo a similar experience to real astronomers, using the same tools and methods to figure out the phenomena they are trying to understand. Specifically, they use images and spectra of stars, and organize the data via tables and plots to find trends that will then enable them to answer their questions. The inquiry also includes a "thinking tool" to help connect the trends students observe to the larger picture of stellar evolution. We include a description of the goals of the inquiry, the activity description, the motivations and thoughts that went into the design of the inquiry, and reflections on how the inquiry activity worked in practice.

preprint2009arXiv

Deep Keck u-band imaging of the Hubble Ultra Deep Field: A catalog of z~3 Lyman Break Galaxies

We present a sample of 407 z~3 Lyman break galaxies (LBGs) to a limiting isophotal u-band magnitude of 27.6 mag in the Hubble Ultra Deep Field (UDF). The LBGs are selected using a combination of photometric redshifts and the u-band drop-out technique enabled by the introduction of an extremely deep u-band image obtained with the Keck I telescope and the blue channel of the LRIS spectrometer. The Keck u-band image, totaling 9 hours of integration time, has a one sigma depth of 30.7 mag arcsec^-2, making it one of the most sensitive u-band images ever obtained. The u-band image also substantially improves the accuracy of photometric redshift measurements of ~50% of the z~3 Lyman break galaxies, significantly reducing the traditional degeneracy of colors between z~3 and z~0.2 galaxies. This sample provides the most sensitive, high-resolution multi-filter imaging of reliably identified z~3 LBGs for morphological studies of galaxy formation and evolution and the star formation efficiency of gas at high redshift.

preprint2006arXiv

The Chandra View of the Supernova Remnant 0506-68.0 in the Large Magellanic Cloud

A new Chandra observation of SNR 0506-68.0 (also called N23) reveals a complex, highly structured morphology in the low energy X-ray band and an isolated compact central object in the high energy band. Spectral analysis indicates that the X-ray emission overall is dominated by thermal gas whose composition is consistent with swept-up ambient material. There is a strong gradient in ambient density across the diameter of the remnant. Toward the southeast, near a prominent star cluster, the emitting density is 10 - 23 cm^{-3} while toward the northwest it has dropped to a value of only 1 cm^{-3}. The total extent of the X-ray remnant is 100" by 120" (24 pc x 29 pc for a distance of 50 kpc), somewhat larger than previously known. The remnant's age is estimated to be ~4600 yr. One part of the remnant shows evidence for enhanced O, Ne, and perhaps Mg abundances, which is interpreted as evidence for ejecta from a massive star core collapse supernova. The compact central object has a luminosity of a few times 10^{33} ergs/s and no obvious radio or optical counterpart. It does not show an extended nebula or pulsed emission as expected from a young energetic pulsar, but resembles the compact central objects seen in other core collapse SNe, such as Cas A.