Source author record

Alex Hagen

Alex Hagen 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

23works
10topics
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

23 published item(s)

preprint2022arXiv

DBCal: Density Based Calibration of classifier predictions for uncertainty quantification

Measurement of uncertainty of predictions from machine learning methods is important across scientific domains and applications. We present, to our knowledge, the first such technique that quantifies the uncertainty of predictions from a classifier and accounts for both the classifier's belief and performance. We prove that our method provides an accurate estimate of the probability that the outputs of two neural networks are correct by showing an expected calibration error of less than 0.2% on a binary classifier, and less than 3% on a semantic segmentation network with extreme class imbalance. We empirically show that the uncertainty returned by our method is an accurate measurement of the probability that the classifier's prediction is correct and, therefore has broad utility in uncertainty propagation.

preprint2022arXiv

Reduction of detection limit and quantification uncertainty due to interferent by neural classification with abstention

Many measurements in the physical sciences can be cast as counting experiments, where the number of occurrences of a physical phenomenon informs the prevalence of the phenomenon's source. Often, detection of the physical phenomenon (termed signal) is difficult to distinguish from naturally occurring phenomena (termed background). In this case, the discrimination of signal events from background can be performed using classifiers, and they may range from simple, threshold-based classifiers to sophisticated neural networks. These classifiers are often trained and validated to obtain optimal accuracy, however we show that the optimal accuracy classifier does not generally coincide with a classifier that provides the lowest detection limit, nor the lowest quantification uncertainty. We present a derivation of the detection limit and quantification uncertainty in the classifier-based counting experiment case. We also present a novel abstention mechanism to minimize the detection limit or quantification uncertainty \emph{a posteriori}. We illustrate the method on two data sets from the physical sciences, discriminating Ar-37 and Ar-39 radioactive decay from non-radioactive events in a gas proportional counter, and discriminating neutrons from photons in an inorganic scintillator and report results therefrom.

preprint2020arXiv

Scaling the training of particle classification on simulated MicroBooNE events to multiple GPUs

Measurements in Liquid Argon Time Projection Chamber (LArTPC) neutrino detectors, such as the MicroBooNE detector at Fermilab, feature large, high fidelity event images. Deep learning techniques have been extremely successful in classification tasks of photographs, but their application to LArTPC event images is challenging, due to the large size of the events. Events in these detectors are typically two orders of magnitude larger than images found in classical challenges, like recognition of handwritten digits contained in the MNIST database or object recognition in the ImageNet database. Ideally, training would occur on many instances of the entire event data, instead of many instances of cropped regions of interest from the event data. However, such efforts lead to extremely long training cycles, which slow down the exploration of new network architectures and hyperparameter scans to improve the classification performance. We present studies of scaling a LArTPC classification problem on multiple architectures, spanning multiple nodes. The studies are carried out on simulated events in the MicroBooNE detector. We emphasize that it is beyond the scope of this study to optimize networks or extract the physics from any results here. Institutional computing at Pacific Northwest National Laboratory and the SummitDev machine at Oak Ridge National Laboratory's Leadership Computing Facility have been used. To our knowledge, this is the first use of state-of-the-art Convolutional Neural Networks for particle physics and their attendant compute techniques onto the DOE Leadership Class Facilities. We expect benefits to accrue particularly to the Deep Underground Neutrino Experiment (DUNE) LArTPC program, the flagship US High Energy Physics (HEP) program for the coming decades.

preprint2016arXiv

On the Classification of UGC1382 as a Giant Low Surface Brightness Galaxy

We provide evidence that UGC1382, long believed to be a passive elliptical galaxy, is actually a giant low surface brightness (GLSB) galaxy which rivals the archetypical GLSB Malin 1 in size. Like other GLSB galaxies, it has two components: a high surface brightness disk galaxy surrounded by an extended low surface brightness (LSB) disk. For UGC1382, the central component is a lenticular system with an effective radius of 6 kpc. Beyond this, the LSB disk has an effective radius of ~38 kpc and an extrapolated central surface brightness of ~26 mag/arcsec^2. Both components have a combined stellar mass of ~8x10^10 M_sun, and are embedded in a massive (10^10 M_sun) low-density (<3 M_sun/pc^2) HI disk with a radius of 110 kpc, making this one of the largest isolated disk galaxies known. The system resides in a massive dark matter halo of at least 2x10^12 M_sun. Although possibly part of a small group, its low density environment likely plays a role in the formation and retention of the giant LSB and HI disks. We model the spectral energy distributions and find that the LSB disk is likely older than the lenticular component. UGC1382 has UV-optical colors typical of galaxies transitioning through the green valley. Within the LSB disk are spiral arms forming stars at extremely low efficiencies. The gas depletion time scale of ~10^11 yr suggests that UGC1382 may be a very long term resident of the green valley. We find that the formation and evolution of the LSB disk is best explained by the accretion of gas-rich LSB dwarf galaxies.

preprint2016arXiv

The Astropy Problem

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

preprint2015arXiv

HST Emission Line Galaxies at z ~ 2: Comparing Physical Properties of Lyman Alpha and Optical Emission Line Selected Galaxies

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

preprint2015arXiv

Learning from FITS: Limitations in use in modern astronomical research

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard, however, is showing its age. Developed in the late 1970s, the FITS authors made a number of implementation choices that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not anticipate the challenges which we are facing today in astronomical computing. Difficulties we now face include, but are not limited to, addressing the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets, and capturing significantly more complex metadata and data relationships. There are members of the community today who find some or all of these limitations unworkable, and have decided to move ahead with storing data in other formats. If this fragmentation continues, we risk abandoning the advantages of broad interoperability, and ready archivability, that the FITS format provides for astronomy. In this paper we detail some selected important problems which exist within the FITS standard today. These problems may provide insight into deeper underlying issues which reside in the format and we provide a discussion of some lessons learned. It is not our intention here to prescribe specific remedies to these issues; rather, it is to call attention of the FITS and greater astronomical computing communities to these problems in the hope that it will spur action to address them.

preprint2015arXiv

Surveying Galaxy Proto-clusters in Emission: A Large-scale Structure at z=2.44 and the Outlook for HETDEX

Galaxy proto-clusters at z >~ 2 provide a direct probe of the rapid mass assembly and galaxy growth of present day massive clusters. Because of the need of precise galaxy redshifts for density mapping and the prevalence of star formation before quenching, nearly all the proto-clusters known to date were confirmed by spectroscopy of galaxies with strong emission lines. Therefore, large emission-line galaxy surveys provide an efficient way to identify proto-clusters directly. Here we report the discovery of a large-scale structure at z = 2.44 in the HETDEX Pilot Survey. On a scale of a few tens of Mpc comoving, this structure shows a complex overdensity of Lya emitters (LAE), which coincides with broad-band selected galaxies in the COSMOS/UltraVISTA photometric and zCOSMOS spectroscopic catalogs, as well as overdensities of intergalactic gas revealed in the Lya absorption maps of Lee et al. (2014). We construct mock LAE catalogs to predict the cosmic evolution of this structure. We find that such an overdensity should have already broken away from the Hubble flow, and part of the structure will collapse to form a galaxy cluster with 10^14.5 +- 0.4 M_sun by z = 0. The structure contains a higher median stellar mass of broad-band selected galaxies, a boost of extended Lya nebulae, and a marginal excess of active galactic nuclei relative to the field, supporting a scenario of accelerated galaxy evolution in cluster progenitors. Based on the correlation between galaxy overdensity and the z = 0 descendant halo mass calibrated in the simulation, we predict that several hundred 1.9 < z < 3.5 proto-clusters with z = 0 mass of > 10^14.5 M_sun will be discovered in the 8.5 Gpc^3 of space surveyed by the Hobby Eberly Telescope Dark Energy Experiment.

preprint2015arXiv

The Biases of Optical Line-Ratio Selection for Active Galactic Nuclei, and the Intrinsic Relationship between Black Hole Accretion and Galaxy Star Formation

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

preprint2015arXiv

The Dust Attenuation Curve versus Stellar Mass for Emission Line Galaxies at z ~ 2

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

preprint2015arXiv

The Evolution of the Far-UV Luminosity Function and Star Formation Rate Density of the Chandra Deep Field South from z=0.2-1.2 with Swift/UVOT

We use deep Swift UV/Optical Telescope (UVOT) near-ultraviolet (1600A to 4000A) imaging of the Chandra Deep Field South to measure the rest-frame far-UV (FUV; 1500A) luminosity function (LF) in four redshift bins between z=0.2 and 1.2. Our sample includes 730 galaxies with u < 24.1 mag. We use two methods to construct and fit the LFs: the traditional V_max method with bootstrap errors and a maximum likelihood estimator. We observe luminosity evolution such that M* fades by ~2 magnitudes from z~1 to z~0.3 implying that star formation activity was substantially higher at z~1 than today. We integrate our LFs to determine the FUV luminosity densities and star formation rate densities from z=0.2 to 1.2. We find evolution consistent with an increase proportional to (1+z)^1.9 out to z~1. Our luminosity densities and star formation rates are consistent with those found in the literature, but are, on average, a factor of ~2 higher than previous FUV measurements. In addition, we combine our UVOT data with the MUSYC survey to model the galaxies' ultraviolet-to-infrared spectral energy distributions and estimate the rest-frame FUV attenuation. We find that accounting for the attenuation increases the star formation rate densities by ~1 dex across all four redshift bins.

preprint2015arXiv

Young, Star-forming Galaxies and their local Counterparts: the Evolving Relationship of Mass-SFR-Metallicity since z ~ 2.1

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

preprint2014arXiv

3D-HST Emission Line Galaxies at z ~ 2: Discrepancies in the Optical/UV Star Formation Rates

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

preprint2014arXiv

HST Emission Line Galaxies at z ~ 2: The Ly-alpha Escape Fraction

We compare the H-beta line strengths of 1.90 < z < 2.35 star-forming galaxies observed with the near-IR grism of the Hubble Space Telescope with ground-based measurements of Ly-alpha from the HETDEX Pilot Survey and narrow-band imaging. By examining the line ratios of 73 galaxies, we show that most star-forming systems at this epoch have a Ly-alpha escape fraction below ~6%. We confirm this result by using stellar reddening to estimate the effective logarithmic extinction of the H-beta emission line (c_Hbeta = 0.5) and measuring both the H-beta and Ly-alpha luminosity functions in a ~ 100,000 cubic Mpc volume of space. We show that in our redshift window, the volumetric Ly-alpha escape fraction is at most 4.4+/-2.1(1.2)%, with an additional systematic ~25% uncertainty associated with our estimate of extinction. Finally, we demonstrate that the bulk of the epoch's star-forming galaxies have Ly-alpha emission line optical depths that are significantly greater than that for the underlying UV continuum. In our predominantly [O~III] 5007-selected sample of galaxies, resonant scattering must be important for the escape of Ly-alpha photons.

preprint2014arXiv

HST Emission Line Galaxies at z ~ 2: The Mystery of Neon

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

preprint2014arXiv

Physical and Morphological Properties of [O II] Emitting Galaxies in the HETDEX Pilot Survey

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

preprint2014arXiv

Spectral Energy Distribution Fitting of HETDEX Pilot Survey Lyman-alpha Emitters in COSMOS and GOODS-N

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

preprint2014arXiv

Studying Gender in Conference Talks -- data from the 223rd meeting of the American Astronomical Society

We present a study on the gender balance, in speakers and attendees, at the recent major astronomical conference, the American Astronomical Society meeting 223, in Washington, DC. We conducted an informal survey, yielding over 300 responses by volunteers at the meeting. Each response included gender data about a single talk given at the meeting, recording the gender of the speaker and all question-askers. In total, 225 individual AAS talks were sampled. We analyze basic statistical properties of this sample. We find that the gender ratio of the speakers closely matched the gender ratio of the conference attendees. The audience asked an average of 2.8 questions per talk. Talks given by women had a slightly higher number of questions asked (3.2$\pm$0.2) than talks given by men (2.6$\pm$0.1). The most significant result from this study is that while the gender ratio of speakers very closely mirrors that of conference attendees, women are under-represented in the question-asker category. We interpret this to be an age-effect, as senior scientists may be more likely to ask questions, and are more commonly men. A strong dependence on the gender of session chairs is found, whereby women ask disproportionately fewer questions in sessions chaired by men. While our results point to laudable progress in gender-balanced speaker selection, we believe future surveys of this kind would help ensure that collaboration at such meetings is as inclusive as possible.

preprint2014arXiv

The HETDEX Pilot Survey V: The Physical Origin of Lyman-alpha Emitters Probed by Near-infrared Spectroscopy

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

preprint2013arXiv

Galaxy redshift surveys with sparse sampling

Survey observations of the three-dimensional locations of galaxies are a powerful approach to measure the distribution of matter in the universe, which can be used to learn about the nature of dark energy, physics of inflation, neutrino masses, etc. A competitive survey, however, requires a large volume (e.g., Vsurvey is roughly 10 Gpc3) to be covered, and thus tends to be expensive. A "sparse sampling" method offers a more affordable solution to this problem: within a survey footprint covering a given survey volume, Vsurvey, we observe only a fraction of the volume. The distribution of observed regions should be chosen such that their separation is smaller than the length scale corresponding to the wavenumber of interest. Then one can recover the power spectrum of galaxies with precision expected for a survey covering a volume of Vsurvey (rather than the volume of the sum of observed regions) with the number density of galaxies given by the total number of observed galaxies divided by Vsurvey (rather than the number density of galaxies within an observed region). We find that regularly-spaced sampling yields an unbiased power spectrum with no window function effect, and deviations from regularly-spaced sampling, which are unavoidable in realistic surveys, introduce calculable window function effects and increase the uncertainties of the recovered power spectrum. While we discuss the sparse sampling method within the context of the forthcoming Hobby-Eberly Telescope Dark Energy Experiment, the method is general and can be applied to other galaxy surveys.

preprint2013arXiv

Searching for Neutral Hydrogen Halos around z ~ 2.1 and z ~ 3.1 Ly-alpha Emitting Galaxies

We search for evidence of diffuse Ly-alpha emission from extended neutral hydrogen surrounding Ly-alpha emitting galaxies (LAEs) using deep narrow-band images of the Extended Chandra Deep Field South. By stacking the profiles of 187 LAEs at z = 2.06, 241 LAEs at z = 3.10, and 179 LAEs at z = 3.12, and carefully performing low-surface brightness photometry, we obtain mean surface brightness maps that reach 9.9, 8.7, and 6.2 * 10^{-19} ergs cm^{-2} s^{-1} arcsec^{-2} in the emission line. We undertake a thorough investigation of systematic uncertainties in our surface brightness measurements, and find that our limits are 5--10 times larger than would be expected from Poisson background fluctuations; these uncertainties are often underestimated in the literature. At z ~ 3.1, we find evidence for extended halos with small scale lengths of 5--8 kpc in some, but not all of our sub-samples. We demonstrate that sub-samples of LAEs with low equivalent widths and brighter continuum magnitudes are more likely to possess such halos. At z ~ 2.1, we find no evidence of extended Ly-alpha emission down to our detection limits. Through Monte-Carlo simulations, we also show that we would have detected large diffuse LAE halos if they were present in our data sets. We compare these findings to other measurements in the literature, and discuss possible instrumental and astrophysical reasons for the discrepancies.

preprint2013arXiv

The Spectrally Resolved Lyman-alpha Emission of Three Lyman-alpha Selected Field Galaxies at z~2.4 from the HETDEX Pilot Survey

We present new results on the spectrally resolved Lyman-alpha (LyA) emission of three LyA emitting field galaxies at z~2.4 with high LyA equivalent width (>100 Angstroms) and LyA luminosity (~10^43 erg/s). At 120 km/s (FWHM) spectral resolution, the prominent double-peaked LyA profile straddles the systemic velocity, where the velocity zero-point is determined from spectroscopy of the galaxies' rest-frame optical nebular emission lines. The average velocity offset from systemic of the stronger redshifted emission component for our sample is 176 km/s while the average total separation between the redshifted and main blueshifted emission components is 380 km/s. These measurements are a factor of ~2 smaller than for UV continuum-selected galaxies that show LyA in emission with lower LyA equivalent width. We compare our LyA spectra to the predicted line profiles of a spherical "expanding shell" LyA radiative transfer grid that models large-scale galaxy outflows. Specifically blueward of the systemic velocity where two galaxies show a weak, highly blueshifted (by ~1000 km/s) tertiary emission peak, the model line profiles are a relatively poor representation of the observed spectra. Since the neutral gas column density has a dominant influence over the shape of the LyA line profile, we caution against equating the observed LyA velocity offset with a physical outflow velocity, especially at lower spectral resolution where the unresolved LyA velocity offset is a convoluted function of several degenerate parameters. Referring to rest-frame ultraviolet and optical Hubble Space Telescope imaging, we find that galaxy-galaxy interactions may play an important role in inducing a starburst that results in copious LyA emission, as well as perturbing the gas distribution and velocity field which have strong influence over the LyA emission line profile.

preprint2013arXiv

To Stack or Not to Stack: Spectral Energy Distribution Properties of Lya-Emitting Galaxies at z=2.1

We use the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) GOODS-S multi-wavelength catalog to identify counterparts for 20 Lya Emitting (LAE) galaxies at z=2.1. We build several types of stacked Spectral Energy Distributions (SEDs) of these objects. We combine photometry to form average and median flux-stacked SEDs, and postage stamp images to form average and median image-stacked SEDs. We also introduce scaled flux stacks that eliminate the influence of variation in overall brightness. We use the SED fitting code SpeedyMC to constrain the physical properties of individual objects and stacks. Our LAEs at z = 2.1 have stellar masses ranging from 2x10^7 Msun - 8x10^9 Msun (median = 3x10^8 Msun), ages ranging from 4 Myr to 500 Myr (median =100 Myr), and E(B-V) between 0.02 and 0.24 (median = 0.12). We do not observe strong correlations between Lya equivalent width (EW) and stellar mass, age, or E(B-V). The Lya radiative transfer (q) factors of our sample are predominantly close to one and do not correlate strongly with EW or E(B-V), implying that Lya radiative transfer prevents Lya photons from resonantly scattering in dusty regions. The SED parameters of the flux stacks match the average and median values of the individual objects, with the flux-scaled median SED performing best with reduced uncertainties. Median image-stacked SEDs provide a poor representation of the median individual object, and none of the stacking methods captures the large dispersion of LAE properties.