Formation History Of Galaxies Research Articles

The XXL Survey

Exploiting a sample of galaxies drawn from the XXL-North multiwavelength survey, we present an analysis of the stellar population properties of galaxies at 0.1 ≤ z ≤ 0.5, by studying galaxy fractions and the star formation rate (SFR)–stellar mass (M⋆) relation. Furthermore, we exploit and compare two parametrisations of environment. When adopting a definition of “global” environment, we consider separately cluster virial (r ≤ 1r200) and outer (1r200 < r ≤ 3r200) members and field galaxies. We also distinguish between galaxies that belong or do not belong to superclusters, but never find systematic differences between the two subgroups. When considering the “local” environment, we take into account the projected number density of galaxies in a fixed aperture of 1 Mpc in the sky. We find that regardless of the environmental definition adopted, the fraction of blue or star-forming galaxies is the highest in the field or least dense regions and the lowest in the virial regions of clusters or highest densities. Furthermore, the fraction of star-forming galaxies is higher than the fraction of blue galaxies, regardless of the environment. This result is particularly evident in the virial cluster regions, most likely reflecting the different star formation histories of galaxies in different environments. Also the overall SFR–M⋆ relation does not seem to depend on the parametrisation adopted. Nonetheless, the two definitions of environment lead to different results as far as the fraction of galaxies in transition between the star-forming main sequence and the quenched regime is concerned. In fact, using the local environment the fraction of galaxies below the main sequence is similar at low and high densities, whereas in clusters (and especially within the virial radii) a population with reduced SFR with respect to the field is observed. Our results show that the two parametrisations adopted to describe the environment have different physical meanings, i.e. are intrinsically related to different physical processes acting on galaxy populations and are able to probe different physical scales.

Combining stellar populations with orbit-superposition dynamical modelling: the formation history of the lenticular galaxy NGC 3115

Abstract We present a combination of the Schwarzschild orbit-superposition dynamical modelling technique with the spatially resolved mean stellar age and metallicity maps to uncover the formation history of galaxies. We apply this new approach to a remarkable five-pointing mosaic of VLT/MUSE observations obtained by Guérou et al. (2016) extending to a maximum galactocentric distance of ${\sim } {120}{\, {\rm arcsec}}\ \left({5.6}\, {\rm kpc}\right)$ along the major axis, corresponding to ∼2.5Re. Our method first identifies ‘families’ of orbits from the dynamical model that represent dynamically distinct structures of the galaxy. Individual ages and metallicities of these components are then fit for using the stellar-population information. Our results highlight components of the galaxy that are distinct in the combined stellar dynamics/populations space, which implies distinct formation paths. We find evidence for a dynamically cold, metal-rich disc, consistent with a gradual in situ formation. This disc is embedded in a generally old population of stars, with kinematics ranging from dispersion dominated in the centre to an old, diffuse, metal-poor stellar halo at the extremities. We find also a direct correlation between the dominant dynamical support of these components, and their associated age, akin to the relation observed in the Milky Way. This approach not only provides a powerful model for inferring the formation history of external galaxies but also paves the way to a complete population-dynamical model.

Horizon-AGN virtual observatory – 1. SED-fitting performance and forecasts for future imaging surveys

Abstract Using the light-cone from the cosmological hydrodynamical simulation horizon-AGN, we produced a photometric catalogue over 0 < z < 4 with apparent magnitudes in COSMOS, Dark Energy Survey, Large Synoptic Survey Telescope (LSST)-like, and Euclid-like filters at depths comparable to these surveys. The virtual photometry accounts for the complex star formation history (SFH) and metal enrichment of horizon-AGN galaxies, and consistently includes magnitude errors, dust attenuation, and absorption by intergalactic medium. The COSMOS-like photometry is fitted in the same configuration as the COSMOS2015 catalogue. We then quantify random and systematic errors of photometric redshifts, stellar masses, and star formation rates (SFR). Photometric redshifts and redshift errors capture the same dependencies on magnitude and redshift as found in COSMOS2015, excluding the impact of source extraction. COSMOS-like stellar masses are well recovered with a dispersion typically lower than 0.1 dex. The simple SFHs and metallicities of the templates induce a systematic underestimation of stellar masses at z < 1.5 by at most 0.12 dex. SFR estimates exhibit a dust-induced bimodality combined with a larger scatter (typically between 0.2 and 0.6 dex). We also use our mock catalogue to predict photometric redshifts and stellar masses in future imaging surveys. We stress that adding Euclid near-infrared photometry to the LSST-like baseline improves redshift accuracy especially at the faint end and decreases the outlier fraction by a factor ∼2. It also considerably improves stellar masses, reducing the scatter up to a factor 3. It would therefore be mutually beneficial for LSST and Euclid to work in synergy.

The time delay between star formation quenching and morphological transformation of galaxies in clusters: a phase–space view of EDisCS

AbstractWe explore the possible effect of cluster environments on the structure and star formation histories of galaxies by analysing the projected phase–space (PPS) of intermediate-redshift clusters (0.4 ≤ z ≤ 0.8). HST I−band imaging data from the ESO Distant Cluster Survey (EDisCS) allow us to measure deviations of the galaxies’ light distributions from symmetric and smooth profiles using two parameters, Ares (‘asymmetry’) and RFF (residual flux fraction or ‘roughness’). Combining these structural parameters with age-sensitive spectral indicators ($H_{\delta \rm {A}}$, $H_{\gamma \rm {A}}$, and Dn4000), we establish that in all environments younger star-forming galaxies of all morphologies are ‘rougher’ and more asymmetric than older, more quiescent ones. Combining a subset of the EDisCS clusters, we construct a stacked PPS diagram and find a significant correlation between the position of the galaxies on the PPS and their stellar ages, irrespective of their morphology. We also observe an increasing fraction of galaxies with older stellar populations towards the cluster core, while the galaxies’ structural parameters (Ares and RFF) do not seem to segregate strongly with PPS. These results may imply that, under the possible influence of their immediate cluster environment, galaxies have their star formation suppressed earlier, while their structural transformation happens on a longer time-scale as they accumulate and age in the cluster cores.

FIGS: spectral fitting constraints on the star formation history of massive galaxies since the cosmic noon

Abstract We constrain the stellar population properties of a sample of 52 massive galaxies – with stellar mass log (Ms/M⊙) ≳ 10.5 – over the redshift range 0.5 < z < 2 by use of observer-frame optical and near-infrared slitless spectra from Hubble Space Telescope’s ACS and WFC3 grisms. The deep exposures (∼100 ks) allow us to target individual spectra of massive galaxies to F160W = 22.5 AB. Our spectral fitting approach uses a set of six base models adapted to the redshift and spectral resolution of each observation, and fits the weights of the base models, including potential dust attenuation, via a Markov Chain Monte Carlo method. Our sample comprises a mixed distribution of quiescent (19) and star-forming galaxies (33). We quantify the width of the age distribution (Δt) that is found to dominate the variance of the retrieved parameters according to principal component analysis. The population parameters follow the expected trend towards older ages with increasing mass, and Δt appears to weakly anticorrelate with stellar mass, suggesting a more efficient star formation at the massive end. As expected, the redshift dependence of the relative stellar age (measured in units of the age of the Universe at the source) in the quiescent sample rejects the hypothesis of a single burst (aka monolithic collapse). Radial colour gradients within each galaxy are also explored, finding a wider scatter in the star-forming subsample, but no conclusive trend with respect to the population parameters.

The star formation histories of dwarf galaxies in Local Group cosmological simulations

We use the APOSTLE and Auriga cosmological simulations to study the star formation histories (SFHs) of field and satellite dwarf galaxies. Despite sizeable galaxy-to-galaxy scatter, the SFHs of APOSTLE and Auriga dwarfs exhibit robust average trends with galaxy stellar mass: faint field dwarfs ($10^5<M_{\rm star}/M_\odot<10^{6.5}$) have, on average, steadily declining SFHs, whereas brighter dwarfs ($10^{7.5}<M_{\rm star}/M_\odot<10^{9}$) show the opposite trend. Intermediate-mass dwarfs have roughly constant SFHs. Satellites exhibit similar average trends, but with substantially suppressed star formation in the most recent $\sim 5$ Gyr, likely as a result of gas loss due to tidal and ram-pressure stripping after entering the haloes of their primaries. These simple mass and environmental trends are in good agreement with the derived SFHs of Local Group (LG) dwarfs whose photometry reaches the oldest main sequence turnoff. SFHs of galaxies with less deep data show deviations from these trends, but this may be explained, at least in part, by the large galaxy-to-galaxy scatter, the limited sample size, and the large uncertainties of the inferred SFHs. Confirming the predicted mass and environmental trends will require deeper photometric data than currently available, especially for isolated dwarfs.

Recovering the star formation history of galaxies through spectral fitting: Current challenges

AbstractWith the exception of some nearby galaxies, we cannot resolve stars individually. To recover the galaxies star formation history (SFH), the challenge is to extract information from their integrated spectrum. A widely used tool is the full spectral fitting technique. This consists of combining simple stellar populations (SSPs) of different ages and metallicities to match the integrated spectrum. This technique works well for optical spectra, for metallicities near solar and chemical histories not much different from our Galaxy. For everything else there is room for improvement. With telescopes being able to explore further and further away, and beyond the optical, the improvement of this type of tool is crucial. SSPs use as ingredients isochrones, an initial mass function, and a library of stellar spectra. My focus are the stellar libraries, key ingredient for SSPs. Here I talk about the latest developments of stellar libraries, how they influence the SSPs and how to improve them.

Late Bloomer Galaxies: Growing Up in Cosmic Autumn

Abstract Late bloomers (LBs) are massive ( ) galaxies at z &lt; 1 that formed the majority of their stars within ∼2 Gyr of the epoch of observation. Our improved methodology for deriving star formation histories (SFHs) of galaxies at redshifts 0.45 &lt; z &lt; 0.75 from the Carnegie-Spitzer-IMACS Survey includes confidence intervals that robustly distinguish LBs from “old” galaxies. We use simulated SFHs to test for “false positives” and contamination from old galaxies to demonstrate that the late-bloomer population is not an artifact of our template modeling technique. We show that LBs account for ∼20% of z ∼ 0.6 galaxies with masses of the modern Milky Way, with a moderate dependence on mass. We take advantage of a 1% overlap of our sample with HST (CANDELS) imaging to construct a “gold standard” catalog of 74 galaxies with high-confidence SFHs, SEDs, basic data, and HST images to facilitate comparison with future studies by others. This small subset suggests that galaxies with both old and young SFHs cover the full range of morphology and environment (excluding rich groups or clusters), albeit with a mild but suggestive correlation with the local environment. We begin the investigation of whether LBs of sufficient mass and frequency are produced in current-generation ΛCDM-based semianalytic models of galaxy formation. In terms of halo growth, we find a late-assembling halo fraction within a factor of two of our late bloomer fraction. However, sufficiently delaying star formation in such halos may be a challenge for the baryon component of such models.

The Mass-to-light Ratios and the Star Formation Histories of Disk Galaxies

We combine new data from the main sequence (M_* versus SFR) of star-forming galaxies and galaxy colors (from GALEX to Spitzer) with a flexible stellar population scheme to deduce the mass-to-light ratio (\Upsilon_*) of star-forming galaxies from the SPARC and S^4G samples. We find that the main sequence for galaxies, particular the low-mass end, combined with the locus of galaxy colors, constrains the possible star formation histories of disk and dwarf galaxies to a similar shape found by Speagle et al. (2014). Combining the deduced star formation history with stellar population models in the literature produces reliable \Upsilon_* values as a function of galaxy color with an uncertainty of only 0.05 dex. We provide prescriptions to deduce \Upsilon_* for optical and near-IR bandpasses, with near-IR bandpasses having the least uncertainty (\Upsilon_* from 0.40 to 0.55). We also provide the community with a webtool, with flexible stellar population parameters, to generate their own \Upsilon_* values over the wavelength range for most galaxy surveys.

Optical integral field spectroscopy observations applied to simulated galaxies: testing the fossil record method

By means of post-processed cosmological Hydrodynamics simulations we explore the ability of the fossil record method to recover the stellar mass, age gradients, and global/radial star formation and mass growth histories of galaxies observed with an optical Integral Field Unity. We use two simulations of representative Milky Way-sized galaxies, one disk dominated and another bulge-dominated. We generate data cubes emulating MaNGA/SDSS-IV instrumental and observational conditions, and analyze them using the Pipe3D pipeline. For optimal MaNGA instrumental/observational setups, the global masses, ages, and dust extinctions of both galaxies are well recovered within the uncertainties, while for the radial distributions, the trend is to get slightly flatter age gradients and weaker inside-out growth modes than the true ones. The main bias is in recovering younger stellar populations in the inner, older regions, and slightly older stellar populations in the outer, younger regions, as well as lower global stellar masses. These trends are enhanced as lower is the spatial sampling/resolution (less number of fibers), more inclined and dust-attenuated is the observed galaxy, lower is the signal-to-noise ratio, and worse is the seeing. The recovered properties and histories are more affected for the disk-dominated galaxy with a prominent inside-out growth mode than for the bulge-dominated galaxy assembled earlier and more coherently. We make publicly available the code and the mock IFU datacubes used in this paper.

The main sequence of star-forming galaxies – I. The local relation and its bending

By using a set of different SFR indicators, including WISE mid-infrared and Halpha emission, we study the slope of the Main Sequence (MS) of local star forming galaxies at stellar masses larger than 10^{10} M_{\odot}. The slope of the relation strongly depends on the SFR indicator used. In all cases, the local MS shows a bending at high stellar masses with respect to the slope obtained in the low mass regime. While the distribution of galaxies in the upper envelope of the MS is consistent with a log-normal distribution, the lower envelope shows an excess of galaxies, which increases as a function of the stellar mass but varies as a function of the SFR indicator used. The scatter of the best log-normal distribution increases with stellar mass from ~0.3 dex at 10^{10} M_{\odot} to ~0.45 at 10^{11} M_{\odot}. The MS high-mass end is dominated by central galaxies of group sized halos with a red bulge and a disk redder than the lower mass counterparts. We argue that the MS bending in this region is due to two processes: i) the formation of a bulge component as a consequence of the increased merger activity in groups, and ii) the cold gas starvation induced by the hot halo environment, which cuts off the gas inflow onto the disk. Similarly, the increase of the MS scatter at high stellar masses would be explained by the larger spread of star formation histories of central group and cluster galaxies with respect to lower mass systems.

The MOSDEF Survey: Stellar Continuum Spectra and Star Formation Histories of Active, Transitional, and Quiescent Galaxies at 1.4 &lt; z &lt; 2.6

Abstract Using the MOSFIRE Deep Evolution Field (MOSDEF) rest-frame optical spectroscopic survey, we investigate the star formation histories (SFHs) of different galaxy types, ranging from actively star-forming to quiescent at 1.4 ≤ z ≤ 2.6. SFHs are constrained utilizing stellar continuum spectroscopy, specifically through a combination of Balmer absorption lines, the 4000 Å break, and the equivalent width of the Hα emission line. To attain a sufficiently high signal-to-noise ratio (S/N) to conduct these measurements we stack spectra of galaxies with similar spectral types, as determined from their rest-frame U − V and V − J colors. We bin the MOSDEF sample into five spectral types, subdividing the quiescent and star-forming bins to better explore galaxies transitioning between the two. We constrain the average SFHs for each type, finding that quiescent and transitional galaxies in the MOSDEF sample are dominated by an SFH with an average star formation timescale of τ ∼ 0.1–0.2 Gyr. These findings contrast with measurements from the low-redshift Universe where, on average, galaxies form their stars over a more extended time period (τ &gt; 1 Gyr). Furthermore, our spectral index measurements correlate with mass surface density for all spectral types. Finally, we compare the average properties of the galaxies in our transitional bins to investigate possible paths to quiescence, and speculate on the viability of a dusty post-starburst phase.

Molecular gas in two companion cluster galaxies at z = 1.2

Context. Probing both star formation history and evolution of distant cluster galaxies is essential to evaluate the effect of dense environment on shaping the galaxy properties we observe today. Aims. We investigate the effect of cluster environment on the processing of the molecular gas in distant cluster galaxies. We study the molecular gas properties of two star-forming galaxies separated by 6 kpc in the projected space and belonging to a galaxy cluster selected from the Irac Shallow Cluster Survey, at a redshift z = 1.2, that is, ~ 2 Gyr after the cosmic star formation density peak. This work describes the first CO detection from 1 &lt; z &lt; 1.4 star-forming cluster galaxies with no clear reported evidence of active galactic nuclei. Methods. We exploit observations taken with the NOEMA interferometer at ~3 mm to detect CO(2−1) line emission from the two selected galaxies, unresolved by our observations. Results. Based on the CO(2−1) spectrum, we estimate a total molecular gas mass M(H2) = (2.2+0.50.4) × 1010 M⊙, where fully excited gas is assumed, and a dust mass Mdust &lt; 4.2 × 108 M⊙ for the two blended sources. The two galaxies have similar stellar masses and Hα-based star formation rates (SFRs) found in previous work, as well as a large relative velocity of ~400 km s−1 estimated from the CO(2−1) line width. These findings tend to privilege a scenario where both sources contribute to the observed CO(2−1). Using the archival Spitzer MIPS flux at 24 μm we estimate an SFR (24μm) = (28+12−8) M⊙/yr for each of the two galaxies. Assuming that the two sources contribute equally to the observed CO(2−1), our analysis yields a depletion timescale of τdep = (3.9+1.4−1.8) × 108 yr, and a molecular gas to stellar mass ratio of 0.17 ± 0.13 for each of two sources, separately. We also provide a new, more precise measurement of an unknown weighted mean of the redshifts of the two galaxies, z = 1.163 ± 0.001. Conclusions. Our results are in overall agreement with those of other distant cluster galaxies and with model predictions for main sequence (MS) field galaxies at similar redshifts. The two target galaxies have molecular gas mass and depletion times that are marginally compatible with, but smaller than those of MS field galaxies, suggesting that the molecular gas has not been sufficiently refueled. We speculate that the cluster environment might have played a role in preventing the refueling via environmental mechanisms such as galaxy harassment, strangulation, ram-pressure, or tidal stripping. Higher-resolution and higher-frequency observations will enable us to spatially resolve the two sources and possibly distinguish between different gas processing mechanisms.

A Systematic Analysis of Stellar Populations in the Host Galaxies of SDSS Type I QSOs

Abstract We investigate the relationship between host galaxies’ stellar content and active galactic nuclei (AGNs) for optically selected QSOs with z &lt; 0.5. There is a total of 82 QSOs that we select from the Sloan Digital Sky Survey. These 82 QSOs have both Wide-field Infrared Survey Explorer data and measurable stellar content. With the help of the stellar population synthesis code STARLIGHT, we determine the luminosity fractions of AGNs, stellar population ages, and star formation histories (SFHs) of host galaxies. We find that there is a correlation between the SFH and AGN properties, which suggests a possible delay from star formation to AGN. This probably indicates that the AGN activity correlates with the star formation activity, which consistent with a coevolution scheme for black hole and host galaxies.

Inferring the star formation histories of massive quiescent galaxies with bagpipes: evidence for multiple quenching mechanisms

We present Bayesian Analysis of Galaxies for Physical Inference and Parameter EStimation, or BAGPIPES, a new Python tool which can be used to rapidly generate complex model galaxy spectra and to fit these to arbitrary combinations of spectroscopic and photometric data using the MultiNest nested sampling algorithm. We extensively test our ability to recover realistic star-formation histories (SFHs) by fitting mock observations of quiescent galaxies from the MUFASA simulation. We then perform a detailed analysis of the SFHs of a sample of 9289 quiescent galaxies from UltraVISTA with stellar masses, $M_* > 10^{10}\ \mathrm{M_\odot}$ and redshifts $0.25 < z < 3.75$. The majority of our sample exhibit SFHs which rise gradually then quench relatively rapidly, over $1{-}2$ Gyr. This behaviour is consistent with recent cosmological hydrodynamic simulations, where AGN-driven feedback in the low-accretion (jet) mode is the dominant quenching mechanism. At $z > 1$ we also find a class of objects with SFHs which rise and fall very rapidly, with quenching timescales of $< 1$ Gyr, consistent with quasar-mode AGN feedback. Finally, at $z < 1$ we find a population with SFHs which quench more slowly than they rise, over $>3$ Gyr, which we speculate to be the result of diminishing overall cosmic gas supply. We confirm the mass-accelerated evolution (downsizing) trend, and a trend towards more rapid quenching at higher stellar masses. However, our results suggest that the latter is a natural consequence of mass-accelerated evolution, rather than a change in quenching physics with stellar mass. We find $61\pm8$ per cent of $z > 1.5$ massive quenched galaxies undergo significant further evolution by $z = 0.5$. BAGPIPES is available at https://bagpipes.readthedocs.io

Local Group Dwarf Galaxies

AbstractLocal Group dwarf galaxies are a unique astrophysical laboratory because they are the only objects in which we can reliably and precisely characterize the star formation histories of low-mass galaxies going back to the epoch of reionization. There are of order 100 known galaxies less massive than the Small Magellanic Cloud within ~1 Megaparsec of the Milky Way, with a vide variety of star formation history, gas content, and mass to light ratios. In this overview the current understanding of the formation and evolution of low-mass galaxies across cosmic time will be presented, and the possibility of drawing links between the properties of individual systems and the broader Local Group and cosmological context will be discussed. Local Group dwarfs will remain a uniquely powerful testbed to constrain the properties of dark matter and to evaluate the performance of simulations for the foreseeable future.

Planetary nebulae in the (extra)-galactic context: Probing chemical evolution in star-forming galaxies

AbstractThe populations of planetary nebulae (PNe) probe metallicity and chemical content (and its evolution) of the parent galaxy, giving clues to galaxy formation and evolution. This sub-field of extra-galactic PN research has been particularly active in the recent years. Comparison of data and models yielded estimates of global cosmic enrichment and provided constraints to galaxy formation history. In external spiral galaxies, the chemical contents of PNe and H II regions can be compared to disclose possible evolution of the radial metallicity gradient, which is, in turn, a powerful constraint to galactic chemical evolutionary models. In the Milky Way, recent PN progenitor dating and new chemical abundances offer an updated look into our own Galaxy. Collectively, Galactic and extra-galactic radial metallicity gradients from emission-line probes (PNe and H II regions) can be compared to have a cosmological outlook on galactic evolution.

Resolving Quiescent Galaxies at z ≳ 2. II. Direct Measures of Rotational Support

Abstract Stellar kinematics provide insights into the masses and formation histories of galaxies. At high redshifts, spatially resolving the stellar kinematics of quiescent galaxies is challenging due to their compact sizes. Using deep near-infrared spectroscopy, we have measured the resolved stellar kinematics of four quiescent galaxies at z = 1.95–2.64, introduced in Paper I, that are gravitationally lensed by galaxy clusters. Analyses of two of these have previously been reported individually by Newman et al. and Toft et al., and for the latter, we present new observations. All four galaxies show significant rotation and can be classified as “fast rotators.” In the three systems for which the lensing constraints permit a reconstruction of the source, we find that all are likely to be highly flattened (intrinsic ellipticities of ≈0.75–0.85) disk-dominated galaxies with rapid rotation speeds of V max = 290–352 km s−1 and predominantly rotational support, as indicated by the ratio Compared to coeval star-forming galaxies of similar mass, the quiescent galaxies have smaller V/σ. Given their high masses, , we argue that these galaxies are likely to evolve into “slow rotator” elliptical galaxies whose specific angular momentum is reduced by a factor of 5–10. This provides strong evidence for merger-driven evolution of massive galaxies after quenching. Consistent with indirect evidence from earlier morphological studies, our small but unique sample suggests that the kinematic transformations that produced round, dispersion-supported elliptical galaxies were not generally coincident with quenching. Such galaxies probably emerged later via mergers that increased their masses and sizes while also eroding their rotational support.

The formation of solar-neighbourhood stars in two generations separated by 5 billion years.

The chemical compositions of stars encode those of the gas from which they formed, providing important clues regarding the formation histories of galaxies. A powerful diagnostic is the abundance of α elements (O, Mg, Si, S, Ca and Ti) relative to iron, [α/Fe]. The α elements are synthesized and injected into the interstellar medium by type II supernovae, which occur about ten million years after their originating stars form; by contrast, iron is returned to the interstellar medium by type Ia supernovae, which occur after a much longer timescale of roughly one billion years1. Periods of rapid star formation therefore tend to produce high-[α/Fe] stellar populations (because only type II supernovae have time to contribute to interstellar-medium enrichment as the stellar population forms), whereas low-[α/Fe] stars require periods of star formation that last more than a few billion years (over which timescales type Ia supernovaebegin to affect the elemental composition of the interstellar medium more strongly than type II supernovae). The existence of two distinct groups of stars in the solar neighbourhood2-7, one with high [α/Fe] and the other with low [α/Fe], therefore suggests two different origins, but the mechanism by which this bimodal distribution arose remains unknown. Here we use a model of disk-galaxy evolution to show that the two episodes of star formation8 predicted by the 'cold flow' theory of galactic gas accretion9,10 also explain the observed chemical bimodality. In this scenario, the high-[α/Fe] stars form early, during an initial phase of accretion that involves infalling streams of cold primordial gas. There is then a hiatus of around two billion years until the shock-heated gas in the galactic dark-matter halo has cooled as a result of radiation and can itself commence accretion. The low-[α/Fe] stars form during this second phase. The peaks in these two star-formation episodes are separated by around five billion years. In addition, the large-scale variation in the abundance patterns of these two stellar populations that has been observed for the Milky Way5,7 is partially explained by the spatial variation in this gas-accretion history.

Stellar populations of HII galaxies

Aims. We present a UV to mid-IR spectral energy distribution (SED) study of a large sample of SDSS DR13 HII galaxies. These galaxies are selected as starbursts (EW(Hα) &gt; 50Å) and for their high-excitation locus in the upper left region of the BPT diagram. Their photometry was derived from the cross-matched GALEX, SDSS, UKDISS, and WISE catalogs. Methods. We used CIGALE modeling and a SED fitting routine with the parametrization of a three-burst star formation history, and a comprehensive analysis of all other model parameters. We were able to estimate the contribution of the underlying old stellar population to the observed equivalent width of Hβ, and allow for more accurate burst age determination. Results. We found that the star formation histories of HII Galaxies can be reproduced remarkably well by three major eras of star formation. In addition, the SED fitting results indicate that in all cases the current burst produces a small percent of the total stellar mass, i.e., the bulk of stellar mass in HII galaxies has been produced by the past episodes of star formation, and also indicate that at a given age the Hβ luminosity depends only on the mass of young stars favoring a universal IMF for massive stars. Most importantly, the current star formation episodes are maximum starbursts that produce stars at the highest possible rate.