Star Formation History Research Articles

Likelihood of white dwarf binaries to dominate the astrophysical gravitational wave background in the mHz band

Context. The astrophysical gravitational wave background (AGWB) is a collective signal of astrophysical gravitational wave sources dominated by compact binaries. One key science goal of current and future gravitational wave detectors is to obtain its measurement. Aims. We aim to determine the population of compact binaries dominating the AGWB in the mHz band. We revisit and update an earlier work by Farmer & Phinney (2003, MNRAS, 346, 1197) to model the astrophysical gravitational wave background sourced by extragalactic white dwarf binaries in the mHz frequency band. Methods. We calculated the signal using a single-metallicity model for the white dwarf population in the Universe using the global star formation history. Results. We estimate the white dwarf AGWB amplitude to be ∼60% higher than the earlier estimate. We also find that the overall shape of the white dwarf AGWB shows a good fit with a broken power law combined with an exponential cut-off. Conclusions. We compare our results to present-day best estimates for the background due to black hole and neutron star binaries, finding that the white dwarf component is likely to dominate in the mHz band. We provide an order-of-magnitude estimate that explains this hierarchy and we comment on the implications for future missions that aim to detect the AGWB. We also note that the black hole AGWB may only be detectable at high frequency. We outline several improvements that can be made to our estimate, however, these points are unlikely to change our main conclusion, which posits that the white dwarf AGWB dominates the mHz band.

Constraints on cosmic star formation history via a new modeling of the radio luminosity function of star-forming galaxies

Context. Radio wavelengths offer a unique possibility to trace the total star-formation rate (SFR) in galaxies, both obscured and unobscured. To probe the dust-unbiased star-formation history, an accurate measurement of the radio luminosity function (LF) for star-forming galaxies (SFGs) is crucial. Aims. We make use of an SFG sample (5900 sources) from the Very Large Array (VLA) COSMOS 3 GHz data to perform a new modeling of the radio LF. By integrating the analytical LF, we aim to calculate the history of the cosmic SFR density (SFRD) from z ∼ 5 onwards. Methods. For the first time, we use both models of the pure luminosity evolution (PLE) and joint luminosity+density evolution (LADE) to fit the LFs directly to the radio data using a full maximum-likelihood analysis, considering the sample completeness correction. We also incorporate updated observations of local radio LFs and radio source counts into the fitting process to obtain additional constraints. Results. We find that the PLE model cannot be used to describe the evolution of the radio LF at high redshift (z > 2). By construct, our LADE models can successfully fit a large amount of data on radio LFs and source counts of SFGs from recent observations. The Akaike information criterion (AIC) also demonstrates that the LADE model is superior to the PLE model. We therefore conclude that density evolution is genuinely indispensable in modeling the evolution of SFG radio LFs. Our SFRD curve shows a good fit to the SFRD points derived by previous radio estimates. In view of the fact that our radio LFs are not biased, as opposed those of previous studies performed by fitting the 1/Vmax LF points, our SFRD results should be an improvement on these previous estimates. Below z ∼ 1.5, our SFRD matches a published multiwavelength compilation, while our SFRD turns over at a slightly higher redshift (2 < z < 2.5) and falls more rapidly out to high redshift.

Quantifying the Effects of Known Unknowns on Inferred High-redshift Galaxy Properties: Burstiness, IMF, and Nebular Physics

The era of the James Webb Space Telescope ushers stellar population models into uncharted territories, particularly at the high-redshift frontier. In a companion paper, we apply the Prospector Bayesian framework to jointly infer galaxy redshifts and stellar population properties from broadband photometry as part of the UNCOVER survey. Here we present a comprehensive error budget in spectral energy distribution (SED) modeling. Using a sample selected to have photometric redshifts higher than 9, we quantify the systematic shifts stemming from various model choices in inferred stellar mass, star formation rate (SFR), and age. These choices encompass different timescales for changes in the star formation history (SFH), nonuniversal stellar initial mass functions (IMF), and the inclusion of variable nebular abundances, gas density, and ionizing photon budget. We find that the IMF exerts the strongest influence on the inferred properties: the systematic uncertainties can be as much as 1 dex, 2–5 times larger than the formal reported uncertainties in mass and SFR, and importantly, exceed the scatter seen when using different SED fitting codes. Although the assumptions on the lower end of the IMF induce degeneracy, our findings suggest that a common practice in the literature of assessing uncertainties in SED-fitting processes by comparing multiple codes is substantively underestimating the true systematic uncertainty. Highly stochastic SFHs change the inferred SFH by much larger than the formal uncertainties, and introduce ∼0.8 dex systematics in SFR averaged over a short timescale and ∼0.3 dex systematics in average age. Finally, employing a flexible nebular emission model causes ∼0.2 dex systematic increase in mass and SFR, comparable to the formal uncertainty. This paper constitutes an initial step toward a complete uncertainty estimate in SED modeling.

A Massive Quiescent Galaxy in a Group Environment at z = 4.53

We report on the spectroscopic confirmation of a massive quiescent galaxy at z spec = 4.53 in the COSMOS field. The object was first identified as a galaxy with suppressed star formation at z phot ∼ 4.65 from the COSMOS2020 catalog. The follow-up spectroscopy with Keck/MOSFIRE in the K band reveals faint [O ii] emission and the Balmer break, indicative of evolved stellar populations. We fit the spectral energy distribution using photometry and a spectrum to infer physical properties. The obtained stellar mass is high (M * ∼ 1010.8 M ⊙) and the current star formation rate is more than 1 dex below that of main-sequence galaxies at z = 4.5. Its star formation history suggests that this galaxy experienced rapid quenching from z ∼ 5. The galaxy is among the youngest quiescent galaxies confirmed so far at z spec > 3 with z form ∼ 5.2 (200 Myr ago), which is the epoch when 50% of the total stellar mass was formed. A unique aspect of the galaxy is that it is in an extremely dense region; there are four massive star-forming galaxies at 4.4 < z phot < 4.7 located within 150 physical kpc from the galaxy. Interestingly, three of them have virial radii that strongly overlap with that of the central quiescent galaxy (∼70 kpc), suggesting that the overdensity region is likely the highest-redshift candidate of a dense group with a spectroscopically confirmed quiescent galaxy at the center. The group provides us with a unique opportunity to gain insights into the role of the group environment in quenching at z ∼ 5, which corresponds to the formation epoch of massive elliptical galaxies in the local Universe.

Constraining the stellar populations of ultra-diffuse galaxies in the MATLAS survey using spectral energy distribution fitting

ABSTRACT We use spectral energy distribution fitting to place constraints on the stellar populations of 59 ultra-diffuse galaxies (UDGs) in the low-to-moderate density fields of the MATLAS survey. We use the routine prospector, coupled with archival data in the optical from the Dark Energy Camera Legacy Survey, and near- and mid-infrared imaging from the Wide-field Infrared Survey Explorer, to recover the stellar masses, ages, metallicities, and star formation time-scales of the UDGs. We find that a subsample of the UDGs lies within the scatter of the mass–metallicity relation (MZR) for local classical dwarfs. However, another subsample is more metal-poor, being consistent with the evolving MZR at high redshift. We investigate UDG positioning trends in the mass–metallicity plane as a function of surface brightness, effective radius, axis ratio, local volume density, mass-weighted age, star formation time-scale, globular cluster (GC) counts, and GC specific frequency. We find that our sample of UDGs can be separated into two main classes: Class A: comprised of UDGs with lower stellar masses, prolonged star formation histories (SFHs), more elongated, inhabiting less dense environments, hosting fewer GCs, younger, consistent with the classical dwarf MZR, and fainter. Class B: UDGs with higher stellar masses, rapid SFHs, rounder, inhabiting the densest of our probed environments, hosting on average the most numerous GC systems, older, consistent with the high-redshift MZR (i.e. consistent with early-quenching), and brighter. The combination of these properties suggests that UDGs of Class A are consistent with a ‘puffed-up dwarf’ formation scenario, while UDGs of Class B seem to be better explained by ‘failed galaxy’ scenarios.

Hints of a disrupted binary dwarf galaxy in the Sagittarius stream

ABSTRACT In this work, we look for evidence of a non-unity mass ratio binary dwarf galaxy merger in the Sagittarius stream. Simulations of such a merger show that, upon merging with a host, particles from the less massive galaxy will often mostly be found in the extended stream and less so in the central remnant. Motivated by these simulations, we use the Apache Point Observatory Galactic Evolution Experiment Data Release 17 chemical data from approximately 1100 stars in both the Sagittarius remnant and stream to look for evidence of contamination from a second dwarf galaxy. We separate the Sagittarius data into its remnant and stream and compare the [Mg/Fe] content of the two populations. In particular, we select [Mg/Fe] to search for hints of unique star formation histories among our sample stars. Comparing the stream and remnant populations, we find regions that have distinct [Mg/Fe] distributions for fixed [Fe/H], in addition to distinct chemical tracks in [Mg/Fe]–[Fe/H] abundance space. We show that there are large regions of the tracks for which the probability of the two samples being drawn from the same distribution is very low (p &amp;lt; 0.05). Furthermore, we show that the two tracks can be fit with unique star formation histories using simple, one-zone galactic chemical evolution models. While more work must be done to discern whether the hypothesis presented here is true, our work hints at the possibility that Sagittarius may consist of two dwarf galaxy progenitors.

Effect of cosmic rays and ionizing radiation on observational ultraviolet plasma diagnostics in the circumgalactic medium

ABSTRACT The relevance of some galactic feedback mechanisms, in particular cosmic-ray (CR) feedback and the hydrogen ionizing radiation field, has been challenging to definitively describe in a galactic context, especially far outside the galaxy in the circumgalactic medium (CGM). Theoretical and observational uncertainties prevent conclusive interpretations of multiphase CGM properties derived from ultraviolet (UV) diagnostics. We conduct three-dimensional magnetohydrodynamic simulations of a section of a galactic disc with star formation and feedback, including radiative heating from stars, a UV background, and CR feedback. We utilize the temperature phases present in our simulations to generate Cloudy models to derive spatially and temporally varying synthetic UV diagnostics. We find that radiative effects without additional heating mechanisms are not able to produce synthetic diagnostics in the observed ranges. For low CR diffusivity $\kappa _{\rm {cr}}=10^{28} \rm {cm}^2 \rm {s}^{-1}$, CR streaming heating in the outflow helps our synthetic line ratios roughly match observed ranges by producing transitional temperature gas (T ∼ 105–106 K). High CR diffusivity $\kappa _{\rm {cr}}=10^{29} \rm {cm}^2 \rm {s}^{-1}$, with or without CR streaming heating, produced transitional temperature gas. The key parameter controlling the production of this gas phase remains unclear, as the different star formation history and outflow evolution itself influences these diagnostics. Our work demonstrates the use of UV plasma diagnostics to differentiate between galactic/circumgalactic feedback models.

A population of faint, old, and massive quiescent galaxies at 3

Here we present a sample of 12 massive quiescent galaxy candidates at z∼3-4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z\\sim 3-4$$\\end{document} observed with the James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec). These galaxies were pre-selected from the Hubble Space Telescope imaging and 10 of our sources were unable to be spectroscopically confirmed by ground based spectroscopy. By combining spectroscopic data from NIRSpec with multi-wavelength imaging data from the JWST Near Infrared Camera (NIRCam), we analyse their stellar populations and their formation histories. We find that all of our galaxies classify as quiescent based on the reconstruction of their star formation histories but show a variety of quenching timescales and ages. All our galaxies are massive (∼0.1-1.2×1011\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 0.1-1.2\ imes 10^{11}$$\\end{document} M⊙\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_\\odot$$\\end{document}), with masses comparable to massive galaxies in the local Universe. We find that the oldest galaxy in our sample formed ∼1.0×1011\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 1.0\ imes 10^{11}$$\\end{document} M⊙\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_\\odot$$\\end{document} of mass within the first few hundred million years of the Universe and has been quenched for more than a billion years by the time of observation at z∼3.2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z\\sim 3.2$$\\end{document} (∼2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 2$$\\end{document} billion years after the Big Bang). Our results point to very early formation of massive galaxies requiring a high conversion rate of baryons to stars in the early Universe.

On the Significance of the Thick Disks of Disk Galaxies

Thick disks are a prevalent feature observed in numerous disk galaxies, including our own Milky Way. Their significance has been reported to vary widely, ranging from a few percent to 100% of the disk mass, depending on the galaxy and the measurement method. We use the NewHorizon simulation, which has high spatial and stellar mass resolutions, to investigate the issue of the thick-disk mass fraction. We also use the NewHorizon2 simulation, which was run on the same initial conditions, but additionally traced nine chemical elements. Based on a sample of 27 massive disk galaxies with M * > 1010 M ⊙ in NewHorizon, the contribution of the thick disk was found to be 20% ± 11% in r-band luminosity or 35% ± 15% in mass to the overall galactic disk, which seems in agreement with observational data. The vertical profiles of 0, 22, and 5 galaxies are best fitted by 1, 2, or 3 sech2 components, respectively. The NewHorizon2 data show that the selection of thick-disk stars based on a single [α/Fe] cut is contaminated by stars of different kinematic properties, while missing the bulk of kinematically thick disk stars. Vertical luminosity profile fits recover the key properties of thick disks reasonably well. The majority of stars are born near the galactic midplane with high circularity and get heated with time via fluctuations in the force field. Depending on the star formation and merger histories, galaxies may naturally develop thick disks with significantly different properties.

The PHANGS-AstroSat Atlas of Nearby Star-forming Galaxies

We present the Physics at High Angular resolution in Nearby GalaxieS (PHANGS)-AstroSat atlas, which contains UV imaging of 31 nearby star-forming galaxies captured by the Ultraviolet Imaging Telescope on the AstroSat satellite. The atlas provides a homogeneous data set of far-UV and near-UV maps of galaxies within a distance of 22 Mpc and a median angular resolution of 1.″4 (corresponding to a physical scale between 25 and 160 pc). After subtracting a uniform UV background and accounting for Milky Way extinction, we compare our estimated flux densities to GALEX observations, finding good agreement. We find candidate extended UV disks around the galaxies NGC 6744 and IC 5332. We present the first statistical measurements of the clumping of the UV emission and compare it to the clumping of molecular gas traced with the Atacama Large Millimeter/submillimeter Array (ALMA). We find that bars and spiral arms exhibit the highest degree of clumping, and the molecular gas is even more clumped than the far-UV (FUV) emission in galaxies. We investigate the variation of the ratio of observed FUV to Hα in different galactic environments and kiloparsec-sized apertures. We report that ∼65% of the variation of the log10(FUV/Hα) can be described through a combination of dust attenuation with star formation history parameters. The PHANGS-AstroSat atlas enhances the multiwavelength coverage of our sample, offering a detailed perspective on star formation. When integrated with PHANGS data sets from ALMA, the Very Large Telescope-MUSE, the Hubble Space Telescope, and JWST, it develops our comprehensive understanding of attenuation curves and dust attenuation in star-forming galaxies.

Exploring the faintest end of mid-infrared luminosity functions up to z ≃ 5 with the JWST CEERS survey

ABSTRACT Mid-infrared (MIR) light from galaxies is sensitive to dust-obscured star formation activities because it traces the characteristic emission of dust heated by young, massive stars. By constructing the MIR luminosity functions (LFs), we are able to quantify the overall dusty star formation history and the evolution of galaxies over cosmic time. In this work, we report the first rest-frame MIR LFs at 7.7, 10, 12.8, 15, 18, and 21 μm as well as the total IR LF from the JWST Cosmic Evolution Early Release Science (CEERS) survey. We identify 506 galaxies at z = 0–5.1 in the CEERS survey that also have optical photometry from the Hubble Space Telescope. With the unprecedented sensitivity of the JWST, we probe the faintest end of the LFs at z = 0–1 down to L* ∼ 107L⊙, ∼2 orders of magnitude fainter than those from the previous generation of IR space telescopes. Our findings connect well with and continue the faint end of the MIR LFs from the deepest observations in past works. As a proxy of star formation history, we present the MIR-based luminosity density up to z ≃ 4.0, marking the first probe of the early Universe by JWST Mid-Infrared Instrument.

Beyond mirkwood: Enhancing SED Modeling with Conformal Predictions

Traditional spectral energy distribution (SED) fitting techniques face uncertainties due to assumptions in star formation histories and dust attenuation curves. We propose an advanced machine learning-based approach that enhances flexibility and uncertainty quantification in SED fitting. Unlike the fixed NGBoost model used in mirkwood, our approach allows for any scikit-learn-compatible model, including deterministic models. We incorporate conformalized quantile regression to convert point predictions into error bars, enhancing interpretability and reliability. Using CatBoost as the base predictor, we compare results with and without conformal prediction, demonstrating improved performance using metrics such as coverage and interval width. Our method offers a more versatile and accurate tool for deriving galaxy physical properties from observational data.

Resolving cosmic star formation histories of present-day bulges, discs, and spheroids with ProFuse

ABSTRACT We present the first look at star formation histories of galaxy components using ProFuse, a new technique to model the 2D distribution of light across multiple wavelengths using simultaneous spectral and spatial fitting of purely imaging data. We present a number of methods to classify galaxies structurally/morphologically, showing the similarities and discrepancies between these schemes. We show the variation in componentwise mass functions that can occur simply due to the use of a different classification method, which is most dramatic in separating bulges and spheroids. Rather than identifying the best-performing scheme, we use the spread of classifications to quantify uncertainty in our results. We study the cosmic star formation history (CSFH), forensically derived using ProFuse with a sample of ∼7000 galaxies from the Galaxy And Mass Assembly survey. Remarkably, the forensic CSFH recovered via both our method (ProFuse) and traditional Spectral Energy Distribution (SED) fitting (ProSpect) are not only exactly consistent with each other over the past 8 Gyr, but also with the in situ CSFH measured using ProSpect. Furthermore, we separate the CSFH by contributions from spheroids, bulges, and discs. While the vast majority (70 per cent) of present-day star formation takes place in the disc population, we show that 50 per cent of the stars that formed at cosmic noon (8–12 Gyr ago) now reside in spheroids, and present-day bulges are composed of stars that were primarily formed in the very early Universe, with half their stars already formed ∼12 Gyr ago.

Dust enrichment in the circum-galactic medium

ABSTRACT To understand the origin of dust in the circum-galactic medium (CGM), we develop a dust enrichment model. We describe each of the central galaxy and its CGM as a single zone, and consider the mass exchange between them through galactic inflows and outflows. We calculate the evolution of the gas, metal, and dust masses in the galaxy and the CGM. In the galaxy, we include stellar dust production and interstellar dust processing following our previous models. The dust in the galaxy is transported to the CGM via galactic outflows, and it is further processed by dust destruction (sputtering) in the CGM. We parametrize the time-scale or efficiency of each process and investigate the effect on the dust abundance in the CGM. We find that the resulting dust mass is sensitive to the dust destruction in the CGM, and the dust supply from galactic outflows, both of which directly regulate the dust abundance in the CGM. The inflow time-scale also affects the dust abundance in the CGM because it determines the gas mass evolution (thus, the star formation history) in the galaxy. The dust abundance in the CGM, however, is insensitive to stellar dust formation in the galaxy at later epochs because the dust production is dominated by dust growth in the interstellar medium. We also find that the resulting dust mass in the CGM is consistent with the value derived from a large sample of SDSS galaxies.

Strategies for Obtaining Robust Spectral Energy Distribution Fitting Parameters for Galaxies at z ∼ 1 and z ∼ 2 in the Absence of Infrared Data

Robust estimation of star formation rates (SFRs) at higher redshifts (z ≳ 1) using UV–optical–near-infrared (NIR) photometry is contingent on the ability of spectral energy distribution (SED) fitting to constrain the dust attenuation, stellar metallicity, and star formation history (SFH) simultaneously. IR-derived dust luminosities can help break the degeneracy between these parameters, but IR data are often not available. Here, we explore strategies for SED fitting at z ≳ 1 in the absence of IR data using a sample of log M * > 10.2 star-forming galaxies from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) for which 24 μm data are available. We adopt the total IR luminosity (L TIR) obtained from 24 μm as the “ground truth,” which allows us to assess how well it can be recovered (as L dust) from UV–optical–NIR SED fitting. We test a variety of dust attenuation models, stellar population synthesis models, metallicity assumptions, and SFHs separately to identify which assumptions maximize the agreement (correlation and linearity) between L TIR and L dust. We find that a flexible dust attenuation law performs best. For stellar populations, we find that Bruzual & Charlot models are favored over those of Eldridge et al. Fixing the stellar metallicity at solar value is preferred to other fixed values or leaving it as a free parameter. For SFHs, we find that minimizing the variability in the recent (<100 Myr) SFH improves the agreement with L TIR. Finally, we provide a catalog of galaxy parameters (including M * and SFR) for CANDELS galaxies with logM*>8 and 0.7 < z < 1.3, obtained using the models we found to be the most robust.

Revised Constraints on the Fast Radio Burst Population from the First CHIME/FRB Catalog

In this paper, we investigate the fast radio burst (FRB) population using the first CHIME/FRB catalog. We first reconstruct the extragalactic dispersion measure–redshift relation (DME–z relation) from well-localized FRBs, then use it to infer the redshift and isotropic energy of the first CHIME/FRB catalog. The intrinsic energy distribution is modeled by the power law with an exponential cutoff, and the selection effect of the CHIME telescope is modeled by a two-parametric function of specific fluence. For the intrinsic redshift distribution, the star formation history (SFH) model and five other SFH-related models are considered. We construct the joint likelihood of fluence, energy, and redshift, and all the free parameters are constrained simultaneously using the Bayesian inference method. The Bayesian information criterion (BIC) is used to choose the model that best matches the observational data. For comparison, we fit our models with two data samples, i.e., the Full sample and the Gold sample. The power-law index and cutoff energy are tightly constrained to be 1.8 ≲ α ≲ 1.9 and log(Ec/erg)≈42 , which are almost independent of the redshift distribution model and the data sample we choose. The parameters involving the selection effect strongly depend on the data sample but are insensitive to the redshift distribution model. According to BIC, the pure SFH model is strongly disfavored by both the Full sample and Gold sample. For the other five SFH-related redshift distribution models, most of them can match the data well if the parameters are properly chosen. Therefore, with the present data, it is still premature to draw a conclusion on the FRB population.

Stellar Population Properties in the Stellar Streams around SPRC047

We have investigated the properties (e.g., age, metallicity) of the stellar populations of a ringlike tidal stellar stream (or streams) around the edge-on galaxy SPRC047 (z = 0.031) using spectral energy distribution (SED) fits to integrated broadband aperture flux densities. We used visual images in six different bands and Spitzer/IRAC 3.6 μm data. We have attempted to derive best-fit stellar population parameters (metallicity, age) in three noncontiguous segments of the stream. Due to the very low surface brightness of the stream, we have performed a deconvolution with a Richardson–Lucy–type algorithm of the low spatial resolution 3.6 μm IRAC image, thereby reducing the effect of the point-spread function aliased emission from the bright edge-on central galaxy at the locations of our three stream segments. Our SED fits that used several different star formation (SF) history priors, from an exponentially decaying SF burst to continuous SF, indicate that the age–metallicity–dust degeneracy is not resolved, most likely because of inadequate wavelength coverage and low signal-to-noise ratios of the low surface brightness features. We also discuss how future deep visual–near-infrared observations, combined with absolute flux calibration uncertainties at or below the 1% level, complemented by equally well absolute flux-calibrated observations in ultraviolet and mid-infrared bands, would improve the accuracy of broadband SED fitting results for low surface brightness targets, such as stellar streams around nearby galaxies that are not resolved into stars.

On the Evolutionary History of a Simulated Disk Galaxy as Seen by Phylogenetic Trees

Phylogenetic methods have long been used in biology and more recently have been extended to other fields—for example, linguistics and technology—to study evolutionary histories. Galaxies also have an evolutionary history and fall within this broad phylogenetic framework. Under the hypothesis that chemical abundances can be used as a proxy for the interstellar medium’s DNA, phylogenetic methods allow us to reconstruct hierarchical similarities and differences among stars—essentially, a tree of evolutionary relationships and thus history. In this work, we apply phylogenetic methods to a simulated disk galaxy obtained with a chemodynamical code to test the approach. We found that at least 100 stellar particles are required to reliably portray the evolutionary history of a selected stellar population in this simulation, and that the overall evolutionary history is reliably preserved when the typical uncertainties in the chemical abundances are smaller than 0.08 dex. The results show that the shapes of the trees are strongly affected by the age–metallicity relation, as well as the star formation history of the galaxy. We found that regions with low star formation rates produce shorter trees than regions with high star formation rates. Our analysis demonstrates that phylogenetic methods can shed light on the process of galaxy evolution.

Resolving Clumpy versus Extended Lyα in Strongly Lensed, High-redshift Lyα Emitters

We present six strongly gravitationally lensed Lyα emitters (LAEs) at z ∼ 4–5 with Hubble Space Telescope (HST) narrowband imaging isolating Lyα. Through complex radiative transfer Lyα encodes information about the spatial distribution and kinematics of the neutral hydrogen upon which it scatters. We investigate the galaxy properties and Lyα morphologies of our sample. Many previous studies of high-redshift LAEs have been limited in Lyα spatial resolution. In this work we take advantage of high-resolution Lyα imaging boosted by lensing magnification, allowing us to probe subgalactic scales that are otherwise inaccessible at these redshifts. We use broadband imaging from HST (rest-frame UV) and Spitzer (rest-frame optical) in spectral energy distribution fitting, providing estimates of the stellar masses (∼108–109 M ⊙), stellar population ages (t 50 < 40 Myr), and amounts of dust (A V ∼ 0.1–0.6, statistically consistent with zero). We employ nonparametric star formation histories to probe the young stellar populations which create the Lyα. We also examine the offsets between the Lyα and stellar continuum, finding small upper limits of offsets (<0.″1) consistent with studies of low-redshift LAEs, indicating our galaxies are not interacting or merging. Finally, we find a bimodality in our sample’s Lyα morphologies: clumpy and extended. We find a suggestive trend: our LAEs with clumpy Lyα are generally younger than the LAEs with extended Lyα, suggesting a possible correlation with age.

A Global Semianalytic Model of the First Stars and Galaxies Including Dark Matter Halo Merger Histories

We present a new self-consistent semianalytic model of the first stars and galaxies to explore the high-redshift (z ≥ 15) Population III (PopIII) and metal-enriched star formation histories. Our model includes the detailed merger history of dark matter halos generated with Monte Carlo merger trees. We calibrate the minimum halo mass for PopIII star formation from recent hydrodynamical cosmological simulations that simultaneously include the baryon–dark matter streaming velocity, Lyman–Werner (LW) feedback, and molecular hydrogen self-shielding. We find an overall increase in the resulting star formation rate density (SFRD) compared to calibrations based on previous simulations (e.g., the PopIII SFRD is over an order of magnitude higher at z = 35−15). We evaluate the effect of the halo-to-halo scatter in this critical mass and find that it increases the PopIII stellar mass density by a factor ∼1.5 at z ≥ 15. Additionally, we assess the impact of various semianalytic/analytic prescriptions for halo assembly and star formation previously adopted in the literature. For example, we find that models assuming smooth halo growth computed via abundance matching predict SFRDs similar to the merger tree model for our fiducial model parameters, but that they may underestimate the PopIII SFRD in cases of strong LW feedback. Finally, we simulate subvolumes of the Universe with our model both to quantify the reduction in total star formation in numerical simulations due to a lack of density fluctuations on spatial scales larger than the simulation box, and to determine spatial fluctuations in SFRD due to the diversity in halo abundances and merger histories.