ABSTRACT
We present a method to flexibly and self-consistently determine individual galaxies’ star formation rates (SFRs) from their host haloes’ potential well depths, assembly histories, and ...redshifts. The method is constrained by galaxies’ observed stellar mass functions, SFRs (specific and cosmic), quenched fractions, ultraviolet (UV) luminosity functions, UV–stellar mass relations, IRX–UV relations, auto- and cross-correlation functions (including quenched and star-forming subsamples), and quenching dependence on environment; each observable is reproduced over the full redshift range available, up to 0 < z < 10. Key findings include the following: galaxy assembly correlates strongly with halo assembly; quenching correlates strongly with halo mass; quenched fractions at fixed halo mass decrease with increasing redshift; massive quenched galaxies reside in higher-mass haloes than star-forming galaxies at fixed galaxy mass; star-forming and quenched galaxies’ star formation histories at fixed mass differ most at z < 0.5; satellites have large scatter in quenching time-scales after infall, and have modestly higher quenched fractions than central galaxies; Planck cosmologies result in up to 0.3 dex lower stellar – halo mass ratios at early times; and, none the less, stellar mass–halo mass ratios rise at z > 5. Also presented are revised stellar mass – halo mass relations for all, quenched, star-forming, central, and satellite galaxies; the dependence of star formation histories on halo mass, stellar mass, and galaxy SSFR; quenched fractions and quenching time-scale distributions for satellites; and predictions for higher-redshift galaxy correlation functions and weak lensing surface densities. The public data release (DR1) includes the massively parallel (>105 cores) implementation (the UniverseMachine), the newly compiled and remeasured observational data, derived galaxy formation constraints, and mock catalogues including lightcones.
The spectral absorption lines in early-type galaxies contain a wealth of information regarding the detailed abundance pattern, star formation history, and stellar initial mass function (IMF) of the ...underlying stellar population. Using our new population synthesis model that accounts for the effect of variable abundance ratios of 11 elements, we analyze very high quality absorption line spectra of 38 early-type galaxies and the nuclear bulge of M31. These data extend to 1 mu m and they therefore include the IMF-sensitive spectral features Na I, Ca II, and FeH at 0.82 mu m, 0.86 mu m, and 0.99 mu m, respectively. The models fit the data well, with typical rms residuals <, ~1%. Strong constraints on the IMF and therefore the stellar mass-to-light ratio, (M/L) sub(stars), are derived for individual galaxies. We find that the IMF becomes increasingly bottom-heavy with increasing velocity dispersion and Mg/Fe, At the lowest dispersions and Mg/Fe values the derived IMF is consistent with the Milky Way (MW) IMF, while at the highest dispersions and Mg/Fe values the derived IMF contains more low-mass stars (is more bottom-heavy) than even a Salpeter IMF. Our best-fit (M/L) sub(stars) values do not exceed dynamically based M/L values. We also apply our models to stacked spectra of four metal-rich globular clusters in M31 and find an (M/L) sub(stars) that implies fewer low-mass stars than a MW IMF, again agreeing with dynamical constraints. We discuss other possible explanations for the observed trends and conclude that variation in the IMF is the simplest and most plausible.
Stellar population synthesis (SPS) provides the link between the stellar and dust content of galaxies and their observed spectral energy distributions. In the present work, we perform a comprehensive ...calibration of our own flexible SPS (FSPS) model against a suite of data. These data include ultraviolet, optical, and near-IR photometry, surface brightness fluctuations, and integrated spectra of star clusters in the Magellanic Clouds (MCs), M87, M31, and the Milky Way (MW), and photometry and spectral indices of both quiescent and post-starburst galaxies at z {approx} 0. Several public SPS models are intercompared, including the models of Bruzual and Charlot (BC03), Maraston (M05), and FSPS. The relative strengths and weaknesses of these models are evaluated, with the following conclusions: (1) the FSPS and BC03 models compare favorably with MC data at all ages, whereas M05 colors are too red and the age dependence is incorrect; (2) all models yield similar optical and near-IR colors for old metal-poor systems, and yet they all provide poor fits to the integrated J - K and V - K colors of both MW and M31 star clusters; (3) FSPS is able to fit all of the ultraviolet data because both the post-asymptotic giant branch (post-AGB) and horizontal branch evolutionary phases are handled flexibly, while the BC03 and M05 models fail in the far-UV, and both far- and near-UV, respectively; (4) all models predict ugr colors too red, D{sub n}4000 strengths too strong, and Hdelta{sub A} strengths too weak compared to massive red sequence galaxies, under the assumption that such galaxies are composed solely of old metal-rich stars; and (5) FSPS and, to a lesser extent, BC03 can reproduce the optical and near-IR colors of post-starburst galaxies, while M05 cannot. Reasons for these discrepancies are explored. The failure at predicting the ugr colors, D{sub n}4000, and Hdelta{sub A} strengths can be explained by some combination of a minority population of metal-poor stars, young stars, blue straggler and/or blue horizontal branch (HB) stars, but not by appealing to inadequacies in either theoretical stellar atmospheres or canonical evolutionary phases (e.g., the main-sequence turnoff). The different model predictions in the near-IR for intermediate age systems are due to different treatments of the thermally pulsating asymptotic giant branch stellar evolutionary phase. We emphasize that due to a lack of calibrating star cluster data in regions of the metallicity-age plane relevant for galaxies, all of these models continue to suffer from serious uncertainties that are difficult to quantify.
Accounting for nebular emission when modeling galaxy spectral energy distributions (SEDs) is important, as both line and continuum emissions can contribute significantly to the total observed flux. ...In this work, we present a new nebular emission model integrated within the Flexible Stellar Population Synthesis code that computes the line and continuum emission for complex stellar populations using the photoionization code Cloudy. The self-consistent coupling of the nebular emission to the matched ionizing spectrum produces emission line intensities that correctly scale with the stellar population as a function of age and metallicity. This more complete model of galaxy SEDs will improve estimates of global gas properties derived with diagnostic diagrams, star formation rates based on H , and physical properties derived from broadband photometry. Our models agree well with results from other photoionization models and are able to reproduce observed emission from H ii regions and star-forming galaxies. Our models show improved agreement with the observed H ii regions in the Ne iii/O ii plane and show satisfactory agreement with He ii emission from z = 2 galaxies, when including rotating stellar models. Models including post-asymptotic giant branch stars are able to reproduce line ratios consistent with low-ionization emission regions. The models are integrated into current versions of FSPS and include self-consistent nebular emission predictions for MIST and Padova+Geneva evolutionary tracks.
We investigate the integrated properties of massive ( ) rotating single-star stellar populations for a variety of initial rotation rates ( , 0.2, 0.4, 0.5, and 0.6). We couple the new MESA Isochrone ...and Stellar Tracks (MIST) models to the Flexible Stellar Population Synthesis (FSPS) package, extending the stellar population synthesis models to include the contributions from very massive stars ( ), which can be significant in the first ∼4 Myr after a starburst. These models predict ionizing luminosities that are consistent with recent observations of young nuclear star clusters. We also construct composite stellar populations assuming a distribution of initial rotation rates. Even in low-metallicity environments where rotation has a significant effect on the evolution of massive stars, we find that stellar population models require a significant contribution from fast-rotating ( ) stars in order to sustain the production of ionizing photons beyond a few Myr following a starburst. These results have potentially important implications for cosmic reionization by massive stars and the interpretation of nebular emission lines in high-redshift star-forming galaxies.
Nonparametric star formation histories (SFHs) have long promised to be the "gold standard" for galaxy spectral energy distribution (SED) modeling as they are flexible enough to describe the full ...diversity of SFH shapes, whereas parametric models rule out a significant fraction of these shapes a priori. However, this flexibility is not fully constrained even with high-quality observations, making it critical to choose a well-motivated prior. Here, we use the SED-fitting code Prospector to explore the effect of different nonparametric priors by fitting SFHs to mock UV-IR photometry generated from a diverse set of input SFHs. First, we confirm that nonparametric SFHs recover input SFHs with less bias and return more accurate errors than do parametric SFHs. We further find that, while nonparametric SFHs robustly recover the overall shape of the input SFH, the primary determinant of the size and shape of the posterior star formation rate as a function of time (SFR(t)) is the choice of prior, rather than the photometric noise. As a practical demonstration, we fit the UV-IR photometry of ∼6000 galaxies from the Galaxy and Mass Assembly survey and measure scatters between priors to be 0.1 dex in mass, 0.8 dex in SFR100 Myr, and 0.2 dex in mass-weighted ages, with the bluest star-forming galaxies showing the most sensitivity. An important distinguishing characteristic for nonparametric models is the characteristic timescale for changes in SFR(t). This difference controls whether galaxies are assembled in bursts or in steady-state star formation, corresponding respectively to (feedback-dominated/accretion-dominated) models of galaxy formation and to (larger/smaller) confidence intervals derived from SED fitting. High-quality spectroscopy has the potential to further distinguish between these proposed models of SFR(t).
We discuss the approach of searching the lowest mass dwarf galaxies, , in the general field, using integrated light surveys. By exploring the limiting surface brightness-spatial resolution ( ...eff,lim−θ) parameter space, we suggest that faint field dwarfs in the Local Volume, between 3 and 10 Mpc, are expected to be detected very effectively and in large numbers using integrated light photometric surveys, complementary to the classical star counts method. We use a sample of dwarf galaxies in the Local Group to construct relations between their photometric and structural parameters, M*- eff,V and M*-Reff. We use these relations, along with assumed functional forms for the halo mass function and the stellar mass-halo mass (SMHM) relation, to calculate the lowest detectable stellar masses in the Local Volume and the expected number of galaxies as a function of the limiting surface brightness and spatial resolution. The number of detected galaxies depends mostly on the limiting surface brightness for distances >3 Mpc, while spatial resolution starts to play a role for galaxies at distances >8 Mpc. Surveys with eff,lim ∼ 30 mag arcsec−2 should be able to detect galaxies with stellar masses down to ∼104 M in the Local Volume. Depending on the form of the SMHM relation, the expected number of dwarf galaxies with distances between 3 and 10 Mpc is 0.04-0.35 per square degree, assuming a limiting surface brightness of ∼29-30 mag arcsec−2 and a spatial resolution <4″. We plan to search for a population of low-mass dwarf galaxies in the field by performing a blank wide field photometric survey with the Dragonfly Telephoto Array, an imaging system optimized for the detection of extended ultra low surface brightness structures.
We use the first Gaia data release, combined with the RAVE and APOGEE spectroscopic surveys, to investigate the origin of halo stars within kpc from the Sun. We identify halo stars kinematically as ...moving at a relative speed of at least 220 km s−1 with respect to the local standard of rest. These stars are generally less metal-rich than the disk, but surprisingly, half of our halo sample is comprised of stars with . The orbital directions of these metal-rich halo stars are preferentially aligned with the disk rotation, in sharp contrast with the intrinsically isotropic orbital distribution of the metal-poor halo stars. We find similar properties in the Latte cosmological zoom-in simulation of a Milky Way-like galaxy from the FIRE project. In Latte, metal-rich halo stars formed primarily inside of the solar circle, whereas lower-metallicity halo stars preferentially formed at larger distances (extending beyond the virial radius). This suggests that metal-rich halo stars in the solar neighborhood actually formed in situ within the Galactic disk, rather than having been accreted from satellite systems. These stars, currently on halo-like orbits, therefore have likely undergone substantial radial migration/heating.
ABSTRACT 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 ( ), masses ( ), and metallicities ( ). The models are self-consistently and continuously evolved from the pre-main sequence (PMS) 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 PMS to the end of core helium burning for . 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/.