BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Benjamin D. Johnson

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.SRastro-ph.IMhep-ph

Catalog footprint

What is connected

45works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

45 published item(s)

preprint2024arXiv

Inside the bubble: exploring the environments of reionisation-era Lyman-$α$ emitting galaxies with JADES and FRESCO

We present a study of the environments of 17 Lyman-$α$ (Ly$α$) emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Ly$α$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Ly$α$ velocity offsets ($Δv_\text{Ly$α$}<300\,\mathrm{km\,s^{-1}}$) with moderately high Ly$α$ escape fractions ($f_\mathrm{esc,\,Lyα}>5\%$) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical `bubble' sizes of the order of $R_\text{ion}\sim0.1$-$1\,\mathrm{pMpc}$ are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z\sim7.3$, suggesting Ly$α$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV}\gtrsim-20\,\mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Ly$α$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of JWST multi-object and slitless spectroscopy in acquiring a unique view of the early Universe during cosmic reionisation via the most distant LAEs.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Molecular Gas Reservoirs in Massive Quiescent Galaxies at $\mathrm{z\sim0.7}$ Linked to Late Time Star Formation

We explore how the presence of detectable molecular gas depends on the inferred star formation histories (SFHs) in 8 massive, quiescent galaxies at $\mathrm{z\sim0.7}$. Half of the sample have clear detections of molecular gas, traced by CO(2-1). We find that the molecular gas content is unrelated to the rate of star formation decline prior to the most recent 1 Gyr, suggesting that the gas reservoirs are not leftover from their primary star formation epoch. However, the recent SFHs of CO-detected galaxies demonstrate evidence for secondary bursts of star formation in their last Gyr. The fraction of stellar mass formed in these secondary bursts ranges from $\mathrm{f_{burst}\approx0.3-6\%}$, and ended between $\mathrm{t_{end\mbox{-}burst}\approx0-330~Myr}$ ago. The CO-detected galaxies form a higher fraction of mass in the last Gyr ($\mathrm{f_{M_{1Gyr}}=2.6\pm1.8\%}$) compared to the CO-undetected galaxies ($\mathrm{f_{M_{1Gyr}}=0.2\pm0.1\%}$). The galaxies with gas reservoirs have enhanced late-time star formation, highlighting this as a contributing factor to the observed heterogeneity in the gas reservoirs in high-redshift quiescent galaxies. We find that the amount of gas and star formation driven by these secondary bursts are inconsistent with that expected from dry minor mergers, and instead are likely driven by recently-accreted gas i.e., gas-rich minor mergers. This conclusion would not have been made based on $\mathrm{SFR_{UV+IR}}$ measurements alone, highlighting the power of detailed SFH modeling in the interpretation of gas reservoirs. Larger samples are needed to understand the frequency of low-level rejuvenation among quiescent galaxies at intermediate redshifts, and to what extent this drives the diversity of molecular gas reservoirs.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Tilt in the Dark Matter Halo of the Galaxy

Recent observations of the stellar halo have uncovered the debris of an ancient merger, Gaia-Sausage-Enceladus, estimated to have occurred ~8 Gyr ago. Follow-up studies have associated GSE with a large-scale tilt in the stellar halo that links two well-known stellar over-densities in diagonally opposing octants of the Galaxy (the Hercules-Aquila Cloud and Virgo Overdensity; HAC and VOD). In this paper, we study the plausibility of such unmixed merger debris persisting over several Gyr in the Galactic halo. We employ the simulated stellar halo from Naidu et al. (2021), which reproduces several key properties of the merger remnant, including the large-scale tilt. By integrating the orbits of these simulated stellar halo particles, we show that adoption of a spherical halo potential results in rapid phase mixing of the asymmetry. However, adopting a tilted halo potential preserves the initial asymmetry in the stellar halo for many Gyr. The asymmetry is preserved even when a realistic growing disk is added to the potential. These results suggest that HAC and VOD are long-lived structures that are associated with GSE and that the dark matter halo of the Galaxy is tilted with respect to the disk and aligned in the direction of HAC-VOD. Such halo-disk misalignment is common in modern cosmological simulations. Lastly, we study the relationship between the local and global stellar halo in light of a tilted global halo comprised of highly radial orbits. We find that the local halo offers a dynamically biased view of the global halo due to its displacement from the Galactic Center.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8

We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $τ_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($τ_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langleτ_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and "late bloomers" are uncommon. In summary, our analysis supports the "grow & quench" framework and is consistent with a wide and continuously-populated diversity of quenching timescales.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

How Well Can We Measure Galaxy Dust Attenuation Curves? The Impact of the Assumed Star-Dust Geometry Model in SED Fitting

One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the Simba cosmological simulation with the Prospector SED fitting code. We consider the impact of the commonly-assumed uniform screen model and introduce a new non-uniform screen model parameterized by the fraction of unobscured stellar light. This non-uniform screen model allows for a non-zero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a non-uniform screen model, with the median offset in UV attenuation decreasing from $-0.30$ dex with a uniform screen model to $-0.17$ dex with the non-uniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the non-uniform screen model, decreasing the SFR offset on average by $0.12$ dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Live Fast, Die $α$-Enhanced: The Mass-Metallicity-$α$ Relation of the Milky Way's Disrupted Dwarf Galaxies

The Milky Way's satellite galaxies ("surviving dwarfs") have been studied for decades as unique probes of chemical evolution in the low-mass regime. Here we extend such studies to the "disrupted dwarfs", whose debris constitutes the stellar halo. We present abundances ([Fe/H], [$α$/Fe]) and stellar masses for nine disrupted dwarfs with $M_{\star}\approx10^{6}-10^{9}M_{\odot}$ from the H3 Survey (Sagittarius, $Gaia$-Sausage-Enceladus, Helmi Streams, Sequoia, Wukong/LMS-1, Cetus, Thamnos, I'itoi, Orphan/Chenab). The surviving and disrupted dwarfs are chemically distinct: at fixed mass, the disrupted dwarfs are systematically metal-poor and $α$-enhanced. The disrupted dwarfs define a mass-metallicity relation (MZR) with a similar slope as the $z=0$ MZR followed by the surviving dwarfs, but offset to lower metallicities by $Δ$[Fe/H]$\approx0.3-0.4$ dex. Dwarfs with larger offsets from the $z=0$ MZR are more $α$-enhanced. In simulations as well as observations, galaxies with higher $Δ$[Fe/H] formed at higher redshifts -- exploiting this, we infer the disrupted dwarfs have typical star-formation truncation redshifts of $z_{\rm{trunc}}{\sim}1-2$. We compare the chemically inferred $z_{\rm{trunc}}$ with dynamically inferred accretion redshifts and find almost all dwarfs are quenched only after accretion. The differences between disrupted and surviving dwarfs are likely because the disrupted dwarfs assembled their mass rapidly, at higher redshifts, and within denser dark matter halos that formed closer to the Galaxy. Our results place novel archaeological constraints on low-mass galaxies inaccessible to direct high-$z$ studies: (i) the redshift evolution of the MZR along parallel tracks but offset to lower metallicities extends to $M_{\star}\approx10^{6}-10^{9}M_{\odot}$; (ii) galaxies at $z\approx2-3$ are $α$-enhanced with [$α$/Fe]$\approx0.4$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Recovering the star formation histories of recently-quenched galaxies: the impact of model and prior choices

Accurate models of the star formation histories (SFHs) of recently-quenched galaxies can provide constraints on when and how galaxies shut down their star formation. The recent development of "non-parametric" SFH models promises the flexibility required to make these measurements. However, model and prior choices significantly affect derived SFHs, particularly for post-starburst galaxies (PSBs) which have sharp changes in their recent SFH. In this paper, we create mock PSBs, then use the Prospector SED fitting software to test how well four different SFH models recover key properties. We find that a two-component parametric model performs well for our simple mock galaxies, but is sensitive to model mismatches. The fixed- and flexible-bin non-parametric models included in Prospector are able to rapidly quench a major burst of star formation, but systematically underestimate the post-burst age by up to 200 Myr. We develop a custom SFH model that allows for additional flexibility in the recent SFH. Our flexible non-parametric model is able to constrain post-burst ages with no significant offset and just ~90 Myr of scatter. Our results suggest that while standard non-parametric models are able to recover first-order quantities of the SFH (mass, SFR, average age), accurately recovering higher-order quantities (burst fraction, quenching time) requires careful consideration of model flexibility. These mock recovery tests are a critical part of future SFH studies. Finally, we show that our new, public SFH model is able to accurately recover the properties of mock star-forming and quiescent galaxies and is suitable for broader use in the SED fitting community.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Schrodinger's Galaxy Candidate: Puzzlingly Luminous at $z\approx17$, or Dusty/Quenched at $z\approx5$?

$JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5μ$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $Λ$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $<2.5$" have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749's photometry. The emission lines at $z\approx5$ conspire to boost the $>2μ$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($Δz\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

SQuIGGLE: Studying Quenching in Intermediate-z Galaxies -- Gas, AnguLar Momentum, and Evolution

We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.5<z<~0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be useful analogs for high-redshift quenching galaxies, but low enough redshift that multi-wavelength follow-up observations are feasible with modest telescope investments. We use the Prospector code to infer the stellar population properties and non-parametric star formation histories of all galaxies in the sample. We find that SQuIGGLE galaxies are both very massive (M* ~ 10^11.25 Msun) and quenched, with inferred star formation rates <~1Msun/yr, more than an order of magnitude below the star-forming main sequence. The best-fit star formation histories confirm that these galaxies recently quenched a major burst of star formation: >75% of SQuIGGLE galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200Myr before observation. We find that SQuIGGLE galaxies are on average younger and more burst-dominated than most other z<~1 post-starburst samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full SQuIGGLE survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, AGN incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Diverse Molecular Gas Content of Massive Galaxies Undergoing Quenching at z~1

We present a detailed study of the molecular gas content and stellar population properties of three massive galaxies at 1 < z < 1.3 that are in different stages of quenching. The galaxies were selected to have a quiescent optical/near-infrared spectral energy distribution and a relatively bright emission at 24 micron, and show remarkably diverse properties. CO emission from each of the three galaxies is detected in deep NOEMA observations, allowing us to derive molecular gas fractions Mgas/Mstar of 13-23%. We also reconstruct the star formation histories by fitting models to the observed photometry and optical spectroscopy, finding evidence for recent rejuvenation in one object, slow quenching in another, and rapid quenching in the third system. To better constrain the quenching mechanism we explore the depletion times for our sample and other similar samples at z~0.7 from the literature. We find that the depletion times are highly dependent on the method adopted to measure the star formation rate: using the UV+IR luminosity we obtain depletion times about 6 times shorter than those derived using dust-corrected [OII] emission. When adopting the star formation rates from spectral fitting, which are arguably more robust, we find that recently quenched galaxies and star-forming galaxies have similar depletion times, while older quiescent systems have longer depletion times. These results offer new, important constraints for physical models of galaxy quenching.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Mass of the Milky Way from the H3 Survey

The mass of the Milky Way is a critical quantity which, despite decades of research, remains uncertain within a factor of two. Until recently, most studies have used dynamical tracers in the inner regions of the halo, relying on extrapolations to estimate the mass of the Milky Way. In this paper, we extend the hierarchical Bayesian model applied in Eadie & Jurić (2019) to study the mass distribution of the Milky Way halo; the new model allows for the use of all available 6D phase-space measurements. We use kinematic data of halo stars out to $142~{\rm kpc}$, obtained from the H3 Survey and $\textit{Gaia}$ EDR3, to infer the mass of the Galaxy. Inference is carried out with the No-U-Turn sampler, a fast and scalable extension of Hamiltonian Monte Carlo. We report a median mass enclosed within $100~{\rm kpc}$ of $\rm M(<100 \; kpc) = 0.69_{-0.04}^{+0.05} \times 10^{12} \; M_\odot$ (68% Bayesian credible interval), or a virial mass of $\rm M_{200} = M(<216.2_{-7.5}^{+7.5} \; kpc) = 1.08_{-0.11}^{+0.12} \times 10^{12} \; M_\odot$, in good agreement with other recent estimates. We analyze our results using posterior predictive checks and find limitations in the model's ability to describe the data. In particular, we find sensitivity with respect to substructure in the halo, which limits the precision of our mass estimates to $\sim 15\%$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Diffuse Metal-Poor Component of the Sagittarius Stream Revealed by the H3 Survey

The tidal disruption of the Sagittarius dwarf galaxy has generated a spectacular stream of stars wrapping around the entire Galaxy. We use data from $Gaia$ and the H3 Stellar Spectroscopic Survey to identify 823 high-quality Sagittarius members based on their angular momenta. The H3 Survey is largely unbiased in metallicity, and so our sample of Sagittarius members is similarly unbiased. Stream stars span a wide range in [Fe/H] from $-0.2$ to $\approx -3.0$, with a mean overall metallicity of $\langle$[Fe/H]$\rangle=-0.99$. We identify a strong metallicity-dependence to the kinematics of the stream members. At [Fe/H]$\gt -0.8$ nearly all members belong to the well-known cold ($σ_v \lt 20$ km/s) leading and trailing arms. At intermediate metallicities ($-1.9 \lt$[Fe/H]$\lt -0.8$) a significant population (24$\%$) emerges of stars that are kinematically offset from the cold arms. These stars also appear to have hotter kinematics. At the lowest metallicities ([Fe/H]$\lesssim-2$), the majority of stars (69$\%$) belong to this kinematically-offset diffuse population. Comparison to simulations suggests that the diffuse component was stripped from the Sagittarius progenitor at earlier epochs, and therefore resided at larger radius on average, compared to the colder metal-rich component. We speculate that this kinematically diffuse, low metallicity, population is the stellar halo of the Sagittarius progenitor system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A New Census of the 0.2 < z < 3.0 Universe, Part I: The Stellar Mass Function

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Brackett-$γ$ as a Gold-standard Test of Star Formation Rates Derived from SED Fitting

Using a local reference sample of 21 galaxies, we compare observations of the $λ$2.16 $μ$m Brackett-$γ$ (Br$γ$) hydrogen recombination line with predictions from the Prospector Bayesian inference framework, which was used to fit the broadband photometry of these systems. This is a clean test of the spectral-energy-distribution-derived star formation rates (SFRs), as dust is expected to be optically thin at this wavelength in nearly all galaxies; thus, the internal conversion of SFR to predicted line luminosity does not depend strongly on the adopted dust model and posterior dust parameters, as is the case for shorter-wavelength lines such as H$α$. We find that Prospector predicts Br$γ$ luminosities and equivalent widths with small offsets ($\sim$0.05 dex), and scatter ($\sim$0.2 dex), consistent with measurement uncertainties, though we caution that the derived offset is dependent on the choice of stellar isochrones. We demonstrate that even when the Prospector-derived dust attenuation does not well describe, e.g., H$α$ line properties or observed reddening between H$α$ and Br$γ$, the underlying SFRs are accurate, as verified by the dust-free Br$γ$ comparison. Finally, we discuss in what ways Br$γ$ might be able to help constrain model parameters when treated as an input to the model, and comment on its potential as an accurate monochromatic SFR indicator in the era of JWST multiobject near-IR spectroscopy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High-resolution spectroscopy of the GD-1 stellar stream localizes the perturber near the orbital plane of Sagittarius

The $100^\circ$-long thin stellar stream in the Milky Way halo, GD-1, has an ensemble of features that may be due to dynamical interactions. Using high-resolution MMT/Hectochelle spectroscopy we show that a spur of GD-1-like stars outside of the main stream are kinematically and chemically consistent with the main stream. In the spur, as in the main stream, GD-1 has a low intrinsic radial velocity dispersion, $σ_{V_r}\lesssim1\,\rm km\,s^{-1}$, is metal-poor, $\rm [Fe/H]\approx-2.3$, with little $\rm [Fe/H]$ spread and some variation in $\rm [α/Fe]$ abundances, which point to a common globular cluster progenitor. At a fixed location along the stream, the median radial velocity offset between the spur and the main stream is smaller than $0.5\,\rm km\,s^{-1}$, comparable to the measurement uncertainty. A flyby of a massive, compact object can change orbits of stars in a stellar stream and produce features like the spur observed in GD-1. In this scenario, the radial velocity of the GD-1 spur relative to the stream constrains the orbit of the perturber and its current on-sky position to $\approx5,000\,\rm deg^2$. The family of acceptable perturber orbits overlaps the stellar and dark-matter debris of the Sagittarius dwarf galaxy in present-day position and velocity. This suggests that GD-1 may have been perturbed by a globular cluster or an extremely compact dark-matter subhalo formerly associated with Sagittarius.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

MINESweeper: Spectrophotometric Modeling of Stars in the Gaia Era

We present MINESweeper, a tool to measure stellar parameters by jointly fitting observed spectra and broadband photometry to model isochrones and spectral libraries. This approach enables the measurement of spectrophotometric distances, in addition to stellar parameters such as Teff, log(g), [Fe/H], [a/Fe], and radial velocity. MINESweeper employs a Bayesian framework and can easily incorporate a variety of priors, including Gaia parallaxes. Mock data are fit in order to demonstrate how the precision of derived parameters depends on evolutionary phase and SNR. We then fit a selection of data in order to validate the model outputs. Fits to a variety of benchmark stars including Procyon, Arcturus, and the Sun result in derived stellar parameters that are in good agreement with the literature. We then fit combined spectra and photometry of stars in the open and globular clusters M92, M13, M3, M107, M71, and M67. Derived distances, [Fe/H], [a/Fe], and log(g)-Teff, relations are in overall good agreement with literature values, although there are trends between metallicity and log(g), within clusters that point to systematic uncertainties at the ~0.1 dex level. Finally, we fit a large sample of stars from the H3 Spectroscopic Survey in which high quality Gaia parallaxes are also available. These stars are fit without the Gaia parallaxes so that the geometric parallaxes can serve as an independent test of the spectrophotometric distances. Comparison between the two reveals good agreement within their formal uncertainties after accounting for the Gaia zero point uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Powderday: Dust Radiative Transfer for Galaxy Simulations

We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Timing the Early Assembly of the Milky Way with the H3 Survey

The archaeological record of stars in the Milky Way opens a uniquely detailed window into the early formation and assembly of galaxies. Here we use 11,000 main-sequence turn-off stars with well-measured ages, [Fe/H], [$α$/Fe], and orbits from the H3 Survey and Gaia to time the major events in the early Galaxy. Located beyond the Galactic plane, $1\lesssim |Z|/\rm kpc \lesssim4$, this sample contains three chemically distinct groups: a low metallicity population, and low-$α$ and high-$α$ groups at higher metallicity. The age and orbit distributions of these populations show that: 1) the high-$α$ group, which includes both disk stars and the in-situ halo, has a star-formation history independent of eccentricity that abruptly truncated $8.3\pm0.1$ Gyr ago ($z\simeq1$); 2) the low metallicity population, which we identify as the accreted stellar halo, is on eccentric orbits and its star formation truncated $10.2.^{+0.2}_{-0.1}$ Gyr ago ($z\simeq2$); 3) the low-$α$ population is primarily on low eccentricity orbits and the bulk of its stars formed less than 8 Gyr ago. These results suggest a scenario in which the Milky Way accreted a satellite galaxy at $z\approx2$ that merged with the early disk by $z\approx1$. This merger truncated star formation in the early high-$α$ disk and perturbed a fraction of that disk onto halo-like orbits. The merger enabled the formation of a chemically distinct, low-$α$ disk at $z\lesssim1$. The lack of any stars on halo-like orbits at younger ages indicates that this event was the last significant disturbance to the Milky Way disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Towards Studying Hierarchical Assembly in Real Time: A Milky Way Progenitor Galaxy at z = 2.36 under the Microscope

We use Hubble Space Telescope (HST) imaging and near-infrared spectroscopy from Keck/MOSFIRE to study the sub-structure around the progenitor of a Milky Way-mass galaxy in the Hubble Frontier Fields (HFF). Specifically, we study an $r_e = 40^{+70}_{-30}$pc, $M_{\star} \sim 10^{8.2} M_{\odot}$ rest-frame ultra-violet luminous "clump" at a projected distance of $\sim$100~pc from a $M_{\star} \sim 10^{9.8}$M$_{\odot}$ galaxy at $z = 2.36$ with a magnification $μ= 5.21$. We measure the star formation history of the clump and galaxy by jointly modeling the broadband spectral energy distribution from HST photometry and H$α$ from MOSFIRE spectroscopy. Given our inferred properties (e.g., mass, metallicity, dust) of the clump and galaxy, we explore scenarios in which the clump formed \emph{in-situ} (e.g., a star forming complex) or \emph{ex-situ} (e.g., a dwarf galaxy being accreted). If it formed \emph{in-situ}, we conclude that the clump is likely a single entity as opposed to a aggregation of smaller star clusters, making it one of the most dense star clusters cataloged. If it formed \emph{ex-situ}, then we are witnessing an accretion event with a 1:40 stellar mass ratio. However, our data alone are not informative enough to distinguish between \emph{in-situ} and \emph{ex-situ} scenarios to a high level of significance. We posit that the addition of high-fidelity metallicity information, such as [OIII]4363Å, which can be detected at modest S/N with only a few hours of JWST/NIRSpec time, may be a powerful discriminant. We suggest that studying larger samples of moderately lensed sub-structures across cosmic time can provide unique insight into the hierarchical formation of galaxies like the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Lower Limit on the Mass of Our Galaxy from the H3 Survey

The timing argument provides a lower limit on the mass of the Milky Way. We find, using a sample of 32 stars at $R > 60$ kpc drawn from the H3 Spectroscopic Survey and mock catalogs created from published numerical simulations, that M$_{200} > 0.91\times 10^{12}$ M$_\odot$ with 90% confidence. We recommend using this limit to refine the allowed prior mass range in more complex and sophisticated statistical treatments of Milky Way dynamics. The use of such a prior would have significantly reduced many previously published uncertainty ranges. Our analysis suggests that the most likely value of M$_{200}$ is $\sim 1.4 \times 10^{12}$ M$_\odot$, but establishing this as the Milky Way mass requires a larger sample of outer halo stars and a more complete analysis of the inner halo stars in H3. The imminent growth in the sample of outer halo stars due to ongoing and planned surveys will make this possible.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Resolving the Metallicity Distribution of the Stellar Halo with the H3 Survey

The Galactic stellar halo is predicted to have formed at least partially from the tidal disruption of accreted dwarf galaxies. This assembly history should be detectable in the orbital and chemical properties of stars. The H3 Survey is obtaining spectra for 200,000 stars, and, when combined with Gaia data, is providing detailed orbital and chemical properties of Galactic halo stars. Unlike previous surveys of the halo, the H3 target selection is based solely on magnitude and Gaia parallax; the survey therefore provides a nearly unbiased view of the entire stellar halo at high latitudes. In this paper we present the distribution of stellar metallicities as a function of Galactocentric distance and orbital properties for a sample of 4232 kinematically-selected halo giants to 100 kpc. The stellar halo is relatively metal-rich, [Fe/H]=-1.2, and there is no discernable metallicity gradient over the range $6<R_{\rm gal}<100$ kpc. However, the halo metallicity distribution is highly structured including distinct metal-rich and metal-poor components at $R_{\rm gal}<10$ kpc and $R_{\rm gal}>30$ kpc, respectively. Metal-poor stars with [Fe/H]$<-2$ are a small population of the halo at all distances and orbital categories. We associate the "in-situ" stellar halo with stars displaying thick-disk chemistry on halo-like orbits; such stars are confined to $|z|<10$ kpc. The majority of the stellar halo is resolved into discrete features in orbital-chemical space, suggesting that the bulk of the stellar halo formed from the accretion and tidal disruption of dwarf galaxies. (ABRIDGED)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models

This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ($5 \leq \rm \log(Age)\;[yr] \leq 10.3$), masses ($0.1 \leq M/M_{\odot} \leq 300$), and metallicities ($-2.0 \leq \rm [Z/H] \leq 0.5$). The models are self-consistently and continuously evolved from the pre-main sequence to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the pre-main sequence to the end of core helium burning for $-4.0 \leq \rm [Z/H] < -2.0$. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.

0 citations0 reviewsNo code linkedOpen access
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]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Local Group as a time machine: studying the high-redshift Universe with nearby galaxies

We infer the UV luminosities of Local Group galaxies at early cosmic times ($z \sim 2$ and $z \sim 7$) by combining stellar population synthesis modeling with star formation histories derived from deep color-magnitude diagrams constructed from Hubble Space Telescope (HST) observations. Our analysis provides a basis for understanding high-$z$ galaxies - including those that may be unobservable even with the James Webb Space Telescope (JWST) - in the context of familiar, well-studied objects in the very low-$z$ Universe. We find that, at the epoch of reionization, all Local Group dwarfs were less luminous than the faintest galaxies detectable in deep HST observations of blank fields. We predict that JWST will observe $z \sim 7$ progenitors of galaxies similar to the Large Magellanic Cloud today; however, the HST Frontier Fields initiative may already be observing such galaxies, highlighting the power of gravitational lensing. Consensus reionization models require an extrapolation of the observed blank-field luminosity function at $z \approx 7$ by at least two orders of magnitude in order to maintain reionization. This scenario requires the progenitors of the Fornax and Sagittarius dwarf spheroidal galaxies to be contributors to the ionizing background at $z \sim 7$. Combined with numerical simulations, our results argue for a break in the UV luminosity function from a faint-end slope of $α\sim -2$ at $M_{\rm UV} < -13$ to $α\sim -1.2$ at lower luminosities. Applied to photometric samples at lower redshifts, our analysis suggests that HST observations in lensing fields at $z \sim 2$ are capable of probing galaxies with luminosities comparable to the expected progenitor of Fornax.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A pilot study using deep infrared imaging to constrain the star formation history of the XUV stellar populations in NGC 4625

In a LCDM universe, disk galaxies' outer regions are the last to form. Characterizing their contents is critical for understanding the ongoing process of disk formation, but observing outer disk stellar populations is challenging due to their low surface brightness. We present extremely deep 3.6 micron observations (Spitzer IRAC) of NGC 4625, a galaxy known for its radially extended ultraviolet emitting stellar population. We combine the new imaging with archival UV imaging from the GALEX mission to derive multiwavelength radial profiles for NGC 4625 and compare them to stellar populations models. The colors can be explained by the young stellar population that is responsible for the UV emission and indicate that the current star formation rates in the outermost disk are recent. Extended star formation in NGC 4625 may have been initiated by an interaction with neighboring galaxies NGC 4618 and NGC 4625a, supporting speculation that minor interactions are a common trigger for outer disk star formation and late stage disk growth.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Empirical ugri-UBVRc Transformations for Galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spitzer Local Volume Legacy (LVL) SEDs and Physical Properties

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Mice at play in the CALIFA survey: A case study of a gas-rich major merger between first passage and coalescence

We present optical integral field spectroscopy (IFS) observations of the Mice, a major merger between two massive (>10^11Msol) gas-rich spirals NGC4676A and B, observed between first passage and final coalescence. The spectra provide stellar and gas kinematics, ionised gas properties and stellar population diagnostics, over the full optical extent of both galaxies. The Mice provide a perfect case study highlighting the importance of IFS data for improving our understanding of local galaxies. The impact of first passage on the kinematics of the stars and gas has been significant, with strong bars likely induced in both galaxies. The barred spiral NGC4676B exhibits a strong twist in both its stellar and ionised gas disk. On the other hand, the impact of the merger on the stellar populations has been minimal thus far: star formation induced by the recent close passage has not contributed significantly to the global star formation rate or stellar mass of the galaxies. Both galaxies show bicones of high ionisation gas extending along their minor axes. In NGC4676A the high gas velocity dispersion and Seyfert-like line ratios at large scaleheight indicate a powerful outflow. Fast shocks extend to ~6.6kpc above the disk plane. The measured ram pressure and mass outflow rate (~8-20Msol/yr) are similar to superwinds from local ULIRGs, although NGC4676A has only a moderate infrared luminosity of 3x10^10Lsol. Energy beyond that provided by the mechanical energy of the starburst appears to be required to drive the outflow. We compare the observations to mock kinematic and stellar population maps from a merger simulation. The models show little enhancement in star formation during and following first passage, in agreement with the observations. We highlight areas where IFS data could help further constrain the models.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Spitzer Local Volume Legacy (LVL) Global Optical Photometry

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Very Faint End of the UV Luminosity Function over Cosmic Time: Constraints from the Local Group Fossil Record

We present a new technique to estimate the evolution of the very faint end of the UV luminosity function (LF) out to $z\sim5$. Measured star formation histories (SFHs) from the fossil record of Local Group galaxies are used to reconstruct the LF down to M$_{UV}\sim-5$ at $z\sim5$ and M$_{UV}\sim-1.5$ at $z<1$. Such faint limits are well beyond the current observational limits and are likely to remain beyond the limits of next generation facilities. The reconstructed LFs, when combined with direct measurements of the LFs at higher luminosity, are well-fit by a standard Schechter function with no evidence of a break to the faintest limits probed by this technique. The derived faint end slope, $α$, steepens from $\approx-1.2$ at $z<1$ to $\approx-1.6$ at $4<z<5$. We test the effects of burstiness in the SFHs and find the recovered LFs to be only modestly affected. Incompleteness corrections for the faintest Local Group galaxies and the (unlikely) possibility of significant luminosity-dependent destruction of dwarf galaxies between high redshift and the present epoch are important uncertainties. These and other uncertainties can be mitigated with more detailed modeling and future observations. The reconstructed faint end LF from the fossil record can therefore be a powerful and complementary probe of the high redshift faint galaxies believed to play a key role in the reionization of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Encoding the infrared excess (IRX) in the NUVrK color diagram for star-forming galaxies

We present an empirical method of assessing the star formation rate (SFR) of star-forming galaxies based on their locations in the rest-frame color-color diagram (NUV-r) vs (r-K). By using the Spitzer 24 micron sample in the COSMOS field (~16400 galaxies with 0.2 < z < 1.3) and a local GALEX-SDSS-SWIRE sample (~700 galaxies with z < 0.2), we show that the mean infrared excess <IRX>= < L_IR / L_UV > can be described by a single vector, NRK, that combines the two colors. The calibration between <IRX> and NRK allows us to recover the IR luminosity, L_IR, with an accuracy of ~0.21 dex for the COSMOS sample and ~0.27 dex for the local one. The SFRs derived with this method agree with the ones based on the observed (UV+IR) luminosities and on the spectral energy distribution fitting for the vast majority (~85 %) of the star-forming population. Thanks to a library of model galaxy SEDs with realistic prescriptions for the star formation history, we show that we need to include a two-component dust model (i.e., birth clouds and diffuse ISM) and a full distribution of galaxy inclinations in order to reproduce the behavior of the <IRX> stripes in the NUVrK diagram. In conclusion, the NRK method, based only on rest-frame UV and optical colors available in most of the extragalactic fields, offers a simple alternative of assessing the SFR of star-forming galaxies in the absence of far-IR or spectral diagnostic observations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dust-Corrected Star Formation Rates of Galaxies. II. Combinations of Ultraviolet and Infrared Tracers

We present new calibrations of far-ultraviolet (FUV) attenuation as derived from the total infrared to FUV luminosity ratio (IRX) and the FUV-NUV color. We find that the IRX-corrected FUV luminosities are tightly and linearly correlated with the attenuation-corrected Hα luminosities (as measured from the Balmer decrement), with a rms scatter of $\pm 0.09$ dex. The ratios of these attenuation-corrected FUV to Hα luminosities are consistent with evolutionary synthesis model predictions, assuming a constant star formation rate over 100 Myr, solar metallicity and either a Salpeter or a Kroupa IMF with lower and upper mass limits of 0.1 and 100\msun. The IRX-corrected FUV to Balmer-corrected Hα luminosity ratios do not show any trend with other galactic properties over the ranges covered by our sample objects. In contrast, FUV attenuation derived from the FUV-NUV color (UV spectral slope) show much larger random and systematic uncertainties. When compared to either Balmer-corrected Hα luminosities or IRX-corrected FUV luminosities the color-corrected FUV luminosities show $\sim 2.5$ times larger rms scatter, and systematic nonlinear deviations as functions of luminosity and other parameters. Linear combinations of 25um and 1.4GHz radio continuum luminosities with the observed FUV luminosities are also well correlated with the Balmer-corrected Hα luminosities. These results provide useful prescriptions for deriving attenuation-corrected star formation rates of galaxies based on linear combinations of UV and IR or radio luminosities, which are presented in convenient tabular form. Comparisons of our calibrations with attenuation corrections in the literature and with dust attenuation laws are also made.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

SSGSS: The Spitzer-SDSS-GALEX Spectroscopic Survey

The Spitzer-SDSS-GALEX Spectroscopic Survey (SSGSS) provides a new sample of 101 star-forming galaxies at z < 0.2 with unprecedented multi-wavelength coverage. New mid- to far-infrared spectroscopy from the Spitzer Space Telescope is added to a rich suite of previous imaging and spectroscopy, including ROSAT, Galaxy Evolution Explorer, Sloan Digital Sky Survey, Two Micron All Sky Survey, and Spitzer/SWIRE. Sample selection ensures an even coverage of the full range of normal galaxy properties, spanning two orders of magnitude in stellar mass, color, and dust attenuation. In this paper we present the SSGSS data set, describe the science drivers, and detail the sample selection, observations, data reduction, and quality assessment. Also in this paper, we compare the shape of the thermal continuum and the degree of silicate absorption of these typical, star-forming galaxies to those of starburst galaxies. We investigate the link between star formation rate, infrared luminosity, and total polycyclic aromatic hydrocarbon luminosity, with a view to calibrating the latter for spectral energy distribution models in photometric samples and at high redshift. Last, we take advantage of the 5-40 micron spectroscopic and far-infrared photometric coverage of this sample to perform detailed fitting of the Draine et al. dust models, and investigate the link between dust mass and star formation history and active galactic nucleus properties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Emission by Dust and Stars of Nearby Galaxies in the Herschel KINGFISH Survey

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large scale outflows from z ~ 0.7 starburst galaxies identified via ultra-strong MgII quasar absorption lines

(Abridged) Star formation-driven outflows are a critical phenomenon in theoretical treatments of galaxy evolution, despite the limited ability of observations to trace them across cosmological timescales. If the strongest MgII absorption-line systems detected in the spectra of background quasars arise in such outflows, "ultra-strong" MgII (USMgII) absorbers would identify significant numbers of galactic winds over a huge baseline in cosmic time, in a manner independent of the luminous properties of the galaxy. To this end, we present the first detailed imaging and spectroscopic study of the fields of two USMgII absorber systems culled from a statistical absorber catalog, with the goal of understanding the physical processes leading to the large velocity spreads that define such systems. Each field contains two bright emission-line galaxies at similar redshift (dv < 300 km/s) to that of the absorption. Lower-limits on their instantaneous star formation rates (SFR) from the observed OII and Hb line fluxes, and stellar masses from spectral template fitting indicate specific SFRs among the highest for their masses at z~0.7. Additionally, their 4000A break and Balmer absorption strengths imply they have undergone recent (~0.01 - 1 Gyr) starbursts. The concomitant presence of two rare phenomena - starbursts and USMgII absorbers - strongly implies a causal connection. We consider these data and USMgII absorbers in general in the context of various popular models, and conclude that galactic outflows are generally necessary to account for the velocity extent of the absorption. We favour starburst driven outflows over tidally-stripped gas from a major interaction which triggered the starburst as the energy source for the majority of systems. Finally, we discuss the implications of these results and speculate on the overall contribution of such systems to the global SFR density at z~0.7.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The GALEX Arecibo SDSS Survey II: The Star Formation Efficiency of Massive Galaxies

We use measurements of the HI content, stellar mass and star formation rates in ~190 massive galaxies with stellar masses greater than 10^10 Msun, obtained from the Galex Arecibo SDSS Survey (GASS) described in Paper I (Catinella et al. 2010) to explore the global scaling relations associated with the bin-averaged ratio of the star formation rate over the HI mass, which we call the HI-based star formation efficiency (SFE). Unlike the mean specific star formation rate, which decreases with stellar mass and stellar mass surface density, the star formation efficiency remains relatively constant across the sample with a value close to SFE = 10^-9.5 yr^-1 (or an equivalent gas consumption timescale of ~3 Gyr). Specifically, we find little variation in SFE with stellar mass, stellar mass surface density, NUV-r color and concentration. We interpret these results as an indication that external processes or feedback mechanisms that control the gas supply are important for regulating star formation in massive galaxies. An investigation into the detailed distribution of SFEs reveals that approximately 5% of the sample shows high efficiencies with SFE > 10^-9 yr^-1, and we suggest that this is very likely due to a deficiency of cold gas rather than an excess star formation rate. Conversely, we also find a similar fraction of galaxies that appear to be gas-rich for their given specific star-formation rate, although these galaxies show both a higher than average gas fraction and lower than average specific star formation rate. Both of these populations are plausible candidates for "transition" galaxies, showing potential for a change (either decrease or increase) in their specific star formation rate in the near future. We also find that 36+/-5% of the total HI mass density and 47+/-5% of the total SFR density is found in galaxies with stellar mass greater than 10^10 Msun. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Polycyclic Aromatic Hydrocarbons in Galaxies at z~0.1: the Effect of Star Formation and AGN

We present the analysis of the Polycyclic Aromatic Hydrocarbon (PAH) spectra of a sample of 92 typical star forming galaxies at 0.03 < z < 0.2 observed with the Spitzer IRS. We compare the relative strengths of PAH emission features with SDSS optical diagnostics to probe the relationship between PAH grain properties and star formation and AGN activity. Short-to-long wavelength PAH ratios, and in particular the 7.7-to-11.3 micron feature ratio, are strongly correlated with the star formation diagnostics D_n(4000) and H-alpha equivalent width, increasing with younger stellar populations. This ratio also shows a significant difference between active and non-active galaxies, with the active galaxies exhibiting weaker 7.7 micron emission. A hard radiation field as measured by [OIII]/H-beta and [NeIII]_15.6/[NeII]_12.8 affects PAH ratios differently depending on whether this field results from starburst activity or an AGN. Our results are consistent with a picture in which larger PAH molecules grow more efficiently in richer media and in which smaller PAH molecules are preferentially destroyed by AGN.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The GALEX Arecibo SDSS Survey. I. Gas Fraction Scaling Relations of Massive Galaxies and First Data Release

We introduce the GALEX Arecibo SDSS Survey (GASS), an on-going large program that is gathering high quality HI-line spectra using the Arecibo radio telescope for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05, selected from the SDSS spectroscopic and GALEX imaging surveys. The galaxies are observed until detected or until a low gas mass fraction limit (1.5-5%) is reached. This paper presents the first Data Release, consisting of ~20% of the final GASS sample. We use this data set to explore the main scaling relations of HI gas fraction with galaxy structure and NUV-r colour. A large fraction (~60%) of the galaxies in our sample are detected in HI. We find that the atomic gas fraction decreases strongly with stellar mass, stellar surface mass density and NUV-r colour, but is only weakly correlated with galaxy bulge-to-disk ratio (as measured by the concentration index of the r-band light). We also find that the fraction of galaxies with significant (more than a few percent) HI decreases sharply above a characteristic stellar surface mass density of 10^8.5 Msun kpc^-2. The fraction of gas-rich galaxies decreases much more smoothly with stellar mass. One of the key goals of GASS is to identify and quantify the incidence of galaxies that are transitioning between the blue, star-forming cloud and the red sequence of passively-evolving galaxies. Likely transition candidates can be identified as outliers from the mean scaling relations between gas fraction and other galaxy properties. [abridged]

0 citations0 reviewsNo code linkedOpen access