The distributions of the initial main-sequence binary parameters are one of the key ingredients in obtaining evolutionary predictions for compact binary (BH–BH/BH–NS/NS–NS) merger rates. Until now, ...such calculations were done under the assumption that initial binary parameter distributions were independent. For the first time, we implement empirically derived inter-correlated distributions of initial binary parameters primary mass (M1), mass ratio (q), orbital period (P), and eccentricity (e). Unexpectedly, the introduction of inter-correlated initial binary parameters leads to only a small decrease in the predicted merger rates by a factor of ≲2–3 relative to the previously used non-correlated initial distributions. The formation of compact object mergers in the isolated classical binary evolution favours initial binaries with stars of comparable masses (q ≈ 0.5–1) at intermediate orbital periods (log P (days) = 2–4). New distributions slightly shift the mass ratios towards lower values with respect to the previously used flat q distribution, which is the dominant effect decreasing the rates. New orbital periods (∼1.3 more initial systems within log P (days) = 2–4), together with new eccentricities (higher), only negligibly increase the number of progenitors of compact binary mergers. Additionally, we discuss the uncertainty of merger rate predictions associated with possible variations of the massive-star initial mass function (IMF). We argue that evolutionary calculations should be normalized to a star formation rate (SFR) that is obtained from the observed amount of UV light at wavelength 1500 Å (an SFR indicator). In this case, contrary to recent reports, the uncertainty of the IMF does not affect the rates by more than a factor of ∼2. Any change to the IMF slope for massive stars requires a change of SFR in a way that counteracts the impact of IMF variations on compact object merger rates. In contrast, we suggest that the uncertainty in cosmic SFR at low metallicity can be a significant factor at play.
Abstract
In this paper, we study the impact of different galaxy statistics and empirical metallicity scaling relations on the merging rates and properties of compact object binaries. Firstly, we ...analyze the similarities and differences of using the star formation rate functions versus stellar mass functions as galaxy statistics for the computation of cosmic star formation rate density. We then investigate the effects of adopting the Fundamental Metallicity Relation versus a classic Mass Metallicity Relation to assign metallicity to galaxies with given properties. We find that when the Fundamental Metallicity Relation is exploited, the bulk of the star formation occurs at relatively high metallicities, even at high redshift; the opposite holds when the Mass Metallicity Relation is employed, since in this case the metallicity at which most of the star formation takes place strongly decreases with redshift. We discuss the various reasons and possible biases giving rise to this discrepancy. Finally, we show the impact of these different astrophysical prescriptions on the merging rates and properties of compact object binaries; specifically, we present results for the redshift-dependent merging rates and for the chirp mass and time delay distributions of the merging binaries.
The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in ...an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters, and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10−2 yr−1 from binary stars, 5 × 10−5 yr−1 from globular clusters, and 10−5 yr−1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc3). These models are thus in tension with the detection of GW170817 with an observed rate of 1.5−1.2+3.2$1.5^{+3.2}_{-1.2}$1.5−1.2+3.2 yr−1 (per 100 Mpc3; LIGO/Virgo 90% credible limits). Our results imply that either the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; that physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play.
All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported following the O1/O2 runs have near-zero effective spins. There are only three potential explanations for this. If the BH spin ...magnitudes are large, then: (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Then there is also the possibility that (iii) the BH spin magnitudes are small. We consider the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, including revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can simultaneously match the observed BH-BH merger rate density and BH masses and BH-BH effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be key for obtaining an improved reproduction of the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely to be overestimated if the merger GW170729 hosts a BH more massive than 50
M
⊙
. We also estimate rates of black hole-neutron star (BH-NS) mergers from recent LIGO/Virgo observations. If, in fact. angular momentum transport in massive stars is efficient, then any (electromagnetic or gravitational wave) observation of a rapidly spinning BH would indicate either a very effective tidal spin up of the progenitor star (homogeneous evolution, high-mass X-ray binary formation through case A mass transfer, or a spin- up of a Wolf-Rayet star in a close binary by a close companion), significant mass accretion by the hole, or a BH formation through the merger of two or more BHs (in a dense stellar cluster).
Recent observational and theoretical studies indicate that the stellar initial mass function (IMF) varies systematically with the environment (star formation rate – SFR, metallicity). Although the ...exact dependence of the IMF on those properties is likely to change with improving observational constraints, the reported trend in the shape of the IMF appears robust. We present the first study aiming to evaluate the effect of the IMF variations on the measured cosmic SFR density (SFRD) as a function of metallicity and redshift,
f
SFR
(
Z
,
z
). We also study the expected number and metallicity of white dwarf, neutron star, and black hole progenitors under different IMF assumptions. Applying the empirically driven IMF variations described by the integrated galactic IMF (IGIMF) theory, we revise
f
SFR
(
Z
,
z
) obtained in our previous study that assumed a universal IMF. We find a lower SFRD at high redshifts and a higher fraction of metal-poor stars being formed than previously determined. In the local Universe, our calculation applying the IGIMF theory suggests more white dwarf and neutron star progenitors in comparison with the universal IMF scenario, while the number of black hole progenitors remains unaffected.
Abstract
New observational facilities are probing astrophysical transients such as stellar explosions and gravitational-wave sources at ever-increasing redshifts, while also revealing new features in ...source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation history (
(
Z
,
z
)
) as an input. Large uncertainties remain in the shape and evolution of this function. In this work, we propose a simple analytical function for
(
Z
,
z
)
. Variations of this function can be easily interpreted because the parameters link to its shape in an intuitive way. We fit our analytical function to the star-forming gas of the cosmological TNG100 simulation and find that it is able to capture the main behavior well. As an example application, we investigate the effect of systematic variations in the
(
Z
,
z
)
parameters on the predicted mass distribution of locally merging binary black holes. Our main findings are that (i) the locations of features are remarkably robust against variations in the metallicity-dependent cosmic star formation history, and (ii) the low-mass end is least affected by these variations. This is promising as it increases our chances of constraining the physics that govern the formation of these objects (
https://github.com/LiekeVanSon/SFRD_fit/tree/7348a1ad0d2ed6b78c70d5100fb3cd2515493f02/
).
Our current knowledge of the star-forming metallicity of galaxies relies primarily on gas-phase oxygen abundance measurements. However, these do not always allow an accurate description of ...differences in stellar evolution and feedback, which are driven by variations in iron abundance. alpha -elements (such as oxygen) and iron are produced by sources that operate on different timescales and the link between them is not straightforward.
We explore the origin of the O/Fe --specific SFR (sSFR) relation, linking chemical abundances to galaxy formation timescales. This relation is adhered to by star-forming galaxies across redshifts according to cosmological simulations and basic theoretical expectations. Its apparent universality makes it suitable for trading the readily available oxygen for iron abundance. We show that the relation is determined by the relative iron production efficiency of core-collapse and type Ia supernovae and the delay-time distribution of the latter ---uncertain factors that could be constrained empirically with the O/Fe ---sSFR relation.
We compile and homogenise a literature sample of star-forming galaxies with observational iron abundance determinations to place first constraints on the O/Fe --sSFR relation over a wide range of sSFR. The relation shows a clear evolution towards lower O/Fe with decreasing sSFR and a flattening above log$_ $(sSFR/yr)>-9. These results are broadly consistent with expectations, but better constraints are needed to inform the models. We independently derive the relation from old Milky Way stars and find remarkable agreement between the two, as long as the recombination-line absolute oxygen abundance scale is used in conjunction with stellar metallicity measurements.
Trading oxygen for iron Chruślińska, M; Pakmor, R; Matthee, J ...
Astronomy and astrophysics (Berlin),
06/2024, Volume:
686
Journal Article
Peer reviewed
Our current knowledge of the star-forming metallicity of galaxies relies primarily on gas-phase oxygen abundance measurements. However, these do not always allow an accurate description of ...differences in stellar evolution and feedback, which are driven by variations in iron abundance. α-elements (such as oxygen) and iron are produced by sources that operate on different timescales and the link between them is not straightforward. We explore the origin of the O/Fe–specific SFR (sSFR) relation, linking chemical abundances to galaxy formation timescales. This relation is adhered to by star-forming galaxies across redshifts according to cosmological simulations and basic theoretical expectations. Its apparent universality makes it suitable for trading the readily available oxygen for iron abundance. We show that the relation is determined by the relative iron production efficiency of core-collapse and type Ia supernovae and the delay-time distribution of the latter – uncertain factors that could be constrained empirically with the O/Fe–sSFR relation. We compile and homogenise a literature sample of star-forming galaxies with observational iron abundance determinations to place first constraints on the O/Fe–sSFR relation over a wide range of sSFR. The relation shows a clear evolution towards lower O/Fe with decreasing sSFR and a flattening above log10(sSFR/yr) > − 9. These results are broadly consistent with expectations, but better constraints are needed to inform the models. We independently derive the relation from old Milky Way stars and find remarkable agreement between the two, as long as the recombination-line absolute oxygen abundance scale is used in conjunction with stellar metallicity measurements.
Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the
Gaia
space ...mission. Its light curve exhibited five distinct brightening episodes, reaching up to
I
= 12 mag, and it was covered in great detail with almost 25 000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and
Gaia
around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and
Gaia
located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57 ± 0.05
M
⊙
and 0.36 ± 0.03
M
⊙
at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by
Gaia
as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.
New observational facilities are probing astrophysical transients such as stellar explosions and gravitational wave (GW) sources at ever increasing redshifts, while also revealing new features in ...source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation history (\(\mathcal{S}(Z,z)\)) as an input. Large uncertainties remain in the shape and evolution of this function. In this work, we propose a simple analytical function for \(\mathcal{S}(Z,z)\). Variations of this function can be easily interpreted, because the parameters link to its shape in an intuitive way. We fit our analytical function to the star-forming gas of the cosmological TNG100 simulation and find that it is able to capture the main behaviour well. As an example application, we investigate the effect of systematic variations in the \(\mathcal{S}(Z,z)\) parameters on the predicted mass distribution of locally merging binary black holes (BBH). Our main findings are: I) the locations of features are remarkably robust against variations in the metallicity-dependent cosmic star formation history, and II) the low mass end is least affected by these variations. This is promising as it increases our chances to constrain the physics that governs the formation of these objects.