Realistic Star Formation Histories Research Articles

Predicting the binary black hole population of the Milky Way with cosmological simulations

Binary black holes are the primary endpoint of massive stellar evolution. Their properties provide a unique opportunity to constrain binary evolution, which is still poorly understood. In this paper, we predict the inventory of binary black holes and their merger products in/around the Milky Way, and detail their main properties. We present the first combination of a high-resolution cosmological simulation of a Milky Way-mass galaxy with a binary population synthesis model. The hydrodynamic simulation, taken from the FIRE project, provides a cosmologically realistic star formation history for the galaxy and its stellar halo and satellites. We apply a metallicity-dependent evolutionary model to the star particles to produce individual binary black holes. We find that a million binary black holes have merged in the model Milky Way, and 3 million binaries are still present, with an average mass of 28 Msun per binary. Because the black hole progenitors are biased towards low metallicity stars, half reside in the stellar halo and satellites and 40 per cent of the binaries were formed outside the main galaxy. This trend increases with the masses of the black holes. The numbers and mass distribution of the merged systems is compatible with the LIGO/Virgo detections. Observations of these black holes will be challenging, both with electromagnetic methods and LISA. We find that a cosmologically realistic star formation history, with self-consistent metal enrichment and Galactic accretion history, are key ingredients for determining binary black hole rates that can be compared with observations to constrain massive binary evolution.

Analysis of the Intrinsic Mid-infrared L band to Visible–Near-infrared Flux Ratios in Spectral Synthesis Models of Composite Stellar Populations

Abstract We analyze the intrinsic flux ratios of various visible–near-infrared filters with respect to 3.5 μm for simple and composite stellar populations (CSPs), and their dependence on age, metallicity, and star formation history (SFH). UV/optical light from stars is reddened and attenuated by dust, where different sightlines across a galaxy suffer varying amounts of extinction. Tamura et al. (2009) developed an approximate method to correct for dust extinction on a pixel-by-pixel basis, dubbed the “ ” method, by comparing the observed flux ratio to an empirical estimate of the intrinsic ratio of visible and ∼3.5 μm data. Through extensive modeling, we aim to validate the “ ” method for various filters spanning the visible through near-infrared wavelength range, for a wide variety of simple and CSPs. Combining Starburst99 and BC03 models, we built spectral energy distributions (SEDs) of simple (SSP) and composite (CSP) stellar populations for various realistic SFHs, while taking metallicity evolution into account. We convolve various 0.44–1.65 μm filter throughput curves with each model SED to obtain intrinsic flux ratios . When unconstrained in redshift, the total allowed range of is 0.6–4.7, or almost a factor of eight. At known redshifts, and in particular at low redshifts (z ≲ 0.01), is predicted to span a narrow range of 0.6–1.9, especially for early-type galaxies (0.6–0.7), and is consistent with observed values. The method can therefore serve as a first-order dust-correction method for large galaxy surveys that combine JWST (rest-frame 3.5 μm) and HST (rest-frame visible–near-IR) data.

The Spectral Evolution of the First Galaxies. III. Simulated James Webb Space Telescope Spectra of Reionization-epoch Galaxies with Lyman-continuum Leakage

Abstract Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman-continuum escape fractions (f esc) at redshifts –9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hβ)–β diagram (previously proposed as a tool for identifying galaxies with very high f esc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman-continuum leakage. Based on our models, we expect that essentially all galaxies that exhibit rest-frame Å to have . Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of galaxies can in principle bias the selection, but substantial model deficiencies of this type should at the same time be evident from offsets in the observed distribution of galaxies in the EW(Hβ)–β diagram compared to the simulated distribution. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.

SIMULTANEOUS ESTIMATION OF PHOTOMETRIC REDSHIFTS AND SED PARAMETERS: IMPROVED TECHNIQUES AND A REALISTIC ERROR BUDGET

We seek to improve the accuracy of joint galaxy photometric redshift estimation and spectral energy distribution (SED) fitting. By simulating different sources of uncorrected systematic errors, we demonstrate that if the uncertainties on the photometric redshifts are estimated correctly, so are those on the other SED fitting parameters, such as stellar mass, stellar age, and dust reddening. Furthermore, we find that if the redshift uncertainties are over(under)-estimated, the uncertainties in SED parameters tend to be over(under)-estimated by similar amounts. These results hold even in the presence of severe systematics and provide, for the first time, a mechanism to validate the uncertainties on these parameters via comparison with spectroscopic redshifts. We propose a new technique (annealing) to re-calibrate the joint uncertainties in the photo-z and SED fitting parameters without compromising the performance of the SED fitting + photo-z estimation. This procedure provides a consistent estimation of the multidimensional probability distribution function in SED fitting + z parameter space, including all correlations. While the performance of joint SED fitting and photo-z estimation might be hindered by template incompleteness, we demonstrate that the latter is "flagged" by a large fraction of outliers in redshift, and that significant improvements can be achieved by using flexible stellar populations synthesis models and more realistic star formation histories. In all cases, we find that the median stellar age is better recovered than the time elapsed from the onset of star formation [abridged].

Constraining the properties of AGN host galaxies with spectral energy distribution modelling

[abridged] We use the latest release of CIGALE, a galaxy SED fitting model relying on energy balance, to study the influence of an AGN in estimating both the SFR and stellar mass in galaxies, as well as the contribution of the AGN to the power output of the host. Using the galaxy formation SAM GALFORM, we create mock galaxy SEDs using realistic star formation histories (SFH) and add an AGN of Type 1, Type 2, or intermediate type whose contribution to the bolometric luminosity can be variable. We perform an SED fitting of these catalogues with CIGALE assuming three different SFHs: a single- and double-exponentially-decreasing, and a delayed SFH. Constraining thecontribution of an AGN to the LIR (fracAGN) is very challenging for fracAGN<20%, with uncertainties of ~5-30% for higher fractions depending on the AGN type, while FIR and sub-mm are essential. The AGN power has an impact on the estimation of $M_*$ in Type 1 and intermediate type AGNs but has no effect for galaxies hosting Type 2 AGNs. We find that in the absence of AGN emission, the best estimates of $M_*$ are obtained using the double-exponentially-decreasing model but at the expense of realistic ages of the stellar population. The delayed SFH model provides good estimates of $M_*$ and SFR, with a maximum offset of 10% as well as better estimates of the age. Our analysis shows that the underestimation of the SFR increases with fracAGN for Type 1 systems, as well as for low contributions of an intermediate AGN type, but it is quite insensitive to the emission of Type 2 AGNs up to fracAGN~45%. Similarly the UV emission is critical in accurately retrieving the $M_*$ for Type 1 and intermediate type AGN, and the SFR of all of the three AGN types. We show that the presence of AGN emission introduces a scatter to the SFR-$M_*$ main sequence relation derived from SED fitting, which is driven by the uncertainties on $M_*$.

Binary White Dwarfs in the Galactic Halo

We use the stellar population synthesis code SeBa (Portegies Zwart &amp; Verbunt (1996), Toonen, Nelemans &amp; Portegies Zwart (2012)) to study the halo white dwarf population. Here we assume a Kroupa initial mass function and compare 4 models, varying two parameters: the star formation (SF) history of the halo (either continuous SF during 2.5 Gyr, which started 13.2 Gyr ago, or a SF burst during 360 Myr, which started 12.9 Gyr ago – see the left panel of the figure) and the binary fraction of the halo (either 100% single stars, or 100% binaries). White dwarf cooling models (Althaus et al. (2009) and Renedo et al. (2010)) allow us to plot the halo white dwarf luminosity function for these 4 models, as is done in the right panel of the figure. Combined with an assumption about the density distribution of halo stars, we will study which of these white dwarfs Gaia can see, and what that can tell us about the initial parameter distributions in the halo. In the near future, we plan to use the Munich-Groningen semi-analytical galaxy formation model (Starkenburg et al. (2013)), to obtain key ingredients for the population synthesis modeling, such as a realistic star formation history (see the left panel of the figure).

QUANTIFYING NON-STAR-FORMATION-ASSOCIATED 8 μm DUST EMISSION IN NGC 628

Combining Ha and IRAC images of the nearby spiral galaxy NGC 628, we find that between 30-43% of its 8um dust emission is not related to recent star formation. Contributions from dust heated by young stars are separated by identifying HII regions in the Ha map and using these areas as a mask to determine the 8um dust emission that must be due to heating by older stars. Corrections are made for sub-detection-threshold HII regions, photons escaping from HII regions and for young stars not directly associated to HII regions (i.e. 10-100 Myr old stars). A simple model confirms this amount of 8um emission can be expected given dust and PAH absorption cross-sections, a realistic star-formation history, and the observed optical extinction values. A Fourier power spectrum analysis indicates that the 8um dust emission is more diffuse than the Ha emission (and similar to observed HI), supporting our analysis that much of the 8um-emitting dust is heated by older stars. The 8um dust-to-Ha emission ratio declines with galactocentric radius both within and outside of HII regions, probably due to a radial increase in disk transparency. In the course of this work, we have also found that intrinsic diffuse Ha fractions may be lower than previously thought in galaxies, if the differential extinction between HII regions and diffuse regions is taken into account.

Star formation history and the SED of galaxies: insights from resolved stars

AbstractThe star formation history (SFH) of galaxies is a principle uncertainty in SED modeling, and simple parameterizations of the SFH in typical SED fitting techniques may introduce biases in the resulting derived parameters. It is possible to constrain the SFH of galaxies more tightly through the observations of resolved stellar colour-magnitude diagrams with HST. This work is a first attempt to combine constraints on galaxy SFH from resolved stars with broadband SED modeling from the UV to the IR. This combination allows for the effects of different realistic SFHs on the SED to be quantified.

INTEGRATED STELLAR POPULATIONS: CONFRONTING PHOTOMETRY WITH SPECTROSCOPY

We investigate the ability of spectroscopic techniques to yield realistic star formation histories (SFHs) for the bulges of spiral galaxies based on a comparison with their observed broadband colors. Full spectrum fitting to optical spectra indicates that recent (within ~1 Gyr) star formation activity can contribute significantly to the V-band flux, whilst accounting for only a minor fraction of the stellar mass budget which is made up primarily of old stars. Furthermore, recent implementations of stellar population (SP) models reveal that the inclusion of a more complete treatment of the thermally pulsating asymptotic giant branch (TP-AGB) phase to SP models greatly increases the NIR flux for SPs of ages 0.2-2 Gyr. Comparing the optical--NIR colors predicted from population synthesis fitting, using models which do not include all stages of the TP-AGB phase, to the observed colors reveals that observed optical--NIR colors are too red compared to the model predictions. However, when a 1 Gyr SP from models including a full treatment the TP-AGB phase is used, the observed and predicted colors are in good agreement. This has strong implications for the interpretation of stellar populations, dust content, and SFHs derived from colors alone.

Cosmological Gamma-Ray Bursts and Evolution of Galaxies

Evolution of the rate density of cosmological gamma-ray bursts (GRBs) is calculated and compared to the BATSE brightness distribution in the context of binary neutron-star mergers as the source of GRBs, taking account of the realistic star formation history in the universe and evolution of compact binary systems. We tried two models of the evolution of cosmic star formation rate (SFR): one is based on recent observations of SFRs at high redshifts, while the other is based on a galaxy evolution model of stellar population synthesis that reproduces the present-day colors of galaxies. It is shown that the binary merger scenario of GRBs naturally results in the comoving rate-density evolution of \propto (1+z)^{2-2.5} up to z ~ 1, that has been suggested independently from the compatibility between the number-brightness distribution and duration-brightness correlation. If the cosmic SFR has its peak at z ~ 1--2 as suggested by recent observations, the effective power-index of GRB photon spectrum, \alpha >~ 1.5$ is favored, that is softer than the recent observational determination of \alpha = 1.1 \pm 0.3. However, high redshift starbursts (z >~ 5) in elliptical galaxies, that have not yet been detected, can alleviate this discrepancy. The redshift of GRB970508 is likely about 2, just below the upper limit that is recently determined, and the absorption system at z = 0.835 seems not to be the site of the GRB.