Population Synthesis Models Research Articles

Constructing binary neutron star initial data with high spins, high compactnesses, and high mass ratios

The construction of accurate and consistent initial data for various binary parameters is a critical ingredient for numerical relativity simulations of the compact binary coalescence. In this article, we present an upgrade of the pseudospectral SGRID code, which enables us to access even larger regions of the binary neutron star parameter space. As a proof of principle, we present a selected set of first simulations based on initial configurations computed with the new code version. In particular, we simulate two millisecond pulsars close to their breakup spin, highly compact neutron stars with masses at about $98\%$ of the maximum supported mass of the employed equation of state, and an unequal mass systems with mass ratios even outside the range predicted by population synthesis models ($q = 2.03$). The discussed code extension will help us to simulate previously unexplored binary configurations. This is a necessary step to construct and test new gravitational wave approximants and to interpret upcoming binary neutron star merger observations. When we construct initial data, one has to specify various parameters, such as a rotation parameter for each star. Some of these parameters do not have direct physical meaning, which makes comparisons with other methods or models difficult. To facilitate this, we introduce simple estimates for the initial spin, momentum, mass, and center of mass of each individual star.

The Exoplanet Population Observation Simulator. II. Population Synthesis in the Era of Kepler

Abstract The collection of planetary system properties derived from large surveys such as Kepler provides critical constraints on planet formation and evolution. These constraints can only be applied to planet formation models, however, if the observational biases and selection effects are properly accounted for. Here we show how epos, the Exoplanet Population Observation Simulator, can be used to constrain planet formation models by comparing the Bern planet population synthesis models to the Kepler exoplanetary systems. We compile a series of diagnostics, based on occurrence rates of different classes of planets and the architectures of multiplanet systems within 1 au, that can be used as benchmarks for future and current modeling efforts. Overall, we find that a model with 100-seed planetary cores per protoplanetary disk provides a reasonable match to most diagnostics. Based on these diagnostics we identify physical properties and processes that would result in the Bern model more closely matching the known planetary systems. These are as follows: moving the planet trap at the inner disk edge outward; increasing the formation efficiency of mini-Neptunes; and reducing the fraction of stars that form observable planets. We conclude with an outlook on the composition of planets in the habitable zone, and highlight that the majority of simulated planets smaller than 1.7 Earth radii in this zone are predicted to have substantial hydrogen atmospheres. The software used in this paper is available online for public scrutiny at https://github.com/GijsMulders/epos.

Binary population synthesis models for core-collapse gamma-ray burst progenitors

ABSTRACT Long-duration gamma-ray bursts (GRBs) are understood to be the final fate for a subset of massive, stripped envelope, rapidly rotating stars. Beyond this, our knowledge of the progenitor systems is limited. Using the Binary Population and Spectral Synthesis (bpass) stellar evolution models, we investigate the possibility that some massive stars in binaries can maintain the angular momentum required for jet production, while still loosing their outer envelope through winds or binary interactions. We find that a total hydrogen mass of MH < 5 × 10−4 M⊙ and a helium ejecta mass fraction of FHe < 0.20 provide the best thresholds for the supernova type II/Ibc and Ib/Ic divisions, respectively. Tidal interactions in binaries are accounted for by applying a tidal algorithm to post-process the stellar evolution models output by bpass. We show that the observed volumetric GRB rate evolution can be recreated using two distinct pathways and plausible distributions for burst parameters. In the first pathway, stars are spun up by mass accretion into a quasi-homogeneous state. In the second, tides maintain rotation where otherwise the star would spin-down. Both lead to type Ic supernova progenitors, and a metallicity distribution consistent with the GRB host galaxy population. The inferred core angular momentum threshold for jet production is consistent with theoretical requirements for collapsars, given the assumptions made in our model. We can therefore reproduce several aspects of core-collapse supernova/GRB observation and theory simultaneously. We discuss the predicted observable properties of GRB progenitors and their surviving companions.

Asymptotic kinematics of globular clusters: The emergence of a Tully–Fisher relation

ABSTRACT Using a recent homogeneous sample of 40 high-quality velocity dispersion profiles for Galactic globular clusters, we study the regime of low gravitational acceleration relevant to the outskirts of these systems. We find that a simple empirical profile having a central Gaussian component and a constant large-radius asymptote, σ∞, accurately describes the variety of observed velocity dispersion profiles. We use published population synthesis models, carefully tailored to each individual cluster, to estimate mass-to-light ratios from which total stellar masses, M, are inferred. We obtain a clear scaling, reminiscent of the galactic Tully–Fisher relation of $\sigma _{\infty }[\, \mathrm{km \, s}^{-1}]= 0.084^{+0.075}_{-0.040} (\mathrm{{\it M}/M}_{\odot })^{0.3 \pm 0.051}$, which is interesting to compare to the deep modified Newtonian dynamics (MOND) limit of $\sigma _{\infty } [\mathrm{km \, s}^{-1}]=0.2(\mathrm{{\it M}/M}_{\odot })^{0.25}$. Under a Newtonian interpretation, our results constitute a further restriction on models where initial conditions are crafted to yield the outer flattening observed today. Within a modified gravity scheme, because the globular clusters studied are not isolated objects in the deep MOND regime, the results obtained point towards a modified gravity where the external field effect of MOND does not appear, or is significantly suppressed.

Evolutionary Population Synthesis model with binary stars – Yunnan-II model

AbstractBy considering a modified version of the evolutionary population synthesis (EPS) model for stellar populations (SPs) comprising binary stars, the retrieved galaxy and HII-region parameters/properties differ from the case of neglecting binary stars. The retrieved age, stellar metallicity and mass of galaxies increase (e.g. ∼ 0.2 dex when using spectral fitting algorithm), whilst the star formation rate decreases (∼0.2 dex). The radiation fields from intermediate-age SPs with binary stars could be potentially important ionizing sources in HII regions. Under this possibility, the theoretical division between star forming galaxy and AGN on the diagnostic diagrams would move towards the up-right corner and the retrieved gaseous metallicity would decrease.Our prediction for the birth rate of binary neutron stars in SPs ranges from 10−9 to 10−6${\M {^\minus 1_\odot}} $ yr−1 when the kick velocity is from 0 to 190 km s−1.

Mentari: A pipeline to model the galaxy SED using semi analytic models

AbstractWe build a theoretical picture of how the light from galaxies evolves across cosmic time. In particular, we predict the evolution of the galaxy spectral energy distribution (SED) by carefully integrating the star formation and metal enrichment histories of semi-analytic model (SAM) galaxies and combining these with stellar population synthesis models which we call mentari. Our SAM combines prescriptions to model the interplay between gas accretion, star formation, feedback process, and chemical enrichment in galaxy evolution. From this, the SED of any simulated galaxy at any point in its history can be constructed and compared with telescope data to reverse engineer the various physical processes that may have led to a particular set of observations. The synthetic SEDs of millions of simulated galaxies from mentari can cover wavelengths from the far UV to infrared, and thus can tell a near complete story of the history of galaxy evolution.

The MOSDEF Survey: The Metallicity Dependence of X-Ray Binary Populations at z ∼ 2

Abstract Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low-metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity (L X) per star formation rate (SFR) with redshift due to the decrease of metallicity (Z) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 z ∼ 2 star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for ∼1500 galaxies in the CANDELS fields at 1.37 < z < 3.80. Using Chandra data from the Chandra AEGIS-X Deep, Chandra Deep Field North, and Chandra Deep Field South surveys, we stack the X-ray data at the galaxy locations in bins of redshift and Z because the galaxies are too faint to be individually detected. In agreement with previous studies, the average L X/SFR of our z ∼ 2 galaxy sample is enhanced by ≈0.4–0.8 dex relative to local HMXB L X–SFR scaling relations. Splitting our sample by Z, we find that L X/SFR and Z are anticorrelated with 97% confidence. This observed Z dependence for HMXB-dominated galaxies is consistent with both the local L X–SFR–Z relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak owing to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of L X/SFR and the Z dependence of HMXBs.

To use or not to use synthetic stellar spectra in population synthesis models?

ABSTRACT Stellar population synthesis (SPS) models are invaluable to study star clusters and galaxies. They provide means to extract stellar masses, stellar ages, star formation histories, chemical enrichment, and dust content of galaxies from their integrated spectral energy distributions, colours, or spectra. As most models, they contain uncertainties that can hamper our ability to model and interpret observed spectra. This work aims at studying a specific source of model uncertainty: the choice of an empirical versus a synthetic stellar spectral library. Empirical libraries suffer from limited coverage of parameter space, while synthetic libraries suffer from modelling inaccuracies. Given our current inability to have both ideal stellar-parameter coverage with ideal stellar spectra, what should one favour: better coverage of the parameters (synthetic library) or better spectra on a star-by-star basis (empirical library)? To study this question, we build a synthetic stellar library mimicking the coverage of an empirical library, and SPS models with different choices of stellar library tailored to these investigations. Through the comparison of model predictions and the spectral fitting of a sample of nearby galaxies, we learned that predicted colours are more affected by the coverage effect than the choice of a synthetic versus empirical library; the effects on predicted spectral indices are multiple and defy simple conclusions; derived galaxy ages are virtually unaffected by the choice of the library, but are underestimated when SPS models with limited parameter coverage are used; metallicities are robust against limited HRD coverage, but are underestimated when using synthetic libraries.

Recalibrating the cosmic star formation history

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 α.

The 50–100 pc scale parent stellar populations of Type II supernovae and limitations of single star evolution models

ABSTRACT There is observational evidence of a dearth in core-collapse supernova (ccSN) explosions from stars with zero-age main-sequence (ZAMS) mass M0 ≈ 17–30M⊙, referred to as the ‘red supergiant problem’. However, simulations now predict that above 20 M⊙ we should indeed only expect stars within certain pockets of M0 to produce a visible SN explosion. Validating these predictions requires large numbers of ccSNe of different types with measured M0, which is challenging. In this paper, we explore the reliability of using host galaxy emission lines and the H α equivalent width to constrain the age, and thus the M0 of ccSNe progenitors. We use Binary Population and Spectral Synthesis models to infer a stellar population age from MUSE observations of the ionized gas properties and H α EW at the location of eleven ccSNe with reliable M0 measurements. Comparing our results to published M0 values, we find that models that do not consider binary systems yield stellar ages that are systematically too young (thus M0 too large), whereas accounting for binary system interactions typically overpredict the stellar age (thus underpredict M0). Taking into account the effects of photon leakage bring our M0 estimates in much closer agreement with expectations. These results highlight the need for careful modelling of diffuse environments, such as are present in the vicinity of Type II SNe, before ionized emission line spectra can be used as reliable tracers of progenitor stellar age.

The K2-HERMES Survey: age and metallicity of the thick disc

ABSTRACT Asteroseismology is a promising tool to study Galactic structure and evolution because it can probe the ages of stars. Earlier attempts comparing seismic data from the Kepler satellite with predictions from Galaxy models found that the models predicted more low-mass stars compared to the observed distribution of masses. It was unclear if the mismatch was due to inaccuracies in the Galactic models, or the unknown aspects of the selection function of the stars. Using new data from the K2 mission, which has a well-defined selection function, we find that an old metal-poor thick disc, as used in previous Galactic models, is incompatible with the asteroseismic information. We use an importance-sampling framework, which takes the selection function into account, to fit for the metallicities of a population synthesis model using spectroscopic data. We show that spectroscopic measurements of [Fe/H] and [α/Fe] elemental abundances from the GALAH survey indicate a mean metallicity of log (Z/Z⊙) = −0.16 for the thick disc. Here Z is the effective solar-scaled metallicity, which is a function of [Fe/H] and [α/Fe]. With the revised disc metallicities, for the first time, the theoretically predicted distribution of seismic masses show excellent agreement with the observed distribution of masses. This indirectly verifies that the asteroseismic mass scaling relation is good to within five per cent. Assuming the asteroseismic scaling relations are correct, we estimate the mean age of the thick disc to be about 10 Gyr, in agreement with the traditional idea of an old α-enhanced thick disc.

Constraining disk evolution prescriptions of planet population synthesis models with observed disk masses and accretion rates

While planets are commonly discovered around main-sequence stars, the processes leading to their formation are still far from being understood. Current planet population synthesis models, which aim to describe the planet formation process from the protoplanetary disk phase to the time exoplanets are observed, rely on prescriptions for the underlying properties of protoplanetary disks where planets form and evolve. The recent development in measuring disk masses and disk-star interaction properties, i.e., mass accretion rates, in large samples of young stellar objects demand a more careful comparison between the models and the data. We performed an initial critical assessment of the assumptions made by planet synthesis population models by looking at the relation between mass accretion rates and disk masses in the models and in the currently available data. We find that the currently used disk models predict mass accretion rate in line with what is measured, but with a much lower spread of values than observed. This difference is mainly because the models have a smaller spread of viscous timescales than what is needed to reproduce the observations. We also find an overabundance of weakly accreting disks in the models where giant planets have formed with respect to observations of typical disks. We suggest that either fewer giant planets have formed in reality or that the prescription for planet accretion predicts accretion on the planets that is too high. Finally, the comparison of the properties of transition disks with large cavities confirms that in many of these objects the observed accretion rates are higher than those predicted by the models. On the other hand, PDS70, a transition disk with two detected giant planets in the cavity, shows mass accretion rates well in line with model predictions.

ROME/REA: A Gravitational Microlensing Search for Exoplanets Beyond the Snow Line on a Global Network of Robotic Telescopes

Planet population synthesis models predict an abundance of planets with semimajor axes between 1 and 10 au, yet they lie at the edge of the detection limits of most planet finding techniques. Discovering these planets and studying their distribution is critical to understanding the physical processes that drive planet formation. ROME/REA is a gravitational microlensing project whose main science driver is to discover exoplanets in the cold outer regions of planetary systems. To achieve this, it uses a novel approach combining a multiband survey with reactive follow-up observations, exploiting the unique capabilities of the Las Cumbres Observatory global network of robotic telescopes combined with a Target and Observation Manager system. We present the main science objectives and a technical overview of the project, including initial results.

The Star Formation History and Dynamics of the Ultra-diffuse Galaxy Dragonfly 44 in MOND and MOG

Abstract The observed line-of-sight velocity dispersion σ los of the ultra-diffuse galaxy Dragonfly 44 (DF44) requires a Newtonian dynamical mass-to-light ratio of solar units. This is well outside the acceptable limits of our stellar population synthesis (SPS) models, which we construct using the integrated galactic initial mass function (IGIMF) theory. Assuming DF44 is in isolation and using Jeans analysis, we calculate σ los profiles of DF44 in Milgromian dynamics (MOND) and modified gravity (MOG) theories without invoking dark matter. Comparing with the observed kinematics, the best-fitting MOND model has and a constant orbital anisotropy of . In MOG, we first fix its two theoretical parameters α and μ based on previous fits to the observed rotation curve data of The HI Nearby Galaxy Survey (THINGS). The DF44 σ los profile is best fit with , larger than plausible SPS values. MOG produces a σ los profile for DF44 with acceptable M dyn/L I and isotropic orbits if α and μ are allowed to vary. MOND with the canonical a 0 can explain DF44 at the 2.40σ confidence level (1.66%) if considering both its observed kinematics and typical star formation histories in an IGIMF context. However, MOG is ruled out at 5.49σ (P-value of 4.07 × 10−8) if its free parameters are fixed at the highest values consistent with THINGS data.

Predicting the LISA white dwarf binary population in the Milky Way with cosmological simulations

ABSTRACT White dwarf binaries with orbital periods below 1 h will be the most numerous sources for the space-based gravitational wave detector Laser Interferometer Space Antenna (LISA). Based on thousands of individually resolved systems, we will be able to constrain binary evolution and provide a new map of the Milky Way and its close surroundings. In this paper we predict the main properties of populations of different types of detached white dwarf binaries detected by LISA over time. For the first time, we combine a high-resolution cosmological simulation of a Milky Way-mass galaxy (taken from the FIRE project) with a binary population synthesis model for low- and intermediate-mass stars. Our Galaxy model therefore provides a cosmologically realistic star formation and metallicity history for the Galaxy and naturally produces its different components such as the thin and thick disc, the bulge, the stellar halo, and satellite galaxies and streams. Thanks to the simulation, we show how different Galactic components contribute differently to the gravitational wave signal, mostly due to their typical age and distance distributions. We find that the dominant LISA sources will be He–He double white dwarfs (DWDs) and He–CO DWDs with important contributions from the thick disc and bulge. The resulting sky map of the sources is different from previous models, with important consequences for the searches for electromagnetic counterparts and data analysis. We also emphasize that much of the science-enabling information regarding white dwarf binaries, such as the chirp mass and the sky localization, becomes increasingly rich with long observations, including an extended mission up to 8 yr.

The rate of planet–star coalescences due to tides and stellar evolution

ABSTRACT Orbits of close-in planets can shrink significantly due to dissipation of tidal energy in a host star. This process can result in star–planet coalescence within the Galactic lifetime. In some cases, such events can be accompanied by an optical or/and UV/X-ray transient. Potentially, these outbursts can be observed in near future with new facilities such as LSST from distances about few Mpc. We use a population synthesis model to study this process and derive the rate of star–planet mergers of different types. Mostly, planets are absorbed by red giants. However, these events, happening with the rate about 3 per year, mostly do not produce detectable transients. The rate of mergers with main sequence stars depends on the effectiveness of tidal dissipation; for reasonable values of stellar tidal quality factor, such events happen in a Milky Way-like galaxy approximately once in 70 yr or more rarely. This rate is dominated by planets with low masses. Such events do not produce bright transients having maximum luminosities ≲ 1036.5 erg s−1. Brighter events, related to massive planets, with maximum luminosity ∼1037.5–1038 erg s−1, have the rate nearly five times smaller.

GABE: Galaxy Assembly with Binary Evolution

We developed a new semi-analytic galaxy formation model: Galaxy Assembly with Binary Evolution (GABE). For the first time, we introduce binary evolution into semi-analytic models of galaxy formation by using the Yunnan-II stellar population synthesis model, which includes various binary interactions. When implementing our galaxy formation model onto the merger trees extracted from the Millennium simulation, it can reproduce a large body of observational results. We find that in the local universe, the model including binary evolution reduces the luminosity at optical and infrared wavelengths slightly, but it increases the luminosity at ultraviolet wavelengths significantly, especially in FUV band. The resulting luminosity function does not change very much over SDSS optical bands and infrared band, but the predicted colors are bluer, especially when the FUV band is under consideration. The new model allows us to explore the physics of various high energy events related to the remnants of binary stars, such as type Ia supernovae, short gamma-ray bursts and gravitational wave events, and their relation with host galaxies in a cosmological context.

Stroopwafel: simulating rare outcomes from astrophysical populations, with application to gravitational-wave sources

ABSTRACT Gravitational-wave observations of double compact object (DCO) mergers are providing new insights into the physics of massive stars and the evolution of binary systems. Making the most of expected near-future observations for understanding stellar physics will rely on comparisons with binary population synthesis models. However, the vast majority of simulated binaries never produce DCOs, which makes calculating such populations computationally inefficient. We present an importance sampling algorithm, stroopwafel, that improves the computational efficiency of population studies of rare events, by focusing the simulation around regions of the initial parameter space found to produce outputs of interest. We implement the algorithm in the binary population synthesis code compas, and compare the efficiency of our implementation to the standard method of Monte Carlo sampling from the birth probability distributions. stroopwafel finds ∼25–200 times more DCO mergers than the standard sampling method with the same simulation size, and so speeds up simulations by up to two orders of magnitude. Finding more DCO mergers automatically maps the parameter space with far higher resolution than when using the traditional sampling. This increase in efficiency also leads to a decrease of a factor of ∼3–10 in statistical sampling uncertainty for the predictions from the simulations. This is particularly notable for the distribution functions of observable quantities such as the black hole and neutron star chirp mass distribution, including in the tails of the distribution functions where predictions using standard sampling can be dominated by sampling noise.

Resolved and Integrated Stellar Masses in the SDSS-IV/MaNGA Survey. II. Applications of PCA-based Stellar Mass Estimates

Abstract A galaxy’s stellar mass is one of its most fundamental properties, but it remains challenging to measure reliably. With the advent of very large optical spectroscopic surveys, efficient methods that can make use of low signal-to-noise spectra are needed. With this in mind, we created a new software package for estimating effective stellar mass-to-light ratios that uses a principal component analysis (PCA) basis set to optimize the comparison between observed spectra and a large library of stellar population synthesis models. In Paper I, we showed that with a set of six PCA basis vectors we could faithfully represent most optical spectra from the Mapping Nearby Galaxies at APO (MaNGA) survey, and we tested the accuracy of our M/L estimates using synthetic spectra. Here, we explore sources of systematic error in our mass measurements by comparing our new measurements to data from the literature. We compare our stellar mass surface density estimates to kinematics-derived dynamical mass surface density measurements from the DiskMass Survey and find some tension between the two that could be resolved if the disk scale heights used in the kinematic analysis were overestimated by a factor of ∼1.5. We formulate an aperture-corrected stellar mass catalog for the MaNGA survey, and compare to previous stellar mass estimates based on multiband optical photometry, finding typical discrepancies of 0.1 dex. Using the spatially resolved MaNGA data, we evaluate the impact of estimating total stellar masses from spatially unresolved spectra, and we explore how the biases that result from unresolved spectra depend upon the galaxy’s dust extinction and star formation rate. Finally, we describe an SDSS Value-Added Catalog that will include both spatially resolved and total (aperture-corrected) stellar masses for MaNGA galaxies.

Constraints on the production and escape of ionizing radiation from the emission-line spectra of metal-poor star-forming galaxies

ABSTRACT We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and density-bounded H ii regions, active-galactic-nucleus (AGN) narrow-line regions, and radiative shocks computed all using the same physically consistent description of element abundances and depletion on to dust grains down to very low metallicities. We compare these models with a reference sample of metal-poor star-forming galaxies and Lyman-continuum (LyC) leakers at various redshifts, which allows the simultaneous exploration of more spectral diagnostics than typically available at once for individual subsamples. We confirm that current single- and binary-star population synthesis models do not produce hard-enough radiation to account for the high-ionization emission of the most metal-poor galaxies. Introducing either an AGN or radiative-shock component brings models into agreement with observations. A published model including X-ray binaries is an attractive alternative to reproduce the observed rise in He iiλ4686/H β ratio with decreasing oxygen abundance in metal-poor star-forming galaxies, but not the high observed He iiλ4686/Hβ ratios of galaxies with large EW(Hβ). A source of harder ionizing radiation appears to be required in these extreme objects, such as an AGN or radiative-shock component, perhaps linked to an initial-mass-function bias towards massive stars at low metallicity. This would also account for the surprisingly high [O i]/[O iii] ratios of confirmed LyC leakers relative to ionization-bounded models. We find no simple by-eye diagnostic of the nature of ionizing sources and the escape of LyC photon, which require proper simultaneous fits of several lines to be discriminated against.