Source author record

Evan D. Skillman

Evan D. Skillman 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

101works
5topics
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

101 published item(s)

preprint2023arXiv

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come.

preprint2022arXiv

CLASSY V: The impact of aperture effects on the inferred nebular properties of local star-forming galaxies

Strong nebular emission lines are an important diagnostic tool for tracing the evolution of star-forming galaxies across cosmic time. However, different observational setups can affect these lines, and the derivation of the physical nebular properties. We analyze 12 local star-forming galaxies from the COS Legacy Spectroscopy SurveY (CLASSY) to assess the impact of using different aperture combinations on the determination of the physical conditions and gas-phase metallicity. We compare optical spectra observed with the SDSS aperture, which has a 3" of diameter similar to COS, to IFU and longslit spectra, including new LBT/MODS observations of five CLASSY galaxies. We calculate the reddening, electron densities and temperatures, metallicities, star formation rates, and equivalent widths (EWs). We find that measurements of the electron densities and temperatures, and metallicity remained roughly constant with aperture size, indicating that the gas conditions are relatively uniform for this sample. However, using the IFU observations of 3 galaxies, we find that the E(B-V) values derived from the Balmer ratios decrease ( by up to 53%) with increasing aperture size. The values change most significantly in the center of the galaxies, and level out near the COS aperture diameter of 2.5". We examine the relative contributions from the gas and stars using the H$α$ and [OIII] $λ$5007 EWs as a function of aperture light fraction, but find little to no variations within a given galaxy. These results imply that the optical spectra provide nebular properties appropriate for the FUV CLASSY spectra, even when narrow 1.0" long-slit observations are used.

preprint2022arXiv

The COS Legacy Archive Spectroscopy SurveY (CLASSY) Treasury Atlas

Far-ultraviolet (FUV; ~1200-2000 angstroms) spectra are fundamental to our understanding of star-forming galaxies, providing a unique window on massive stellar populations, chemical evolution, feedback processes, and reionization. The launch of JWST will soon usher in a new era, pushing the UV spectroscopic frontier to higher redshifts than ever before, however, its success hinges on a comprehensive understanding of the massive star populations and gas conditions that power the observed UV spectral features. This requires a level of detail that is only possible with a combination of ample wavelength coverage, signal-to-noise, spectral-resolution, and sample diversity that has not yet been achieved by any FUV spectral database. We present the COS Legacy Spectroscopic SurveY (CLASSY) treasury and its first high level science product, the CLASSY atlas. CLASSY builds on the HST archive to construct the first high-quality (S/N_1500 >~ 5/resel), high-resolution (R~15,000) FUV spectral database of 45 nearby (0.002 < z < 0.182) star-forming galaxies. The CLASSY atlas, available to the public via the CLASSY website, is the result of optimally extracting and coadding 170 archival+new spectra from 312 orbits of HST observations. The CLASSY sample covers a broad range of properties including stellar mass (6.2 < logM_star(M_sol) < 10.1), star formation rate (-2.0 < log SFR (M_sol/yr) < +1.6), direct gas-phase metallicity (7.0 < 12+log(O/H) < 8.8), ionization (0.5 < O_32 < 38.0), reddening (0.02 < E(B-V < 0.67), and nebular density (10 < n_e (cm^-3) < 1120). CLASSY is biased to UV-bright star-forming galaxies, resulting in a sample that is consistent with z~0 mass-metallicity relationship, but is offset to higher SFRs by roughly 2 dex, similar to z >~2 galaxies. This unique set of properties makes the CLASSY atlas the benchmark training set for star-forming galaxies across cosmic time.

preprint2022arXiv

The JWST Resolved Stellar Populations Early Release Science Program I.: NIRCam Flux Calibration

We use globular cluster data from the Resolved Stellar Populations Early Release Science (ERS) program to validate the flux calibration for the Near Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST). We find a significant flux offset between the eight short wavelength detectors, ranging from 1-23% (about 0.01-0.2 mag) that affects all NIRCam imaging observations. We deliver improved zeropoints for the ERS filters and show that alternate zeropoints derived by the community also improve the calibration significantly. We also find that the detector offsets appear to be time variable by up to at least 0.1 mag.

preprint2020arXiv

A Demonstration of Extremely Low Latency $γ$-ray, X-Ray & UV Follow-Up of a Millisecond Radio Transient

We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5σ$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3σ$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science.

preprint2020arXiv

CHAOS IV: Gas-Phase Abundance Trends From The First Four CHAOS Galaxies

The chemical abundances of spiral galaxies, as probed by HII regions across their disks, are key to understanding the evolution of galaxies over a wide range of environments. We present LBT/MODS spectra of 52 HII regions in NGC3184 as part of the CHemical Abundances Of Spirals (CHAOS) project. We explore the direct-method gas-phase abundance trends for the first four CHAOS galaxies, using temperature measurements from one or more auroral line detections in 190 individual HII regions. We find the dispersion in Te-Te relationships is dependent on ionization, as characterized by F_5007/F_3727, and so recommend ionization-based temperature priorities for abundance calculations. We confirm our previous results that [NII] and [SIII] provide the most robust measures of electron temperature in low-ionization zones, while [OIII] provides reliable electron temperatures in high-ionization nebula. We measure relative and absolute abundances for O, N, S, Ar, and Ne. The four CHAOS galaxies marginally conform with a universal O/H gradient, as found by empirical IFU studies when plotted relative to effective radius. However, after adjusting for vertical offsets, we find a tight universal N/O gradient of alpha_N/O = -0.33 dex/Re with sigma_tot. = 0.08 for Rg/Re < 2.0, where N is dominated by secondary production. Despite this tight universal N/O gradient, the scatter in the N/O-O/H relationship is significant. Interestingly, the scatter is similar when N/O is plotted relative to O/H or S/H. The observable ionic states of S probe lower ionization and excitation energies than O, which might be more appropriate for characterizing abundances in metal-rich HII regions.

preprint2020arXiv

CHAOS V: Recombination Line Carbon Abundances in M101

The CHAOS project is building a large database of LBT H II region spectra in nearby spiral galaxies to use direct abundances to better determine the dispersion in metallicity as a function of galactic radius. Here, we present CHAOS LBT observations of C II $λ$4267 emission detected in 10 H II regions in M 101, and, using a new photoionization model based ionization correction factor, we convert these measurements into total carbon abundances. A comparison with M101 C II recombination line observations from the literature shows excellent agreement, and we measure a relatively steep gradient in log(C/H) of -0.37 +/- 0.06 dex/R_e. The C/N observations are consistent with a constant value of log(C/N) = 0.84 with a dispersion of only 0.09 dex, which, given the different nucleosynthetic sources of C and N, is challenging to understand. We also note that when plotting N/O versus O/H, all of the H II regions with detections of CII $λ$4267 present N/O abundances at the minimum of the scatter in N/O at a given value of O/H. If the high surface brightness necessary for the detection of the faint recombination lines is interpreted as an indicator of H II region youth, then this may point to a lack of nitrogen pollution in the youngest H II regions. In the future, we anticipate that the CHAOS project will significantly increase the total number of C II $λ$4267 measurements in extragalactic H II regions.

preprint2020arXiv

First Co-spatial Comparison of Stellar, Neutral-, and Ionized-gas Metallicities in a metal-rich galaxy: M83

We carry out a comparative analysis of the metallicities from the stellar, neutral-gas, and ionized-gas components in the metal-rich spiral galaxy M83. We analyze spectroscopic observations taken with the Hubble Space Telescope (HST), the Large Binocular Telescope (LBT) and the Very Large Telescope (VLT). We detect a clear depletion of the HI gas, as observed from the HI column densities in the nuclear region of this spiral galaxy. We find column densities of log[$N$(HI) cm$^{-2}$] $<$ 20.0 at galactocentric distances of $<$ 0.18 kpc, in contrast to column densities of log[$N$(HI) cm$^{-2}$] $\sim$ 21.0 in the galactic disk, a trend observed in other nearby spiral galaxies. We measure a metallicity gradient of $-$0.03 $\pm$ 0.01 dex kpc$^{-1}$ for the ionized gas, comparable to the metallicity gradient of a local benchmark of 49 nearby star-forming galaxies of $-$0.026 $\pm$ 0.002 dex kpc$^{-1}$. Our co-spatial metallicity comparison of the multi-phase gas and stellar populations shows excellent agreement outside of the nucleus of the galaxy hinting at a scenario where the mixing of newly synthesized metals from the most massive stars in the star clusters takes longer than their lifetimes ($\sim$10 Myr). Finally, our work shows that caution must be taken when studying the metallicity gradient of the neutral-gas component in star-forming galaxies, since this can be strongly biased, as these environments can be dominated by molecular gas. In these regions the typical metallicity tracers can provide inaccurate abundances as they may trace both the neutral- and molecular-gas components.

preprint2020arXiv

PHAT XX. AGB stars and other cool giants in M31 star clusters

The presence of AGB stars in clusters provides key constraints for stellar models, as has been demonstrated with historical data from the Magellanic Clouds. In this work, we look for candidate AGB stars in M31 star clusters from the Panchromatic Hubble Andromeda Treasury (PHAT) survey. Our photometric criteria selects stars brighter than the tip of the red giant branch, which includes the bulk of the thermally-pulsing AGB stars as well as early-AGB stars and other luminous cool giants expected in young stellar populations (e.g. massive red supergiants, and intermediate-mass red helium-burning stars). The AGB stars can be differentiated, a posteriori, using the ages already estimated for our cluster sample. 937 candidates are found within the cluster aperture radii, half (450) of which are very likely cluster members. Cross-matching with additional databases reveals two carbon stars and ten secure variables among them. The field-corrected age distribution reveals the presence of young supergiants peaking at ages smaller than 100 Myr, followed by a long tail of AGB stars extending up to the oldest possible ages. This long tail reveals the general decrease in the numbers of AGB stars from initial values of 50e-6/Msun at 100 Myr down to 5e-6/Msun at 10 Gyr. Theoretical models of near-solar metallicity reproduce this general trend, although with localized discrepancies over some age intervals, whose origin is not yet identified. The entire catalogue is released together with finding charts to facilitate follow-up studies.

preprint2020arXiv

The Leoncino Dwarf Galaxy: Exploring the Low-Metallicity End of the Luminosity-Metallicity and Mass-Metallcity Relations

Extremely metal-poor (XMP) galaxies are low-mass, star-forming galaxies with gas-phase oxygen abundances below 12+log(O/H) = 7.35 (~1/20 Zsun). Galaxy evolution scenarios suggest three pathways to form an XMP: (1) secular evolution at low galaxy masses, (2) slow evolution in voids, or (3) dilution of measured abundances from infall of pristine gas. The recently discovered XMP galaxy Leoncino, with an oxygen abundance below 3% Zsun, provides an opportunity to explore these different scenarios. Using Hubble Space Telescope imaging of the resolved stellar populations of Leoncino, we measure the distance to the galaxy to be D=12.1 (+1.7/-3.4) Mpc and find that Leoncino is located in an under-dense environment. Leoncino has a compact morphology, hosts a population of young, massive stars, has a high gas-to-star mass ratio, and shows signs of interaction with a galaxy nearby on the sky, UGC 5186. Similar to nearly all XMP galaxies known in the nearby universe, Leoncino is offset from the Luminosity-Metallicity (LZ) relation. Yet, Leoncino is consistent with the stellar Mass-Metallicity (MZ) relation defined by Local Volume galaxies. Thus, our results suggest that the offset from the LZ relation is due to higher recent star formation, likely triggered by a minor interaction, while the low oxygen abundance is consistent with the expectation that low-mass galaxies will undergo secular evolution marked by inefficient star formation and metal-loss via galactic winds. This is in contrast to XMP galaxies that are outliers in both the LZ and MZ relations; in such cases, the low oxygen abundances are best explained by dilution due to the infall of pristine gas. We also discuss why quiescent XMP galaxies are underrepresented in current surveys.

preprint2019arXiv

Comparing the Quenching Times of Faint M31 and Milky Way Satellite Galaxies

We present the star formation histories (SFHs) of 20 faint M31 satellites ($-12 \lesssim M_V \lesssim -6$) that were measured by modeling sub-horizontal branch (HB) depth color-magnitude diagrams constructed from Hubble Space Telescope (HST) imaging. Reinforcing previous results, we find that virtually all galaxies quenched between 3 and 9 Gyr ago, independent of luminosity, with a notable concentration $3-6$ Gyr ago. This is in contrast to the Milky Way (MW) satellites, which are generally either faint with ancient quenching times or luminous with recent ($<3$ Gyr) quenching times. We suggest that systematic differences in the quenching times of M31 and MW satellites may be a reflection of the varying accretion histories of M31 and the MW. This result implies that the formation histories of low-mass satellites may not be broadly representative of low-mass galaxies in general. Among the M31 satellite population we identify two distinct groups based on their SFHs: one with exponentially declining SFHs ($τ\sim 2$ Gyr) and one with rising SFHs with abrupt quenching. We speculate how these two groups could be related to scenarios for a recent major merger involving M31. The Cycle 27 HST Treasury survey of M31 satellites will provide well-constrained ancient SFHs to go along with the quenching times we measure here. The discovery and characterization of M31 satellites with $M_V \gtrsim -6$ would help quantify the relative contributions of reionization and environment to quenching of the lowest-mass satellites.

preprint2019arXiv

Galactic Winds in Low-Mass Galaxies

Stellar-feedback driven outflows are predicted to play a fundamental role in the baryon cycle of low-mass galaxies. However, observational constraints of winds in nearby dwarf galaxies are limited as outflows are transient, intrinsically low-surface brightness features, and, thus, difficult to detect. Using deep Hapha observations, we search for winds in a sample of twelve nearby dwarfs (M_* ~ 10^7 - 10^9.3 Msun) which host on-going or recent starbursts. We detect features which we classify as winds in 6 galaxies, fountain candidates in 5 galaxies, and diffuse ISM in 1 system. Winds are found preferentially in galaxies with centrally concentrated star formation, while fountains are found in galaxies with spatially distributed star formation. We suggest that the concentration of star formation is a predictor for whether a low-mass galaxy will develop a wind. The spatial extent of all detected ionized gas is limited (<1/10 virial radius) and would still be considered the ISM by cosmological simulations. Our observations suggest that the majority of material expelled from dwarfs does not escape to the intergalactic medium but remains in the halo and may be recycled to the galaxies. Derived mass-loading factors range from 0.2-7 (with only a weak dependency on circular velocity or stellar mass), in tension with higher values in simulations needed to reproduce realistic low-mass galaxies and resolve discrepancies with LambdaCDM. The sample is part of the panchromatic STARBurst IRegular Dwarf Survey - STARBIRDS - designed to characterize the starburst phenomenon in dwarf galaxies. We also report a previously uncatalogued nearby galaxy (J1118+7913).

preprint2016arXiv

ALFALFA Discovery of the Most Metal-Poor Gas-Rich Galaxy Known: AGC 198691

We present spectroscopic observations of the nearby dwarf galaxy AGC 198691. This object is part of the Survey of HI in Extremely Low-Mass Dwarfs (SHIELD) project, which is a multi-wavelength study of galaxies with HI masses in the range of 10$^{6}$-10$^{7.2}$~M$_{\odot}$ discovered by the ALFALFA survey. We have obtained spectra of the lone HII region in AGC 198691 with the new high-throughput KPNO Ohio State Multi-Object Spectrograph (KOSMOS) on the Mayall 4-m as well as with the Blue Channel spectrograph on the MMT 6.5-m telescope. These observations enable the measurement of the temperature-sensitive [OIII]$λ$4363 line and hence the determination of a "direct" oxygen abundance for AGC 198691. We find this system to be an extremely metal-deficient (XMD) system with an oxygen abundance of 12+log(O/H) = 7.02 $\pm$ 0.03, making AGC 198691 the lowest-abundance star-forming galaxy known in the local universe. Two of the five lowest-abundance galaxies known have been discovered by the ALFALFA blind HI survey; this high yield of XMD galaxies represents a paradigm shift in the search for extremely metal-poor galaxies.

preprint2016arXiv

Carbon and Oxygen Abundances in Low Metallicity Dwarf Galaxies

The study of carbon and oxygen abundances yields information on the time evolution and nucleosynthetic origins of these elements, yet remains relatively unexplored. At low metallicities (12+log(O/H) < 8.0), nebular carbon measurements are limited to rest-frame UV collisionally excited emission lines. Therefore, we present UV spectrophotometry of 12 nearby, low-metallicity, high-ionization HII regions in dwarf galaxies obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. We present the first analysis of the C/O ratio in local galaxies based solely on simultaneous significant detections of the UV O^+2 and C^+2 collisionally excited lines in seven of our targets and five objects from the literature, to create a final sample of 12 significant detections. Our sample is complemented by optical SDSS spectra, from which we measured the nebular physical conditions and oxygen abundances using the direct method. At low metallicity (12+log(O/H) < 8.0), no clear trend is evident in C/O vs. O/H for the present sample given the large dispersion observed. When combined with recombination line observations at higher values of O/H, a general trend of increasing C/O with increasing O/H is also viable, but with some significant outliers. Additionally, we find the C/N ratio appears to be constant (but with significant scatter) over a large range in oxygen abundance, indicating carbon is predominantly produced by similar nucleosynthetic mechanisms as nitrogen. If true, and our current understanding of nitrogen production is correct, this would indicate that primary production of carbon (a flat trend) dominates at low metallicity, but quasi-secondary production (an increasing trend) becomes prominent at higher metallicities. A larger sample will be needed to determine the true nature and dispersion of the relation.

preprint2016arXiv

CHAOS III: Gas-Phase Abundances in NGC5457

The CHemical Abundances of Spirals (CHAOS) project leverages the combined power of the Large Binocular Telescope with the broad spectral range and sensitivity of the Multi Object Double Spectrograph (MODS) to measure direct abundances in large samples of HII regions in spiral galaxies. We present LBT MODS observations of 109 Hii regions in NGC5457, of which 74 have robust measurements of key auroral lines, a factor of 3 larger than all previous published detections of auroral lines in the HII regions of NGC5457. Comparing the temperatures derived from the different ionic species we find: (1) strong correlations of T[NII] with T[SIII] and T[OIII], consistent with little or no intrinsic scatter; (2) a correlation of T[SIII] with T[OIII], but with significant intrinsic dispersion; (3) overall agreement between T[NII], T[SII], and T[OII], as expected, but with significant outliers; (4) the correlations of T[NII] with T[SIII] and T[OIII] match the predictions of photoionization modeling while the correlation of T[SIII] with T[OIII] is offset from the prediction of photoionization modeling. Based on these observations, which include significantly more observations of lower excitation HII regions, missing in many analyses, we inspect the commonly used ionization correction factors (ICFs) for unobserved ionic species and propose new empirical ICFs for S and Ar. We have discovered an unexpected population of HII regions with a significant offset to low values in Ne/O, which defies explanation. We derive radial gradients in O/H and N/O which agree with previous studies. Our large observational database allows us to examine the dispersion in abundances, and we find intrinsic dispersions of 0.074 in O/H and 0.095 in N/O (at a given radius). We stress that this measurement of the intrinsic dispersion comes exclusively from direct measurements of HII regions in NGC5457.

preprint2016arXiv

Chemistry and Kinematics of the Late-Forming Dwarf Irregular Galaxies Leo A, Aquarius, and Sagittarius DIG

We present Keck/DEIMOS spectroscopy of individual stars in the relatively isolated Local Group dwarf galaxies Leo A, Aquarius, and the Sagittarius dwarf irregular galaxy. The three galaxies--but especially Leo A and Aquarius--share in common delayed star formation histories relative to many other isolated dwarf galaxies. The stars in all three galaxies are supported by dispersion. We found no evidence of stellar velocity structure, even for Aquarius, which has rotating HI gas. The velocity dispersions indicate that all three galaxies are dark matter-dominated, with dark-to-baryonic mass ratios ranging from $4.4^{+1.1}_{-0.8}$ (SagDIG) to $9.6^{+2.5}_{-1.8}$ (Aquarius). Leo A and SagDIG have lower stellar metallicities than Aquarius, and they also have higher gas fractions, both of which would be expected if Aquarius were farther along in its chemical evolution. The metallicity distribution of Leo A is inconsistent with a Closed or Leaky Box model of chemical evolution, suggesting that the galaxy was pre-enriched or acquired external gas during star formation. The metallicities of stars increased steadily for all three galaxies, but possibly at different rates. The [$α$/Fe] ratios at a given [Fe/H] are lower than that of the Sculptor dwarf spheroidal galaxy, which indicates more extended star formation histories than Sculptor, consistent with photometrically derived star formation histories. Overall, the bulk kinematic and chemical properties for the late-forming dwarf galaxies do not diverge significantly from those of less delayed dwarf galaxies, including dwarf spheroidal galaxies.

preprint2016arXiv

Exploring Systematic Effects in the Relation Between Stellar Mass, Gas Phase Metallicity, and Star Formation Rate

There is evidence that the well-established mass-metallicity relation in galaxies is correlated with a third parameter: star formation rate (SFR). The strength of this correlation may be used to disentangle the relative importance of different physical processes (e.g., infall of pristine gas, metal-enriched outflows) in governing chemical evolution. However, all three parameters are susceptible to biases that might affect the observed strength of the relation between them. We analyze possible sources of systematic error, including sample bias, application of signal-to-noise ratio cuts on emission lines, choice of metallicity calibration, uncertainty in stellar mass determination, aperture effects, and dust. We present the first analysis of the relation between stellar mass, gas phase metallicity, and SFR using strong line abundance diagnostics from Dopita et al. (2013) for ~130,000 star-forming galaxies in the Sloan Digital Sky Survey and provide a detailed comparison of these diagnostics in an appendix. Using these abundance diagnostics yields a 30-55% weaker anti-correlation between metallicity and SFR at fixed stellar mass than that reported by Mannucci et al. (2010). We find that, for all abundance diagnostics, the anti-correlation with SFR is stronger for the relatively few galaxies whose current SFRs are elevated above their past average SFRs. This is also true for the new abundance diagnostic of Dopita et al. (2016), which gives anti-correlation between metallicity and SFR only in the high specific star formation rate (sSFR) regime, in contrast to the recent results of Kashino et al. (2016). The poorly constrained strength of the relation between stellar mass, metallicity, and SFR must be carefully accounted for in theoretical studies of chemical evolution.

preprint2016arXiv

Isolating the young stellar population in the outer disk of NGC 300

The recent star formation history (SFH) in the outer disk of NGC 300 is presented through the analysis of color magnitude diagrams (CMDs). We analyze resolved stellar photometry by creating CMDs from four Hubble Space Telescope fields containing a combination of images from the Advanced Camera for Surveys and the UVIS imager aboard the Wide Field Camera 3. From the best models of these CMDs, we derive the SFH in order to extract the young stellar component for the past 200 Myrs. We find that the young stellar disk of NGC 300 is unbroken out to at least ${\sim}$8 scale lengths (including an upper limit out to ${\sim}$10 scale lengths) with $r{_s}$ = 1.4 ${\pm}$ 0.1 kpc, which is similar to the total stellar surface brightness profile. This unbroken profile suggests that NGC 300 is undisturbed, similar to the isolated disk galaxy NGC 2403. We compare the environments of NGC 300, NGC 2403, and M33 along with the properties of the gas and stellar disks. We find that the disturbed HI outer disk morphology is not accompanied by a break in the young stellar disk. This may indicate that processes which affect the outer HI morphology may not leave an imprint on the young stellar disk.

preprint2016arXiv

Panchromatic Hubble Andromeda Treasury XVI. Star Cluster Formation Efficiency and the Clustered Fraction of Young Stars

We use the Panchromatic Hubble Andromeda Treasury (PHAT) survey dataset to perform spatially resolved measurements of star cluster formation efficiency ($Γ$), the fraction of stellar mass formed in long-lived star clusters. We use robust star formation history and cluster parameter constraints, obtained through color-magnitude diagram analysis of resolved stellar populations, to study Andromeda's cluster and field populations over the last $\sim$300 Myr. We measure $Γ$ of 4-8% for young, 10-100 Myr old populations in M31. We find that cluster formation efficiency varies systematically across the M31 disk, consistent with variations in mid-plane pressure. These $Γ$ measurements expand the range of well-studied galactic environments, providing precise constraints in an HI-dominated, low intensity star formation environment. Spatially resolved results from M31 are broadly consistent with previous trends observed on galaxy-integrated scales, where $Γ$ increases with increasing star formation rate surface density ($Σ_{\mathrm{SFR}}$). However, we can explain observed scatter in the relation and attain better agreement between observations and theoretical models if we account for environmental variations in gas depletion time ($τ_{\mathrm{dep}}$) when modeling $Γ$, accounting for the qualitative shift in star formation behavior when transitioning from a H$_2$-dominated to a HI-dominated interstellar medium. We also demonstrate that $Γ$ measurements in high $Σ_{\mathrm{SFR}}$ starburst systems are well-explained by $τ_{\mathrm{dep}}$-dependent fiducial $Γ$ models.

preprint2016arXiv

The Distance to M104

M104 (NGC 4594; the Sombrero galaxy) is a nearby, well-studied elliptical galaxy included in scores of surveys focused on understanding the details of galaxy evolution. Despite the importance of observations of M104, a consensus distance has not yet been established. Here, we use newly obtained Hubble Space Telescope optical imaging to measure the distance to M104 based on the tip of the red giant branch method. Our measurement yields the distance to M104 to be 9.55 +/- 0.13 +/- 0.31 Mpc equivalent to a distance modulus of 29.90 +/- 0.03 +/- 0.07 mag. Our distance is an improvement over previous results as we use a well-calibrated, stable distance indicator, precision photometry in a optimally selected field of view, and a Bayesian Maximum Likelihood technique that reduces measurement uncertainties. The most discrepant previous results are due to Tully-Fisher method distances, which are likely inappropriate for M104 given its peculiar morphology and structure. Our results are part of a larger program to measure accurate distances to a sample of well-known spiral galaxies (including M51, M74, and M63) using the tip of the red giant branch method.

preprint2016arXiv

The Distance to M51

Great investments of observing time have been dedicated to the study of nearby spiral galaxies with diverse goals ranging from understanding the star formation process to characterizing their dark matter distributions. Accurate distances are fundamental to interpreting observations of these galaxies, yet many of the best studied nearby galaxies have distances based on methods with relatively large uncertainties. We have started a program to derive accurate distances to these galaxies. Here we measure the distance to M51 - the Whirlpool galaxy - from newly obtained Hubble Space Telescope optical imaging using the tip of the red giant branch method. We measure the distance modulus to be 8.58+/-0.10 Mpc (statistical), corresponding to a distance modulus of 29.67+/-0.02 mag. Our distance is an improvement over previous results as we use a well-calibrated, stable distance indicator, precision photometry in a optimally selected field of view, and a Bayesian Maximum Likelihood technique that reduces measurement uncertainties.

preprint2016arXiv

The ISLANDS project I: Andromeda XVI, An Extremely Low Mass Galaxy not Quenched by Reionization

Based on data aquired in 13 orbits of HST time, we present a detailed evolutionary history of the M31 dSph satellite Andromeda XVI, including its life-time star formation history, the spatial distribution of its stellar populations, and the properties of its variable stars. And XVI is characterized by prolonged star formation activity from the oldest epochs until star formation was quenched ~6 Gyr ago, and, notably, only half of the mass in stars of And XVI was in place 10 Gyr ago. And XVI appears to be a low mass galaxy for which the early quenching by either reionization or starburst feedback seems highly unlikely, and thus, is most likely due to an environmental effect (e.g., an interaction), possibly connected to a late infall in the densest regions of the Local Group. Studying the star formation history as a function of galactocentric radius, we detect a mild gradient in the star formation history: the star formation activity between 6 and 8 Gyr ago is significantly stronger in the central regions than in the external regions, although the quenching age appears to be the same, within 1 Gyr. We also report the discovery of 9 RR Lyrae stars, 8 of which belong to And XVI. The RR Lyrae stars allow a new estimate of the distance, (m-M)0= 23.72+/-0.09 mag, which is marginally larger than previous estimates based on the tip of the red giant branch.

preprint2016arXiv

The Panchromatic Hubble Andromeda Treasury XV. The BEAST: Bayesian Extinction and Stellar Tool

We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys.

preprint2015arXiv

A clear age-velocity dispersion correlation in Andromeda's stellar disk

The stellar kinematics of galactic disks are key to constraining disk formation and evolution processes. In this paper, for the first time, we measure the stellar age-velocity dispersion correlation in the inner 20 kpc (3.5 disk scale lengths) of M31 and show that it is dramatically different from that in the Milky Way. We use optical Hubble Space Telescope/Advanced Camera for Surveys photometry of 5800 individual stars from the Panchromatic Hubble Andromeda Treasury (PHAT) survey and Keck/DEIMOS radial velocity measurements of the same stars from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. We show that the average line-of-sight velocity dispersion is a steadily increasing function of stellar age exterior to R=10 kpc, increasing from 30 km/s for the young upper main sequence stars to 90 km/s for the old red giant branch stars. This monotonic increase implies that a continuous or recurring process contributed to the evolution of the disk. Both the slope and normalization of the dispersion vs. age relation are significantly larger than in the Milky Way, allowing for the possibility that the disk of M31 has had a more violent history than the disk of the Milky Way, more in line with cosmological predictions. We also find evidence for an inhomogeneous distribution of stars from a second kinematical component in addition to the dominant disk component. One of the largest and hottest high-dispersion patches is present in all age bins, and may be the signature of the end of the long bar.

preprint2015arXiv

Calibrating UV Star Formation Rates for Dwarf Galaxies from STARBIRDS

Integrating our knowledge of star formation traced by observations at different wavelengths is essential for correctly interpreting and comparing star formation activity in a variety of systems and environments. This study compares extinction corrected integrated ultraviolet (UV) emission from resolved galaxies with color-magnitude diagram (CMD) based star formation rates (SFRs) derived from resolved stellar populations and CMD fitting techniques in 19 nearby starburst and post-starburst dwarf galaxies. The datasets are from the panchromatic STARBurst IRregular Dwarf Survey (STARBIRDS) and include deep legacy GALEX UV imaging, HST optical imaging, and Spitzer MIPS imaging. For the majority of the sample, the integrated near UV fluxes predicted from the CMD-based SFRs - using four different models - agree with the measured, extinction corrected, integrated near UV fluxes from GALEX images, but the far UV predicted fluxes do not. Further, we find a systematic deviation between the SFRs based on integrated far UV luminosities and existing scaling relations, and the SFRs based on the resolved stellar populations. This offset is not driven by different star formation timescales, variations in SFRs, UV attenuation, nor stochastic effects. This first comparison between CMD-based SFRs and an integrated FUV emission SFR indicator suggests that the most likely cause of the discrepancy is the theoretical FUV-SFR calibration from stellar evolutionary libraries and/or stellar atmospheric models. We present an empirical calibration of the FUV-based SFR relation for dwarf galaxies, with uncertainties, which is ~53% larger than previous relations.

preprint2015arXiv

CARMA CO Observations of Three Extremely Metal-Poor, Star-Forming Galaxies

We present sensitive CO (J = 1 - 0) emission line observations of three metal-poor dwarf irregular galaxies Leo P (Z ~ 3% Z_Solar), Sextans A (Z ~ 7.5% Z_Solar), and Sextans B (Z ~ 7.5% Z_Solar), all obtained with the Combined Array for Millimeter-wave Astronomy (CARMA) interferometer. While no CO emission was detected, the proximity of the three systems allows us to place very stringent (4 sigma) upper limits on the CO luminosity (L_CO) in these metal-poor galaxies. We find the CO luminosities to be L_CO < 2900 K km/s pc^2 for Leo P, L_CO < 12400 K km/s pc^2 for Sextans A, and L_CO < 9700 K km/s pc^2 for Sextans B. Comparison of our results with recent observational estimates of the factor for converting between L_CO and the mass of molecular hydrogen, as well as theoretical models, provides further evidence that either the CO-to-H_2 conversion factor increases sharply as metallicity decreases, or that stars are forming in these three galaxies very efficiently, requiring little molecular hydrogen.

preprint2015arXiv

CHAOS I: Direct Chemical Abundances for HII Regions in NGC 628

The CHemical Abundances of Spirals (CHAOS) project leverages the combined power of the Large Binocular Telescope (LBT) with the broad spectral range and sensitivity of the Multi Object Double Spectrograph (MODS) to measure "direct" abundances in large samples of HII regions in spiral galaxies. We present LBT MODS observations of 62 HII regions in the nearby NGC628. We measure one or more auroral lines ([OIII] 4363, [NII] 5755, [SIII] 6312, or [OII] 7320,7330) in a large number of HII regions (40). Comparing derived temperatures from multiple auroral line measurements, we find: a strong correlation between temperatures based on [SIII] and [NII]; and large discrepancies for some temperatures based on [OII] and [OIII]. These trends are consistent with other observations in the literature, yet, given the widespread use and acceptance of [OIII] as a temperature determinant, the magnitude of the T[OIII] discrepancies still came as a surprise. Based on these results, we conduct a uniform abundance analysis using the temperatures derived from [SIII] and [NII], and report the gas-phase abundance gradients for NGC628. Relative abundances of S/O, Ne/O, and Ar/O are constant across the galaxy, consistent with no systematic change in the upper IMF over the sampled range in metallicity. These alpha-element ratios, along with N/O, all show small dispersions consistent with no intrinsic dispersion (0.05) over 70% of the azimuthally averaged radius. We interpret these results as an indication that, at a given radius, the interstellar medium in NGC628 is chemically well-mixed. Unlike the nearly temperature-independent gradients, O/H abundances have a larger intrinsic dispersion of ~0.13 dex. We posit that this dispersion represents an upper limit to the true dispersion in O/H at a given radius and that some of that dispersion is due to systematic uncertainties arising from temperature measurements.

preprint2015arXiv

CHAOS II: Gas-Phase Abundances in NGC 5194

We have observed NGC5194 (M51a) as part of the CHemical Abundances of Spirals (CHAOS) project. Using the Multi Object Double Spectrographs (MODS) on the Large Binocular Telescope (LBT) we are able to measure one or more of the temperature-sensitive auroral lines ([O III] 4363, [N II] 5755, [S III] 6312) and thus measure "direct" gas-phase abundances in 29 individual HII regions. [O III] 4363 is only detected in two HII regions both of which show indications of excitation by shocks. We compare our data to previous direct abundances measured in NGC5194 and find excellent agreement for all but one region (Delta[log(O/H)] ~ 0.04). We find no evidence of trends in Ar/O, Ne/O, or S/O within NGC5194 or compared to other galaxies. We find modest negative gradients in both O/H and N/O with very little scatter (sigma < 0.08 dex), most of which can be attributed to random error and not to intrinsic dispersion. The gas-phase abundance gradient is consistent with the gradients observed in other interacting galaxies, which tend to be shallower than gradients measured in isolated galaxies. The N/O ratio (<log(N/O)> = -0.62) suggests secondary nitrogen production is responsible for a significantly larger fraction of nitrogen (e.g., factor of 8-10) relative to primary production mechanisms than predicted by theoretical models.

preprint2015arXiv

Characterizing the Star Formation of the Low-Mass SHIELD Galaxies from Hubble Space Telescope Imaging

The Survey of HI in Extremely Low-mass Dwarfs (SHIELD) is an on-going multi-wavelength program to characterize the gas, star formation, and evolution in gas-rich, very low-mass galaxies that populate the faint end of the galaxy luminosity function. The galaxies were selected from the first ~10% of the HI ALFALFA survey based on their low HI mass and low baryonic mass. Here, we measure the star-formation properties from optically resolved stellar populations for 12 galaxies using a color-magnitude diagram fitting technique. We derive lifetime average star-formation rates (SFRs), recent SFRs, stellar masses, and gas fractions. Overall, the recent SFRs are comparable to the lifetime SFRs with mean birthrate parameter of 1.4, with a surprisingly narrow standard deviation of 0.7. Two galaxies are classified as dwarf transition galaxies (dTrans). These dTrans systems have star-formation and gas properties consistent with the rest of the sample, in agreement with previous results that some dTrans galaxies may simply be low-luminosity dIrrs. We do not find a correlation between the recent star-formation activity and the distance to the nearest neighboring galaxy, suggesting that the star-formation process is not driven by gravitational interactions, but regulated internally. Further, we find a broadening in the star-formation and gas properties (i.e., specific SFRs, stellar masses, and gas fractions) compared to the generally tight correlation found in more massive galaxies. Overall, the star-formation and gas properties indicate these very low-mass galaxies host a fluctuating, non-deterministic, and inefficient star-formation process.

preprint2015arXiv

Leo P: An Unquenched Very Low-Mass Galaxy

Leo P is a low-luminosity dwarf galaxy discovered through the blind HI Arecibo Legacy Fast ALFA (ALFALFA) survey. The HI and follow-up optical observations have shown that Leo P is a gas-rich dwarf galaxy with active star formation, an underlying older population, and an extremely low oxygen abundance. We have obtained optical imaging with the Hubble Space Telescope to two magnitudes below the red clump in order to study the evolution of Leo P. We refine the distance measurement to Leo P to be 1.62+/-0.15 Mpc, based on the luminosity of the horizontal branch stars and 10 newly identified RR Lyrae candidates. This places the galaxy at the edge of the Local Group, ~0.4 Mpc from Sextans B, the nearest galaxy in the NGC 3109 association of dwarf galaxies of which Leo P is clearly a member. The star responsible for ionizing the HII region is most likely an O7V or O8V spectral type, with a stellar mass >25 Msun. The presence of this star provides observational evidence that massive stars at the upper-end of the initial mass function are capable of being formed at star formation rates as low as ~10^-5 Msun/yr. The best-fitting star formation history derived from the resolved stellar populations of Leo P using the latest PARSEC models shows a relatively constant star formation rate over the lifetime of the galaxy. The modeled luminosity characteristics of Leo P at early times are consistent with low-luminosity dSph Milky Way satellites, suggesting that Leo P is what a low-mass dSph would look like if it evolved in isolation and retained its gas. Despite the very low mass of Leo P, the imprint of reionization on its star formation history is subtle at best, and consistent with being totally negligible. The isolation of Leo P, and the total quenching of star formation of Milky Way satellites of similar mass, implies that local environment dominates the quenching of the Milky Way satellites.

preprint2015arXiv

Leo P: How Many Metals can a Very Low-Mass, Isolated Galaxy Retain?

Leo P is a gas-rich dwarf galaxy with an extremely low gas-phase oxygen abundance (3% solar). The isolated nature of Leo P enables a quantitative measurement of metals lost solely due to star formation feedback. We present an inventory of the oxygen atoms in Leo P based on the gas-phase oxygen abundance measurement, the star formation history, and the chemical enrichment evolution derived from resolved stellar populations. The star formation history also provides the total amount of oxygen produced. Overall, Leo P has retained 5 % of its oxygen; 25% of the retained oxygen is in the stars while 75% is in the gas phase. This is considerably lower than the 20-25% calculated for massive galaxies, supporting the trend for less efficient metal retention for lower mass galaxies. The retention fraction is higher than that calculated for other alpha elements (Mg, Si, Ca) in dSph Milky Way satellites of similar stellar mass and metallicity. Accounting only for the oxygen retained in stars, our results are consistent with those derived for the alpha elements in dSph galaxies. Thus, under the assumption that the dSph galaxies lost the bulk of their gas mass through an environmental process such as tidal stripping, the estimates of retained metal fractions represent underestimates by roughly a factor of four. Because of its isolation, Leo P provides an important datum for the fraction of metals lost as a function of galaxy mass due to star formation.

preprint2015arXiv

Panchromatic Hubble Andromeda Treasury XII. Mapping Stellar Metallicity Distributions in M31

We present a study of spatial variations in the metallicity of old red giant branch stars in the Andromeda galaxy. Photometric metallicity estimates are derived by interpolating isochrones for over seven million stars in the Panchromatic Hubble Andromeda Treasury (PHAT) survey. This is the first systematic study of stellar metallicities over the inner 20 kpc of Andromeda's galactic disk. We see a clear metallicity gradient of $-0.020\pm0.004$ dex/kpc from $\sim4-20$ kpc assuming a constant RGB age. This metallicity gradient is derived after correcting for the effects of photometric bias and completeness and dust extinction and is quite insensitive to these effects. The unknown age gradient in M31's disk creates the dominant systematic uncertainty in our derived metallicity gradient. However, spectroscopic analyses of galaxies similar to M31 show that they typically have small age gradients that make this systematic error comparable to the 1$σ$ error on our metallicity gradient measurement. In addition to the metallicity gradient, we observe an asymmetric local enhancement in metallicity at radii of 3-6 kpc that appears to be associated with Andromeda's elongated bar. This same region also appears to have an enhanced stellar density and velocity dispersion.

preprint2015arXiv

PHAT Stellar Cluster Survey. II. Andromeda Project Cluster Catalog

We construct a stellar cluster catalog for the Panchromatic Hubble Andromeda Treasury (PHAT) survey using image classifications collected from the Andromeda Project citizen science website. We identify 2,753 clusters and 2,270 background galaxies within ~0.5 deg$^2$ of PHAT imaging searched, or ~400 kpc$^2$ in deprojected area at the distance of the Andromeda galaxy (M31). These identifications result from 1.82 million classifications of ~20,000 individual images (totaling ~7 gigapixels) by tens of thousands of volunteers. We show that our crowd-sourced approach, which collects >80 classifications per image, provides a robust, repeatable method of cluster identification. The high spatial resolution Hubble Space Telescope images resolve individual stars in each cluster and are instrumental in the factor of ~6 increase in the number of clusters known within the survey footprint. We measure integrated photometry in six filter passbands, ranging from the near-UV to the near-IR. PHAT clusters span a range of ~8 magnitudes in F475W (g-band) luminosity, equivalent to ~4 decades in cluster mass. We perform catalog completeness analysis using >3000 synthetic cluster simulations to determine robust detection limits and demonstrate that the catalog is 50% complete down to ~500 solar masses for ages <100 Myr. We include catalogs of clusters, background galaxies, remaining unselected candidates, and synthetic cluster simulations, making all information publicly available to the community. The catalog published here serves as the definitive base data product for PHAT cluster science, providing a census of star clusters in an L$^*$ spiral galaxy with unmatched sensitivity and quality.

preprint2015arXiv

Testing Density Wave Theory with Resolved Stellar Populations around Spiral Arms in M81

Stationary density waves rotating at a constant pattern speed $Ω_{\rm P}$ would produce age gradients across spiral arms. We test whether such age gradients are present in M81 by deriving the recent star formation histories (SFHs) of 20 regions around one of M81's grand-design spiral arms. For each region, we use resolved stellar populations to determine the SFH by modeling the observed color-magnitude diagram (CMD) constructed from archival Hubble Space Telescope (HST) F435W and F606W imaging. Although we should be able to detect systematic time delays in our spatially-resolved SFHs, we find no evidence of star formation propagation across the spiral arm. Our data therefore provide no convincing evidence for a stationary density wave with a single pattern speed in M81, and instead favor the scenario of kinematic spiral patterns that are likely driven by tidal interactions with the companion galaxies M82 and NGC 3077.

preprint2015arXiv

The effects of He I 10830 on helium abundance determinations

Observations of helium and hydrogen emission lines from metal-poor extragalactic H II regions provide an independent method for determining the primordial helium abundance, Y_p. Traditionally, the emission lines employed are in the visible wavelength range, and the number of suitable lines is limited. Furthermore, when using these lines, large systematic uncertainties in helium abundance determinations arise due to the degeneracy of physical parameters, such as temperature and density. Recently, Izotov, Thuan, & Guseva (2014) have pioneered adding the He 10830 infrared emission line in helium abundance determinations. The strong electron density dependence of He 10830 makes it ideal for better constraining density, potentially breaking the degeneracy with temperature. We revisit our analysis of the dataset published by Izotov, Thuan, & Stasinska (2007) and incorporate the newly available observations of He 10830 by scaling them using the observed-to-theoretical Paschen-gamma ratio. The solutions are better constrained, in particular for electron density, temperature, and the neutral hydrogen fraction, improving the model fit to data, with the result that more spectra now pass screening for quality and reliability, in addition to a standard 95% confidence level cut. Furthermore, the addition of He 10830 decreases the uncertainty on the helium abundance for all galaxies, with reductions in the uncertainty ranging from 10-80%. Overall, we find a reduction in the uncertainty on Y_p by over 50%. From a regression to zero metallicity, we determine Y_p = 0.2449 +/- 0.0040, consistent with the BBN result, Y_p = 0.2470 +/- 0.0002, based on the Planck determination of the baryon density. The dramatic improvement in the uncertainty from incorporating He 10830 strongly supports the case for simultaneous (thus not requiring scaling) observations of visible and infrared helium emission line spectra.

preprint2015arXiv

The High-Mass Stellar Initial Mass Function in M31 Clusters

We have undertaken the largest systematic study of the high-mass stellar initial mass function (IMF) to date using the optical color-magnitude diagrams (CMDs) of 85 resolved, young (4 Myr < t < 25 Myr), intermediate mass star clusters (10^3-10^4 Msun), observed as part of the Panchromatic Hubble Andromeda Treasury (PHAT) program. We fit each cluster's CMD to measure its mass function (MF) slope for stars >2 Msun. For the ensemble of clusters, the distribution of stellar MF slopes is best described by $Γ=+1.45^{+0.03}_{-0.06}$ with a very small intrinsic scatter. The data also imply no significant dependencies of the MF slope on cluster age, mass, and size, providing direct observational evidence that the measured MF represents the IMF. This analysis implies that the high-mass IMF slope in M31 clusters is universal with a slope ($Γ=+1.45^{+0.03}_{-0.06}$) that is steeper than the canonical Kroupa (+1.30) and Salpeter (+1.35) values. Using our inference model on select Milky Way (MW) and LMC high-mass IMF studies from the literature, we find $Γ_{\rm MW} \sim+1.15\pm0.1$ and $Γ_{\rm LMC} \sim+1.3\pm0.1$, both with intrinsic scatter of ~0.3-0.4 dex. Thus, while the high-mass IMF in the Local Group may be universal, systematics in literature IMF studies preclude any definitive conclusions; homogenous investigations of the high-mass IMF in the local universe are needed to overcome this limitation. Consequently, the present study represents the most robust measurement of the high-mass IMF slope to date. We have grafted the M31 high-mass IMF slope onto widely used sub-solar mass Kroupa and Chabrier IMFs and show that commonly used UV- and Halpha-based star formation rates should be increased by a factor of ~1.3-1.5 and the number of stars with masses >8 Msun are ~25% fewer than expected for a Salpeter/Kroupa IMF. [abridged]

preprint2015arXiv

The link between mass distribution and starbursts in dwarf galaxies

Recent studies have shown that starburst dwarf galaxies have steeply rising rotation curves in their inner parts, pointing to a close link between the intense star formation and a centrally concentrated mass distribution (baryons and dark matter). More quiescent dwarf irregulars typically have slowly rising rotation curves, although some "compact" irregulars with steep, inner rotation curves exist. We analyze archival Hubble Space Telescope images of two nearby "compact" irregular galaxies (NGC 4190 and NGC 5204), which were selected solely on the basis of their dynamical properties and their proximity. We derive their recent star-formation histories by fitting color-magnitude diagrams of resolved stellar populations, and find that the star-formation properties of both galaxies are consistent with those of known starburst dwarfs. Despite the small sample, this strongly reinforces the notion that the starburst activity is closely related to the inner shape of the potential well.

preprint2015arXiv

The Panchromatic Hubble Andromeda Treasury VIII: A Wide-Area, High-Resolution Map of Dust Extinction in M31

We map the distribution of dust in M31 at 25pc resolution, using stellar photometry from the Panchromatic Hubble Andromeda Treasury. We develop a new mapping technique that models the NIR color-magnitude diagram (CMD) of red giant branch (RGB) stars. The model CMDs combine an unreddened foreground of RGB stars with a reddened background population viewed through a log-normal column density distribution of dust. Fits to the model constrain the median extinction, the width of the extinction distribution, and the fraction of reddened stars. The resulting extinction map has >4 times better resolution than maps of dust emission, while providing a more direct measurement of the dust column. There is superb morphological agreement between the new map and maps of the extinction inferred from dust emission by Draine et al. 2014. However, the widely-used Draine & Li (2007) dust models overpredict the observed extinction by a factor of ~2.5, suggesting that M31's true dust mass is lower and that dust grains are significantly more emissive than assumed in Draine et al. (2014). The discrepancy we identify is consistent with similar findings in the Milky Way by the Planck Collaboration (2015), but has a more complex dependence on parameters from the Draine & Li (2007) dust models. We also show that the discrepancy with the Draine et al. (2014) map is lowest where the interstellar radiation field has a harder spectrum than average. We discuss possible improvements to the CMD dust mapping technique, and explore further applications.

preprint2015arXiv

The Panchromatic Hubble Andromeda Treasury XI: The Spatially-Resolved Recent Star Formation History of M31

We measure the recent star formation history (SFH) across M31 using optical images taken with the \texit{Hubble Space Telescope} as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit the color-magnitude diagrams in ~9000 regions that are ~100 pc $\times$ 100 pc in projected size, covering a 0.5 square degree area (~380 kpc$^2$, deprojected) in the NE quadrant of M31. We show that the SFHs vary significantly on these small spatial scales but that there are also coherent galaxy-wide fluctuations in the SFH back to ~500 Myr, most notably in M31's 10-kpc star-forming ring. We find that the 10-kpc ring is at least 400 Myr old, showing ongoing star formation over the past ~500 Myr. This indicates the presence of molecular gas in the ring over at least 2 dynamical times at this radius. We also find that the ring's position is constant throughout this time, and is stationary at the level of 1 km/s, although there is evidence for broadening of the ring due to diffusion of stars into the disk. Based on existing models of M31's ring features, the lack of evolution in the ring's position makes a purely collisional ring origin highly unlikely. We find that the global SFR has been fairly constant over the last ~500 Myr, though it does show a small increase at 50 Myr that is 1.3 times the average SFR over the past 100 Myr. During the last ~500 Myr, ~60% of all SF occurs in the 10-kpc ring. Finally, we find that in the past 100 Myr, the average SFR over the PHAT survey area is $0.28\pm0.03$ M$_\odot$ yr$^{-1}$ with an average deprojected intensity of $7.3 \times 10^{-4}$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$, which yields a total SFR of ~0.7 M$_\odot$ yr$^{-1}$ when extrapolated to the entire area of M31's disk. This SFR is consistent with measurements from broadband estimates. [abridged]

preprint2015arXiv

The Star Formation Histories of Local Group Dwarf Galaxies III. Characterizing Quenching in Low-Mass Galaxies

We explore the quenching of low-mass galaxies (10^4 < Mstar < 10^8 Msun) as a function of lookback time using the star formation histories (SFHs) of 38 Local Group dwarf galaxies. The SFHs were derived from analyzing color-magnitude diagrams of resolved stellar populations in archival Hubble Space Telescope/Wide Field Planetary Camera 2 imaging. We find: (1) Lower mass galaxies quench earlier than higher mass galaxies; (2) Inside of virial radius there is no correlation between a satellite's current proximity to a massive host and its quenching epoch; (3) There are hints of systematic differences in quenching times of M31 and Milky Way (MW) satellites, although the sample sample size and uncertainties in the SFHs of M31 dwarfs prohibit definitive conclusions. Combined with literature results, we qualitatively consider the redshift evolution (z=0-1) of the quenched galaxy fraction over ~7 dex in stellar mass (10^4 < Mstar < 10^11.5 Msun). The quenched fraction of all galaxies generally increases toward the present, with both the lowest and highest mass systems exhibiting the largest quenched fractions at all redshifts. In contrast, galaxies between Mstar ~ 10^8-10^10 Msun have the lowest quenched fractions. We suggest that such intermediate-mass galaxies are the least efficient at quenching. Finally, we compare our quenching times with predictions for infall times of low-mass galaxies associated with the MW. We find that some of the lowest-mass satellites (e.g., CVn II, Leo IV) may have been quenched before infall while higher mass satellites (e.g., Leo I, Fornax) typically quench ~1-4 Gyr after infall.

preprint2014arXiv

ALFALFA Discovery of the Nearby Gas-rich Dwarf Galaxy Leo P. V. Neutral Gas Dynamics and Kinematics

We present new HI spectral line imaging of the extremely metal-poor, star-forming dwarf irregular galaxy Leo P. Our HI images probe the global neutral gas properties and the local conditions of the interstellar medium (ISM). The HI morphology is slightly elongated along the optical major-axis. We do not find obvious signatures of interaction or infalling gas at large spatial scales. The neutral gas disk shows obvious rotation, although the velocity dispersion is comparable to the rotation velocity. The rotation amplitude is estimated to be V_c = 15 +/- 5 km/s. Within the HI radius probed by these observations, the mass ratio of gas to stars is roughly 2:1, while the ratio of the total mass to the baryonic mass is ~15:1. We use this information to place Leo P on the baryonic Tully-Fisher relation, testing the baryonic content of cosmic structures in a sparsely populated portion of parameter space that has hitherto been occupied primarily by dwarf spheroidal galaxies. We detect the signature of two temperature components in the neutral ISM of Leo P; the cold and warm components have characteristic velocity widths of 4.2 +/- 0.9 km/s and 10.1 +/- 1.2 km/s, corresponding to kinetic temperatures of ~1100 K and ~6200 K, respectively. The cold HI component is unresolved at a physical resolution of 200 pc. The highest HI surface densities are observed in close physical proximity to the single HII region. A comparison of the neutral gas properties of Leo P with other extremely metal-deficient (XMD) galaxies reveals that Leo P has the lowest neutral gas mass of any known XMD, and that the dynamical mass of Leo P is more than two orders of magnitude smaller than any known XMD with comparable metallicity.

preprint2014arXiv

Comparing M31 and Milky Way Satellites: The Extended Star Formation Histories of Andromeda II and Andromeda XVI

We present the first comparison between the lifetime star formation histories (SFHs) of M31 and Milky Way (MW) satellites. Using the Advanced Camera for Surveys aboard the Hubble Space Telescope, we obtained deep optical imaging of Andromeda II (M$_{V} = -$12.0; log(M$_{\star}$/M$_{\odot}$) $\sim$ 6.7) and Andromeda XVI (M$_{V} = -$7.5; log(M$_{\star}$/M$_{\odot}$) $\sim$ 4.9) yielding color-magnitude diagrams that extend at least 1 magnitude below the oldest main sequence turnoff, and are similar in quality to those available for the MW companions. And II and And XVI show strikingly similar SFHs: both formed 50-70% of their total stellar mass between 12.5 and 5 Gyr ago (z$\sim$5-0.5) and both were abruptly quenched $\sim$ 5 Gyr ago (z$\sim$0.5). The predominance of intermediate age populations in And XVI makes it qualitatively different from faint companions of the MW and clearly not a pre-reionization fossil. Neither And II nor And XVI appears to have a clear analog among MW companions, and the degree of similarity in the SFHs of And II and And XVI is not seen among comparably faint-luminous pairs of MW satellites. These findings provide hints that satellite galaxy evolution may vary substantially among hosts of similar stellar mass. Although comparably deep observations of more M31 satellites are needed to further explore this hypothesis, our results underline the need for caution when interpreting satellite galaxies of an individual system in a broader cosmological context.

preprint2014arXiv

Delayed Star Formation in Isolated Dwarf Galaxies: HST Star Formation History of the Aquarius Dwarf Irregular

We have obtained deep images of the highly isolated (d = 1 Mpc) Aquarius dwarf irregular galaxy (DDO 210) with the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). The resulting color-magnitude diagram (CMD) reaches more than a magnitude below the oldest main-sequence turnoff, allowing us to derive the star formation history (SFH) over the entire lifetime of the galaxy with a timing precision of ~10% of the lookback time. Using a maximum likelihood fit to the CMD we find that only ~10% of all star formation in Aquarius took place more than 10 Gyr ago (lookback time equivalent to redshift z ~2). The star formation rate increased dramatically ~6-8 Gyr ago (z ~ 0.7-1.1) and then declined until the present time. The only known galaxy with a more extreme confirmed delay in star formation is Leo A, a galaxy of similar M(HI)/M(stellar), dynamical mass, mean metallicity, and degree of isolation. The delayed stellar mass growth in these galaxies does not track the mean dark matter accretion rate from CDM simulations. The similarities between Leo A and Aquarius suggest that if gas is not removed from dwarf galaxies by interactions or feedback, it can linger for several gigayears without cooling in sufficient quantity to form stars efficiently. We discuss possible causes for the delay in star formation including suppression by reionization and late-time mergers. We find reasonable agreement between our measured SFHs and select cosmological simulations of isolated dwarfs. Because star formation and merger processes are both stochastic in nature, delayed star formation in various degees is predicted to be a characteristic (but not a universal) feature of isolated small galaxies.

preprint2014arXiv

Discovery of a Gas-Rich Companion to the Extremely Metal-Poor Galaxy DDO 68

We present HI spectral-line imaging of the extremely metal-poor galaxy DDO 68. This system has a nebular oxygen abundance of only 3% Z$_{\odot}$, making it one of the most metal-deficient galaxies known in the local volume. Surprisingly, DDO 68 is a relatively massive and luminous galaxy for its metal content, making it a significant outlier in the mass-metallicity and luminosity-metallicity relationships. The origin of such a low oxygen abundance in DDO 68 presents a challenge for models of the chemical evolution of galaxies. One possible solution to this problem is the infall of pristine neutral gas, potentially initiated during a gravitational interaction. Using archival HI spectral-line imaging obtained with the Karl G. Jansky Very Large Array, we have discovered a previously unknown companion of DDO 68. This low-mass (M$_{\rm HI}$ $=$ 2.8$\times$10$^{7}$ M$_{\odot}$), recently star-forming (SFR$_{\rm FUV}$ $=$ 1.4$\times$10$^{-3}$ M$_{\odot}$ yr$^{-1}$, SFR$_{\rm Hα}$ $<$ 7$\times$10$^{-5}$ M$_{\odot}$ yr$^{-1}$) companion has the same systemic velocity as DDO 68 (V$_{\rm sys}$ $=$ 506 km s$^{-1}$; D $=$ 12.74$\pm$0.27 Mpc) and is located at a projected distance of 42 kpc. New HI maps obtained with the 100m Robert C. Byrd Green Bank Telescope provide evidence that DDO 68 and this companion are gravitationally interacting at the present time. Low surface brightness HI gas forms a bridge between these objects.

preprint2014arXiv

Distance Determinations to SHIELD Galaxies from HST Imaging

The Survey of HI in Extremely Low-mass Dwarf galaxies (SHIELD) is an on-going multi-wavelength program to characterize the gas, star formation, and evolution in gas-rich, very low-mass galaxies. The galaxies were selected from the first ~10% of the HI ALFALFA survey based on their inferred low HI mass and low baryonic mass, and all systems have recent star formation. Thus, the SHIELD sample probes the faint end of the galaxy luminosity function for star-forming galaxies. Here, we measure the distances to the 12 SHIELD galaxies to be between 5-12 Mpc by applying the tip of the red giant method to the resolved stellar populations imaged by the Hubble Space Telescope. Based on these distances, the HI masses in the sample range from $4\times10^6$ to $6\times10^7$ M$_{\odot}$, with a median HI mass of $1\times10^7$ M$_{\odot}$. The TRGB distances are up to 73% farther than flow-model estimates in the ALFALFA catalog. Because of the relatively large uncertainties of flow model distances, we are biased towards selecting galaxies from the ALFALFA catalog where the flow model underestimates the true distances. The measured distances allow for an assessment of the native environments around the sample members. Five of the galaxies are part of the NGC 672 and NGC 784 groups, which together constitute a single structure. One galaxy is part of a larger linear ensemble of 9 systems that stretches 1.6 Mpc from end to end. Two galaxies reside in regions with 1-4 neighbors, and four galaxies are truly isolated with no known system identified within a radius of 1 Mpc.

preprint2014arXiv

The ACS LCID project. X. The Star Formation History of IC 1613: Revisiting the Over-Cooling Problem

We present an analysis of the star formation history (SFH) of a field near the half light radius in the Local Group dwarf irregular galaxy IC 1613 based on deep Hubble Space Telescope Advanced Camera for Surveys imaging. Our observations reach the oldest main sequence turn-off, allowing a time resolution at the oldest ages of ~1 Gyr. Our analysis shows that the SFH of the observed field in IC 1613 is consistent with being constant over the entire lifetime of the galaxy. These observations rule out an early dominant episode of star formation in IC 1613. We compare the SFH of IC 1613 with expectations from cosmological models. Since most of the mass is in place at early times for low mass halos, a naive expectation is that most of the star formation should have taken place at early times. Models in which star formation follows mass accretion result in too many stars formed early and gas mass fractions which are too low today (the "over-cooling problem"). The depth of the present photometry of IC 1613 shows that, at a resolution of ~1 Gyr, the star formation rate is consistent with being constant, at even the earliest times, which is difficult to achieve in models where star formation follows mass assembly.

preprint2014arXiv

The Panchromatic Hubble Andromeda Treasury V: Ages and Masses of the Year 1 Stellar Clusters

We present ages and masses for 601 star clusters in M31 from the analysis of the six filter integrated light measurements from near ultraviolet to near infrared wavelengths, made as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We derive the ages and masses using a probabilistic technique, which accounts for the effects of stochastic sampling of the stellar initial mass function. Tests on synthetic data show that this method, in conjunction with the exquisite sensitivity of the PHAT observations and their broad wavelength baseline, provides robust age and mass recovery for clusters ranging from $\sim 10^2 - 2 \times 10^6 M_\odot$. We find that the cluster age distribution is consistent with being uniform over the past $100$ Myr, which suggests a weak effect of cluster disruption within M31. The age distribution of older ($>100$ Myr) clusters fall towards old ages, consistent with a power-law decline of index $-1$, likely from a combination of fading and disruption of the clusters. We find that the mass distribution of the whole sample can be well-described by a single power-law with a spectral index of $-1.9 \pm 0.1$ over the range of $10^3-3 \times 10^5 M_\odot$. However, if we subdivide the sample by galactocentric radius, we find that the age distributions remain unchanged. However, the mass spectral index varies significantly, showing best fit values between $-2.2$ and $-1.8$, with the shallower slope in the highest star formation intensity regions. We explore the robustness of our study to potential systematics and conclude that the cluster mass function may vary with respect to environment.

preprint2014arXiv

The Panchromatic Hubble Andromeda Treasury. VI. The reliability of far-ultraviolet flux as a star formation tracer on sub-kpc scales

We have used optical observations of resolved stars from the Panchromatic Hubble Andromeda Treasury (PHAT) to measure the recent (< 500 Myr) star formation histories (SFHs) of 33 FUV-bright regions in M31. The region areas ranged from ~$10^4$ to $10^6$ pc$^2$, which allowed us to test the reliability of FUV flux as a tracer of recent star formation on sub-kpc scales. The star formation rates (SFRs) derived from the extinction-corrected observed FUV fluxes were, on average, consistent with the 100-Myr mean SFRs of the SFHs to within the 1$σ$ scatter. Overall, the scatter was larger than the uncertainties in the SFRs and particularly evident among the smallest regions. The scatter was consistent with an even combination of discrete sampling of the initial mass function and high variability in the SFHs. This result demonstrates the importance of satisfying both the full-IMF and the constant-SFR assumptions for obtaining precise SFR estimates from FUV flux. Assuming a robust FUV extinction correction, we estimate that a factor of 2.5 uncertainty can be expected in FUV-based SFRs for regions smaller than $10^5$ pc$^2$, or a few hundred pc. We also examined ages and masses derived from UV flux under the common assumption that the regions are simple stellar populations (SSPs). The SFHs showed that most of the regions are not SSPs, and the age and mass estimates were correspondingly discrepant from the SFHs. For those regions with SSP-like SFHs, we found mean discrepancies of 10 Myr in age and a factor of 3 to 4 in mass. It was not possible to distinguish the SSP-like regions from the others based on integrated FUV flux.

preprint2014arXiv

The Star Formation Histories of Local Group Dwarf Galaxies I. Hubble Space Telescope / Wide Field Planetary Camera 2 Observations

We present uniformly measured star formation histories (SFHs) of 40 Local Group dwarf galaxies based on color-magnitude diagram (CMD) analysis from archival Hubble Space Telescope imaging. We demonstrate that accurate SFHs can be recovered from CMDs that do not reach the oldest main sequence turn-off (MSTO), but emphasize that the oldest MSTO is critical for precisely constraining the earliest epochs of star formation. We find that: (1) the average lifetime SFHs of dwarf spheroidals (dSphs) can be approximated by an exponentially declining SFH with $τ$ $\sim$ 5 Gyr; (2) lower luminosity dSphs are less likely to have extended SFHs than more luminous dSphs; (3) the average SFHs of dwarf irregulars (dIrrs), transition dwarfs (dTrans), and dwarf ellipticals (dEs) can be approximated by the combination of an exponentially declining SFH ($τ$ $\sim$ 3-4 Gyr) for lookback ages $>$ 10-12 Gyr ago and a constant SFH thereafter; (4) the observed fraction of stellar mass formed prior to z=2 ranges considerably (80\% for galaxies with M $<$ 10$^5$ M$_{\odot}$ to 30\% for galaxies with M$>$10$^7$ M$_{\odot}$) and is largely explained by environment; (5) the distinction between "ultra-faint" and "classical" dSphs is arbitrary; (6) LG dIrrs formed a significantly higher fraction of stellar mass prior to z=2 than the SDSS galaxies from Leiter 2012 and the SFHs from the abundance matching models of Behroozi et al. 2013. This may indicate higher than expected star-formation efficiencies at early times in low mass galaxies. Finally, we provide all the SFHs in tabulated electronic format for use by the community.

preprint2014arXiv

The Star Formation Histories of Local Group Dwarf Galaxies II. Searching For Signatures of Reionization

We search for signatures of reionization in the star formation histories (SFHs) of 38 Local Group dwarf galaxies (10$^4$ $<$ M$_{\star}$ $<$ 10$^9$ M$_{\odot}$). The SFHs are derived from color-magnitude diagrams using archival Hubble Space Telescope/Wide Field Planetary Camera 2 imaging. Only five quenched galaxies (And V, And VI, And XIII, Leo IV, Hercules) are consistent with forming the bulk of their stars before reionization, when full uncertainties are considered. Observations of 13 of the predicted `true fossils' identified by Bovill & Ricotti show that only two (Hercules and Leo IV) indicate star formation quenched by reionization. However, both are within the virial radius of the Milky Way and evidence of tidal disturbance complicates this interpretation. We argue that the late-time gas capture scenario posited by Ricotti for the low mass, gas-rich, and star-forming fossil candidate Leo T is observationally indistinguishable from simple gas retention. Given the ambiguity between environmental effects and reionization, the best reionization fossil candidates are quenched low mass field galaxies (e.g., KKR 25).

preprint2013arXiv

A New Approach to Detailed Structural Decomposition from the SPLASH and PHAT Surveys: Kicked-up Disk Stars in the Andromeda Galaxy?

We characterize the bulge, disk, and halo subcomponents in the Andromeda galaxy (M31) over the radial range 4 < R_proj < 225 kpc. The cospatial nature of these subcomponents renders them difficult to disentangle using surface brightness (SB) information alone, especially interior to ~20 kpc. Our new decomposition technique combines information from the luminosity function (LF) of over 1.5 million bright (20 < m_814W < 22) stars from the Panchromatic Hubble Andromeda Treasury (PHAT) survey, radial velocities of over 5000 red giant branch stars in the same magnitude range from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey, and integrated I-band SB profiles from various sources. We use an affine-invariant Markov chain Monte Carlo algorithm to fit an appropriate toy model to these three data sets. The bulge, disk, and halo SB profiles are modeled as a Sersic, exponential, and cored power-law, respectively, and the LFs are modeled as broken power-laws. We present probability distributions for each of 32 parameters describing the SB profiles and LFs of the three subcomponents. We find that the number of stars with a disk-like LF is ~5% larger than the the number with disk-like (dynamically cold) kinematics, suggesting that some stars born in the disk have been dynamically heated to the point that they are kinematically indistinguishable from halo members. This is the first kinematical evidence for a "kicked-up disk" halo population in M31. The fraction of kicked-up disk stars is consistent with that found in simulations. We also find evidence for a radially varying disk LF, consistent with a negative metallicity gradient in the stellar disk.

preprint2013arXiv

ALFALFA Discovery of the Nearby Gas-Rich Dwarf Galaxy Leo P. I. HI Observations

The discovery of a previously unknown 21cm HI line source identified as an ultra-compact high velocity cloud in the ALFALFA survey is reported. The HI detection is barely resolved by the Arecibo 305m telescope ~4' beam and has a narrow HI linewidth (HPFW of 24 km/s). Further HI observations at Arecibo and with the VLA corroborate the ALFALFA HI detection, provide an estimate of the HI radius, ~1' at the 5 x 10^19 cm^-2 isophote, and show the cloud to exhibit rotation with an amplitude of ~9.0 +/- 1.5 km/s. In other papers, Rhode et al. (2013) show the HI source to have a resolved stellar counterpart and ongoing star forming activity, while Skillman et al. (2013) reveal it as having extremely low metallicity: 12 + log(O/H) = 7.16 +/- 0.04. The HI mass to stellar mass ratio of the object is found to be 2.6. We use the Tully-Fisher template relation in its baryonic form (McGaugh 2012) to obtain a distance estimate D = 1.3 (+0.9,-0.5) Mpc. Additional constraints on the distance are also provided by the optical data of Rhode et al. (2013) and McQuinn et al. (private communication), both indicating a distance in the range of 1.5 to 2.0 Mpc. The three estimates are compatible within their errors. The object appears to be located beyond the dynamical boundaries of, but still in close proximity to the Local Group. Its pristine properties are consistent with the sedate environment of its location. At a nominal distance of 1.75 Mpc, it would have an HI mass of ~1.0 x 10^6 Msun, a stellar mass of ~3.6 x 10^5 Msun, and a dynamical mass within the HI radius of ~1.5 x 10^7 Msun. This discovery supports the idea that optically faint - or altogether dark - low mass halos may be detectable through their non-stellar baryons.

preprint2013arXiv

ALFALFA Discovery of the Nearby Gas-Rich Dwarf Galaxy Leo P. IV. Distance Measurement from LBT Optical Imaging

Leo P is a low-luminosity dwarf galaxy discovered through the blind HI Arecibo Legacy Fast ALFA (ALFALFA) survey. The HI and follow-up optical observations have shown that Leo P is a gas-rich dwarf galaxy with both active star formation and an underlying older population, as well as an extremely low oxygen abundance. Here, we measure the distance to Leo P by applying the tip of the red giant branch (TRGB) distance method to photometry of the resolved stellar population from new Large Binocular Telescope (LBT) V and I band imaging. We measure a distance modulus of 26.19 (+0.17/-0.50) mag corresponding to a distance of 1.72 (+0.14/-0.40) Mpc. Although our photometry reaches 3 magnitudes below the TRGB, the sparseness of the red giant branch (RGB) yields higher uncertainties on the lower limit of the distance. Leo P is outside the Local Group with a distance and velocity consistent with the local Hubble flow. While located in a very low-density environment, Leo P lies within ~0.5 Mpc of a loose association of dwarf galaxies which include NGC 3109, Antlia, Sextans A, and Sextans B, and 1.1 Mpc away from its next nearest neighbor, Leo A. Leo P is one of the lowest metallicity star-forming galaxies known in the nearby universe, comparable in metallicity to I Zw 18 and DDO 68, but with stellar characteristics similar to dwarf spheriodals (dSphs) in the Local Volume such as Carina, Sextans, and Leo II. Given its physical properties and isolation, Leo P may provide an evolutionary link between gas-rich dwarf irregular galaxies and dSphs that have fallen into a Local Group environment and been stripped of their gas.

preprint2013arXiv

ALFALFA Discovery of the Nearby Gas-Rich Dwarf Galaxy Leo~P. III. An Extremely Metal Deficient Galaxy

We present KPNO 4-m and LBT/MODS spectroscopic observations of an HII region in the nearby dwarf irregular galaxy Leo P discovered recently in the Arecibo ALFALFA survey. In both observations, we are able to accurately measure the temperature sensitive [O III] 4363 Angstrom line and determine a "direct" oxygen abundance of 12 + log(O/H) = 7.17 +/- 0.04. Thus, Leo P is an extremely metal deficient (XMD) galaxy, and, indeed, one of the most metal deficient star-forming galaxies ever observed. For its estimated luminosity, Leo P is consistent with the relationship between luminosity and oxygen abundance seen in nearby dwarf galaxies. Leo P shows normal alpha element abundance ratios (Ne/O, S/O, and Ar/O) when compared to other XMD galaxies, but elevated N/O, consistent with the "delayed release" hypothesis for N/O abundances. We derive a helium mass fraction of 0.2509 +0.0184 -0.0123 which compares well with the WMAP + BBN prediction of 0.2483 +/- 0.0002 for the primordial helium abundance. We suggest that surveys of very low mass galaxies compete well with emission line galaxy surveys for finding XMD galaxies. It is possible that XMD galaxies may be divided into two classes: the relatively rare XMD emission line galaxies which are associated with starbursts triggered by infall of low-metallicity gas and the more common, relatively quiescent XMD galaxies like Leo P, with very low chemical abundances due to their intrinsically small masses.

preprint2013arXiv

Comparing the Ancient Star Formation Histories of the Magellanic Clouds

We present preliminary results from a new HST archival program aimed at tightly constraining the ancient (>4 Gyr ago) star formation histories (SFHs) of the field populations of the SMC and LMC. We demonstrate the quality of the archival data by constructing HST/WFPC2-based color-magnitude diagrams (CMDs; M_{F555W} ~ +8) for 7 spatially diverse fields in the SMC and 8 fields in the LMC. The HST-based CMDs are >2 magnitudes deeper than any from ground based observations, and are particularly superior in high surface brightness regions, e.g., the LMC bar, which contain a significant fraction of star formation and are crowding limited from ground based observations. To minimize systematic uncertainties, we derive the SFH of each field using an identical maximum likelihood CMD fitting technique. We then compute an approximate mass weighted average SFH for each galaxy. We find that both galaxies lack a dominant burst of early star formation, which suggests either a suppression or an under-fueling of early star formation. From 10-12 Gyr ago, the LMC experienced a period of enhanced stellar mass growth relative to the SMC. Similar to some previous studies, we find two notable peaks in the SFH of the SMC at ~4.5 and 9 Gyr ago, which could be due to repeated close passages with the LMC, implying an interaction history that has persisted for at least 9 Gyr. We find little evidence for strong periodic behavior in the lifetime SFHs of both MCs, suggesting that repeated encounters with the Milky Way are unlikely. Beginning ~3.5 Gyr ago, both galaxies show increases in their SFHs, in agreement with previous studies, and thereafter, track each other remarkably well. (abridged)

preprint2013arXiv

Drivers of HI Turbulence in Dwarf Galaxies

Neutral hydrogen (HI) velocity dispersions are believed to be set by turbulence in the interstellar medium (ISM). Although turbulence is widely believed to be driven by star formation (SF), recent studies have shown that this driving mechanism may not be dominant in regions of low SF rate surface density (SFRSD), such as found in dwarf galaxies or the outer regions of spirals. We have generated average HI line profiles in a number of nearby dwarfs and low-mass spirals by co-adding HI spectra in regions with either a common radius or SFRSD. We find that the spatially-resolved superprofiles are composed of a central narrow peak (5-15 km/s) with higher velocity wings to either side. With the assumption that the central peak reflects the turbulent velocity dispersion, we compare HI kinematics to local ISM properties, including surface mass densities and measures of SF. The HI velocity dispersion is correlated most strongly with surface mass density, which points at a gravitational origin for turbulence, but it is unclear which instabilities can operate efficiently in these systems. SF energy is produced at a level sufficient to drive HI turbulent motions where SFRSD > 10^-4 Msun yr^-1 kpc^-2. At low SF intensities, SF does not supply enough energy for turbulence, nor does it uniquely determine the velocity dispersion. Nevertheless, SF appears to provide a lower threshold for HI velocity dispersions. We find that coupling efficiency decreases with increasing SFRSD, consistent with a picture where SF couples to the ISM with constant efficiency, but that less of that energy is found in HI at higher SFRSD. We examine a number of potential drivers of HI turbulence, including SF, gravitational instabilities, the magnetorotational instability, and accretion, and find that no single mechanism can drive the observed levels of turbulence at low SFRSD. We discuss possible solutions to this conundrum.

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

New Radial Abundance Gradients for NGC 628 and NGC 2403

Motived by recent ISM studies, we present high quality MMT and Gemini spectroscopic observations of H II regions in the nearby spiral galaxies NGC 628 and NGC 2403 in order to measure their chemical abundance gradients. Using long-slit and multi-object mask optical spectroscopy, we obtained measurements of the temperature sensitive auroral lines [O III] λ4363 and/or [N II] λ5755 at a strength of 4σ or greater in 11 H II regions in NGC 628 and 7 regions in NGC 2403. These observations allow us, for the first time, to derive an oxygen abundance gradient in NGC 628 based solely on "direct" oxygen abundances of H II regions: 12 + log(O/H) = (8.43+/-0.03) + (-0.017+/-0.002) x Rg (dex/kpc), with a dispersion in log(O/H) of σ = 0.10 dex, from 14 regions with a radial coverage of ~2-19 kpc. This is a significantly shallower slope than found by previous "strong-line" abundance studies. In NGC 2403, we derive an oxygen abundance gradient of 12 + log(O/H) = (8.48+/-0.04) + (-0.032+/-0.007) x Rg (dex/kpc), with a dispersion in log(O/H) of σ = 0.07 dex, from 7 H II with a radial coverage of ~1-10 kpc. Additionally, we measure the N, S, Ne, and Ar abundances. We find the N/O ratio decreases with increasing radius for the inner disk, but reaches a plateau past R25 in NGC 628. NGC 2403 also has a negative N/O gradient with radius, but we do not sample the outer disk of the galaxy past R25 and so do not see evidence for a plateau. This bi-modal pattern measured for NGC 628 indicates dominant contributions from secondary nitrogen inside of the R25 transition and dominantly primary nitrogen farther out. As expected for α-process elements, S/O, Ne/O, and Ar/O are consistent with constant values over a range in oxygen abundance.

preprint2013arXiv

The ACS LCID project. IX. Imprints of the early Universe in the radial variation of the star formation history of dwarf galaxies

Based on Hubble Space Telescope observations from the Local Cosmology from Isolated Dwarfs project, we present the star formation histories, as a function of galactocentric radius, of four isolated Local Group dwarf galaxies: two dSph galaxies, Cetus and Tucana, and two transition galaxies (dTrs), LGS-3 and Phoenix. The oldest stellar populations of the dSphs and dTrs are, within the uncertainties, coeval ($\sim 13 Gyr$) at all galactocentric radii. We find that there are no significative differences between the four galaxies in the fundamental properties (such as the normalized star formation rate or age-metallicity relation) of their outer regions (radii greater than four exponential scale lengths); at large radii, these galaxies consist exclusively of old ($\geq 10.5 Gyr$) metal-poor stars. The duration of star formation in the inner regions vary from galaxy to galaxy, and the extended central star formation in the dTrs produces the dichotomy between dSph and dTr galaxy types. The dTr galaxies show prominent radial stellar population gradients: the centers of these galaxies host young ($\leq 1 Gyr$) populations while the age of the last formation event increases smoothly with increasing radius. This contrasts with the two dSph galaxies. Tucana shows a similar, but milder, gradient, but no gradient in age is detected Cetus. For the three galaxies with significant stellar population gradients, the exponential scale length decreases with time. These results are in agreement with outside-in scenarios of dwarf galaxy evolution, in which a quenching of the star formation toward the center occurs as the galaxy runs out of gas in the outskirts.

preprint2013arXiv

The ACS LCID Project. VIII. The short-period Cepheids of Leo A

We present the results of a new search for variable stars in the Local Group dwarf galaxy Leo A, based on deep photometry from the Advanced Camera for Surveys onboard the Hubble Space Telescope. We detected 166 bona fide variables in our field, of which about 60 percent are new discoveries, and 33 candidate variables. Of the confirmed variables, we found 156 Cepheids, but only 10 RR Lyrae stars despite nearly 100 percent completeness at the magnitude of the horizontal branch. The RR Lyrae stars include 7 fundamental and 3 first-overtone pulsators, with mean periods of 0.636 and 0.366 day, respectively. From their position on the period-luminosity (PL) diagram and light-curve morphology, we classify 91, 58, and 4 Cepheids as fundamental, first-overtone, and second-overtone mode Classical Cepheids (CC), respectively, and two as population II Cepheids. However, due to the low metallicity of Leo A, about 90 percent of the detected Cepheids have periods shorter than 1.5 days. Comparison with theoretical models indicate that some of the fainter stars classified as CC could be Anomalous Cepheids. We estimate the distance to Leo A using the tip of the RGB (TRGB) and various methods based on the photometric and pulsational properties of the Cepheids and RR Lyrae stars. The distances obtained with the TRGB and RR Lyrae stars agree well with each other while that from the Cepheid PL relations is somewhat larger, which may indicate a mild metallicity effect on the luminosity of the short-period Cepheids. Due to its very low metallicity, Leo A thus serves as a valuable calibrator of the metallicity dependencies of the variable star luminosities.

preprint2013arXiv

The ACS Nearby Galaxy Survey Treasury XI. The Remarkably Undisturbed NGC 2403 Disk

We present detailed analysis of color-magnitude diagrams of NGC2403, obtained from a deep (m<28) Hubble Space Telescope (HST) Wide Field Planetary Camera 2 observation of the outer disk of NGC2403, supplemented by several shallow (m<26) HST Advanced Camera for Surveys fields. We derive the spatially resolved star formation history of NGC2403 out to 11 disk scale lengths. In the inner portions of the galaxy, we compare the recent star formation rates (SFRs) we derive from the resolved stars with those measured using GALEX FUV + Spitzer 24-micron fluxes, finding excellent agreement between the methods. Our measurements also show that the radial gradient in recent SFR mirror s the disk exponential profile to 11 scale lengths with no break, extending to SFR densities a factor of 100 lower than those that can be measured with GALEX and Spitzer (2x10^{-6} M_{\sun} yr^{-1} kpc^{-2}). Furthermore, we find that the cumulative stellar mass of the disk was formed at similar times at all radii. We compare these characteristics of NGC2403 to those of its "morphological twins," NGC300 and M33, showing that the structure and age distributions of the NGC2403 disk are more similar to those of the relatively isolated system NGC300 than to those of the Local Group analog M33. We also discuss the environments and HI morphologies of these three nearby galaxies, comparing them to integrated light studies of larger samples of more distant galaxy disks. Taken together, the physical properties and evolutionary history of NGC2403 suggest that the galaxy has had no close encounters with other M81 group members and may be falling into the group for the first time.

preprint2013arXiv

The primordial helium abundance from updated emissivities

Observations of metal-poor extragalactic H II regions allow the determination of the primordial helium abundance, Y_p. The He I emissivities are the foundation of the model of the H II region's emission. Porter, Ferland, Storey, & Detisch (2012) have recently published updated He I emissivities based on improved photoionization cross-sections. We incorporate these new atomic data and update our recent Markov Chain Monte Carlo analysis of the dataset published by Izotov, Thuan, & Stasinska (2007). As before, cuts are made to promote quality and reliability, and only solutions which fit the data within 95% confidence level are used to determine the primordial He abundance. The previously qualifying dataset is almost entirely retained and with strong concordance between the physical parameters. Overall, an upward bias from the new emissivities leads to a decrease in Y_p. In addition, we find a general trend to larger uncertainties in individual objects (due to changes in the emissivities) and an increased variance (due to additional objects included). From a regression to zero metallicity, we determine Y_p = 0.2465 +/- 0.0097, in good agreement with the Planck result of Y_p = 0.2485 +/- 0.0002. In the future, a better understanding of why a large fraction of spectra are not well fit by the model will be crucial to achieving an increase in the precision of the primordial helium abundance determination.

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

Observational Constraints on the Molecular Gas Content in Nearby Starburst Dwarf Galaxies

Using star formation histories derived from optically resolved stellar populations in nineteen nearby starburst dwarf galaxies observed with the Hubble Space Telescope, we measure the stellar mass surface densities of stars newly formed in the bursts. By assuming a star formation efficiency (SFE), we then calculate the inferred gas surface densities present at the onset of the starbursts. Assuming a SFE of 1%, as is often assumed in normal star-forming galaxies, and assuming that the gas was purely atomic, translates to very high HI surface densities (~10^2-10^3 Msun pc^-2), which are much higher than have been observed in dwarf galaxies. This implies either higher values of SFE in these dwarf starburst galaxies or the presence of significant amounts of H_2 in dwarfs (or both). Raising the assumed SFEs to 10% or greater (in line with observations of more massive starbursts associated with merging galaxies), still results in HI surface densities higher than observed in 10 galaxies. Thus, these observations appear to require that a significant fraction of the gas in these dwarf starbursts galaxies was in the molecular form at the onset of the bursts. Our results imply molecular gas column densities in the range 10^19-10^21 cm^-2 for the sample. In those galaxies where CO observations have been made, these densities correspond to values of the CO-H_2 conversion factor (X_CO) in the range >3-80x10^20 cm^-2 (K km s^-1)^-1, or up to 40x greater than Galactic X_CO values.

preprint2012arXiv

On the Origin of the Supergiant HI Shell and Putative Companion in NGC 6822

We present new Hubble Space Telescope Advanced Camera for Surveys imaging of six positions spanning 5.8 kpc of the HI major axis of the Local Group dIrr NGC 6822, including both the putative companion galaxy and the large HI hole. The resulting deep color magnitude diagrams show that NGC 6822 has formed >50% of its stars in the last ~5 Gyr. The star formation histories of all six positions are similar over the most recent 500 Myr, including low-level star formation throughout this interval and a weak increase in star formation rate during the most recent 50 Myr. Stellar feedback can create the giant HI hole, assuming that the lifetime of the structure is longer than 500 Myr; such long-lived structures have now been observed in multiple systems and may be the norm in galaxies with solid-body rotation. The old stellar populations (red giants and red clump stars) of the putative companion are consistent with those of the extended halo of NGC 6822; this argues against the interpretation of this structure as a bona fide interacting companion galaxy and against its being linked to the formation of the HI hole via an interaction. Since there is no evidence in the stellar population of a companion galaxy, the most likely explanation of the extended HI structure in NGC 6822 is a warped disk inclined to the line of sight.

preprint2012arXiv

PHAT Stellar Cluster Survey I. Year 1 Catalog and Integrated Photometry

The Panchromatic Hubble Andromeda Treasury (PHAT) survey is an on-going Hubble Space Telescope (HST) multi-cycle program to obtain high spatial resolution imaging of one-third of the M31 disk at ultraviolet through near-infrared wavelengths. In this paper, we present the first installment of the PHAT stellar cluster catalog. When completed, the PHAT cluster catalog will be among the largest and most comprehensive surveys of resolved star clusters in any galaxy. The exquisite spatial resolution achieved with HST has allowed us to identify hundreds of new clusters that were previously inaccessible with existing ground-based surveys. We identify 601 clusters in the Year 1 sample, representing more than a factor of four increase over previous catalogs within the current survey area (390 arcmin^2). This work presents results derived from the first \sim25% of the survey data; we estimate that the final sample will include \sim2500 clusters. For the Year 1 objects, we present a catalog with positions, radii, and six-band integrated photometry. Along with a general characterization of the cluster luminosities and colors, we discuss the cluster luminosity function, the cluster size distributions, and highlight a number of individually interesting clusters found in the Year 1 search.

preprint2012arXiv

The Chandra Local Volume Survey: The X-ray Point Source Catalog of NGC 300

We present the source catalog of a new Chandra ACIS-I observation of NGC 300 obtained as part of the Chandra Local Volume Survey. Our 63 ks exposure covers ~88% of the D25 isophote (R~6.3 kpc) and yields a catalog of 95 X-ray point sources detected at high significance to a limiting unabsorbed 0.35-8 keV luminosity of ~10^36 erg s^-1. Sources were cross-correlated with a previous XMM-Newton catalog, and we find 75 "X-ray transient candidate" sources that were detected by one observatory, but not the other. We derive an X-ray scale length of 1.7+/-0.2 kpc and a recent star formation rate of 0.12 Msun yr^-1, in excellent agreement with optical observations. Deep, multi-color imaging from the Hubble Space Telescope, covering ~32% of our Chandra field, was used to search for optical counterparts to the X-ray sources, and we have developed a new source classification scheme to determine which sources are likely X-ray binaries, supernova remnants, and background AGN candidates. Finally, we present the X-ray luminosity functions (XLFs) at different X-ray energies, and we find the total NGC 300 X-ray point source population to be consistent with other late-type galaxies hosting young stellar populations (<50 Myr). We find XLF of sources associated with older stellar populations has a steeper slope than the XLF of X-ray sources coinciding with young stellar populations, consistent with theoretical predictions.

preprint2012arXiv

The Nature of Starbursts: III. The Spatial Distribution of Star Formation

We map the spatial distribution of recent star formation over a few x 100 Myr timescales in fifteen starburst dwarf galaxies using the location of young blue helium burning stars identified from optically resolved stellar populations in archival Hubble Space Telescope observations. By comparing the star formation histories from both the high surface brightness central regions and the diffuse outer regions, we measure the degree to which the star formation has been centrally concentrated during the galaxies' starbursts, using three different metrics for the spatial concentration. We find that the galaxies span a full range in spatial concentration, from highly centralized to broadly distributed star formation. Since most starbursts have historically been identified by relatively short timescale star formation tracers (e.g., Halpha emission), there could be a strong bias towards classifying only those galaxies with recent, centralized star formation as starbursts, while missing starbursts that are spatially distributed.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury I: Bright UV Stars in the Bulge of M31

As part of the Panchromatic Hubble Andromeda Treasury (PHAT) multi-cycle program, we observed a 12' \times 6.5' area of the bulge of M31 with the WFC3/UVIS filters F275W and F336W. From these data we have assembled a sample of \sim4000 UV-bright, old stars, vastly larger than previously available. We use updated Padova stellar evolutionary tracks to classify these hot stars into three classes: Post-AGB stars (P-AGB), Post-Early AGB (PE-AGB) stars and AGB-manqué stars. P-AGB stars are the end result of the asymptotic giant branch (AGB) phase and are expected in a wide range of stellar populations, whereas PE-AGB and AGB-manqué (together referred to as the hot post-horizontal branch; HP-HB) stars are the result of insufficient envelope masses to allow a full AGB phase, and are expected to be particularly prominent at high helium or α abundances when the mass loss on the RGB is high. Our data support previous claims that most UV-bright sources in the bulge are likely hot (extreme) horizontal branch stars (EHB) and their progeny. We construct the first radial profiles of these stellar populations, and show that they are highly centrally concentrated, even more so than the integrated UV or optical light. However, we find that this UV-bright population does not dominate the total UV luminosity at any radius, as we are detecting only the progeny of the EHB stars that are the likely source of the UVX. We calculate that only a few percent of MS stars in the central bulge can have gone through the HP-HB phase and that this percentage decreases strongly with distance from the center. We also find that the surface density of hot UV-bright stars has the same radial variation as that of low-mass X-ray binaries. We discuss age, metallicity, and abundance variations as possible explanations for the observed radial variation in the UV-bright population.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury III. Measuring Ages and Masses of Partially Resolved Stellar Clusters

The apparent age and mass of a stellar cluster can be strongly affected by stochastic sampling of the stellar initial mass function, when inferred from the integrated color of low mass clusters (less than ~10^4 solar masses). We use simulated star clusters to show that these effects are minimized when the brightest, rapidly evolving stars in a cluster can be resolved, and the light of the fainter, more numerous unresolved stars can be analyzed separately. When comparing the light from the less luminous cluster members to models of unresolved light, more accurate age estimates can be obtained than when analyzing the integrated light from the entire cluster under the assumption that the initial mass function is fully populated. We show the success of this technique first using simulated clusters, and then with a stellar cluster in M31. This method represents one way of accounting for the discrete, stochastic sampling of the stellar initial mass function in less massive clusters and can be leveraged in studies of clusters throughout the Local Group and other nearby galaxies.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury IV. A Probabilistic Approach to Inferring the High Mass Stellar Initial Mass Function and Other Power-law Functions

We present a probabilistic approach for inferring the parameters of the present day power-law stellar mass function (MF) of a resolved young star cluster. This technique (a) fully exploits the information content of a given dataset; (b) accounts for observational uncertainties in a straightforward way; (c) assigns meaningful uncertainties to the inferred parameters; (d) avoids the pitfalls associated with binning data; and (e) is applicable to virtually any resolved young cluster, laying the groundwork for a systematic study of the high mass stellar MF (M > 1 Msun). Using simulated clusters and Markov chain Monte Carlo sampling of the probability distribution functions, we show that estimates of the MF slope, α, are unbiased and that the uncertainty, Δα, depends primarily on the number of observed stars and stellar mass range they span, assuming that the uncertainties on individual masses and the completeness are well-characterized. Using idealized mock data, we compute the lower limit precision on α and provide an analytic approximation for Δα as a function of the observed number of stars and mass range. We find that ~ 3/4 of quoted literature uncertainties are smaller than the theoretical lower limit. By correcting these uncertainties to the theoretical lower limits, we find the literature studies yield <α>=2.46 with a 1-σ dispersion of 0.35 dex. We verify that it is impossible for a power-law MF to obtain meaningful constraints on the upper mass limit of the IMF. We show that avoiding substantial biases in the MF slope requires: (1) including the MF as a prior when deriving individual stellar mass estimates; (2) modeling the uncertainties in the individual stellar masses; and (3) fully characterizing and then explicitly modeling the completeness for stars of a given mass. (abridged)

preprint2012arXiv

Tracing Cold HI Gas in Nearby, Low-Mass Galaxies

We analyze line-of-sight atomic hydrogen (HI) line profiles of 31 nearby, low-mass galaxies selected from the Very Large Array - ACS Nearby Galaxy Survey Treasury (VLA-ANGST) and The HI Nearby Galaxy Survey (THINGS) to trace regions containing cold (T $\lesssim$ 1400 K) HI from observations with a uniform linear scale of 200 pc/beam. Our galaxy sample spans four orders of magnitude in total HI mass and nine magnitudes in M_B. We fit single and multiple component functions to each spectrum to isolate the cold, neutral medium given by a low dispersion (<6 km/s) component of the spectrum. Most HI spectra are adequately fit by a single Gaussian with a dispersion of 8-12 km/s. Cold HI is found in 23 of 27 (~85%) galaxies after a reduction of the sample size due to quality control cuts. The cold HI contributes ~20% of the total line-of-sight flux when found with warm HI. Spectra best fit by a single Gaussian, but dominated by cold HI emission (i.e., have velocity dispersions <6 km/s) are found primarily beyond the optical radius of the host galaxy. The cold HI is typically found in localized regions and is generally not coincident with the very highest surface density peaks of the global HI distribution (which are usually areas of recent star formation). We find a lower limit for the mass fraction of cold-to-total HI gas of only a few percent in each galaxy.

preprint2012arXiv

VLA-ANGST: A high-resolution HI Survey of Nearby Dwarf Galaxies

We present the "Very Large Array survey of Advanced Camera for Surveys Nearby Galaxy Survey Treasury galaxies (VLA-ANGST)." VLA-ANGST is a National Radio Astronomy Observatory Large Program consisting of high spectral (0.6-2.6 km/s) and spatial (~6") resolution observations of neutral, atomic hydrogen (HI) emission toward 35 nearby dwarf galaxies from the ANGST survey. ANGST is a systematic HST survey to establish a legacy of uniform multi-color photometry of resolved stars for a volume-limited sample of nearby galaxies (D\lesssim4 Mpc). VLA-ANGST provides VLA HI observations of the sub-sample of ANGST galaxies with recent star formation that are observable from the northern hemisphere and that were not observed in the "The HI Nearby Galaxy Survey" (THINGS). The overarching scientific goal of VLA-ANGST is to investigate fundamental characteristics of the neutral interstellar medium (ISM) of dwarf galaxies. Here we describe the VLA observations, the data reduction, and the final VLA-ANGST data products. We present an atlas of the integrated HI maps, the intensity-weighted velocity fields, the second moment maps as a measure for the velocity dispersion of the HI, individual channel maps, and integrated HI spectra for each VLA-ANGST galaxy. We closely follow the observational setup and data reduction of THINGS to achieve comparable sensitivity and angular resolution. A major difference, however, is the high velocity resolution of the VLA-ANGST observations (0.65 and 1.3km/s for the majority of the galaxies). The VLA-ANGST data products are made publicly available at: https://science.nrao.edu/science/surveys/vla-angst. With available star formation histories from resolved stellar populations and lower resolution ancillary observations from the FIR to the UV, VLA-ANGST will enable detailed studies of the relationship between the ISM and star formation in dwarf galaxies on a ~100 pc scale.

preprint2011arXiv

A Deep Chandra View of the NGC 404 Central Engine

We present the results of a 100 ks {\it Chandra} observation of the NGC 404 nuclear region. The long exposure and excellent spatial resolution of {\it Chandra} has enabled us to critically examine the nuclear environment of NGC 404, which is known to host a nuclear star cluster and potentially an intermediate-mass black hole (on the order of a few times $10^5$ \Msun). We find two distinct X-ray sources: a hard, central point source coincident with the optical and radio centers of the galaxy, and a soft extended region that is coincident with areas of high H$α$ emission and likely recent star formation. When we fit the 0.3-8 keV spectra of each region separately, we find the hard nuclear point source to be dominated by a power law (\PL = 1.88), while the soft off-nuclear region is best fit by a thermal plasma model ($kT$ = 0.67 keV). We therefore find evidence for both a power law component and hot gas in the nuclear region of NGC 404. We estimate the 2-10 keV luminosity to be 1.3$^{+0.8}_{-0.5}\times10^{37}$ erg s$^{-1}$. A low-level of diffuse X-ray emission is detected out to $\sim$15\as ($\sim$0.2 kpc) from the nucleus. We compare our results to the observed relationships between power law photon index and Eddington ratio for both X-ray binaries and low luminosity active galaxies and find NGC 404 to be consistent with other low luminosity active galaxies. We therefore favor the conclusion that NGC 404 harbors an intermediate-mass black hole accreting at a very low level.

preprint2011arXiv

An MCMC determination of the primordial helium abundance

Spectroscopic observations of the chemical abundances in metal-poor H II regions provide an independent method for estimating the primordial helium abundance. H II regions are described by several physical parameters such as electron density, electron temperature, and reddening, in addition to y, the ratio of helium to hydrogen. It had been customary to estimate or determine self-consistently these parameters to calculate y. Frequentist analyses of the parameter space have been shown to be successful in these determinations, and Markov Chain Monte Carlo (MCMC) techniques have proven to be very efficient in sampling this parameter space. Nevertheless, accurate determination of the primordial helium abundance from observations of H II regions is constrained by both systematic and statistical uncertainties. In an attempt to better reduce the latter, and better characterize the former, we apply MCMC methods to the large dataset recently compiled by Izotov, Thuan, & Stasinska (2007). To improve the reliability of the determination, a high quality dataset is needed. In pursuit of this, a variety of cuts are explored. The efficacy of the He I 4026 emission line as a constraint on the solutions is first examined, revealing the introduction of systematic bias through its absence. As a clear measure of the quality of the physical solution, a χ^2 analysis proves instrumental in the selection of data compatible with the theoretical model. In addition, the method also allows us to exclude systems for which parameter estimations are statistical outliers. As a result, the final selected dataset gains in reliability and exhibits improved consistency. Regression to zero metallicity yields Y_p = 0.2534 \pm 0.0083, in broad agreement with the WMAP result. The inclusion of more observations shows promise for further reducing the uncertainty, but more high quality spectra are required.

preprint2011arXiv

How Typical Are The Local Group Dwarf Galaxies?

We compare the cumulative star formation histories (SFHs) of Local Group (LG) dwarf galaxies with those in the volume-limited ACS Nearby Galaxy Survey Treasury (ANGST) sample (D < 4 Mpc), in order to understand how typical the LG dwarf galaxies are relative to those in the nearby universe. The SFHs were derived in a uniform manner from high quality optical color-magnitude diagrams constructed from Hubble Space Telescope imaging. We find that the {\it mean} cumulative SFHs of the LG dwarfs are comparable to the mean cumulative SFHs of the ANGST sample for the three different morphological types (dwarf spheroidals/ellipticals: dSph/dE; dwarf irregulars: dI; transition dwarfs: dTrans). We also discuss effects such as population gradients and systematic uncertainties in the stellar models that may influence the derived SFHs. Both the ANGST and Local Group dwarf galaxies show a consistent and strong morphology-density relationship, emphasizing the importance of environment in the evolution of dwarf galaxies. Specifically, we confirm that dIs are found at lower densities and higher luminosities than dSphs, within this large sample. We also find that dTrans are located in similar environments to those occupied by dwarf irregular galaxies, but have systematically lower luminosities that are more comparable to those of dwarf spheroidals. The similarity of the SFHs and morphology-density relationships of the LG and ANGST dwarf galaxies suggests that the LG dwarfs are a good representation of dwarf galaxies in the local universe.

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

Observational Constraints on Red and Blue Helium Burning Sequences

We derive the optical luminosity, colors, and ratios of the blue and red helium burning (HeB) stellar populations from archival Hubble Space Telescope observations of nineteen starburst dwarf galaxies and compare them with theoretical isochrones from Padova stellar evolution models across metallicities from Z=0.001 to 0.009. We find that the observational data and the theoretical isochrones for both blue and red HeB populations overlap in optical luminosities and colors and the observed and predicted blue to red HeB ratios agree for stars older than 50 Myr over the time bins studied. These findings confirm the usefulness of applying isochrones to interpret observations of HeB populations. However, there are significant differences, especially for the red HeB population. Specifically we find: (1) offsets in color between the observations and theoretical isochrones of order 0.15 mag (0.5 mag) for the blue (red) HeB populations brighter than M_V ~ -4 mag, which cannot be solely due to differential extinction; (2) blue HeB stars fainter than M_V ~ -3 mag are bluer than predicted; (3) the slope of the red HeB sequence is shallower than predicted by a factor of ~3; and (4) the models overpredict the ratio of the most luminous blue to red HeB stars corresponding to ages <50 Myr. Additionally, we find that for the more metal-rich galaxies in our sample (Z> 0.5 Zsolar) the red HeB stars overlap with the red giant branch stars in the color magnitude diagrams, thus reducing their usefulness as indicators of star formation for ages >100 Myr.

preprint2011arXiv

Re-examining High Abundance SDSS Mass-Metallicity Outliers: High N/O, Evolved Wolf-Rayet Galaxies?

We present new MMT spectroscopic observations of four dwarf galaxies representative of a larger sample observed by the Sloan Digital Sky Survey (SDSS) and identified by Peeples et al. (2008) as low-mass, high oxygen abundance outliers from the mass-metallicity relation. Peeples et al. (2008) showed that these four objects (with metallicity estimates of 8.5 =< 12 + log(O/H) =< 8.8) have oxygen abundance offsets of 0.4-0.6 dex from the M_B luminosity-metallicity relation. Our new observations extend the wavelength coverage to include the [OII] 3726,3729 doublet, which adds leverage in oxygen abundance estimates and allows measurements of N/O ratios. All four spectra are low excitation, with relatively high N/O ratios (N/O >~ 0.10), each of which tend to bias estimates based on strong emission lines toward high oxygen abundances. These spectra all fall in a regime where the "standard" strong line methods for metallicity determinations are not well calibrated either empirically or by photoionization modeling. By comparing our spectra directly to photoionization models, we estimate oxygen abundances in the range of 7.9 =< 12 + log(O/H) =< 8.4, consistent with the scatter of the mass-metallicity relation. We discuss the physical nature of these galaxies that leads to their unusual spectra (and previous classification as outliers), finding their low excitation, elevated N/O, and strong Balmer absorption are consistent with the properties expected from galaxies evolving past the "Wolf-Rayet galaxy" phase. We compare our results to the "main" sample of Peeples et al. (2008) and conclude that they are outliers primarily due to enrichment of nitrogen relative to oxygen, and not due to unusually high oxygen abundances for their masses or luminosities.

preprint2011arXiv

The ACS Nearby Galaxy Survey Treasury VII. The NGC 4214 Starburst and the Effects of Star Formation History on Dwarf Morphology

We present deep Hubble Space Telescope Wide Field Planetary Camera 2 (WFPC2) optical observations obtained as part of the ACS Nearby Galaxy Survey Treasury (ANGST) as well as early release Wide Field Camera 3 (WFC3) ultra-violet and infrared observations of the nearby dwarf starbursting galaxy NGC 4214. Our data provide a detailed example of how covering such a broad range in wavelength provides a powerful tool for constraining the physical properties of stellar populations. The deepest data reach the ancient red clump at M_F814W -0.2. All of the optical data reach the main sequence turnoff for stars younger than ~300 Myr, and the blue He burning sequence for stars younger than 500 Myr. The full CMD-fitting analysis shows that all three fields in our data set are consistent with ~75% of the stellar mass being older than 8 Gyr, in spite of showing a wide range in star formation rates at the present day. Thus, our results suggest that the scale length of NGC 4214 has remained relatively constant for many Gyr. As previously noted by others, we also find the galaxy has recently ramped up production, consistent with its bright UV luminosity and its population of UV-bright massive stars. In the central field we find UV point sources with F336W magnitudes as bright as -9.9. These are as bright as stars with masses of at least 52-56 M_sun and ages near 4 Myr in stellar evolution models. Assuming a standard IMF, our CMD is well-fitted by an increase in star formation rate beginning 100 Myr ago. The stellar populations of this late-type dwarf are compared with those of NGC 404, an early-type dwarf that is also the most massive galaxy in its local environment. The late-type dwarf appears to have a similar high fraction of ancient stars, suggesting that these dominant galaxies may form at early epochs even if they have low total mass and very different present-day morphologies.

preprint2011arXiv

The ACS Nearby Galaxy Survey Treasury VIII. The Global Star Formation Histories of 60 Dwarf Galaxies in the Local Volume

We present uniformly measured star formation histories (SFHs) of 60 nearby (D~4Mpc) dwarf galaxies based on CMDs of resolved stellar populations from images taken with HST as part of the ACS Nearby Galaxy Survey Treasury program (ANGST). This volume-limited sample contains 12 dSph/dE, 5 dwarf spiral, 28 dI, 12 transition, and 3 tidal dwarf galaxies. From the best fit SFHs we find three significant results: (1) the average dwarf galaxy formed >50% of its stars by z~2 and 60% of its stars by z~1, regardless of current morphological type; (2) the mean SFHs of dIs, dTrans, and dSphs are similar over most of cosmic time, with the clearest differences appearing during the most recent 1 Gyr; and (3) the mean values are inconsistent with simple SFH models, e.g., exponentially declining SFRs. The mean SFHs are in general agreement with the cosmic SFH, although we observe offsets near z~1 that could be evidence that low mass systems experienced delayed star formation relative to more massive galaxies. The sample shows a strong density-morphology relationship, such that the dSphs are less isolated than dIs. We find that the transition from a gas-rich to gas-poor galaxy cannot be solely due to internal mechanisms such as stellar feedback, and instead is likely the result of external mechanisms, e.g., ram pressure and tidal forces. The average transition dwarf is slightly less isolated and less gas-rich than the typical dI. Further, the transition dwarfs can be divided into two groups: interacting and gas-poor or isolated and gas-rich, suggesting two possible evolutionary pathways.

preprint2011arXiv

The Formation of Kiloparsec-Scale HI Holes in Dwarf Galaxies

The origin of kpc-scale holes in the atomic hydrogen (H I) distributions of some nearby dwarf irregular galaxies presents an intriguing problem. Star formation histories (SFHs) derived from resolved stars give us the unique opportunity to study past star forming events that may have helped shape the currently visible H I distribution. Our sample of five nearby dwarf irregular galaxies spans over an order of magnitude in both total H I mass and absolute B-band magnitude and is at the low mass end of previously studied systems. We use Very Large Array H I line data to estimate the energy required to create the centrally dominant hole in each galaxy. We compare this energy estimate to the past energy released by the underlying stellar populations computed from SFHs derived from data taken with the Hubble Space Telescope. The inferred integrated stellar energy released within the characteristic ages exceeds our energy estimates for creating the holes in all cases, assuming expected efficiencies. Therefore, it appears that stellar feedback provides sufficient energy to produce the observed holes. However, we find no obvious signature of single star forming events responsible for the observed structures when comparing the global SFHs of each galaxy in our sample to each other or to those of dwarf irregular galaxies reported in the literature. We also fail to find evidence of a central star cluster in FUV or Halpha imaging. We conclude that large H I holes are likely formed from multiple generations of star formation and only under suitable interstellar medium conditions.

preprint2011arXiv

The M81 Group Dwarf Irregular Galaxy DDO 165. I. High Velocity Neutral Gas in a Post-Starburst System

We present new multi-configuration VLA HI spectral line observations of the M81 group dIrr post-starburst galaxy DDO 165. The HI morphology is complex, with multiple column density peaks surrounding a large region of very low HI surface density that is offset from the center of the stellar distribution. The bulk of the neutral gas is associated with the southern section of the galaxy; a secondary peak in the north contains ~15% of the total HI mass. These components appear to be kinematically distinct, suggesting that either tidal processes or large-scale blowout have recently shaped the ISM of DDO 165. Using spatially-resolved position-velocity maps, we find multiple localized high-velocity gas features. Cross-correlating with radius-velocity analyses, we identify eight shell/hole structures in the ISM with a range of sizes (~400-900 pc) and expansion velocities (~7-11 km/s). These structures are compared with narrow- and broad-band imaging from KPNO and HST. Using the latter data, recent works have shown that DDO 165's previous "burst" phase was extended temporally (>1 Gyr). We thus interpret the high-velocity gas features, HI holes, and kinematically distinct components of the galaxy in the context of the immediate effects of "feedback" from recent star formation. In addition to creating HI holes and shells, extended star formation events are capable of creating localized high velocity motion of the surrounding interstellar material. A companion paper connects the energetics from the HI and HST data.

preprint2011arXiv

The M81 Group Dwarf Irregular Galaxy DDO 165. II. Connecting Recent Star Formation with ISM Structures and Kinematics

We compare the stellar populations and complex neutral gas dynamics of the M81 group dIrr galaxy DDO 165 using data from the HST and the VLA. Paper I identified two kinematically distinct HI components, multiple localized high velocity gas features, and eight HI holes and shells (the largest of which spans ~2.2x1.1 kpc). Using the spatial and temporal information from the stellar populations in DDO 165, we compare the patterns of star formation over the past 500 Myr with the HI dynamics. We extract localized star formation histories within 6 of the 8 HI holes identified in Paper I, as well as 23 other regions that sample a range of stellar densities and neutral gas properties. From population synthesis modeling, we derive the energy outputs (from stellar winds and supernovae) of the stellar populations within these regions over the last 100 Myr, and compare with refined estimates of the energies required to create the HI holes. In all cases, we find that "feedback" is energetically capable of creating the observed structures in the ISM. Numerous regions with significant energy inputs from feedback lack coherent HI structures but show prominent localized high velocity gas features; this feedback signature is a natural product of temporally and spatially distributed star formation. In DDO 165, the extended period of heightened star formation activity (lasting more than 1 Gyr) is energetically capable of creating the observed holes and high velocity gas features in the neutral ISM.

preprint2011arXiv

The Survey of HI in Extremely Low-mass Dwarfs (SHIELD)

We present first results from the "Survey of HI in Extremely Low-mass Dwarfs" (SHIELD), a multi-configuration EVLA study of the neutral gas contents and dynamics of galaxies with HI masses in the 10^6-10^7 Solar mass range detected by the Arecibo Legacy Fast ALFA (ALFALFA) survey. We describe the survey motivation and concept demonstration using VLA imaging of 6 low-mass galaxies detected in early ALFALFA data products. We then describe the primary scientific goals of SHIELD and present preliminary EVLA and WIYN 3.5m imaging of the 12 SHIELD galaxies. With only a few exceptions, the neutral gas distributions of these extremely low-mass galaxies are centrally concentrated. In only 1 system have we detected HI column densities higher than 10^21 cm^-2. Despite this, the stellar populations of all of these systems are dominated by blue stars. Further, we find ongoing star formation as traced by H alpha emission in 10 of the 11 galaxies with H alpha imaging obtained to date. Taken together these results suggest that extremely low-mass galaxies are forming stars in conditions different from those found in more massive systems. While detailed dynamical analysis requires the completion of data acquisition, the most well-resolved system is amenable to meaningful position-velocity analysis. For AGC 749237, we find well-ordered rotation of 30 km/s at a distance of ~40 arcseconds from the dynamical center. At the adopted distance of 3.2 Mpc, this implies the presence of a >1x10^8 Solar mass dark matter halo and a baryon fraction < ~0.1.

preprint2010arXiv

A New Approach to Systematic Uncertainties and Self-Consistency in Helium Abundance Determinations

Tests of big bang nucleosynthesis and early universe cosmology require precision measurements for helium abundance determinations. However, efforts to determine the primordial helium abundance via observations of metal poor H II regions have been limited by significant uncertainties. This work builds upon previous work by providing an updated and extended program in evaluating these uncertainties. Procedural consistency is achieved by integrating the hydrogen based reddening correction with the helium based abundance calculation, i.e., all physical parameters are solved for simultaneously. We include new atomic data for helium recombination and collisional emission based upon recent work by Porter et al. and wavelength dependent corrections to underlying absorption are investigated. The set of physical parameters has been expanded here to include the effects of neutral hydrogen collisional emission. Because of a degeneracy between the solutions for density and temperature, the precision of the helium abundance determinations is limited. Also, at lower temperatures (T \lesssim 13,000 K) the neutral hydrogen fraction is poorly constrained resulting in a larger uncertainty in the helium abundances. Thus the derived errors on the helium abundances for individual objects are larger than those typical of previous studies. The updated emissivities and neutral hydrogen correction generally raise the abundance. From a regression to zero metallicity, we find Y_p as 0.2561 \pm 0.0108, in broad agreement with the WMAP result. Tests with synthetic data show a potential for distinct improvement, via removal of underlying absorption, using higher resolution spectra. A small bias in the abundance determination can be reduced significantly and the calculated helium abundance error can be reduced by \sim 25%.

preprint2010arXiv

Mapping systematic errors in helium abundance determinations using Markov Chain Monte Carlo

Monte Carlo techniques have been used to evaluate the statistical and systematic uncertainties in the helium abundances derived from extragalactic H~II regions. The helium abundance is sensitive to several physical parameters associated with the H~II region. In this work, we introduce Markov Chain Monte Carlo (MCMC) methods to efficiently explore the parameter space and determine the helium abundance, the physical parameters, and the uncertainties derived from observations of metal poor nebulae. Experiments with synthetic data show that the MCMC method is superior to previous implementations (based on flux perturbation) in that it is not affected by biases due to non-physical parameter space. The MCMC analysis allows a detailed exploration of degeneracies, and, in particular, a false minimum that occurs at large values of optical depth in the He~I emission lines. We demonstrate that introducing the electron temperature derived from the [O~III] emission lines as a prior, in a very conservative manner, produces negligible bias and effectively eliminates the false minima occurring at large optical depth. We perform a frequentist analysis on data from several "high quality" systems. Likelihood plots illustrate degeneracies, asymmetries, and limits of the determination. In agreement with previous work, we find relatively large systematic errors, limiting the precision of the primordial helium abundance for currently available spectra.

preprint2010arXiv

Spatially Resolved PAH Emission Features in Nearby, Low Metallicity, Star-Forming Galaxies

Low-resolution, mid-infrared Spitzer/IRS spectral maps are presented for three nearby, low-metallicity dwarf galaxies (NGC 55, NGC 3109 and IC 5152) for the purpose of examining the spatial distribution and variation of polycyclic aromatic hydrocarbon (PAH) emission. The sample straddles a metallicity of 12+log(O/H)~8.0, a transition point below which PAH intensity empirically drops and the character of the interstellar medium changes. We derive quantitative radiances of PAH features and atomic lines on both global and spatially-resolved scales. The Spitzer spectra, combined with extensive ancillary data from the UV through the mid-infrared, allow us to examine changes in the physical environments and in PAH feature radiances down to a physical scale of 50 pc. We discuss correlations between various PAH emission feature and atomic line radiances. The (6.2 micron)/(11.3 micron), (7.7 micron)/(11.3 micron), (8.6 micron)/(11.3 micron), (7.7 micron)/(6.2 micron), and (8.6 micron)/(6.2 micron) PAH radiance ratios are found to be independent of position across all three galaxies, although the ratios do vary from galaxy to galaxy. As seen in other galaxies, we find no variation in the grain size distribution as a function of local radiation field strength. Absolute PAH feature intensities as measured by a ratio of PAH/(24 micron) radiances are seen to vary both positionally within a given galaxy, and from one galaxy to another when integrated over the full observed extent of each system. We examine direct comparisons of CC mode PAH ratios (7.7 micron)/(6.2 micron) and (8.6 micron)/(6.2 micron) to the mixed (CC/CH) mode PAH ratio (7.7 micron)/(11.3 micron). We find little variation in either mode, and no difference in trends between modes. While the local conditions change markedly over the observed regions of these galaxies, the properties of PAH emission show a remarkable degree of uniformity.

preprint2010arXiv

The ACS LCID Project. II. Faint Variable Stars in the Isolated Dwarf Irregular Galaxy IC1613

(Abridged) We present the results of a new search for variable stars in the Local Group (LG) isolated dwarf galaxy IC1613, based on 24 orbits of F475W and F814W photometry from the ACS camera onboard the HST. We detected 259 candidate variables in this field, of which only 13 (all of them bright Cepheids) were previously known. Out of the confirmed variables, we found 90 RR Lyrae stars, 49 classical Cepheids (including 36 new discoveries), and 38 eclipsing binary stars for which we could determine a period. The RR Lyrae include 61 fundamental (RRab) and 24 first-overtone (RRc) pulsators, and 5 pulsating in both modes simultaneously (RRd). As for the majority of LG dwarfs, the mean periods of the RRab and RRc (0.611 and 0.334 day, respectively) as well as the fraction of overtone pulsators (f_c=0.28) place this galaxy in the intermediate regime between the Oosterhoff types. From their position on the period-luminosity diagram and light-curve morphology, we can unambiguously classify 25 and 14 Cepheids as fundamental and first-overtone mode pulsators, respectively. Another two are clearly second-overtone Cepheids, the first ones to be discovered beyond the Magellanic Clouds. We estimate the distance to IC1613 using various methods based on the photometric and pulsational properties of the Cepheids and RR Lyrae stars. The values we find are in very good agreement with each other and with previous estimates based on independent methods. When corrected to a common reddening of E(B-V)=0.025 and true LMC distance modulus of (m-M)_{LMC,0}=18.515+-0.085, we find that all the distance determinations from the literature converge to a common value of (m-M)_0=24.400+-0.014 (statistical), or 760 kpc.

preprint2010arXiv

The Nature of Starbursts : II. The Duration of Starbursts in Dwarf Galaxies

The starburst phenomenon can shape the evolution of the host galaxy and the surrounding intergalactic medium. The extent of the evolutionary impact is partly determined by the duration of the starburst, which has a direct correlation with both the amount of stellar feedback and the development of galactic winds, particularly for smaller mass dwarf systems. We measure the duration of starbursts in twenty nearby, ongoing, and "fossil" starbursts in dwarf galaxies based on the recent star formation histories derived from resolved stellar population data obtained with the Hubble Space Telescope. Contrary to the shorter times of 3-10 Myr often cited, the starburst durations we measure range from 450 - 650 Myr in fifteen of the dwarf galaxies and up to 1.3 Gyr in four galaxies; these longer durations are comparable to or longer than the dynamical timescales for each system. The same feedback from massive stars that may quench the flickering SF does not disrupt the overall burst event in our sample of galaxies. While five galaxies present fossil bursts, fifteen galaxies show ongoing bursts and thus the final durations may be longer than we report here for these systems. One galaxy shows a burst that has been ongoing for only 20 Myr; we are likely seeing the beginning of a burst event in this system. Using the duration of the starbursts, we calculate that the bursts deposited 10^(53.9)-10^(57.2) ergs of energy into the interstellar medium through stellar winds and supernovae and produced 3%-26% of the host galaxy's mass.

preprint2010arXiv

The Nature of Starbursts: I. The Star Formation Histories of Eighteen Nearby Starburst Dwarf Galaxies

We use archival HST observations of resolved stellar populations to derive the star formation histories (SFHs) of 18 nearby starburst dwarf galaxies. In this first paper we present the observations, color-magnitude diagrams, and the SFHs of the 18 starburst galaxies, based on a homogeneous approach to the data reduction, differential extinction, and treatment of photometric completeness. We adopt a star formation rate (SFR) threshold normalized to the average SFR of the individual system as a metric for classifying starbursts in SFHs derived from resolved stellar populations. This choice facilitates finding not only currently bursting galaxies but also "fossil" bursts increasing the sample size of starburst galaxies in the nearby (D<8 Mpc) universe. Thirteen of the eighteen galaxies are experiencing ongoing bursts and five galaxies show fossil bursts. From our reconstructed SFHs, it is evident that the elevated SFRs of a burst are sustained for hundreds of Myr with variations on small timescales. A long >100 Myr temporal baseline is thus fundamental to any starburst definition or identification method. The longer lived bursts rule out rapid "self-quenching" of starbursts on global scales. The bursting galaxies' gas consumption timescales are shorter than the Hubble time for all but one galaxy confirming the short-lived nature of starbursts based on fuel limitations. Additionally, we find the strength of the Hα emission usually correlates with the CMD based SFR during the last 4-10 Myr. However, in four cases, the Hα emission is significantly less than what is expected for models of starbursts; the discrepancy is due to the SFR changing on timescales of a few Myr. The inherently short timescale of the Hα emission limits identifying galaxies as starbursts based on the current characteristics which may or may not be representative of the recent SFH of a galaxy.

preprint2009arXiv

A Spitzer Study of Asymptotic Giant Branch Stars. III. Dust Production and Gas Return in Local Group Dwarf Irregular Galaxies

We present the third and final part of a census of Asymptotic Giant Branch (AGB) stars in Local Group dwarf irregular galaxies. Papers I and II presented the results for WLM and IC 1613. Included here are Phoenix, LGS 3, DDO 210, Leo A, Pegasus dIrr, and Sextans A. Spitzer photometry at 3.6, 4.5, 5.8, and 8 are presented, along with a more thorough treatment of background galaxy contamination than was presented in papers I and II. We find that at least a small population of completely optically obscured AGB stars exists in each galaxy, regardless of the galaxy's metallicity, but that higher-metallicity galaxies tend to harbor more stars with slight IR excesses. The optical incompleteness increases for the redder AGB stars, in line with the expectation that some AGB stars are not detected in the optical due to large amounts of extinction associated with in situ dust production. Overall, there is an underrepresentation of 30% - 40% in the optical AGB within the 1 sigma errors for all of the galaxies in our sample. This undetected population is large enough to affect star formation histories derived from optical color-magnitude diagrams. As measured from the [3.6] - [4.5] color excesses, we find average stellar mass-loss rates ranging from 3.1E-7 - 6.6E-6 solar masses per year, and integrated galaxy mass-loss rates ranging from 4.4E-5 - 1.4E-3 solar masses per year. The integrated mass-loss rate is sufficient to sustain the current star formation rate in only LGS 3 and DDO 210, requiring either significant non-dusty mass loss or gas accretion in Phoenix, Leo A, Pegasus dIrr, Sextans A, WLM, and IC 1613 if they are to maintain their status as gas-rich galaxies.

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

Does Stellar Feedback Create HI Holes? An HST/VLA Study of Holmberg II

We use deep HST/ACS F555W and F814W photometry of resolved stars in the M81 Group dwarf irregular galaxy Ho II to study the hypothesis that the holes identified in the neutral ISM (HI) are created by stellar feedback. From the deep photometry, we construct color-magnitude diagrams (CMDs) and measure the star formation histories (SFHs) for stars contained in HI holes from two independent holes catalogs, as well as select control fields, i.e., similar sized regions that span a range of HI column densities. Converting the recent SFHs into stellar feedback energies, we find that enough energy has been generated to have created all holes. However, the required energy is not always produced over a time scale that is less than the estimated kinematic age of the hole. The combination of the CMDs, recent SFHs, and locations of young stars shows that the stellar populations inside HI holes are not coherent, single-aged, stellar clusters, as previously suggested, but rather multi-age populations distributed across each hole. From a comparison of the modeled and observed integrated magnitudes, and the locations and energetics of stars inside of HI holes, we propose a potential new model: a viable mechanism for creating the observed HI holes in Ho II is stellar feedback from multiple generations of SF spread out over tens or hundreds of Myr, and thus, the concept of an age for an HI hole is intrinsically ambiguous. We further find that \halpha and 24 micron emission, tracers of the most recent star formation, do not correlate well with the positions of the HI holes. However, UV emission, which traces star formation over roughly the last 100 Myr, shows a much better correlation with the locations of the HI holes.

preprint2009arXiv

Star Formation in Dwarf Galaxies of the Nearby Centaurus A Group

We present Halpha narrow-band imaging of 17 dwarf irregular galaxies (dIs) in the nearby Centaurus A Group. Although all large galaxies of the group have a current or recent enhanced star formation episode, the dIs have normal star formation rates and do not contain a larger fraction of dwarf starbursts than other nearby groups. Relative distances between dIs and larger galaxies of the group can be computed in 3D since most of them have now fairly accurately known distances. We find that the dI star formation rates do not depend on local environment, and in particular they do not show any correlation with the distance of the dI to the nearest large galaxy of the group. There is a clear morphology-density relation in the Centaurus A Group, similarly to the Sculptor and Local Groups, in the sense that dEs/dSphs tend to be at small distances from the more massive galaxies of the group, while dIs are on average at larger distances. We find four transition dwarfs in the Group, dwarfs that show characteristics of both dE/dSphs and dIs, and which contain cold gas but no current star formation. Interestingly the transition dwarfs have an average distance to the more massive galaxies which is intermediate between those of the dEs/dSphs and dIs, and which is quite large: 0.54 +- 0.31 Mpc. This large distance poses some difficulty for the most popular scenarios proposed for transforming a dI into a dE/dSph (ram-pressure with tidal stripping or galaxy harassment). If the observed transition dwarfs are indeed missing links between dIs and dE/dSphs, their relative isolation makes it less likely to have been produced by these mechanisms. We propose that an inhomogeneous IGM containing higher density clumps would be able to ram-pressure stripped the dIs at such large distances.

preprint2009arXiv

Stellar Metallicities and Kinematics in a Gas-Rich Dwarf Galaxy: First Calcium Triplet Spectroscopy of RGB Stars in WLM

We present the first determination of the radial velocities and metallicities of 78 red giant stars in the isolated dwarf irregular galaxy WLM. Observations of the calcium II triplet in these stars were made with FORS2 at the VLT-UT2 in two separated fields of view in WLM, and the [Fe/H] values were conformed to the Carretta & Gratton (1997) metallicity scale. The mean metallicity is <[Fe/H]> = -1.27 +/- 0.04 dex, with a standard deviation of 0.37. We find that the stars in the inner field are more metal rich by [Fe/H] =0.30 +/- 0.06 dex. These results are in agreement with previous photometric studies that found a radial population gradient, as well as the expectation of higher metallicities in the central star forming regions. Age estimates using Victoria-Regina stellar models show that the youngest stars in the sample (< 6 Gyr) are more metal rich by [Fe/H] = 0.32 +/- 0.08 dex. These stars also show a lower velocity dispersion at all elliptical radii compared to the metal-poor stars. Kinematics for the whole red giant sample suggest a velocity gradient approximately half that of the gas rotation curve, with the stellar component occupying a thicker disk decoupled from the HI rotation plane. Taken together, the kinematics, metallicities, and ages in our sample suggest a young metal-rich, and kinematically cold stellar population in the central gas-rich regions of WLM, surrounded by a separate dynamically hot halo of older, metal poor stars.

preprint2009arXiv

The ACS LCID Project. I. Short-Period Variables in the Isolated Dwarf Spheroidal Galaxies Cetus & Tucana

(abridged) We present the first study of the variable star populations in the isolated dwarf spheroidal galaxies (dSph) Cetus and Tucana. Based on Hubble Space Telescope images obtained with the Advanced Camera for Surveys in the F475W and F814W bands, we identified 180 and 371 variables in Cetus and Tucana, respectively. The vast majority are RR Lyrae stars. In Cetus we also found three anomalous Cepheids, four candidate binaries and one candidate long-period variable (LPV), while six anomalous Cepheids and seven LPV candidates were found in Tucana. Of the RR Lyrae stars, 147 were identified as fundamental mode (RRab) and only eight as first-overtone mode (RRc) in Cetus, with mean periods of 0.614 and 0.363 day, respectively. In Tucana we found 216 RRab and 82 RRc giving mean periods of 0.604 and 0.353 day. These values place both galaxies in the so-called Oosterhoff Gap, as is generally the case for dSph. We calculated the distance modulus to both galaxies using different approaches based on the properties of RRab and RRc, namely the luminosity-metallicity and period-luminosity-metallicity relations, and found values in excellent agreement with previous estimates using independent methods: (m-M)_{0,Cet}=24.46+-0.12 and (m-M)_{0,Tuc}=24.74+-0.12, corresponding to 780+-40 kpc and 890+-50 kpc. We also found numerous RR Lyrae variables pulsating in both modes simultaneously (RRd): 17 in Cetus and 60 in Tucana. Tucana is, after Fornax, the second dSph in which such a large fraction of RRd (~17%) has been observed. We provide the photometry and pulsation parameters for all the variables, and compare the latter with values from the literature for well-studied dSph of the Local Group and Galactic globular clusters.

preprint2009arXiv

The True Durations of Starbursts: HST Observations of Three Nearby Dwarf Starburst Galaxies

The duration of a starburst is a fundamental parameter affecting the evolution of galaxies yet, to date, observational constraints on the durations of starbursts are not well established. Here we study the recent star formation histories (SFHs) of three nearby dwarf galaxies to rigorously quantify the duration of their starburst events using a uniform and consistent approach. We find that the bursts range from ~200 - ~400 Myr in duration resolving the tension between the shorter timescales often derived observationally with the longer timescales derived from dynamical arguments. If these three starbursts are typical of starbursts in dwarf galaxies, then the short timescales (3 - 10 Myr) associated with starbursts in previous studies are best understood as "flickering" events which are simply small components of the larger starburst. In this sample of three nearby dwarfs, the bursts are not localized events. All three systems show bursting levels of star formation in regions of both high and low stellar density. The enhanced star formation moves around the galaxy during the bursts and covers a large fraction of the area of the galaxy. These massive, long duration bursts can significantly affect the structure, dynamics, and chemical evolution of the host galaxy and can be the progenitors of "superwinds" that drive much of the recently chemically enriched material from the galaxy into the intergalactic medium.

preprint2007arXiv

The Aptly Named Phoenix Dwarf Galaxy

The Local Group galaxy Phoenix has the properties of a dwarf spheroidal galaxy, but an adjacent HI cloud has been found to be at the same radial velocity as the stars. The proximity suggests that this cloud is associated with the most recent ($\le 100$ Myr) star formation in Phoenix. We have obtained relatively high sensitivity and high resolution HI imaging with the VLA with the goal of distinguishing between different processes for displacing the gas from the galaxy. Due to the outer curvature of the HI cloud, it appears that expulsion from the galaxy by winds from supernovae is more likely than ram-pressure stripping. The isolation of the galaxy makes tidal stripping highly unlikely. Using a star formation history constructed from HST imaging, we construct a simple kinematic model which implies that the HI cloud is still gravitationally bound to the galaxy. Gas which is expelled from the centers of dwarf galaxies but which remains gravitationally bound may explain the episodic star formation observed in several dwarfs. In the specific case of Phoenix, there may be future star formation in this currently dSph-like galaxy.

preprint2006arXiv

On Extending the Mass-Metallicity Relation of Galaxies by 2.5 Decades in Stellar Mass

We report 4.5 micron luminosities for 27 nearby (D < 5 Mpc) dwarf irregular galaxies measured with the Spitzer Infrared Array Camera. We have constructed the 4.5 micron luminosity-metallicity (L-Z) relation for 25 dwarf galaxies with secure distance and interstellar medium oxygen abundance measurements. The 4.5 micron L-Z relation is 12+log(O/H) = (5.78 +/- 0.21) + (-0.122 +/- 0.012) M_[4.5], where M_[4.5] is the absolute magnitude at 4.5 micron. The dispersion in the near-infrared L-Z relation is smaller than the corresponding dispersion in the optical L-Z relation. The subsequently derived stellar mass-metallicity M-Z relation is 12+log(O/H) = (5.65 +/- 0.23) + (0.298 +/- 0.030) log Mstar. and extends the SDSS M-Z relation to lower mass by about 2.5 dex. We find that the dispersion in the M-Z relation is similar over five orders of magnitude in stellar mass, and that the relationship between stellar mass and interstellar medium metallicity is similarly tight from high-mass to low-mass systems. We find a larger scatter at low mass in the relation between effective yield and total baryonic mass. In fact, there are a few dwarf galaxies with large yields, which is difficult to explain if galactic winds are ubiquitous in dwarf galaxies. The low scatter in the L-Z and M-Z relationships are difficult to understand if galactic superwinds or blowout are responsible for the low metallicities at low mass or luminosity. Naively, one would expect an ever increasing scatter at lower masses, which is not observed.

preprint2004arXiv

Probing The Multiphase Interstellar Medium Of The Dwarf Starburst Galaxy NGC 625 With FUSE Spectroscopy

We present new FUSE spectroscopy of the dwarf starburst galaxy NGC 625. These observations probe multiple phases of the interstellar medium, including the coronal, ionized, neutral and molecular gas. This nearby (D = 3.9 +/- 0.2 Mpc) system shows a clear detection of outflowing coronal gas as traced by OVI 1032 Angstrom absorption. The centroid of the OVI profile is blueshifted with respect to the galaxy systemic velocity by ~ 30 km/sec, suggesting a low-velocity outflow. The implied OVI velocity extent is found to be 100 +/- 20 km/sec, which is fully consistent with the detected HI outflow velocity found in radio synthesis observations. We detect multiple lines of diffuse H2 absorption from the ISM of NGC 625; this is one of only a few extragalactic systems with FUSE detections of H2 lines in the Lyman and Werner bands. We find a potential abundance offset between the neutral and nebular gas that exceeds the errors on the derived column densities. Since such an offset has been found in multiple dwarf galaxies, we discuss the implications of a lower-metallicity halo surrounding the central star forming regions of dwarf galaxies. The apparent offset may be due to saturation of the observed OI line, and higher S/N observations are required to resolve this issue.