Aims.
We present an open-access database that includes a synthetic catalog of black holes (BHs) in the Milky Way, divided by the components disk, bulge, and halo.
Methods.
To calculate the evolution ...of single and binary stars, we used the updated population synthesis code StarTrack. We applied a new model of the star formation history and chemical evolution of Galactic disk, bulge, and halo that was synthesized from observational and theoretical data. This model can be easily employed for other studies of population evolution.
Results.
We find that at the current Milky Way (disk+bulge+halo) contains about 1.2 × 10
8
single BHs with an average mass of about 14
M
⊙
, and 9.3 × 10
6
BHs in binary systems with an average mass of 19
M
⊙
. We present basic statistical properties of the BH population in three Galactic components such as the distributions of BH masses, velocities, or the numbers of BH binary systems in different evolutionary configurations.
Conclusions.
The metallicity of a stellar population has a significant effect on the final BH mass through the stellar winds. The most massive single BH in our simulation of 113
M
⊙
originates from a merger of a BH and a helium star in a low-metallicity stellar environment in the Galactic halo. We constrain that only ∼0.006% of the total Galactic halo mass (including dark matter) can be hidden in the form of stellar origin BHs. These BHs cannot be detected by current observational surveys. We calculated the merger rates for current Galactic double compact objects (DCOs) for two considered common-envelope models: ∼3–81 Myr
−1
for BH-BH, ∼1–9 Myr
−1
for BH-neutron star (NS), and ∼14–59 Myr
−1
for NS-NS systems. We show the evolution of the merger rates of DCOs since the formation of the Milky Way until the current moment with the new star formation model of the Galaxy.
ABSTRACT
Accreting neutron stars (NSs) are one of the main targets for continuous gravitational wave searches, as asymmetric accretion may lead to quadrupolar deformations, or ‘mountains’, on the ...crust of the star, which source gravitational wave (GW) emission at twice the rotation frequency. The GW torque may also impact on the spin evolution of the star, possibly dictating the currently observed spin periods of NSs in low-mass X-ray binaries and leading to the increased spin-down rate observed during accretion in PSR J1023+0038. Previous studies have shown that deformed reaction layers in the crust of the NS lead to thermal and compositional gradients that can lead to GW emission. However, there are no realistic constraints on the level of asymmetry that is expected. In this paper, we consider a natural source of asymmetry, namely the magnetic field, and calculate the density and pressure perturbations that are expected in the crust of accreting NSs. In general, we find that only the outermost reaction layers of the NS are strongly perturbed. The mass quadrupole that we estimate is generally small and cannot explain the increase of spin-down rate of PSR J1023+0038. However, if strong shallow heating sources are present at low densities in the crust, as cooling observations suggest, these layers will be strongly perturbed and the resulting quadrupole could explain the observed spin-down of PSR J1023+0038, and lead to observable GW signals from systems with higher accretion rates.
ABSTRACT
The expected volume of data from the third-generation gravitational waves (GWs) Einstein Telescope (ET) detector would make traditional GWs search methods such as match filtering ...impractical. This is due to the large template bank required and the difficulties in waveforms modelling. In contrast, machine learning (ML) algorithms have shown a promising alternative for GWs data analysis, where ML can be used in developing semi-automatic and automatic tools for the detection and parameter estimation of GWs sources. Compared to second generation detectors, ET will have a wider accessible frequency band but also a lower noise. The ET will have a detection rate for Binary Black Holes (BBHs) and Binary Neutron Stars (BNSs) of the order of 105–106 and 7 × 104 yr−1, respectively. We explored the efficiency of using convolutional neural networks (CNNs) for the detection of BBHs’ mergers in synthetic noisy data that was generated according to ET’s parameters. Without performing data whitening or applying bandpass filtering, we trained four CNN networks with the state-of-the-art performance in computer vision, namely VGG, ResNet, and DenseNet. ResNet has significantly better performance, and was able to detect BBHs sources with SNR of 8 or higher with 98.5 per cent accuracy, and with 92.5 per cent, 85 per cent, 60 per cent, and 62 per cent accuracy for sources with SNR range of 7–8, 6–7, 5–6, and 4–5, respectively. ResNet, in qualitative evaluation, was able to detect a BBH’s merger at 60 Gpc with 4.3 SNR. It was also shown that CNN can be used efficiently for near-real time detection of BBHs.
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).
Context. Mergers of two stellar-origin black holes are a prime source of gravitational waves and are under intensive investigation. One crucial ingredient in their modeling has been neglected: ...pair-instability pulsation supernovae with associated severe mass loss may suppress the formation of massive black holes, decreasing black-hole-merger rates for the highest black-hole masses. Aims. We demonstrate the effects of pair-instability pulsation supernovae on merger rate and mass using populations of double black-hole binaries formed through the isolated binary classical evolution channel. Methods. The mass loss from pair-instability pulsation supernova is estimated based on existing hydrodynamical calculations. This mass loss is incorporated into the StarTrack population synthesis code. StarTrack is used to generate double black-hole populations with and without pair-instability pulsation supernova mass loss. Results. The mass loss associated with pair-instability pulsation supernovae limits the Population I/II stellar-origin black-hole mass to 50 M⊙, in tension with earlier predictions that the maximum black-hole mass could be as high as 100 M⊙. In our model, neutron stars form with mass 1−2 M⊙. We then encounter the first mass gap at 2−5 M⊙ with the compact object absence due to rapid supernova explosions, followed by the formation of black holes with mass 5−50 M⊙, with a second mass gap at 50−135 M⊙ created by pair-instability pulsation supernovae and by pair-instability supernovae. Finally, black holes with masses above 135 M⊙ may potentially form to arbitrarily high mass limited only by the extent of the initial mass function and the strength of stellar winds. Suppression of double black-hole-merger rates by pair-instability pulsation supernovae is negligible for our evolutionary channel. Our standard evolutionary model, with the inclusion of pair-instability pulsation supernovae and pair-instability supernovae, is fully consistent with the Laser Interferometric Gravitational-wave Observatory (LIGO) observations of black-hole mergers: GW150914, GW151226, and LVT151012. The LIGO results are inconsistent with high (≳ 400 km s-1) black hole (BH) natal kicks. We predict the detection of several, and up to as many as ~60, BH-BH mergers with a total mass of 10−150 M⊙ (most likely range: 20−80 M⊙) in the forthcoming ~60 effective days of the LIGO O2 observations, assuming the detectors reach the optimistic target O2 sensitivity.
Context. Current star formation models imply that the binary fraction of Population III stars is non-zero. The evolution of these binaries must have led to the formation of compact object binaries. ...Aims. We estimate the gravitational wave background originating in these binaries and discuss its observability. Methods. The properties of the Population III binaries are investigated using a binary population synthesis code. We numerically model the background and take into account the evolution of eccentric binaries. Results. The gravitational wave background from Population III binaries dominates the spectrum below 100 Hz. If the binary fraction is larger than 10-2, the background will be detectable by Einstein Telescope (ET), Laser Interferometer Space Antenna (LISA), and DECi-Hertz Interferometer Gravitational wave Observatory (DECIGO). Conclusions. The gravitational wave background from Population III binaries will dominate the spectrum below 100 Hz. The instruments LISA, ET, and DECIGO should either see it easily or, in the case of non-detection, provide very strong constraints on the properties of the Population III stars.
Context. Recent studies on stellar evolution have shown that the properties of compact objects strongly depend on the metallicity of the environment in which they were formed. Aims. Using some very ...simple assumptions on the metallicity of the stellar populations, we explore how this property affects the unresolved gravitational-wave background from extragalactic compact binaries. Methods. We obtained a suit of models using population synthesis code, estimated the gravitational-wave background they produce, and discuss its detectability with second- (advanced LIGO, advanced Virgo) and third- (Einstein Telescope) generation detectors. Results. Our results show that the background is dominated by binary black holes for all considered models in the frequency range of terrestrial detectors, and that it could be detected in most cases by advanced LIGO/Virgo, and with Einstein Telescope with a very high signal-to-noise ratio. The observed peak in a gravitational-wave spectrum depends on the metallicity of the stellar population.
Population synthesis of double neutron stars Osłowski, S.; Bulik, T.; Gondek-Rosińska, D. ...
Monthly Notices of the Royal Astronomical Society,
20/May , Volume:
413, Issue:
1
Journal Article
Peer reviewed
Open access
Using the StarTrack binary population synthesis code we model the population of double neutron stars in the Galaxy. We include a detailed treatment of the spin evolution of each pulsar due to ...processes such as spin-down and spin-up during accretion events as well as magnetic field decay. We also model the spatial distribution of double neutron stars by including their natal kicks and subsequent propagation in the Galactic gravitational potential. This synthetic pulsar population is compared to the observed sample of double neutron stars taking into account the selection effects of detection in the radio band, to determine the most likely evolutionary parameters. With these parameters we determine the properties of the double neutron star binaries detectable in gravitational waves by the high-frequency interferometers LIGO and VIRGO. In particular, we discuss the distributions of chirp masses and mass ratios in samples selected by their radio or gravitational wave emission.
Context. The current gravitational wave detectors have reached their operational sensitivity and are nearing detection of compact object binaries. In the coming years, we expect that the Advanced ...LIGO/VIRGO will start taking data. At the same time, there are plans for third generation ground-based detectors such as the Einstein Telescope, and space detectors such as DECIGO. Aims. We discuss the eccentricity distribution of inspiral compact object binaries during they inspiral phase. We analyze the expected distributions of eccentricities at three frequencies that are characteristic of three future detectors: Advanced LIGO/VIRGO (30 Hz), Einstein Telescope (3 Hz), and DECIGO (0.3 Hz). Methods. We use the StarTrack binary population code to investigate the properties of the population of compact binaries in formation. We evolve their orbits until the point that they enter a given detector sensitivity window and analyze the eccentricity distribution at that time. Results. We find that the eccentricities of BH-BH and BH-NS binaries are quite small when entering the Advanced LIGO/VIRGO detector window for all considered models of binary evolution. Even in the case of the DECIGO detector, the typical eccentricities of BH-BH binaries are below 10-4, and the BH-NS eccentricities are smaller than 10-3. Some fraction of NS-NS binaries may have significant eccentricities. Within the range of considered models, we found that a fraction of between 0.2% and 2% NS-NS binaries will have an eccentricity above 0.01 for the Advanced LIGO/VIRGO detectors. For the ET detector, this fraction is between 0.4% and 4%, and for the DECIGO detector it lies between 2% and 27%.