Source author record

Janice C. Lee

Janice C. Lee 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

46works
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

46 published item(s)

preprint2025arXiv

Evolved Supergiants in PHANGS I: Red Supergiants in 19 Galaxies between 5-20 Mpc with HST and JWST

Red supergiants (RSGs) are important for our understanding of supernova progenitors, stellar populations, stellar evolution, mass loss and dust production. Extragalactic surveys of RSGs have a long history in the Local Group, but few studies exist beyond that due to the limited resolution and sensitivity of ground-based and previous space-based infrared observatories. Here we demonstrate the combined power of HST and JWST to push systematic searches of RSGs out to $\sim$20 Mpc. We introduce a catalog of 97057 RSGs -- the largest single-survey release of RSGs -- with masses $\gtrsim$10 M$_{\odot}$ in 19 galaxies from the PHANGS HST+JWST Treasury program. We use HST F814W and JWST F200W photometry to select stars as RSGs based on predicted colors and magnitudes from PARSEC isochrones. The spatial distribution of our recovered RSGs follow the familiar pattern of mostly being concentrated in active star-forming regions such as spiral arms and central starburst rings. The RSG number density on kpc-scales is strongly correlated ($r_s$$\sim$0.82) with local star-formation rate density ($Σ_{SFR}$) traced by extinction-corrected far-ultraviolet (FUV) from GALEX+WISE, and weakly correlated ($r_s$$\sim$0.57) with the total stellar mass density ($Σ_*$), traced by near-infrared emission from WISE+Spitzer. The number of RSGs per mass of stellar populations with ages 6-30 Myr (the likely age range of RSGs $>$10 M$_{\odot}$) is $\sim$1 per 10$^{3.77\pm0.27}$ M$_{\odot}$, assuming constant star-formation rates from FUV+W4. Our sample will be a useful resource for tracking progenitors and feedback sites of future supernovae in PHANGS, age-dating stellar populations, and more.

preprint2024arXiv

ALMA-LEGUS I: The Influence of Galaxy Morphology on Molecular Cloud Properties

We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of the two galaxies with the same sensitivities and resolutions of 13 pc, we directly compare the molecular gas in these two similar galaxies to determine the physical conditions responsible for their large disparity in cluster formation. By fitting size-linewidth relations for the clouds in each galaxy, we find that NGC 1313 has a higher intercept than NGC 7793, implying that its clouds have higher kinetic energies at a given size scale. NGC 1313 also has more clouds near virial equilibrium than NGC 7793, which may be connected to its higher rate of massive cluster formation. However, these virially bound clouds do not show a stronger correlation with young clusters than that of the general cloud population. We find surprisingly small differences between the distributions of molecular cloud populations in the two galaxies, though the largest of those differences are that NGC 1313 has higher surface densities and lower free-fall times.

preprint2024arXiv

ALMA-LEGUS II: The Influence of Sub-Galactic Environment on Molecular Cloud Properties

We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy than there are between the galaxies on a global scale, especially for NGC 1313. There are higher masses, linewidths, pressures, and virial parameters in the arms of NGC 1313 and center of NGC 7793 than in the interarm and outer regions of the galaxies. The massive cluster formation of NGC 1313 may be driven by its greater variation in environments, allowing more clouds with the necessary conditions to arise, although no one parameter seems primarily responsible for the difference in star formation. Meanwhile NGC 7793 has clouds that are as massive and have as much kinetic energy as clouds in the arms of NGC 1313, but have densities and pressures more similar to the interarm regions and so are less inclined to collapse and form stars. The cloud properties in NGC 1313 and NGC 7793 suggest that spiral arms, bars, interarm regions, and flocculent spirals each represent distinct environments with regard to molecular cloud populations. We see surprisingly little difference in surface densities between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions of lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density.

preprint2022arXiv

H$α$ Morphologies of Star Clusters in 16 LEGUS Galaxies: Constraints on HII region evolution timescales

The analysis of star cluster ages in tandem with the morphology of their HII regions can provide insight into the processes that clear a cluster's natal gas, as well as the accuracy of cluster ages and dust reddening derived from Spectral Energy Distribution (SED) fitting. We classify 3757 star clusters in 16 nearby galaxies according to their H$α$ morphology (concentrated, partially exposed, no emission), using Hubble Space Telescope (HST) imaging from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We find: 1) The mean SED ages of clusters with concentrated (1-2 Myr) and partially exposed HII region morphologies (2-3 Myr) indicate a relatively early onset of gas clearing and a short (1-2 Myr) clearing timescale. 2) The reddening of clusters can be overestimated due to the presence of red supergiants, which is a result of stochastic sampling of the IMF in low mass clusters. 3) The age-reddening degeneracy impacts the results of the SED fitting - out of 1408 clusters with M$_*$ $\geq$ 5000 M$_{\odot}$, we find that at least 46 (3%) have SED ages which appear significantly underestimated or overestimated based on H$α$ and their environment, while the total percentage of poor age estimates is expected to be several times larger. 4) Lastly, we examine the dependence of the morphological classifications on spatial resolution. At HST resolution, our conclusions are robust to the distance range spanned by the sample (3-10 Mpc). However, analysis of ground-based H$α$ images shows that compact and partially exposed morphologies frequently cannot be distinguished from each other.

preprint2021arXiv

PHANGS-HST: Star Cluster Spectral Energy Distribution Fitting with CIGALE

The sensitivity and angular resolution of photometric surveys executed by the Hubble Space Telescope (HST) enable studies of individual star clusters in galaxies out to a few tens of megaparsecs. The fitting of spectral energy distributions (SEDs) of star clusters is essential for measuring their physical properties and studying their evolution. We report on the use of the publicly available Code Investigating GALaxy Emission (CIGALE) SED fitting package to derive ages, stellar masses, and reddenings for star clusters identified in the Physics at High Angular resolution in Nearby GalaxieS-HST (PHANGS-HST) survey. Using samples of star clusters in the galaxy NGC 3351, we present results of benchmark analyses performed to validate the code and a comparison to SED fitting results from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We consider procedures for the PHANGS-HST SED fitting pipeline, e.g., the choice of single stellar population models, the treatment of nebular emission and dust, and the use of fluxes versus magnitudes for the SED fitting. We report on the properties of clusters in NGC 3351 and find, on average, the clusters residing in the inner star-forming ring of NGC 3351 are young ($< 10$ Myr) and massive ($10^{5} M_{\odot}$) while clusters in the stellar bulge are significantly older. Cluster mass function fits yield $β$ values around -2, consistent with prior results with a tendency to be shallower at the youngest ages. Finally, we explore a Bayesian analysis with additional physically-motivated priors for the distribution of ages and masses and analyze the resulting cluster distributions.

preprint2021arXiv

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $μ$m imaging has been approved.

preprint2020arXiv

A Comparison of Star-Forming Clumps and Tidal Tails in Local Mergers and High Redshift Galaxies

The Clusters, Clumps, Dust, and Gas in Extreme Star-Forming Galaxies (CCDG) survey with the Hubble Space Telescope includes multi-wavelength imaging of 13 galaxies less than 100 Mpc away spanning a range of morphologies and sizes, from Blue Compact Dwarfs (BCDs) to luminous infrared galaxies (LIRGs), all with star formation rates in excess of hundreds of solar masses per year. Images of 7 merging galaxies in the CCDG survey were artificially redshifted to compare with galaxies at z=0.5, 1, and 2. Most redshifted tails have surface brightnesses that would be visible at z=0.5 or 1 but not at z=2 due to cosmological dimming. Giant star clumps are apparent in these galaxies; the 51 measured have similar sizes, masses and colors as clumps in observed high-z systems in UDF, GEMS, GOODS, and CANDELS surveys. These results suggest that some clumpy high-z galaxies without observable tidal features could be the result of mergers. The local clumps also have the same star formation rate per unit area and stellar surface density as clumps observed at intermediate and high redshift, so they provide insight into the substructure of distant clumps. A total of 1596 star clusters brighter than MV = -9 were identified within the boundaries of the local clumps. The cluster magnitude distribution function is a power law with approximately the same slope (approximately -1 for a number-log luminosity plot) for all the galaxies both inside and outside the clumps and independent of clump surface brightness.

preprint2020arXiv

Candidate LBV stars in galaxy NGC 7793 found via HST photometry + MUSE spectroscopy

Only about 19 Galactic and 25 extra-galactic bona-fide Luminous Blue Variables (LBVs) are known to date. This incomplete census prevents our understanding of this crucial phase of massive star evolution which leads to the formation of heavy binary black holes via the classical channel. With large samples of LBVs one could better determine the duration and maximum stellar luminosity which characterize this phase. We search for candidate LBVs (cLBVs) in a new galaxy, NGC 7793. For this purpose, we combine high spatial resolution images from two Hubble Space Telescope (HST) programs with optical spectroscopy from the Multi Unit Spectroscopic Explorer (MUSE). By combining PSF-fitting photometry measured on F547M, F657N, and F814W images, with restrictions on point-like appearance (at HST resolution) and Hα luminosity, we find 100 potential cLBVs, 36 of which fall in the MUSE fields. Five of the latter 36 sources are promising cLBVs which have M(V) less than or equal to -7 and a combination of: Hα with a P-Cygni profile; no [O I] 6300 emission; weak or no [O III] 5007 emission; large [N II]/Hα relative to H II regions; and [S II] 6716 / [S II] 6731 similar to 1. It is not clear if these five cLBVs are isolated from O-type stars, which would favor the binary formation scenario of LBVs. Our study, which approximately covers one fourth of the optical disc of NGC 7793, demonstrates how by combining the above HST surveys with multi-object spectroscopy from 8-m class telescopes, one can efficiently find large samples of cLBVs in nearby galaxies.

preprint2020arXiv

Deep Transfer Learning for Star Cluster Classification: I. Application to the PHANGS-HST Survey

We present the results of a proof-of-concept experiment which demonstrates that deep learning can successfully be used for production-scale classification of compact star clusters detected in HST UV-optical imaging of nearby spiral galaxies (D < 20 Mpc) in the PHANGS-HST survey. Given the relatively small nature of existing, human-labelled star cluster samples, we transfer the knowledge of state-of-the-art neural network models for real-object recognition to classify star clusters candidates into four morphological classes. We perform a series of experiments to determine the dependence of classification performance on: neural network architecture (ResNet18 and VGG19-BN); training data sets curated by either a single expert or three astronomers; and the size of the images used for training. We find that the overall classification accuracies are not significantly affected by these choices. The networks are used to classify star cluster candidates in the PHANGS-HST galaxy NGC 1559, which was not included in the training samples. The resulting prediction accuracies are 70%, 40%, 40-50%, 50-70% for class 1, 2, 3 star clusters, and class 4 non-clusters respectively. This performance is competitive with consistency achieved in previously published human and automated quantitative classification of star cluster candidate samples (70-80%, 40-50%, 40-50%, and 60-70%). The methods introduced herein lay the foundations to automate classification for star clusters at scale, and exhibit the need to prepare a standardized dataset of human-labelled star cluster classifications, agreed upon by a full range of experts in the field, to further improve the performance of the networks introduced in this study.

preprint2020arXiv

Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ~20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the dataset will provide an unprecedented look into the connections between young stars, HII regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables into physical parameters (e.g., brightness to luminosity, angular to physical sizes of molecular clouds, star clusters and associations). PHANGS-HST has obtained parallel ACS imaging of the galaxy halos in the F606W and F814W bands. Where possible, we use these parallel fields to derive tip of the red giant branch (TRGB) distances to these galaxies. In this paper, we present TRGB distances for 11 galaxies from ~4 to ~15 Mpc, based on the first year of PHANGS-HST observations. Five of these represent the first published TRGB distance measurements (IC 5332, NGC 2835, NGC 4298, NGC 4321, and NGC 4328), and eight of which are the best available distances to these targets. We also provide a compilation of distances for the 118 galaxies in the full PHANGS sample, which have been adopted for the first PHANGS-ALMA public data release.

preprint2020arXiv

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $Σ_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $Σ_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.

preprint2019arXiv

The lifecycle of molecular clouds in nearby star-forming disc galaxies

It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially-resolved ($\sim100$ pc) CO-to-H$α$ flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10-30 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $Σ_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $Σ_{\rm H_2}\leqslant8$M$_{\odot}$pc$^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H$α$ (75-90% of the cloud lifetime), GMCs disperse within just 1-5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.

preprint2016arXiv

A Deeper Look at Faint H$α$ Emission in Nearby Dwarf Galaxies

We present deep H$α$ imaging of three nearby dwarf galaxies, carefully selected to optimize observations with the Maryland-Magellan Tunable Filter (MMTF) on the Magellan 6.5m telescope. An effective bandpass of $\sim$13Å is used, and the images reach 3$σ$ flux limits of $\sim$8$\times10^{-18}$ ergs s$^{-1}$ cm$^{-2}$, which is about an order of magnitude lower than standard narrowband observations obtained by the most recent generation of local H$α$ galaxy surveys. The observations were originally motivated by the finding that the H$α$/FUV flux ratio of galaxies systematically declines as global galactic properties such as the star formation rate and stellar mass decrease. The three dwarf galaxies selected for study have star formation rates, that when calculated from their H$α$ luminosities using standard conversion recipes, are $\sim$50\% of those based on the FUV. Follow-up studies of many of the potential causes for the trends in the H$α$/FUV flux ratio have been performed, but the possibility that previous observations have missed a non-negligible fraction of faint ionized emission in dwarf galaxies has not been investigated. The MMTF observations reveal both diffuse and structured H$α$ emission (filaments, shells, possible single-star HII regions) spanning extents up to 2.5 times larger relative to previous observations. However, only up to an additional $\sim5$\% of H$α$ flux is captured, which does not account for the trends in the H$α$/FUV ratio. Beyond investigation of the H$α$/FUV ratio, the impact of the newly detected extended flux on our understanding of star formation, the properties of HII regions, and the propagation of ionizing photons warrant further investigation.

preprint2016arXiv

Direct Detection of Lyman Continuum Escape from Local Starburst Galaxies with the Cosmic Origins Spectrograph

We report on the detection of Lyman continuum radiation in two nearby starburst galaxies. Tol 0440-381, Tol 1247-232 and Mrk 54 were observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescopes. The three galaxies have radial velocities of ~13,000 km/s, permitting a ~35 A window on the restframe Lyman continuum shortward of the Milky Way Lyman edge at 912 A. The chosen instrument configuration using the G140L grating covers the spectral range from 912 to 2,000 Å. We developed a dedicated background subtraction method to account for temporal and spatial background variations of the detector, which is crucial at the low flux levels around 912 A. This modified pipeline allowed us to significantly improve the statistical and systematic detector noise and will be made available to the community. We detect Lyman continuum in all three galaxies. However, we conservatively interpret the emission in Tol 0440-381 as an upper limit due to possible contamination by geocoronal Lyman series lines. We determined the current star-formation properties from the far-ultraviolet continuum and spectral lines and used synthesis models to predict the Lyman continuum radiation emitted by the current population of hot stars. We discuss the various model uncertainties such as, among others, atmospheres and evolution models. Lyman continuum escape fractions were derived from a comparison between the observed and predicted Lyman continuum fluxes. Tol 1247-232, Mrk 54 and Tol 0440-381 have absolute escape fractions of (4.5 +/- 1.2)%, (2.5 +/- 0.72)% and <(7.1 +/- 1.1)%, respectively.

preprint2016arXiv

GALEX-SDSS-WISE Legacy Catalog (GSWLC): Star Formation Rates, Stellar Masses and Dust Attenuations of 700,000 Low-redshift Galaxies

In this paper, we present GALEX-SDSS-WISE Legacy Catalog (GSWLC), a catalog of physical properties (stellar masses, dust attenuations and star formation rates (SFRs)) of ~700,000 galaxies with SDSS redshifts below 0.3. GSWLC contains galaxies within the GALEX footprint, regardless of a UV detection, covering 90% of SDSS. The physical properties were obtained from UV/optical SED fitting following Bayesian methodology of Salim et al. (2007), with improvements such as blending corrections for low-resolution UV photometry, flexible dust attenuation laws, and emission line corrections. GSWLC includes mid-IR SFRs derived from IR templates based upon 22 micron WISE observations. These estimates are independent of UV/optical SED fitting, in order to separate possible systematics. The paper argues that the comparison of specific SFRs (SSFRs) is more informative and physically motivated than the comparison of SFRs. SSFRs resulting from the UV/optical SED fitting are compared to the mid-IR SSFRs, and to SSFRs from three published catalogs. For "main sequence" galaxies with no AGN contribution all SSFRs are in very good agreement (within 0.1 dex on average). In particular, the widely-used aperture-corrected SFRs from MPA/JHU catalog show no systematic offsets, in contrast to some integral-field spectroscopy results. For galaxies below the main sequence (log SSFR<-11), mid-IR (S)SFRs based on fixed luminosity-SFR conversion are severely biased (up to 2 dex) because the dust is primarily heated by old stars. Furthermore, mid-IR (S)SFRs are overestimated by up to 0.6 dex for galaxies with AGN, presumably due to non-stellar dust heating. UV/optical (S)SFRs are thus preferred to IR-based (S)SFRs for quenched galaxies and those which host AGN.

preprint2016arXiv

The Connection Between Galaxy Environment and the Luminosity Function Slopes of Star-Forming Regions

We present the first study of GALEX far ultra-violet (FUV) luminosity functions of individual star-forming regions within a sample of 258 nearby galaxies spanning a large range in total stellar mass and star formation properties. We identify ~65,000 star-forming regions (i.e., FUV sources), measure each galaxy's luminosity function, and characterize the relationships between the luminosity function slope (alpha) and several global galaxy properties. A final sample of 82 galaxies with reliable luminosity functions are used to define these relationships and represent the largest sample of galaxies with the largest range of galaxy properties used to study the connection between luminosity function properties and galaxy environment. We find that alpha correlates with global star formation properties, where galaxies with higher star formation rates and star formation rate densities (Sigma_SFR) tend to have flatter luminosity function slopes. In addition, we find that neither stochastic sampling of the luminosity function in galaxies with low-number statistics nor the effects of blending due to distance can fully account for these trends. We hypothesize that the flatter slopes in high Sigma_SFR galaxies is due to higher gas densities and higher star formation efficiencies which result in proportionally greater numbers of bright star-forming regions. Finally, we create a composite luminosity function composed of star-forming regions from many galaxies and find a break in the luminosity function at brighter luminosities. However, we find that this break is an artifact of varying detection limits for galaxies at different distances.

preprint2015arXiv

A Transition Mass in the Local Tully-Fisher Relation

We study the stellar mass Tully-Fisher relation (TFR, stellar mass versus rotation velocity) for a morphologically blind selection of emission line galaxies in the field at redshifts 0.1 $<$ z $<$ 0.375. Kinematics ($σ_g$, V$_{rot}$) are measured from emission lines in Keck/DEIMOS spectra and quantitative morphology is measured from V- and I-band Hubble images. We find a transition stellar mass in the TFR, $\log$ M$_*$ = 9.5 M$_{\odot}$. Above this mass, nearly all galaxies are rotation-dominated, on average more morphologically disk-like according to quantitative morphology, and lie on a relatively tight TFR. Below this mass, the TFR has significant scatter to low rotation velocity and galaxies can either be rotation-dominated disks on the TFR or asymmetric or compact galaxies which scatter off. We refer to this transition mass as the "mass of disk formation", M$_{\mathrm{df}}$ because above it all star-forming galaxies form disks (except for a small number of major mergers and highly star-forming systems), whereas below it a galaxy may or may not form a disk.

preprint2015arXiv

Hierarchical Star Formation across the ring galaxy NGC 6503

We present a detailed clustering analysis of the young stellar population across the star-forming ring galaxy NGC 6503, based on the deep HST photometry obtained with the Legacy ExtraGalactic UV Survey (LEGUS). We apply a contour-based map analysis technique and identify in the stellar surface density map 244 distinct star-forming structures at various levels of significance. These stellar complexes are found to be organized in a hierarchical fashion with 95% being members of three dominant super-structures located along the star-forming ring. The size distribution of the identified structures and the correlation between their radii and numbers of stellar members show power-law behaviors, as expected from scale-free processes. The self-similar distribution of young stars is further quantified from their autocorrelation function, with a fractal dimension of ~1.7 for length-scales between ~20 pc and 2.5 kpc. The young stellar radial distribution sets the extent of the star-forming ring at radial distances between 1 and 2.5 kpc. About 60% of the young stars belong to the detected stellar structures, while the remaining stars are distributed among the complexes, still inside the ring of the galaxy. The analysis of the time-dependent clustering of young populations shows a significant change from a more clustered to a more distributed behavior in a time-scale of ~60 Myr. The observed hierarchy in stellar clustering is consistent with star formation being regulated by turbulence across the ring. The rotational velocity difference between the edges of the ring suggests shear as the driving mechanism for this process. Our findings reveal the interesting case of an inner ring forming stars in a hierarchical fashion.

preprint2015arXiv

LEGUS Discovery of a Light Echo Around Supernova 2012aw

We have discovered a luminous light echo around the normal Type II-Plateau Supernova (SN) 2012aw in Messier 95 (M95; NGC 3351), detected in images obtained approximately two years after explosion with the Wide Field Channel 3 on-board the Hubble Space Telescope (HST) by the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). The multi-band observations span from the near-ultraviolet through the optical (F275W, F336W, F438W, F555W, and F814W). The apparent brightness of the echo at the time was ~21--22 mag in all of these bands. The echo appears circular, although less obviously as a ring, with an inhomogeneous surface brightness, in particular, a prominent enhanced brightness to the southeast. The SN itself was still detectable, particularly in the redder bands. We are able to model the light echo as the time-integrated SN light scattered off of diffuse interstellar dust in the SN environment. We have assumed that this dust is analogous to that in the Milky Way with R_V=3.1. The SN light curves that we consider also include models of the unobserved early burst of light from the SN shock breakout. Our analysis of the echo suggests that the distance from the SN to the scattering dust elements along the echo is ~45 pc. The implied visual extinction for the echo-producing dust is consistent with estimates made previously from the SN itself. Finally, our estimate of the SN brightness in F814W is fainter than that measured for the red supergiant star at the precise SN location in pre-SN images, possibly indicating that the star has vanished and confirming it as the likely SN progenitor.

preprint2015arXiv

On the Mass-Metallicity-Star Formation Rate Relation for Galaxies at $z\sim 2$

Recent studies have shown that the local mass-metallicity (M-Z) relation depends on the specific star formation rate (SSFR). Whether such a dependence exists at higher redshifts, and whether the resulting M-Z-SFR relation is redshift invariant, is debated. We re-examine these issues by applying the non-parametric techniques of Salim et al. (2014) to ~130 $z\sim2.3$ galaxies with N2 and O3 measurements from KBSS (Steidel et al. 2014). We find that the KBSS M-Z relation depends on SSFR at intermediate masses, where such dependence exists locally. KBSS and SDSS galaxies of the same mass and SSFR ("local analogs") are similarly offset in the BPT diagram relative to the bulk of local star-forming galaxies, and thus we posit that metallicities can be compared self-consistently at different redshifts as long as the masses and SSFRs of the galaxies are similar. We find that the M-Z-SFR relation of $z\sim2$ galaxies is consistent with the local one at $\log M_*<10$, but is offset up to -0.25 dex at higher masses, so it is altogether not redshift invariant. This high-mass offset could arise from a bias that high-redshift spectroscopic surveys have against high-metallicity galaxies, but additional evidence disfavors this possibility. We identify three causes for the reported discrepancy between N2 and O3N2 metallicities at $z\sim2$: (1) a smaller offset that is also present for SDSS galaxies, which we remove with new N2 calibration, (2) a genuine offset due to differing ISM condition, which is also present in local analogs, (3) an additional offset due to unrecognized AGN contamination.

preprint2015arXiv

Star Cluster Properties in Two LEGUS Galaxies Computed with Stochastic Stellar Population Synthesis Models

We investigate a novel Bayesian analysis method, based on the Stochastically Lighting Up Galaxies (slug) code, to derive the masses, ages, and extinctions of star clusters from integrated light photometry. Unlike many analysis methods, slug correctly accounts for incomplete IMF sampling, and returns full posterior probability distributions rather than simply probability maxima. We apply our technique to 621 visually-confirmed clusters in two nearby galaxies, NGC 628 and NGC 7793, that are part of the Legacy Extragalactic UV Survey (LEGUS). LEGUS provides Hubble Space Telescope photometry in the NUV, U, B, V, and I bands. We analyze the sensitivity of the derived cluster properties to choices of prior probability distribution, evolutionary tracks, IMF, metallicity, treatment of nebular emission, and extinction curve. We find that slug's results for individual clusters are insensitive to most of these choices, but that the posterior probability distributions we derive are often quite broad, and sometimes multi-peaked and quite sensitive to the choice of priors. In contrast, the properties of the cluster population as a whole are relatively robust against all of these choices. We also compare our results from slug to those derived with a conventional non-stochastic fitting code, Yggdrasil. We show that slug's stochastic models are generally a better fit to the observations than the deterministic ones used by Yggdrasil. However, the overall properties of the cluster populations recovered by both codes are qualitatively similar.

preprint2014arXiv

A Critical Look at the Mass-Metallicity-SFR Relation in the Local Universe. I. An Improved Analysis Framework and Confounding Systematics

It has been proposed that the mass-metallicity relation of galaxies exhibits a secondary dependence on star formation rate (SFR), and that the resulting M-Z-SFR relation may be redshift-invariant, i.e., "fundamental." However, conflicting results on the character of the SFR dependence, and whether it exists, have been reported. To gain insight into the origins of the conflicting results, we (a) devise a non-parametric, astrophysically motivated analysis framework based on the offset from the star-forming ("main") sequence at a given stellar mass (relative specific SFR), (b) apply this methodology and perform a comprehensive re-analysis of the local M-Z-SFR relation, based on SDSS, GALEX, and WISE data, and (c) study the impact of sample selection, and of using different metallicity and SFR indicators. We show that metallicity is anti-correlated with specific SFR regardless of the indicators used. We do not find that the relation is spurious due to correlations arising from biased metallicity measurements, or fiber aperture effects. We emphasize that the dependence is weak/absent for massive galaxies ($\log M_*>10.5$), and that the overall scatter in the M-Z-SFR relation does not greatly decrease from the M-Z relation. We find that the dependence is stronger for the highest SSFR galaxies above the star-forming sequence. This two-mode behavior can be described with a broken linear fit in 12+log(O/H) vs. log (SFR$/M_*$), at a given $M_*$. Previous parameterizations used for comparative analysis with higher redshift samples that do not account for the more detailed behavior of the local M-Z-SFR relation may incorrectly lead to the conclusion that those samples follow a different relationship.

preprint2014arXiv

Empirical ugri-UBVRc Transformations for Galaxies

We present empirical color transformations between Sloan Digital Sky Survey ugri and Johnson-Cousins UBVRc photometry for nearby galaxies (D < 11 Mpc). We use the Local Volume Legacy (LVL) galaxy sample where there are 90 galaxies with overlapping observational coverage for these two filter sets. The LVL galaxy sample consists of normal, non-starbursting galaxies. We also examine how well the LVL galaxy colors are described by previous transformations derived from standard calibration stars and model-based galaxy templates. We find significant galaxy color scatter around most of the previous transformation relationships. In addition, the previous transformations show systematic offsets between transformed and observed galaxy colors which are visible in observed color-color trends. The LVL-based $galaxy$ transformations show no systematic color offsets and reproduce the observed color-color galaxy trends.

preprint2014arXiv

Spitzer Local Volume Legacy (LVL) SEDs and Physical Properties

We present the panchromatic spectral energy distributions (SEDs) of the Local Volume Legacy (LVL) survey which consists of 258 nearby galaxies ($D<$11 Mpc). The wavelength coverage spans the ultraviolet to the infrared (1500 $\textrmÅ$ to 24 $μ$m) which is utilized to derive global physical properties (i.e., star formation rate, stellar mass, internal extinction due to dust.). With these data, we find color-color relationships and correlated trends between observed and physical properties (i.e., optical magnitudes and dust properties, optical color and specific star formation rate, and ultraviolet-infrared color and metallicity). The SEDs are binned by different galaxy properties to reveal how each property affects the observed shape of these SEDs. In addition, due to the volume-limited nature of LVL, we utilize the dwarf-dominated galaxy sample to test star formation relationships established with higher-mass galaxy samples. We find good agreement with the star-forming "main-sequence" relationship, but find a systematic deviation in the infrared "main-sequence" at low luminosities. This deviation is attributed to suppressed polycyclic aromatic hydrocarbon (PAH) formation in low metallicity environments and/or the destruction of PAHs in more intense radiation fields occurring near a suggested threshold in sSFR at a value of log($sSFR$) $\sim$ $-$10.2.

preprint2014arXiv

The Relationship between Stellar Mass, Gas Metallicity, and Star Formation Rate for Halpha-selected Galaxies at z~0.8 from the NewHalpha Survey

Using a sample of 299 Ha-selected galaxies at z~0.8, we study the relationship between galaxy stellar mass, gas-phase metallicity, and star formation rate (SFR), and compare to previous results. We use deep optical spectra obtained with the IMACS spectrograph at the Magellan telescope to measure strong oxygen lines. We combine these spectra and metallicities with (1) rest-frame UV-to-optical imaging, which allows us to determine stellar masses and dust attenuation corrections, and (2) Ha narrowband imaging, which provides a robust measure of the instantaneous SFR. Our sample spans stellar masses of 10^9 to 6*10^11 solar masses, SFRs of 0.4 to 270 solar masses per year, and metal abundances of 12+log(O/H)~8.3-9.1 (~0.4-2.6 solar metallicity). The correlations that we find between the Ha-based SFR and stellar mass (i.e., the star-forming "main sequence"), and between the stellar mass and metallicity, are both consistent with previous z~1 studies of star-forming galaxies. We then study the relationship between the three properties using various plane-fitting techniques (Lara-Lopez et al.) and a curve-fitting projection (Mannucci et al.). In all cases, we exclude strong dependence of the M-Z relation on SFR, but are unable to distinguish between moderate and no dependence. Our results are consistent with previous mass-metallicity-SFR studies. We check whether dataset limitations may obscure a strong dependence on the SFR by using mock samples drawn from the SDSS. These experiments reveal that the adopted signal-to-noise cuts may have a significant effect on the measured dependence. Further work is needed to investigate these results, and to test whether a "fundamental metallicity relation" or a "fundamental plane" describes star-forming galaxies across cosmic time.

preprint2014arXiv

The Spitzer Local Volume Legacy (LVL) Global Optical Photometry

We present the global optical photometry of 246 galaxies in the Local Volume Legacy (LVL) survey. The full volume-limited sample consists of 258 nearby (D < 11 Mpc) galaxies whose absolute B-band magnitude span a range of -9.6 < M_B < -20.7 mag. A composite optical (UBVR) data set is constructed from observed UBVR and SDSS ugriz imaging, where the ugriz magnitudes are transformed into UBVR. We present photometry within three galaxy apertures defined at UV, optical, and IR wavelengths. Flux comparisons between these apertures reveal that the traditional optical R25 galaxy apertures do not fully encompass extended sources. Using the larger IR apertures we find color-color relationships where later-type spiral and irregular galaxies tend to be bluer than earlier-type galaxies. These data provide the missing optical emission from which future LVL studies can construct the full panchromatic (UV-optical-IR) spectral energy distributions.

preprint2013arXiv

First Spectroscopic Evidence for High Ionization State and Low Oxygen Abundance in Lya Emitters

We present results from Keck/NIRSPEC and Magellan/MMIRS follow-up spectroscopy of Lya emitters (LAEs) at z=2.2 identified in our Subaru narrowband survey. We successfully detect Ha emission from seven LAEs, and perform a detailed analysis of six LAEs free from AGN activity, two out of which, CDFS-3865 and COSMOS-30679, have [OII] and [OIII] line detections. They are the first [OII]-detected LAEs at high-z, and their [OIII]/[OII] ratios and R23-indices provide the first simultaneous determinations of ionization parameter and oxygen abundance for LAEs. CDFS-3865 has a very high ionization parameter (q_{ion}=2.5^{+1.7}_{-0.8}x10^8 cm s^{-1}) and a low oxygen abundance (12+log(O/H)=7.84^{+0.24}_{-0.25}) in contrast with moderate values of other high-z galaxies such as Lyman-break galaxies (LBGs). COSMOS-30679 also possesses a relatively high ionization parameter (q_{ion}=8^{+10}_{-4}x10^7 cm s^{-1}) and a low oxygen abundance (12+log(O/H)=8.18^{+0.28}_{-0.28}). Both LAEs appear to fall below the mass-metallicity relation of z~2 LBGs. Similarly, a low metallicity of 12+log(O/H)<8.4 is independently indicated for typical LAEs from a composite spectrum and the [NII]/Ha index. Such high ionization parameters and low oxygen abundances can be found in local star-forming galaxies, but this extreme local population occupies only ~0.06% of the SDSS spectroscopic galaxy sample with a number density ~100 times smaller than that of LAEs. With their high ionization parameters and low oxygen abundances, LAEs would represent an early stage of galaxy formation dominated by massive stars in compact star-forming regions. High-q_{ion} galaxies like LAEs would produce ionizing photons efficiently with a high escape fraction achieved by density-bounded HII regions, which would significantly contribute to cosmic reionization at z>6.

preprint2013arXiv

Measuring Galaxy Star Formation Rates From Integrated Photometry: Insights from Color-Magnitude Diagrams of Resolved Stars

We use empirical star formation histories (SFHs), measured from HST-based resolved star color-magnitude diagrams, as input into population synthesis codes to model the broadband spectral energy distributions (SEDs) of ~50 nearby dwarf galaxies (6.5 < log M/M_* < 8.5, with metallicities ~10% solar). In the presence of realistic SFHs, we compare the modeled and observed SEDs from the ultraviolet (UV) through near-infrared (NIR) and assess the reliability of widely used UV-based star formation rate (SFR) indicators. In the FUV through i bands, we find that the observed and modeled SEDs are in excellent agreement. In the Spitzer 3.6micron and 4.5micron bands, we find that modeled SEDs systematically over-predict observed luminosities by up to ~0.2 dex, depending on treatment of the TP-AGB stars in the synthesis models. We assess the reliability of UV luminosity as a SFR indicator, in light of independently constrained SFHs. We find that fluctuations in the SFHs alone can cause factor of ~2 variations in the UV luminosities relative to the assumption of a constant SFH over the past 100 Myr. These variations are not strongly correlated with UV-optical colors, implying that correcting UV-based SFRs for the effects of realistic SFHs is difficult using only the broadband SED. Additionally, for this diverse sample of galaxies, we find that stars older than 100 Myr can contribute from <5% to100% of the present day UV luminosity, highlighting the challenges in defining a characteristic star formation timescale associated with UV emission. We do find a relationship between UV emission timescale and broadband UV-optical color, though it is different than predictions based on exponentially declining SFH models. Our findings have significant implications for the comparison of UV-based SFRs across low-metallicity populations with diverse SFHs.

preprint2013arXiv

The Nature of the Second Parameter in the IRX-beta Relation for Local Galaxies

We present an analysis of 98 galaxies of low-dust content, selected from the Spitzer Local Volume Legacy survey, aimed at examining the relation between the ultraviolet (UV) color, beta, and dust attenuation in normal star-forming galaxies. The IRX-beta diagram relates the total dust attenuation in a galaxy, traced by the far-IR (FIR) to UV ratio, to the observed UV color, indicated by beta. Previous research has indicated that while starburst galaxies exhibit a relatively tight IRX-beta relation, normal star-forming galaxies do not, and have a much larger spread in the total-IR to far-UV (FUV) luminosity for a fixed UV color. We examine the role that the age of the stellar population plays as the "second parameter" responsible for the observed deviation and spread of star-forming galaxies from the starburst relation. We model the FUV to FIR spectral energy distribution (SED) of each galaxy according to two broad bins of star formation history (SFH): constant and instantaneous burst. We find clear trends between stellar population mean age estimators (extinction-corrected FUV/NIR, U-B, and EW(H-alpha)) and the UV color, beta; the trends are mostly driven by the galaxies best-described by instantaneous burst populations. We also find a significant correlation between beta and the mean age directly determined from the best-fit instantaneous models. As already indicated by other authors, the UV attenuation in star-forming galaxies may not be recovered with the UV color alone and is highly influenced by the stellar population's mean age and SFH. Overall, the scatter in the IRX-beta diagram is better correlated with beta than with the perpendicular distance.

preprint2012arXiv

A Dual Narrowband Survey for Hα Emitters at z=2.2: Demonstration of the Technique and Constraints on the Hα Luminosity Function

We present first results from a narrowband imaging program for intermediate redshift emission-line galaxies using the newly commissioned FourStar infrared camera at the 6.5m Magellan telescope. To enable prompt identification of Hα emitters, a pair of custom 1% filters, which sample low-airglow atmospheric windows at 1.19 μm and 2.10 μm, is used to detect both Hα and [OII]λ3727 emission from the same redshift volume at z=2.2. Initial observations are taken over a 130 arcmin^2 area in the CANDELS-COSMOS field. The exquisite image quality resulting from the combination of the instrument, telescope, and standard site conditions (~0.55" FWHM) allows the 1.19 μm and 2.10 μm data to probe 3σ emission-line depths down to 1.0e-17 erg/s/cm^2 and 1.2e-17 erg/s/cm^2 respectively, in less than 10 hours of integration time in each narrowband. For Hα at z=0.8 and z=2.2, these fluxes correspond to observed star formation rates of ~0.3 and ~4 Msun/yr respectively. We find 122 sources with a 1.19 μm excess, and 136 with a 2.10 μm excess, 41 of which show an excess in both bands. The dual narrowband technique, as implemented here, is estimated to identify about >80% of z=2.2 Hα emitters in the narrowband excess population. With the most secure such sample obtained to-date, we compute constraints on the faint-end slope of the z=2.2 Hα luminosity function. These "narrow-deep" FourStar observations have been obtained as part of the larger NewHα Survey, which will combine the data with "wide-shallow" imaging through a similar narrowband filter pair with NEWFIRM at the KPNO/CTIO 4m telescopes, to enable study of both luminous (but rare) and faint emission-line galaxies in the intermediate redshift universe. [Abridged]

preprint2012arXiv

Direct Oxygen Abundances for Low Luminosity LVL Galaxies

We present MMT spectroscopic observations of HII regions in 42 low luminosity galaxies in the LVL. For 31 galaxies, we measured the temperature sensitive [O III] line at a strength of 4 sigma or greater, and thus determine direct oxygen abundances. Our results provide the first direct estimates of oxygen abundance for 19 galaxies. Oxygen abundances were compared to B-band and 4.5 micron luminosities and stellar masses in order to characterize the luminosity-metallicity (L-Z) and mass-metallicity (M-Z) relationships at low-luminosity. We present and analyze a "Combined Select" sample composed of 38 objects (drawn from our parent sample and the literature) with direct oxygen abundances and reliable distance determinations (TRGB or Ceph). Consistent with previous studies, the B-band and 4.5 micron L-Z relationships were found to be 12+log(O/H)=(6.27+/-0.21)+(-0.11+/-0.01)M_B and 12+log(O/H)=(6.10+/-0.21)+(-0.10+/-0.01)M_[4.5] (sigma=0.15 and 0.14). For this sample, we derive a M-Z relationship of 12+log(O/H)=(5.61+/-0.24)+(0.29+/-0.03)log(M*), which agrees with previous studies; however, the dispersion (sigma=0.15) is not significantly lower than that of the L-Z relationships. Because of the low dispersions in these relationships, if an accurate distance is available, the luminosity of a low-luminosity galaxy is often a better indicator of metallicity than that derived using certain strong-line methods, so significant departures from the L-Z relationships may indicate that caution is prudent in such cases. We also revisit the 70/160 micron color metallicity relationship. Additionally, we examine N/O abundance trends with respect to oxygen abundance and B-V color. We find a positive correlation between N/O ratio and B-V color for 0.05\lesssimB-V\lesssim0.75: log(N/O)=(1.18+/-0.9)x(B-V)+(-1.92+/-0.08), with a dispersion of sigma=0.14, that is in agreement with previous studies.

preprint2012arXiv

Nebular Attenuation in Hα-selected Star-forming Galaxies at z=0.8 from the NewHα Survey

We present measurements of the dust attenuation of Hα-selected emission-line galaxies at z=0.8 from the NewHα narrowband survey. The analysis is based on deep follow-up spectroscopy with Magellan/IMACS, which captures the strong rest-frame optical emission lines from [OII] λ3727 to [OIII] λ5007. The spectroscopic sample used in this analysis consists of 341 confirmed Hα emitters. We place constraints on the AGN fraction using diagnostics which can be applied at intermediate redshift. We find that at least 5% of the objects in our spectroscopic sample can be classified as AGN and 2% are composite, i.e. powered by a combination of star-formation and AGN activity. We measure the dust attenuation for individual objects from the ratios of the higher order Balmer lines. The Hβ and Hγ pair of lines is detected with S/N>5 in 55 individual objects and the Hβ and Hδ pair is detected in 50 individual objects. We also create stacked spectra to probe the attenuation in objects without individual detections. The median attenuation at Hα based on the objects with individually detected lines is A(Hα)=0.9+-1.0 magnitudes, in good agreement with the attenuation found in local samples of star-forming galaxies. We find that the z=0.8 galaxies occupy a similar locus of attenuation as a function of magnitude, mass and SFR as a comparison sample drawn from the SDSS DR4. Both the results from the individual z=0.8 galaxies and from the stacked spectra show consistency with the mass -- attenuation and SFR -- attenuation relations found in the local Universe, indicating that these relations are also applicable at intermediate redshift.

preprint2012arXiv

Star Formation Rate Distributions: Inadequacy of the Schechter Function

In this paper we posit that galaxy luminosity functions (LFs) come in two fundamentally different types depending on whether the luminosity traces galaxy stellar mass or its current star formation rate (SFR). Mass function types reflect the older stars and therefore the stellar mass distribution, while SFR function types arise from the young stars and hence the distribution of SFRs. Optical and near-infrared LFs are of the mass function type, and are well fit by a Schechter function (power law with an exponential cutoff at the bright end). In contrast, LFs of the SFR function type are of a different form, one that cannot be adequately described by a Schechter function. We demonstrate this difference by generating SFR distributions for mock samples of galaxies drawn from a Schechter stellar mass distribution along with established empirical relations between the SFR and stellar mass. Compared with the Schechter function, SFR distributions have a shallower decline at the bright end, which can be traced to the large intrinsic scatter of SFRs at any given stellar mass. A superior description of SFR distributions is given by the "Saunders" function, which combines a power law with a Gaussian at the high end. We show that the Schechter-like appearance of UV and H alpha LFs, although they are LFs of SFR function type, results when luminosities are not corrected for dust, or when average statistical corrections are used because individual attenuation measurements are not available. We thus infer that the non-Schechter form of the far-IR LFs is a true reflection of the underlying SFR distribution, rather than the purported artifact of AGN contamination.

preprint2011arXiv

Average Metallicity and Star Formation Rate of Lya Emitters Probed by a Triple Narrow-Band Survey

We present the average metallicity and star-formation rate of Lya emitters (LAEs) measured from our large-area survey with three narrow-band (NB) filters covering the Lya, [OII]3727, and Ha+[NII] lines of LAEs at z=2.2. We select 919 z=2.2 LAEs from Subaru/Suprime-Cam NB data in conjunction with Magellan/IMACS spectroscopy. Of these LAEs, 561 and 105 are observed with KPNO/NEWFIRM near-infrared NB filters whose central wavelengths are matched to redshifted [OII] and Ha nebular lines, respectively. By stacking the near-infrared images of the LAEs, we successfully obtain average nebular-line fluxes of LAEs, the majority of which are too faint to be identified individually by narrow-band imaging or deep spectroscopy. The stacked object has an Ha luminosity of 1.7x10^{42} erg s^{-1} corresponding to a star formation rate (SFR) of 14 M_{sun} yr^{-1}. We place, for the first time, a firm lower limit to the average metallicity of LAEs of Z>~0.09 Z_{sun} (2sigma) based on the [OII]/(Ha+[NII]) index together with photo-ionization models and empirical relations. This lower limit of metallicity rules out the hypothesis that LAEs, so far observed at z~2, are extremely metal poor (Z<2x10^{-2} Z_{sun}) young galaxies at the 4sigma level. This limit is higher than a simple extrapolation of the observed mass-metallicity relation of z~2 UV-selected galaxies toward lower masses (5x10^{8} M_{sun}), but roughly consistent with a recently proposed fundamental mass-metallicity relation when the LAEs' relatively low SFR is taken into account. The Ha and Lya luminosities of our NB-selected LAEs indicate that the escape fraction of Lya photons is ~12-30 %, much higher than the values derived for other galaxy populations at z~2.

preprint2011arXiv

Modeling the Effects of Star Formation Histories on Halpha and Ultra-Violet Fluxes in Nearby Dwarf Galaxies

We consider the effects of non-constant star formation histories (SFHs) on Halpha and GALEX far ultra-violet (FUV) star formation rate (SFR) indicators. Under the assumption of a fully populated Chabrier IMF, we compare the distribution of Halpha-to-FUV flux ratios from ~ 1500 simple, periodic model SFHs with observations of 185 galaxies from the Spitzer Local Volume Legacy survey. We find a set of SFH models that are well matched to the data, such that more massive galaxies are best characterized by nearly constant SFHs, while low mass systems experience bursts amplitudes of ~ 30 (i.e., an increase in the SFR by a factor of 30 over the SFR during the inter-burst period), burst durations of tens of Myr, and periods of ~ 250 Myr; these SFHs are broadly consistent with the increased stochastic star formation expected in systems with lower SFRs. We analyze the predicted temporal evolution of galaxy stellar mass, R-band surface brightness, Halpha-derived SFR, and blue luminosity, and find that they provide a reasonable match to observed flux distributions. We find that our model SFHs are generally able to reproduce both the observed systematic decline and increased scatter in Halpha-to-FUV ratios toward low mass systems, without invoking other physical mechanisms. We also compare our predictions with those from the Integrated Galactic IMF theory with a constant SFR. We find that while both predict a systematic decline in the observed ratios, only the time variable SFH models are capable of producing the observed population of low mass galaxies ($M_{*}$ < 10$^{7}$ Msun) with normal Halpha-to-FUV ratios. These results demonstrate that a variable IMF alone has difficulty explaining the observed scatter in the Halpha-to-FUV ratios. We conclude by considering the limitations of the model SFHs, and discuss the use of additional empirical constraints to improve future SFH modeling efforts.

preprint2011arXiv

The H-alpha Luminosity Function and Star Formation Rate Volume Density at z=0.8 from the NEWFIRM H-alpha Survey

[Abridged] We present new measurements of the H-alpha luminosity function (LF) and SFR volume density for galaxies at z~0.8. Our analysis is based on 1.18$μ$m narrowband data from the NEWFIRM H-alpha Survey, a comprehensive program designed to capture deep samples of intermediate redshift emission-line galaxies using narrowband imaging in the near-infrared. The combination of depth ($\approx1.9\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$ in H-alpha at 3$σ$) and areal coverage (0.82 deg$^2$) complements other recent H-alpha studies at similar redshifts, and enables us to minimize the impact of cosmic variance and place robust constraints on the shape of the LF. The present sample contains 818 NB118 excess objects, 394 of which are selected as H-alpha emitters. Optical spectroscopy has been obtained for 62% of the NB118 excess objects. Empirical optical broadband color classification is used to sort the remainder of the sample. A comparison of the LFs constructed for the four individual fields reveals significant cosmic variance, emphasizing that multiple, widely separated observations are required. The dust-corrected LF is well-described by a Schechter function with L*=10^{43.00\pm0.52} ergs s^{-1}, ϕ*=10^{-3.20\pm0.54} Mpc^{-3}, and α=-1.6\pm0.19. We compare our H-alpha LF and SFR density to those at z<1, and find a rise in the SFR density \propto(1+z)^{3.4}, which we attribute to significant L* evolution. Our H-alpha SFR density of 10^{-1.00\pm0.18} M_sun yr^{-1} Mpc^{-3} is consistent with UV and [O II] measurements at z~1. We discuss how these results compare to other H-alpha surveys at z~0.8, and find that the different methods used to determine survey completeness can lead to inconsistent results. This suggests that future surveys probing fainter luminosities are needed, and more rigorous methods of estimating the completeness should be adopted as standard procedure.

preprint2010arXiv

A GALEX Ultraviolet Imaging Survey of Galaxies in the Local Volume

We present results from a GALEX ultraviolet (UV) survey of a complete sample of 390 galaxies within ~11 Mpc of the Milky Way. The UV data are a key component of the composite Local Volume Legacy (LVL), an ultraviolet-to-infrared imaging program designed to provide an inventory of dust and star formation in nearby spiral and irregular galaxies. The ensemble dataset is an especially valuable resource for studying star formation in dwarf galaxies, which comprise over 80% of the sample. We describe the GALEX survey programs which obtained the data and provide a catalog of far-UV (~1500 Angstroms) and near-UV (~2200 Angstroms) integrated photometry. General UV properties of the sample are briefly discussed. We compute two measures of the global star formation efficiency, the SFR per unit HI gas mass and the SFR per unit stellar mass, to illustrate the significant differences that can arise in our understanding of dwarf galaxies when the FUV is used to measure the SFR instead of H-alpha. We find that dwarf galaxies may not be as drastically inefficient in coverting gas into stars as suggested by prior H-alpha studies. In this context, we also examine the UV properties of late-type dwarf galaxies that appear to be devoid of star formation because they were not detected in previous H-alpha narrowband observations. Nearly all such galaxies in our sample are detected in the FUV, and have FUV SFRs that fall below the limit where the H-alpha flux is robust to Poisson fluctuations in the formation of massive stars. The UV colors and star formation efficiencies of H-alpha-undetected, UV-bright dwarf irregulars appear to be relatively unremarkable with respect to those exhibited by the general population of star-forming galaxies.

preprint2010arXiv

A New Method for Measuring the Upper End of the IMF

A method is presented here for investigating variations in the upper end of the stellar Initial Mass Function (IMF) by probing the production rate of ionizing photons in unresolved, compact star clusters with ages<10 Myr and covering a range of masses. We test this method on the young cluster population in the nearby galaxy M51a, for which multi-wavelength observations from the Hubble Space Telescope are available. Our results indicate that the proposed method can probe the upper end of the IMF in galaxies located out to at least 10 Mpc, i.e., a factor 200 further away than possible by counting individual stars in young compact clusters. Our results for this galaxy show no obvious dependence of the upper mass end of the IMF on the mass of the star cluster, down to ~1000 M_sun, although more extensive analyses involving lower mass clusters and other galaxies are needed to confirm this conclusion.

preprint2010arXiv

The Wyoming Survey for H-alpha. III. A Multi-wavelength Look at Attenuation by Dust in Galaxies out to z~0.4

We report results from the Wyoming Survey for H-alpha (WySH), a comprehensive four-square degree survey to probe the evolution of star-forming galaxies over the latter half of the age of the Universe. We have supplemented the H-alpha data from WySH with infrared data from the Spitzer Wide-area Infrared Extragalactic (SWIRE) Survey and ultraviolet data from the Galaxy Evolution Explorer (GALEX) Deep Imaging Survey. This dataset provides a multi-wavelength look at the evolution of the attenuation by dust, and here we compare a traditional measure of dust attenuation (L(TIR)/L(FUV)) to a diagnostic based on a recently-developed robust star formation rate (SFR) indicator, [H-alpha_obs+24-micron]/H-alpha_obs. With such data over multiple epochs, the evolution in the attenuation by dust with redshift can be assessed. We present results from the ELAIS-N1 and Lockman Hole regions at z~0.16, 0.24, 0.32 and 0.40. While the ensemble averages of both diagnostics are relatively constant from epoch to epoch, each epoch individually exhibits a larger attenuation by dust for higher star formation rates. Hence, an epoch to epoch comparison at a fixed star formation rate suggests a mild decrease in dust attenuation with redshift.

preprint2010arXiv

UV/H-alpha Turmoil

A great deal of our understanding of star formation in the local universe has been built upon an extensive foundation of H-alpha observational studies. However, recent work in the ultraviolet (UV) with GALEX has shown that star formation rates (SFRs) inferred from H-alpha in galactic environments characterized by low stellar and gas densities tend to be less than those based on the UV luminosity. The origin of the discrepancy is actively debated because one possible explanation is that the stellar initial mass function is systematically deficient in high mass stars in such environments. In this contribution, we summarize our work on this topic using a dwarf galaxy dominated sample of ~300 late-type galaxies in the 11 Mpc Local Volume. The sample allows us to examine the discrepancy between H-alpha and UV SFRs using a statistical number of galaxies with activities less than 0.1 Msun/yr. A range of potential causes for such an effect are reviewed. We find that while the IMF hypothesis is not inconsistent with our observations, alternate explanations remain that must be investigated further before a final conclusion can be drawn.

preprint2009arXiv

Chemical Abundances of Seven Irregular and Three Tidal Dwarf Galaxies in the M81 Group

We have derived nebular abundances for 10 dwarf galaxies belonging to the M81 Group, including several galaxies which do not have abundances previously reported in the literature. For each galaxy, multiple H \ii regions were observed with GMOS-N at the Gemini Observatory in order to determine abundances of several elements (oxygen, nitrogen, sulfur, neon, and argon). For seven galaxies, at least one H \ii region had a detection of the temperature sensitive [OIII] $λ$4363 line, allowing a "direct" determination of the oxygen abundance. No abundance gradients were detected in the targeted galaxies and the observed oxygen abundances are typically in agreement with the well known metallicity-luminosity relation. However, three candidate "tidal dwarf" galaxies lie well off this relation, UGC 5336, Garland, and KDG 61. The nature of these systems suggests that UGC 5336 and Garland are indeed recently formed systems, whereas KDG 61 is most likely a dwarf spheroidal galaxy which lies along the same line of sight as the M81 tidal debris field. We propose that these H \ii regions formed from previously enriched gas which was stripped from nearby massive galaxies (e.g., NGC 3077 and M81) during a recent tidal interaction.

preprint2009arXiv

Comparison of H-alpha and UV Star Formation Rates in the Local Volume: Systematic Discrepancies for Dwarf Galaxies

(abridged) Using a complete sample of ~300 star-forming galaxies within 11 Mpc, we evaluate the consistency between star formation rates (SFRs) inferred from the far ultraviolet (FUV) non-ionizing continuum and H-alpha nebular emission, assuming standard conversion recipes in which the SFR scales linearly with luminosity at a given wavelength. Our analysis probes SFRs over 5 orders of magnitude, down to ultra-low activities on the order of ~0.0001 M_sun/yr. The data are drawn from the 11 Mpc H-alpha and Ultraviolet Galaxy Survey (11HUGS), which has obtained H-alpha fluxes from ground-based narrowband imaging, and UV fluxes from imaging with GALEX. For normal spiral galaxies (SFR~1 M_sun/yr), our results are consistent with previous work which has shown that FUV SFRs tend to be lower than H-alpha SFRs before accounting for internal dust attenuation, but that there is relative consistency between the two tracers after proper corrections are applied. However, a puzzle is encountered at the faint end of the luminosity function. As lower luminosity dwarf galaxies, roughly less active than the Small Magellanic Cloud, are examined, H-alpha tends to increasingly under-predict the SFR relative to the FUV. Although past studies have suggested similar trends, this is the first time this effect is probed with a statistical sample for galaxies with SFR~<0.1 M_sun/yr. A range of standard explanations does not appear to be able to account for the magnitude of the systematic. Some recent work has argued for an IMF which is deficient in high mass stars in dwarf and low surface brightness galaxies, and we also consider this scenario.

preprint2009arXiv

Dust-Corrected Star Formation Rates of Galaxies. I. Combinations of H-alpha and Infrared Tracers

We combine H-alpha emission-line and infrared continuum measurements of two samples of nearby galaxies to derive dust attenuation-corrected star formation rates (SFRs). We use a simple energy balance based method that has been applied previously to HII regions in the Spitzer Infrared Nearby Galaxies Survey (SINGS), and extend the methodology to integrated measurements of galaxies. We find that our composite H-alpha+IR based SFRs are in excellent agreement with attenuation-corrected SFRs derived from integrated spectrophotometry, over the full range of SFRs (0.01 -- 80 solar mass per year) and attenuations (0 -- 2.5 mag) studied. We find that the combination of H-alpha and total infrared luminosities provides the most robust SFR measurements, but combinations of H-alpha measurements with monochromatic luminosities at 24 micron and 8 micron perform nearly as well. The calibrations differ significantly from those obtained for HII regions (Calzetti et al. 2007), with the difference attributable to a more evolved population of stars heating the dust. Our results are consistent with a significant component of diffuse dust (the `IR cirrus' component) that is heated by a non-star-forming population. The same methodology can be applied to [OII]$λ$3727 emission-line measurements, and the radio continuum fluxes of galaxies can be applied in place of IR fluxes when the latter are not available. We assess the precision and systematic reliability of all of these composite methods.

preprint2009arXiv

Stellar Populations of Lyman Alpha Emitters at z = 3 - 4 Based on Deep Large Area Surveys in the Subaru-SXDS/UKIDSS-UDS Field

We investigate the stellar populations of Lyman alpha emitters (LAEs) at z=3.1 and 3.7 in 0.65 deg^2 of the Subaru/XMM-Newton Deep Field, based on rest-frame UV-to-optical photometry obtained from the Subaru/XMM-Newton Deep Survey, the UKIDSS/Ultra Deep Survey, and the Spitzer legacy survey of the UKIDSS/UDS. Among a total of 302 LAEs (224 for z=3.1 and 78 for z=3.7), only 11 are detected in the K band, i.e., brighter than K(3sigma)=24.1 mag. Eight of the 11 K-detected LAEs are spectroscopically confirmed. We find that the K-undetected objects, which should closely represent the LAE population as a whole, have low stellar masses of ~ 10^8 - 10^8.5 Msun, modest SFRs of 1 - 100 Msun yr^-1, and modest dust extinction of E(B-V) < 0.2. The K-detected objects are massive, Mstar ~ 10^9 - 10^10.5 Msun, and have significant dust extinction with a median of E(B-V) ~= 0.3. Four K-detected objects with the reddest spectral energy distributions, two of which are spectroscopically confirmed, are heavily obscured with E(B-V) ~ 0.65, and their continua resemble those of some local ULIRGs. Interestingly, they have large Lyman alpha equivalent widths ~= 70 - 250 A. If these four are excluded, our sample has a weak anti-correlation between EW(Lya) and Mstar. We compare the stellar masses and the specific star formation rates (sSFR) of LAEs with those of Lyman-break galaxies, distant red galaxies, submillimetre galaxies, and I- or K-selected galaxies with z_phot ~ 3. We find that the LAE population is the least massive among all the galaxy populations in question, but with relatively high sSFRs, while NIR-detected LAEs have Mstar and sSFR similar to LBGs. Our reddest four LAEs have very high sSFRs in spite of large Mstar, thus occupying a unique region in the Mstar versus sSFR space. (abridged)

preprint2008arXiv

Dwarf Galaxy Starburst Statistics in the Local Volume

An unresolved question in galaxy evolution is whether the star formation histories of low mass systems are preferentially dominated by starbursts or modes that are more quiescent and continuous. Here, we quantify the prevalence of global starbursts in dwarf galaxies at the present epoch, and infer their characteristic durations and amplitudes. The analysis is based on the H-alpha component of the 11 Mpc H-alpha UV Galaxy Survey (11HUGS), which is providing H-alpha and GALEX UV imaging for an approximately volume-limited sample of ~300 star-forming galaxies within 11 Mpc. We first examine the completeness properties of the sample, and then directly tally the number of bursting dwarfs and compute the fraction of star formation that is concentrated in such systems. Our results are consistent with a picture where dwarfs that are currently experiencing massive global bursts are just the ~6% tip of a low-mass galaxy iceberg. Moreover, bursts are only responsible for about a quarter of the total star formation in the overall dwarf population, so the majority of stars in low-mass systems are not formed in this mode today. Spirals and irregulars devoid of H-alpha emission are rare, indicating that the complete cessation of star formation generally does not occur in such galaxies and is not characteristic of the inter-burst state, at least for the more luminous systems with M(B)<-15. The starburst statistics presented here directly constrain the duty cycle and the average burst amplitude under the simplest assumptions where all dwarf irregulars share a common star formation history and undergo similar bursts cycles with equal probability. Uncertainties in such assumptions are discussed in the context of previous work.

preprint2005arXiv

Resonant Readout of a Persistent Current Qubit

We have implemented a resonant circuit that uses a SQUID as a flux-sensitive Josephson inductor for qubit readout. In contrast to the conventional switching current measurement that generates undesired quasi-particles when the SQUID switches to the voltage state, our approach keeps the readout SQUID biased along the supercurrent branch during the measurement. By incorporating the SQUID inductor in a high-Q resonant circuit, we can distinguish the two flux states of a niobium persistent-current (PC) qubit by observing a shift in the resonant frequency of both the magnitude and the phase spectra. The readout circuit was also characterized in the nonlinear regime to investigate its potential use as a nonlinear amplifier.