No Arabic abstract
We present the results of the stellar feedback from Pop III binaries by employing improved, more realistic Pop III evolutionary stellar models. To facilitate a meaningful comparison, we consider a fixed mass of 60 solar masses (Msun) incorporated in Pop III stars, either contained in a single star, or split up in binary stars of 30 Msun each or an asymmetric case of one 45 Msun and one 15 Msun star. Whereas the sizes of the resulting HII regions are comparable across all cases, the HeIII regions around binary stars are significantly smaller than that of the single star. Consequently, the He$^{+}$ 1640 angstrom recombination line is expected to become much weaker. Supernova feedback exhibits great variety due to the uncertainty in possible explosion pathways. If at least one of the component stars dies as a hypernova about ten times more energetic than conventional core-collapse supernovae, the gas inside the host minihalo is effectively blown out, chemically enriching the intergalactic medium (IGM) to an average metallicity of $10^{-4}-10^{-3}$ solar metallicity (Zsun), out to $sim 2$ kpc. The single star, however, is more likely to collapse into a black hole, accompanied by at most very weak explosions. The effectiveness of early chemical enrichment would thus be significantly reduced, in difference from the lower mass binary stars, where at least one component is likely to contribute to heavy element production and dispersal. Important new feedback physics is also introduced if close binaries can form high-mass x-ray binaries, leading to the pre-heating and -ionization of the IGM beyond the extent of the stellar HII regions.
Superluminous supernovae have been proposed to arise from Population III progenitors that explode as pair-instability supernovae. Pop III stars are the first generation of stars in the Universe, and are thought to form as late as $z sim 6$. Future near-infrared imaging facilities such as ULTIMATE-Subaru can potentially detect and identify these PISNe with a dedicated survey. Gravitational lensing by intervening structure in the Universe can aid in the detection of these rare objects by magnifying the high-$z$ source population into detectability. We perform a mock survey with ULTIMATE-Subaru, taking into account lensing by line-of-sight structure to evaluate its impact on the predicted detection rate. We compare a LOS mass reconstruction using observational data from the Hyper Suprime Cam survey to results from cosmological simulations to test their consistency in calculating the magnification distribution in the Universe to high-$z$, but find that the data-based method is still limited by an inability to accurately characterize structure beyond $z sim1.2$. We also evaluate a survey strategy of targeting massive galaxy clusters to take advantage of their large areas of high magnification. We find that targeting clusters can result in a gain of a factor of $sim$two in the predicted number of detected PISNe at $z > 5$, and even higher gains with increasing redshift, given our assumed survey parameters. For the highest-redshift sources at $z sim 7-9$, blank field surveys will not detect any sources, and lensing magnification by massive clusters will be necessary to observe this population.
We explore the hardening of a massive black hole binary embedded in a circum-binary gas disc when the binary and the gas are coplanar and the gas is counter-rotating. The secondary black hole, revolving in the direction opposite to the gas, experiences a drag from gas-dynamical friction and from direct accretion of part of it. Using two-dimensional (2D) hydrodynamical grid simulations we investigate the effect of changing the accretion prescriptions on the dynamics of the secondary black hole which in turn affect the binary hardening and eccentricity evolution. We find that realistic accretion prescriptions lead to results that differ from those inferred assuming accretion of all the gas within the Roche Lobe of the secondary black hole. Different accretion prescriptions result in different discs surface densities which alter the black holes dynamics back. Full 3D SPH realizations of a number of representative cases, run over a shorter interval of time, validate the general trends observed in the less computationally demanding 2D simulations. Initially circular black hole binaries increase only slightly their eccentricity which then oscillates around small values (<0.1) while they harden. By contrast, initially eccentric binaries become more and more eccentric. A semi-analytical model describing the black holes dynamics under accretion only explores the late evolution stages of the binary in an otherwise unperturbed retrograde disc to illustrate how eccentricity evolves with time in relation to the shape of the underlying surface density distribution.
The spatial and velocity distributions of dark matter particles in the Milky Way Halo affect the signals expected to be observed in searches for dark matter. Results from direct detection experiments are often analyzed assuming a simple isothermal distribution of dark matter, the Standard Halo Model (SHM). Yet there has been skepticism regarding the validity of this simple model due to the complicated gravitational collapse and merger history of actual galaxies. In this paper we compare the SHM to the results of cosmological hydrodynamical simulations of galaxy formation to investigate whether or not the SHM is a good representation of the true WIMP distribution in the analysis of direct detection data. We examine two Milky Way-like galaxies from the MaGICC cosmological simulations (a) with dark matter only and (b) with baryonic physics included. The inclusion of baryons drives the shape of the DM halo to become more spherical and makes the velocity distribution of dark matter particles less anisotropic especially at large heliocentric velocities, thereby making the SHM a better fit. We also note that we do not find a significant disk-like rotating dark matter component in either of the two galaxy halos with baryons that we examine, suggesting that dark disks are not a generic prediction of cosmological hydrodynamical simulations. We conclude that in the Solar neighborhood, the SHM is in fact a good approximation to the true dark matter distribution in these cosmological simulations (with baryons) which are reasonable representations of the Milky Way, and hence can also be used for the purpose of dark matter direct detection calculations.
Primordial supermassive stars (SMSs) formed in atomic-cooling halos at z ~ 15 - 20 are leading candidates for the seeds of the first quasars. Past numerical studies of the evolution of SMSs have typically assumed constant accretion rates rather than the highly variable flows in which they form. We model the evolution of SMSs in the cosmological flows that create them using the Kepler stellar evolution and implicit hydrodynamics code. We find that they reach masses of 1 - 2 x $10^5 M_{odot}$ before undergoing direct-collapse to black holes (DCBHs) during or at the end of their main-sequence hydrogen burning, at 1 - 1.5 Myr, regardless of halo mass, spin, or merger history. We also find that realistic, highly-variable accretion histories allow for a much greater diversity of supermassive stellar structures, including in some cases largely thermally relaxed objects, which may provide a significant source of radiative feedback. Our models indicate that the accretion histories predicted for purely atomic-cooling halos may impose a narrow spectrum of masses on the seeds of the first massive quasars, however further studies incorporating realistic feedback will be essential in order to confirm whether or not this holds true in all cases. Our results also indicate that multiple SMSs at disparate stages of evolution can form in these halos, raising the possibility of SMS binaries and supermassive X-ray binaries (SMXBs), as well as DCBH mergers which could be detected by LISA.
Supermassive black holes (SMBHs) are ubiquitous in galaxies with a sizable mass. It is expected that a pair of SMBHs originally in the nuclei of two merging galaxies would form a binary and eventually coalesce via a burst of gravitational waves. So far theoretical models and simulations have been unable to predict directly the SMBH merger timescale from ab-initio galaxy formation theory, focusing only on limited phases of the orbital decay of SMBHs under idealized conditions of the galaxy hosts. The predicted SMBH merger timescales are long, of order Gyrs, which could be problematic for future gravitational wave searches. Here we present the first multi-scale $Lambda$CDM cosmological simulation that follows the orbital decay of a pair of SMBHs in a merger of two typical massive galaxies at $zsim3$, all the way to the final coalescence driven by gravitational wave emission. The two SMBHs, with masses $sim10^{8}$ M$_{odot}$, settle quickly in the nucleus of the merger remnant. The remnant is triaxial and extremely dense due to the dissipative nature of the merger and the intrinsic compactness of galaxies at high redshift. Such properties naturally allow a very efficient hardening of the SMBH binary. The SMBH merger occurs in only $sim10$ Myr after the galactic cores have merged, which is two orders of magnitude smaller than the Hubble time.