Rebounding Cores to Build Star Cluster Multiple Populations

We present a novel approach to the riddle of star cluster multiple populations. Stars form from molecular cores. But not all cores form stars. Following their initial compression, such “failed” cores re-expand, rather than collapsing. We propose that their formation and subsequent dispersal regulate the gas density of cluster-forming clumps and, therefore, their core and star formation rates. Clumps for which failed cores are the dominant core type experience star formation histories with peaks and troughs (i.e., discrete star formation episodes). In contrast, too few failed cores results in smoothly decreasing star formation rates. We identify three main parameters shaping the star formation history of a clump: the star and core formation efficiencies per free-fall time, and the timescale on which failed cores return to the clump gas. The clump mass acts as a scaling factor. We use our model to constrain the density and mass of the Orion Nebula Cluster progenitor clump, and to caution that the star formation histories of starburst clusters may contain close-by peaks concealed by stellar age uncertainties. Our model generates a great variety of star formation histories. Intriguingly, the chromosome maps and O–Na anticorrelations of old globular clusters also present diverse morphologies. This prompts us to discuss our model in the context of globular cluster multiple stellar populations. More massive globular clusters exhibit stronger multiple stellar population patterns, which our model can explain if the formation of the polluting stars requires a given stellar mass threshold.

Predicting the Water Content of Interstellar Objects from Galactic Star Formation Histories

Planetesimals inevitably bear the signatures of their natal environment, preserving in their composition a record of the metallicity of their system’s original gas and dust, albeit one altered by the formation processes. When planetesimals are dispersed from their system of origin, this record is carried with them. As each star is likely to contribute at least 1012 interstellar objects (ISOs), the Galaxy’s drifting population of ISOs provides an overview of the properties of its stellar population through time. Using the EAGLE cosmological simulations and models of protoplanetary formation, our modeling predicts an ISO population with a bimodal distribution in their water mass fraction: objects formed in low-metallicity, typically older, systems have a higher water fraction than their counterparts formed in high-metallicity protoplanetary disks, and these water-rich objects comprise the majority of the population. Both detected ISOs seem to belong to the lower water fraction population; these results suggest they come from recently formed systems. We show that the population of ISOs in galaxies with different star formation histories will have different proportions of objects with high and low water fractions. This work suggests that it is possible that the upcoming Vera C. Rubin Observatory Legacy Survey of Space and Time will detect a large enough population of ISOs to place useful constraints on models of protoplanetary disks, as well as galactic structure and evolution.

Spatially Resolving the Star Formation Histories of Three Nearby Nuclear Star Clusters

The formation of nuclear star clusters (NSCs) remains an open question. In this work, we use spatially-resolved HST/STIS spectroscopic observations of three nearby NSCs (hosted by NGC 5102, NGC 5206, and NGC 205) to constrain their formation histories by exploring radial variations of the stellar populations within each cluster. Utilizing full-spectrum fitting, we find substantial age and metallicity gradients within the central 0.″9 (16 pc) of the NSC in NGC 5102 where populations near the center are young/metal-rich (age ∼400 Myr and [M/H] ∼ −0.4) and become older/metal-poor at larger radii (mean age ∼1 Gyr and mean [M/H] ∼ −1.6 in the radial range [0.″3, 0.″9]). This behavior suggests that the young/metal-rich population at the center was formed from a period of in situ formation, while the older/metal-poor populations were likely formed by inspiraled globular clusters. The two broad populations observed in the NGC 5102 NSC (young/metal-rich and old/metal-poor) appear to be linked to the transition between the two morphological components of the NSC derived from the surface-brightness profile in Nguyen et al. (2018). The radial ranges explored in NGC 5206 and NGC 205 were much smaller due to poor data quality; in NGC 5206 we find a similar metallicity gradient to NGC 5102 (but with much lower significance), while the data for NGC 205 is too poor to reach any conclusions. Overall, this data highlights the links between the morphological and stellar population complexity of NSCs and their formation mechanisms.

Evolution of the Ultraviolet Upturn at 0.3 < z < 1: Exploring Helium-rich Stellar Populations

We measure the near-UV (rest-frame ∼2400 Å) to optical color for early-type galaxies in 12 clusters at 0.3 < z < 1.0. We show that this is a suitable proxy for the more common far-ultraviolet bandpass used to measure the ultraviolet upturn and find that the upturn is detected to z = 0.6 in these data, in agreement with previous work. We find evidence that the strength of the upturn starts to wane beyond this redshift and largely disappears at z = 1. Our data are most consistent with models where early-type galaxies contain minority stellar populations with non-cosmological helium abundances, up to around 46%, formed at z ≥ 3, resembling globular clusters with multiple stellar populations in our Galaxy. This suggests that elliptical galaxies and globular clusters share similar chemical evolution and star formation histories. The vast majority of the stellar mass in these galaxies also must have been in place at z > 3.

On the Properties of Spectroscopically Confirmed Ultra-diffuse Galaxies across Environments

We present new redshift measurements for 19 candidate ultra-diffuse galaxies (UDGs) from the Systematically Measuring Ultra-Diffuse Galaxies (SMUDGes) survey after conducting a long-slit spectroscopic follow-up campaign on 23 candidates with the Large Binocular Telescope. We combine these results with redshift measurements from other sources for 29 SMUDGes and 20 non-SMUDGes candidate UDGs. Together, this sample yields 44 spectroscopically confirmed UDGs (r e ≥ 1.5 kpc and μ g (0) ≥ 24 mag arcsec−2 within uncertainties) and spans cluster and field environments, with all but one projected on the Coma cluster and environs. We find no statistically significant differences in the structural parameters of cluster and noncluster confirmed UDGs, although there are hints of differences among the axis ratio distributions. Similarly, we find no significant structural differences among those in locally dense or sparse environments. However, we observe a significant difference in color with respect to projected clustercentric radius, confirming trends observed previously in statistical UDG samples. This trend strengthens further when considering whether UDGs reside in either cluster or locally dense environments, suggesting starkly different star formation histories for UDGs residing in high- and low-density environments. Of the 16 large (r e ≥ 3.5 kpc) UDGs in our sample, only one is a field galaxy that falls near the early-type galaxy red sequence. No other field UDGs found in low-density environments fall near the red sequence. This finding, in combination with our detection of Galaxy Evolution Explorer NUV flux in nearly half of the UDGs in sparse environments, suggests that field UDGs are a population of slowly evolving galaxies.

On the Impact of Inclination-dependent Attenuation on Derived Star Formation Histories: Results from Disk Galaxies in the Great Observatories Origins Deep Survey Fields

We develop and implement an inclination-dependent attenuation prescription for spectral energy distribution (SED) fitting and study its impact on derived star formation histories. We apply our prescription within the SED fitting code Lightning to a clean sample of 82, z = 0.21–1.35 disk-dominated galaxies in the Great Observatories Origins Deep Survey North and South fields. To compare our inclination-dependent attenuation prescription with more traditional fitting prescriptions, we also fit the SEDs with the inclination-independent Calzetti et al. (2000) attenuation curve. From this comparison, we find that fits to a subset of 58, z < 0.7 galaxies in our sample, utilizing the Calzetti et al. (2000) prescription, recover similar trends with inclination as the inclination-dependent fits for the far-UV-band attenuation and recent star formation rates. However, we find a difference between prescriptions in the optical attenuation (A V ) that is strongly correlated with inclination (p‐value < 10−11). For more face-on galaxies, with i ≲ 50°, (edge-on, i ≈ 90°), the average derived A V is 0.31 ± 0.11 magnitudes lower (0.56 ± 0.16 magnitudes higher) for the inclination-dependent model compared to traditional methods. Further, the ratio of stellar masses between prescriptions also has a significant (p‐value < 10−2) trend with inclination. For i = 0°–65°, stellar masses are systematically consistent between fits, with log 10 ( M ⋆ inc / M ⋆ Calzetti ) = − 0.05 ± 0.03 dex and scatter of 0.11 dex. However, for i ≈ 80°–90°, the derived stellar masses are lower for the Calzetti et al. (2000) fits by an average factor of 0.17 ± 0.03 dex and scatter of 0.13 dex. Therefore, these results suggest that SED fitting assuming the Calzetti et al. (2000) attenuation law potentially underestimates stellar masses in highly inclined disk-dominated galaxies.

Reproducing the UVJ Color Distribution of Star-forming Galaxies at 0.5 < z < 2.5 with a Geometric Model of Dust Attenuation

We analyze the distribution of rest-frame U − V and V − J colors for star-forming galaxies at 0.5 < z < 2.5. Using stellar population synthesis, stochastic star formation histories, and a simple prescription for the dust attenuation that accounts for the shape and inclination of galaxies, we construct a model for the distribution of galaxy colors. With only two free parameters, this model is able to reproduce the observed galaxy colors as a function of redshift and stellar mass remarkably well. Our analysis suggests that the wide range of dust attenuation values measured for star-forming galaxies at a given redshift and stellar mass is almost entirely due to the effect of inclination; if all galaxies at a given stellar mass were observed edge-on, they would show very similar dust attenuation. This result has important implications for the interpretation of dust attenuation measurements, the treatment of UV and IR luminosity, and the comparison between numerical simulations and observations.

Do Current X-Ray Observations Capture Most of the Black-hole Accretion at High Redshifts?

The cosmic black hole accretion density (BHAD) is critical for our understanding of the formation and evolution of supermassive black holes (BHs). However, at high redshifts (z > 3), X-ray observations report BHADs significantly (∼10 times) lower than those predicted by cosmological simulations. It is therefore paramount to constrain the high-z BHAD using independent methods other than direct X-ray detections. The recently established relation between star formation rate and BH accretion rate among bulge-dominated galaxies provides such a chance, as it enables an estimate of the BHAD from the star formation histories (SFHs) of lower-redshift objects. Using the CANDELS Lyα Emission At Reionization (CLEAR) survey, we model the SFHs for a sample of 108 bulge-dominated galaxies at z = 0.7–1.5, and further estimate the BHAD contributed by their high-z progenitors. The predicted BHAD at z ≈ 4–5 is consistent with the simulation-predicted values, but higher than the X-ray measurements (by ≈3–10 times at z = 4–5). Our result suggests that the current X-ray surveys could be missing many heavily obscured Compton-thick active galactic nuclei (AGNs) at high redshifts. However, this BHAD estimation assumes that the high-z progenitors of our z = 0.7–1.5 sample remain bulge-dominated where star formation is correlated with BH cold-gas accretion. Alternatively, our prediction could signify a stark decline in the fraction of bulges in high-z galaxies (with an associated drop in BH accretion). JWST and Origins will resolve the discrepancy between our predicted BHAD and the X-ray results by constraining Compton-thick AGN and bulge evolution at high redshifts.

Constraining Type Ia Supernova Delay Time with Spatially Resolved Star Formation Histories

We present the delay time distribution (DTD) estimates of Type Ia supernovae (SNe Ia) using spatially resolved SN Ia host galaxy spectra from MUSE and MaNGA. By employing a grouping algorithm based on k-means and earth mover’s distances (EMDs), we separated the host galaxy stellar population age distributions (SPADs) into spatially distinct regions and used maximum likelihood method to constrain the DTD of SN Ia progenitors. When a power-law model of the form DTD(t) ∝ t s (t > τ) is used, we find an SN rate decay slope s = − 1.41 − 0.33 + 0.32 and a delay time τ = 120 − 83 + 142 Myr . Moreover, we tested other DTD models, such as a broken power-law model and a two-component power-law model, and found no statistically significant support for these alternative models.