The relationship between stellar populations and the ionizing flux with which they irradiate their surroundings has profound implications for the evolution of the intergalactic medium (IGM). We ...quantify the ionizing flux arising from synthetic stellar populations which incorporate the evolution of interacting binary stars. We determine that these show ionizing flux boosted by 60 per cent at 0.05 ≤ Z ≤ 0.3 Z⊙ and a more modest 10–20 per cent at near-solar metallicities relative to star-forming populations in which stars evolve in isolation. The relation of ionizing flux to observables such as 1500 Å continuum and ultraviolet spectral slope is sensitive to attributes of the stellar population including age, star formation history and initial mass function (IMF). For a galaxy forming 1 M⊙ yr−1, observed at >100 Myr after the onset of star formation, we predict a production rate of photons capable of ionizing hydrogen, N
ion = 1.4 × 1053 s−1 at Z = Z⊙ and 3.5 × 1053 s−1 at 0.1 Z⊙, assuming a Salpeter-like IMF. We evaluate the impact of these issues on the ionization of the IGM, finding that the known galaxy populations can maintain the ionization state of the Universe back to z ∼ 9, assuming that their luminosity functions continue to M
UV = −10, and that constraints on the IGM at z ∼ 2–5 can be satisfied with modest Lyman-continuum photon escape fractions of 4–24 per cent depending on assumed metallicity.
Abstract
We present numerical models of the nebular emission from H ii regions around young stellar populations over a range of compositions and ages. The synthetic stellar populations include both ...single stars and interacting binary stars. We compare these models to the observed emission lines of 254 H ii regions of 13 nearby spiral galaxies and 21 dwarf galaxies drawn from archival data. The models are created using the combination of the bpass (Binary Population and Spectral Synthesis) code with the photoionization code cloudy to study the differences caused by the inclusion of interacting binary stars in the stellar population. We obtain agreement with the observed emission line ratios from the nearby star-forming regions and discuss the effect of binary-star evolution pathways on the nebular ionization of H ii regions. We find that at population ages above 10 Myr, single-star models rapidly decrease in flux and ionization strength, while binary-star models still produce strong flux and high ${\rm \rm O\,{{\small III}}}$/H β ratios. Our models can reproduce the metallicity of H ii regions from spiral galaxies, but we find higher metallicities than previously estimated for the H ii regions from dwarf galaxies. Comparing the equivalent width of H β emission between models and observations, we find that accounting for ionizing photon leakage can affect age estimates for H ii regions. When it is included, the typical age derived for H ii regions is 5 Myr from single-star models, and up to 10 Myr with binary-star models. This is due to the existence of binary-star evolution pathways, which produce more hot Wolf–Rayet and helium stars at older ages. For future reference, we calculate new bpass binary maximal starburst lines as a function of metallicity, and for the total model population, and present these in Appendix A.
ABSTRACT We compare the impacts of uncertainties in both binary population synthesis models and the cosmic star formation history on the predicted rates of gravitational wave (GW) compact binary ...merger events. These uncertainties cause the predicted rates of GW events to vary by up to an order of magnitude. Varying the volume-averaged star formation rate density history of the Universe causes the weakest change to our predictions, while varying the metallicity evolution has the strongest effect. Double neutron star merger rates are more sensitive to assumed neutron star kick velocity than the cosmic star formation history. Varying certain parameters affects merger rates in different ways depending on the mass of the merging compact objects; thus some of the degeneracy may be broken by looking at all the event rates rather than restricting ourselves to one class of mergers.
Recalibrating the cosmic star formation history Wilkins, Stephen M; Lovell, Christopher C; Stanway, Elizabeth R
Monthly Notices of the Royal Astronomical Society,
12/2019, Letnik:
490, Številka:
4
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
The calibrations linking observed luminosities to the star formation rate (SFR) depend on the assumed stellar population synthesis model, initial mass function, star formation and metal ...enrichment history, and whether reprocessing by dust and gas is included. Consequently the shape and normalization of the inferred cosmic star formation history is sensitive to these assumptions. Using v2.2.1 of the Binary Population and Spectral Synthesis (bpass) model we determine a new set of calibration coefficients for the ultraviolet, thermal infrared, and hydrogen recombination lines. These ultraviolet and thermal infrared coefficients are 0.15–0.2 dex higher than those widely utilized in the literature while the H α coefficient is ∼0.35 dex larger. These differences arise in part due to the inclusion binary evolution pathways but predominantly reflect an extension in the IMF to 300 M⊙ and a change in the choice of reference metallicity. We use these new coefficients to recalibrate the cosmic star formation history, and find improved agreement between the integrated cosmic star formation history and the in situ measured stellar mass density as a function of redshift. However, these coefficients produce new tension between SFR densities inferred from the ultraviolet and thermal infrared and those from H α.
ABSTRACT
The detection rate of electromagnetic (EM) and gravitational wave (GW) transients is growing exponentially. As the accuracy of the transient rates will significantly improve over the coming ...decades, so will our understanding of their evolution through cosmic history. To this end, we present predicted rates for EM and GW transients over the age of the universe using Binary Population and Spectral Synthesis (bpass) results combined with four cosmic star formation histories (SFHs). These include a widely used empirical SFH of Madau & Dickinson and those from three cosmological simulations: MilliMillennium, EAGLE, and IllustrisTNG. We find that the choice of SFH changes our predictions: transients with short delay times are most affected by the star formation rate and change up to a factor of 2, while long delay time events tend to depend on the metallicity evolution of star formation and can change the predicted rate up to an order of magnitude. Importantly, we find that the cosmological simulations have very different metallicity evolution that cannot be reproduced by the widely used metallicity model of Langer & Norman, which impacts the binary black hole merger, stripped-envelope supernovae, and LGRBs in the local Universe most acutely. We recommend against using simple prescriptions for the metallicity evolution of the universe when predicting the rates of events that can have long delay times and that are sensitive to metallicity evolution.
The observable characteristics and subsequent evolution of young stellar populations is dominated by their massive stars. As our understanding of those massive stars and the factors affecting their ...evolution improves, so our interpretation of distant, unresolved stellar systems can also advance. As observations increasingly probe the distant Universe, and the rare low-metallicity starbursts nearby, so the opportunity arises for these two fields to complement one another and leads to an improved conception of both stars and galaxies. Here, we review the current state of the art in modeling of massive star-dominated stellar populations and discuss their applications and implications for interpreting the distant Universe. Our principal findings include the following:
Binary evolutionary pathways must be included to understand the stellar populations in early galaxies.
Observations constraining the extreme ultraviolet spectrum of early galaxies are showing that current models are incomplete. The best current guess is that some form of accretion onto compact remnants is required.
The evolution and fates of very massive stars, on the order of 100 M
and above, may be key to fully understand aspects of early galaxies.
Young, massive stars dominate the rest-frame ultraviolet (UV) spectra of star-forming galaxies. At high redshifts (z > 2), these rest-frame UV features are shifted into the observed-frame optical and ...a combination of gravitational lensing, deep spectroscopy and spectral stacking analysis allows the stellar population characteristics of these sources to be investigated. We use our stellar population synthesis code Binary Population and Spectral Synthesis (bpass) to fit two strong rest-frame UV spectral features in published Lyman-break galaxy spectra, taking into account the effects of binary evolution on the stellar spectrum. In particular, we consider the effects of quasi-homogeneous evolution (arising from the rotational mixing of rapidly rotating stars), metallicity and the relative abundance of carbon and oxygen on the observed strengths of He iiλ1640 Å and C ivλ1548, 1551 Å spectral lines. We find that Lyman-break galaxy spectra at z ∼ 2-3 are best fitted with moderately sub-solar metallicities, and with a depleted carbon-to-oxygen ratio. We also find that the spectra of the lowest metallicity sources are best fitted with model spectra in which the He ii emission line is boosted by the inclusion of the effect of massive stars being spun-up during binary mass transfer so these rapidly rotating stars experience quasi-homogeneous evolution.
We have constructed a new code to produce synthetic spectra of stellar populations that includes massive binaries. We have tested this code against the broad-band colours of unresolved young massive ...stellar clusters in nearby galaxies, the equivalent widths of the Red and Blue Wolf–Rayet bumps in star-forming Sloan Digital Sky Survey galaxies and the ultraviolet and optical spectra of the star-forming regions Tol-A and B in NGC 5398. In each case, we find a good agreement between our models and observations. We find that in general binary populations are bluer and have fewer red supergiants, and thus significantly less flux in the I band and at longer wavelengths, than single star populations. Also we find that Wolf–Rayet stars occur over a wider range of ages up to 107 yr in a stellar population including binaries, increasing the ultraviolet flux and Wolf–Rayet spectral features at later times. In addition, we find that nebula emission contributes significantly to these observed properties and must be considered when comparing stellar models with observations of unresolved stellar populations. We conclude that incorporation of massive stellar binaries can improve the agreement between observations and synthetic spectral synthesis codes, particularly for systems with young stellar populations.
ABSTRACT Young, massive stars dominate the rest-frame ultraviolet (UV) spectra of star-forming galaxies. At high redshifts (z > 2), these rest-frame UV features are shifted into the observed-frame ...optical and a combination of gravitational lensing, deep spectroscopy and spectral stacking analysis allows the stellar population characteristics of these sources to be investigated. We use our stellar population synthesis code Binary Population and Spectral Synthesis (bpass) to fit two strong rest-frame UV spectral features in published Lyman-break galaxy spectra, taking into account the effects of binary evolution on the stellar spectrum. In particular, we consider the effects of quasi-homogeneous evolution (arising from the rotational mixing of rapidly rotating stars), metallicity and the relative abundance of carbon and oxygen on the observed strengths of Heiiλ1640Å and Civλ1548, 1551Å spectral lines. We find that Lyman-break galaxy spectra atz 2-3 are best fitted with moderately sub-solar metallicities, and with a depleted carbon-to-oxygen ratio. We also find that the spectra of the lowest metallicity sources are best fitted with model spectra in which the Heii emission line is boosted by the inclusion of the effect of massive stars being spun-up during binary mass transfer so these rapidly rotating stars experience quasi-homogeneous evolution. PUBLICATION ABSTRACT
ABSTRACT
Panchromatic analysis of galaxy spectral energy distributions, spanning from the ultraviolet to the far-infrared, probes not only the stellar population but also the properties of ...interstellar dust through its extinction and long-wavelength re-emission. However, little work has exploited the full power of such fitting to constrain the redshift evolution of dust temperature in galaxies. To do so, we simultaneously fit ultraviolet, optical, and infrared observations of stacked galaxy subsamples at a range of stellar masses and photometric redshifts at 0 < z < 5, using an energy-balance formalism. However, we find UV-emission beyond the Lyman limit in some photometric redshift selected galaxy subsamples, giving rise to the possibility of contaminated observations. We carefully define a robust clean subsample which extends to no further than z ∼ 2. This has consistently lower derived temperatures by $4.0^{+5.0}_{-1.9}$ K, relative to the full sample. We find a linear increase in dust temperature with redshift, with Td (z) = (4.8 ± 1.5) × z + (26.2 ± 1.5) K. Our inferred temperature evolution is consistent with a modest rise in dust temperature with redshift, but inconsistent with some previous analyses. We also find a majority of photometrically selected subsamples at z > 4.5 underpredict the IR emission while giving reasonable fits to the UV-optical. This could be due to a spatial disconnect in the locations of the UV and IR emission peaks, suggesting that an energy-balance formalism may not always be applicable in the distant Universe.