No Arabic abstract
The theory for single stellar evolution predicts a gap in the mass distribution of black holes (BHs) between approximately 45-130M$_{odot}$, the so-called pair-instability mass gap. We examine whether BHs can pollute the gap after accreting from a stellar companion. To this end, we simulate the evolution of isolated binaries using a population synthesis code, where we allow for super-Eddington accretion. Under our most extreme assumptions, we find that at most about 2% of all merging binary BH systems contains a BH with a mass in the pair-instability mass gap, and we find that less than 0.5% of the merging systems has a total mass larger than 90M$_{odot}$. We find no merging binary BH systems with a total mass exceeding 100M$_{odot}$. We compare our results to predictions from several dynamical pathways to pair-instability mass gap events and discuss the distinguishable features. We conclude that the classical isolated binary formation scenario will not significantly contribute to the pollution of the pair-instability mass gap. The robustness of the predicted mass gap for the isolated binary channel is promising for the prospective of placing constraints on (i) the relative contribution of different formation channels, (ii) the physics of the progenitors including nuclear reaction rates, and (iii), tentatively, the Hubble parameter.
Stellar evolution theory predicts a gap in the black hole birth function caused by the pair instability. Presupernova stars that have a core mass below some limiting value, Mlo, after all pulsational activity is finished, collapse to black holes, whereas more massive ones, up to some limiting value, Mhi, explode, promptly and completely, as pair-instability supernovae. Previous work has suggested Mlo is approximately 50 solar masses and Mhi is approximately 130 solar masses. These calculations have been challenged by recent LIGO observations that show many black holes merging with individual masses, Mlo is least some 65 solar masses. Here we explore four factors affecting the theoretical estimates for the boundaries of this mass gap: nuclear reaction rates, evolution in detached binaries, rotation, and hyper-Eddington accretion after black hole birth. Current uncertainties in reaction rates by themselves allow Mlo to rise to 64 solar masses and Mhi as large as 161 solar masses. Rapid rotation could further increase Mlo to about 70 solar masses, depending on the treatment of magnetic torques. Evolution in detached binaries and super-Eddington accretion can, with great uncertainty, increase Mlo still further. Dimensionless Kerr parameters close to unity are allowed for the more massive black holes produced in close binaries, though they are generally smaller.
We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5times 10^8M_{odot}$ black hole with accretion rates varying from $sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases with and without net vertical magnetic flux. The inner regions of all disks have radiation pressure $sim 10^4-10^6$ times the gas pressure. Non-axisymmetric density waves that steepen into spiral shocks form as gas flows towards the black hole. In simulations without net vertical magnetic flux, Reynolds stress generated by the spiral shocks are the dominant mechanism to transfer angular momentum. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when net vertical magnetic flux is sufficiently large. Outflows are formed with speed $sim 0.1-0.4c$. When the accretion rate is smaller than $sim 500 L_{Edd}/c^2$, outflows start around $10$ gravitational radii and the radiative efficiency is $sim 5%-7%$ with both magnetic field configurations. With accretion rate reaching $1500 L_{Edd}/c^2$, most of the funnel region close to the rotation axis becomes optically thick and the outflow only develops beyond $50$ gravitational radii. The radiative efficiency is reduced to $1%$. We always find the kinetic energy luminosity associated with the outflow is only $sim 15%-30%$ of the radiative luminosity. The mass flux lost in the outflow is $sim 15%-50%$ of the net mass accretion rates. We discuss implications of our simulation results on the observational properties of these disks.
By probing the population of binary black hole (BBH) mergers detected by LIGO-Virgo, we can infer properties about the underlying black hole formation channels. A mechanism known as pair-instability (PI) supernova is expected to prevent the formation of black holes from stellar collapse with mass greater than $sim 40-65,M_odot$ and less than $sim 120,M_odot$. Any BBH merger detected by LIGO-Virgo with a component black hole in this gap, known as the PI mass gap, likely originated from an alternative formation channel. Here, we firmly establish GW190521 as an outlier to the stellar-mass BBH population if the PI mass gap begins at or below $65, M_{odot}$. In addition, for a PI lower boundary of $40-50, M_{odot}$, we find it unlikely that the remaining distribution of detected BBH events, excluding GW190521, is consistent with the stellar-mass population.
We perform a binary population synthesis calculation incorporating very massive population (Pop.) III stars up to 1500 $M_odot$, and investigate the nature of binary black hole (BBH) mergers. Above the pair-instability mass gap, we find that the typical primary black hole (BH) mass is 135-340 $M_odot$. The maximum primary BH mass is as massive as 686 $M_odot$. The BBHs with both of their components above the mass gap have low effective inspiral spin $sim$ 0. So far, no conclusive BBH merger beyond the mass gap has been detected, and the upper limit on the merger rate density is obtained. If the initial mass function (IMF) of Pop. III stars is simply expressed as $xi_m(m) propto m^{-alpha}$ (single power law), we find that $alpha gtrsim 2.8$ is needed in order for the merger rate density not to exceed the upper limit. In the future, the gravitational wave detectors such as Einstein telescope and Pre-DECIGO will observe BBH mergers at high redshift. We suggest that we may be able to impose a stringent limit on the Pop. III IMF by comparing the merger rate density obtained from future observations with that derived theoretically.
We discuss properties of the ultra-luminous $X$-ray source in the galaxy M82, NuSTAR J095551+6940.8, containing an accreting neutron star. The neutron star has surface magnetic field $ B_{NS} approx 1.4 times 10^{13 } , {rm G}$ and experiences accretion rate of $9 times 10^{-7} M_odot {rm , yr}^{-1} $. The magnetospheric radius, close to the corotation radius, is $sim 2 times 10^8$ cm. The accretion torque on the neutron star is reduce well below what is expected in a simple magnetospheric accretion due to effective penetration of the stellar magnetic field into the disk beyond the corotation radius. As a result, the radiative force of the surface emission does not lead to strong coronal wind, but pushes plasma along magnetic field lines towards the equatorial disk. The neutron star is nearly an orthogonal rotator, with the angle between the rotation axis and the magnetic moment $geq 80$ degrees. Accretion occurs through optically thick -- geometrically thin and flat accretion curtain, which cuts across the polar cap. High radiation pressure from the neutron star surface is nevertheless smaller than that the ram pressure of the accreting material flowing through the curtain, and thus fails to stop the accretion. At distances below few stellar radii the magnetic suppression of the scattering cross-section becomes important. The $X$-ray luminosity (pulsed and persistent components) comes both from the neutron star surface as a hard $X$-ray component and as a soft component from reprocessing by the accretion disk.